CN105373660A - 基于等效电路的晶体管可靠性表征方法 - Google Patents

基于等效电路的晶体管可靠性表征方法 Download PDF

Info

Publication number
CN105373660A
CN105373660A CN201510772257.3A CN201510772257A CN105373660A CN 105373660 A CN105373660 A CN 105373660A CN 201510772257 A CN201510772257 A CN 201510772257A CN 105373660 A CN105373660 A CN 105373660A
Authority
CN
China
Prior art keywords
transistor
electrical circuit
sensitive parameter
stress
equivalent electrical
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.)
Pending
Application number
CN201510772257.3A
Other languages
English (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.)
Chengdu Gastone Technology Co Ltd
Original Assignee
Chengdu Gastone Technology Co Ltd
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 Chengdu Gastone Technology Co Ltd filed Critical Chengdu Gastone Technology Co Ltd
Priority to CN201510772257.3A priority Critical patent/CN105373660A/zh
Publication of CN105373660A publication Critical patent/CN105373660A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Individual Semiconductor Devices (AREA)
  • Junction Field-Effect Transistors (AREA)

Abstract

本发明提供一种基于等效电路的晶体管可靠性表征方法,包括以下步骤:S1、根据晶体管类型确定晶体管的等效电路;S2、确定用于寿命加速试验的应力的类型和大小,并确定等效电路中的敏感参数退化模型;S3、对晶体管进行寿命加速测试;S4、定时采集晶体管的测试数据,提取等效电路敏感参数,得到等效电路敏感参数随应力和时间的变化关系;S5、对所述敏感参数退化模型中的待定常数进行拟合,得到完整敏感参数退化模型和完整等效电路;S6、对晶体管进行失效机制分析和可靠性分析。本发明可用于模拟器件失效机制对电路性能的影响,可以分析晶体管的失效机制,指导工艺改进,还可以预测晶体管的性能退化量和失效时间,缩短可靠性试验时间。

Description

基于等效电路的晶体管可靠性表征方法
技术领域
本发明属于可靠性评价技术领域,具体涉及一种基于等效电路的晶体管可靠性表征方法。
背景技术
晶体管器件的可靠性是晶体管应用的关键,是产品使用之前必须要解决的问题。目前常用的表征晶体管可靠性的方法,主要是通过可靠性寿命加速测试,根据测试数据随时间的变化趋势,外推得到器件在正常工作条件下的平均失效时间(MTTF)。晶体管器件的失效机制主要包括栅退化、欧姆接触退化、热电子效应、沟道退化等。晶体管工作时,造成器件退化和失效的原因并不只是一种失效机制的作用,而往往是多种失效机制共同作用。这种基于可靠性寿命加速测试的器件可靠性表征方法存在如下一些缺点,首先,它很难用于分析晶体管具体的失效机制;其次,不能反应晶体管性能退化对电路或者系统性能的影响;再者,需要长时间的应力加速测试,效率较低。
发明内容
本发明的目的在于针对现有技术的不足,提供一种基于等效电路的晶体管可靠性表征方法,该方法可以很好地解决现有表征方法难以分析失效机制、效率较低的问题。
为达到上述要求,本发明采取的技术方案是:提供一种基于等效电路的晶体管可靠性表征方法,包括以下步骤:
S1、根据晶体管类型确定晶体管的等效电路;
S2、确定用于寿命加速试验的应力的类型和大小,并根据应力的类型确定等效电路中的敏感参数退化模型;
S3、根据应力的类型和大小对晶体管进行寿命加速测试;
S4、定时采集晶体管的测试数据,根据测试数据提取等效电路敏感参数,得到等效电路敏感参数随应力和时间的变化关系;
S5、根据变化关系对所述敏感参数退化模型中的待定常数进行拟合,得到完整敏感参数退化模型,并将完整敏感参数退化模型带入到所述等效电路中,得到晶体管的完整等效电路;
S6、根据完整敏感参数退化模型对晶体管进行失效机制分析,并根据完整等效电路对晶体管进行电路可靠性分析。
与现有技术相比,本发明具有以下优点:
(1)通过可靠性测量试验中提取的晶体管敏感参数随应力和时间的变化趋势,可以用来分析器件的失效机制,指导工艺改进;
(2)可以根据提取的晶体管敏感参数的变化趋势,建立等效电路外推模型,可预测晶体管的性能退化量和失效时间,从而缩短可靠性试验时间,提高效率;
(3)将该方法建立的等效电路嵌入常用的电路仿真软件,可以模拟晶体管失效机制对电路乃至整个系统的影响,预测电路和系统的性能退化量和失效时间。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,在这些附图中使用相同的参考标号来表示相同或相似的部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1为本发明的流程示意图;
图2为本发明实施例的GaAspHEMT典型的非线性等效电路拓扑。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,以下结合附图及具体实施例,对本申请作进一步地详细说明。为简单起见,以下描述中省略了本领域技术人员公知的某些技术特征。
本发明提供一种基于等效电路的晶体管可靠性表征方法,如图1所示,包括以下步骤:
步骤S1:根据晶体管类型确定晶体管的等效电路。
晶体管类型包括双极晶体管,如BJT、HBT,场效应晶体管,如MOSFET、MESFET、HEMT等,晶体管的材料包括Si、GaAs、GaN、InP等;每种类型的晶体管对应的等效电路拓扑和参数表达式均不同。
本实施例以GaAspHEMT晶体管为例进行阐述,图2为GaAspHEMT典型的非线性等效电路拓扑,忽略了各电极的寄生电感和电容,图2中各个参数的意义如下:
Ids(Vgs,Vds)——受漏、栅电压控制的漏极电流源;
Qgd(Vgs,Vds)——受漏、栅电压控制的漏栅非线性储存电荷;
Qgs(Vgs,Vds)——受漏、栅电压控制的栅源非线性储存电荷;
Qds(Vds)——受漏、源电压控制的非线性储存电荷;
Rin(Vgs,Vds)——受漏、源电压控制的非线性电阻;
Igd、Igs——栅极肖特基二极管电流;
Rs、Rg、Rd——源极、栅极、漏极寄生电阻。
步骤S2:确定用于寿命加速试验的应力的类型和大小,并根据应力的类型确定所述等效电路中的敏感参数退化模型。
晶体管可靠性寿命加速试验所施加的应力类型一般有温度、电压及射频信号等,可以单个应力单独作用或者多种应力同时作用,应力的大小根据试验条件选取。本实施例选取温度为GaAspHEMT寿命加速试验的应力类型,将待测器件分为三组样品,分别施加不同的温度应力,大小分别为170℃、200℃和230℃。
晶体管等效电路中的敏感参数是指对晶体管性能影响最为明显的参数,同时,这些敏感参数也能直观的反应晶体管器件的失效机制。例如,对HEMT器件而言,它的敏感参数包括漏源电流Ids、跨导gm、栅极泄漏电流Ig、栅电阻Rg、源电阻Rs等;这些敏感参数与晶体管失效机制的关系有的是一对一的关系,有的是一对多的关系,需要具体分析,例如,栅极泄漏电流Ig可以反应栅肖特基接触退化的影响,Rs可以反应源极欧姆接触退化的影响,Ids可以反应沟道退化及其它失效机制的综合影响等。
描述晶体管敏感参数退化模型主要有Arrhenius模型、逆幂率模型和Eyring模型,其中,Arrhenius模型只能处理热应力,逆幂率模型只能处理电应力,而Eyring模型及包括热应力也包括电应力,在实际使用时,需要根据寿命加速试验所选取的应力类型来选取适合的敏感参数退化模型。三种敏感参数退化模型的表达式分别为:
Arrhenius模型:
A=a1·eE/KT(1)
式中,a1是一个正常数,E是激活能,K是波尔兹曼常数,T为绝对温度。
逆幂率模型:
A=a2·V-C(2)
式中,a2是一个常数,V是电压应力,C是一个与激活能相关的正常数。
Eyring模型:
A = a 3 T exp ( a 4 K T ) exp ( a 5 + a 6 K T ) - - - ( 3 )
式中,a3,a4,a5,a6是待定常数,由于此方法待定系数较多,所以所需应力类型和测试样品数量也较多。
另外,根据物质的反应速率关系式,可以将等效电路中的敏感参数随时间的退化关系表示为:
B ( t ) = b · e - t τ - - - ( 4 )
式中,b称为退化幅度常数,τ称为退化时间常数。
因此,晶体管等效电路中的敏感参数退化模型表示为:
F(t)=A·B(t)(5)
GaAspHEMT晶体管等效电路中的敏感参数主要有漏极电流Ids、源电阻Rs、栅电阻Rg。由于选取的应力类型为温度应力,因此采用Arrhenius模型来描述GaAspHEMT晶体管敏感参数随应力的退化关系,由式(5)得到GaAspHEMT晶体管敏感参数退化模型可以表示为:
I d s = I d s 0 · e E i K ( 1 T - 1 T 0 ) · e - t τ i - - - ( 6 )
R s = R s 0 · e E s K ( 1 T - 1 T 0 ) · e - t τ s - - - ( 7 )
R g = R g 0 · e E g K ( 1 T - 1 T 0 ) · e - t τ g - - - ( 8 )
式中,T0为常温,Ids0、Rs0、Rg0为各个敏感参数在常温下的初始值,Ei、Es、Eg、τi、τs、τg为参数退化模型的待定常数。
步骤S3:根据应力的类型和大小对晶体管进行寿命加速测试。
步骤S4:定时采集晶体管的测试数据,根据测试数据提取等效电路敏感参数,得到等效电路敏感参数随应力和时间的变化关系。
对待测晶体管施加应力,进行寿命加速试验,定时采集晶体管的DC和RF测试数据,提取等效电路参数,得到等效电路敏感参数随应力和时间的变化关系。
根据步骤S2确定的GaAspHEMT晶体管寿命加速试验的应力类型和大小,对三组待测样品进行寿命加速试验。试验总时长约400~700小时,每24小时采集一次晶体管的DC和RF测试数据,测量GaAspHEMT晶体管的IV曲线,漏极电压Vds范围为0~8V,栅极电压Vgs范围为-1.2~0.5V,并测量不同偏置点下的小信号S参数,频率为0.5~40GHz。图2所示GaAspHEMT非线性等效电路,可以选取EE-HEMT、Angelov等非线性模型来描述等效电路中的非线性受控电流源和非线性存储电荷,然后根据GaAspHEMT晶体管非线性模型参数提取流程,提取图2等效电路中的模型参数值,并可以得到等效电路敏感参数随应力和时间的变化关系。
步骤S5:根据变化关系对所述敏感参数退化模型中的待定常数进行拟合,得到完整敏感参数退化模型,并将完整敏感参数退化模型带入到等效电路中,得到晶体管的完整等效电路。
根据步骤S4得到的GaAspHEMT等效电路敏感参数随应力和时间的变化关系,对等效电路中的敏感参数退化模型中的待定常数进行拟合,确定式(6)~(8)中的Ei、Es、Eg、τi、τs、τg各参数值,得到GaAspHEMT晶体管等效电路中的完整敏感参数退化模型,将该敏感参数退化模型代入到等效电路中,得到GaAspHEMT晶体管的完整等效电路。
步骤S6:根据完整敏感参数退化模型对晶体管进行失效机制分析,并根据完整等效电路对晶体管进行电路可靠性分析。
根据步骤S5中得到的GaAspHEMT晶体管等效电路中的完整敏感参数退化模型,可用来分析GaAspHEMT晶体管的失效机制,栅极泄漏电流Ig可以反应栅肖特基接触退化的影响,Rs可以反应源极欧姆接触退化的影响,Ids可以反应沟道退化及其它失效机制的综合影响等;同时,步骤S5中得到的GaAspHEMT晶体管完整等效电路,可计算得到GaAspHEMT晶体管性能参数在不同应力条件下随时间的退化量,如增益、效率、输出功率等,从而外推晶体管的长期退化规律,预测晶体管的失效时间;另外,将步骤S5得到的GaAspHEMT晶体管完整等效电路嵌入电路设计仿真软件AgilentADS中,用来模拟晶体管失效机制对电路乃至整个系统的影响,预测电路和系统的性能退化量和失效时间。
本发明建立了各种失效机制与器件性能之间的联系,可用于模拟器件失效机制对电路或系统性能的影响;并且可以根据晶体管参数随应力和时间的变化趋势分析晶体管的失效机制,指导工艺改进;另外,该方法还可以预测晶体管的性能退化量和失效时间,从而缩短可靠性试验时间,提高效率。
以上所述实施例仅表示本发明的几种实施方式,其描述较为具体和详细,但并不能理解为对本发明范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明保护范围。因此本发明的保护范围应该以所述权利要求为准。

Claims (3)

1.基于等效电路的晶体管可靠性表征方法,其特征在于,包括以下步骤:
S1、根据晶体管类型确定晶体管的等效电路;
S2、确定用于寿命加速试验的应力的类型和大小,并根据所述应力的类型确定所述等效电路中的敏感参数退化模型;
S3、根据应力的类型和大小对晶体管进行寿命加速测试;
S4、定时采集晶体管的测试数据,根据测试数据提取等效电路敏感参数,得到所述等效电路敏感参数随应力和时间的变化关系;
S5、根据所述变化关系对所述敏感参数退化模型中的待定常数进行拟合,得到完整敏感参数退化模型,并将所述完整敏感参数退化模型带入到所述等效电路中,得到晶体管的完整等效电路;
S6、根据所述完整敏感参数退化模型对晶体管进行失效机制分析,并根据所述完整等效电路对晶体管进行电路可靠性分析。
2.根据权利要求1所述的基于等效电路的晶体管可靠性表征方法,其特征在于,所述步骤S2中应力的类型包括温度和/或电压和/或射频信号。
3.根据权利要求1所述的基于等效电路的晶体管可靠性表征方法,其特征在于,所述步骤S3的敏感参数退化模型包括处理热应力的Arrhenius模型、处理电应力的逆幂率模型和同时处理热应力和电应力的Eyring模型。
CN201510772257.3A 2015-11-12 2015-11-12 基于等效电路的晶体管可靠性表征方法 Pending CN105373660A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510772257.3A CN105373660A (zh) 2015-11-12 2015-11-12 基于等效电路的晶体管可靠性表征方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510772257.3A CN105373660A (zh) 2015-11-12 2015-11-12 基于等效电路的晶体管可靠性表征方法

Publications (1)

Publication Number Publication Date
CN105373660A true CN105373660A (zh) 2016-03-02

Family

ID=55375855

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510772257.3A Pending CN105373660A (zh) 2015-11-12 2015-11-12 基于等效电路的晶体管可靠性表征方法

Country Status (1)

Country Link
CN (1) CN105373660A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105844059A (zh) * 2016-04-19 2016-08-10 成都海威华芯科技有限公司 一种微波大功率晶体管建模方法
CN105956228A (zh) * 2016-04-19 2016-09-21 成都海威华芯科技有限公司 高效率开关类功率放大器的晶体管建模方法
CN108268717A (zh) * 2018-01-18 2018-07-10 上海华虹宏力半导体制造有限公司 一种新型栅诱导漏极泄漏电流模型
CN108445371A (zh) * 2018-01-18 2018-08-24 国网浙江省电力公司舟山供电公司 绝缘栅双极型晶体管使用寿命预分拣方法
CN109190244A (zh) * 2018-08-31 2019-01-11 哈尔滨工业大学 一种n沟道vdmosfet在pbti效应下退化模型的建立方法
CN111832226A (zh) * 2020-07-13 2020-10-27 中国南方电网有限责任公司超高压输电公司柳州局 一种基于辅助粒子滤波的igbt剩余寿命估算方法
CN112782558A (zh) * 2020-12-29 2021-05-11 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) 集成电路失效率获取方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5838947A (en) * 1996-04-02 1998-11-17 Synopsys, Inc. Modeling, characterization and simulation of integrated circuit power behavior
US6112275A (en) * 1990-05-15 2000-08-29 Dallas Semiconductor Corporation Method of communicating over a single wire bus between a host device and a module device which measures thermal accumulation over time
US20120122230A1 (en) * 2010-11-17 2012-05-17 Exxonmobil Research And Engineering Company Methods for identifying high fouling hydrocarbon and for mitigating fouling of process equipment
CN102567560A (zh) * 2010-12-27 2012-07-11 北京圣涛平试验工程技术研究院有限责任公司 Mos器件寿命预计方法及系统
CN103063995A (zh) * 2011-10-21 2013-04-24 北京大学 一种预测soi mosfet器件可靠性寿命的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6112275A (en) * 1990-05-15 2000-08-29 Dallas Semiconductor Corporation Method of communicating over a single wire bus between a host device and a module device which measures thermal accumulation over time
US5838947A (en) * 1996-04-02 1998-11-17 Synopsys, Inc. Modeling, characterization and simulation of integrated circuit power behavior
US20120122230A1 (en) * 2010-11-17 2012-05-17 Exxonmobil Research And Engineering Company Methods for identifying high fouling hydrocarbon and for mitigating fouling of process equipment
CN102567560A (zh) * 2010-12-27 2012-07-11 北京圣涛平试验工程技术研究院有限责任公司 Mos器件寿命预计方法及系统
CN103063995A (zh) * 2011-10-21 2013-04-24 北京大学 一种预测soi mosfet器件可靠性寿命的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
万宁 等: "《AlGaAs/InGaAs PHEMT 栅电流参数退化模型研究》", 《物理学报》 *
白云霞 等: "《基于Arrhenius模型快速评价功率VDMOS可靠性》", 《半导体技术》 *
罗俊 等: "《微纳米CMOS VLSI电路可靠性仿真与设计》", 《微电子学》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105844059A (zh) * 2016-04-19 2016-08-10 成都海威华芯科技有限公司 一种微波大功率晶体管建模方法
CN105956228A (zh) * 2016-04-19 2016-09-21 成都海威华芯科技有限公司 高效率开关类功率放大器的晶体管建模方法
CN105844059B (zh) * 2016-04-19 2019-03-12 成都海威华芯科技有限公司 一种微波大功率晶体管建模方法
CN105956228B (zh) * 2016-04-19 2019-03-12 成都海威华芯科技有限公司 高效率开关类功率放大器的晶体管建模方法
CN108268717A (zh) * 2018-01-18 2018-07-10 上海华虹宏力半导体制造有限公司 一种新型栅诱导漏极泄漏电流模型
CN108445371A (zh) * 2018-01-18 2018-08-24 国网浙江省电力公司舟山供电公司 绝缘栅双极型晶体管使用寿命预分拣方法
CN109190244A (zh) * 2018-08-31 2019-01-11 哈尔滨工业大学 一种n沟道vdmosfet在pbti效应下退化模型的建立方法
CN111832226A (zh) * 2020-07-13 2020-10-27 中国南方电网有限责任公司超高压输电公司柳州局 一种基于辅助粒子滤波的igbt剩余寿命估算方法
CN111832226B (zh) * 2020-07-13 2024-04-26 中国南方电网有限责任公司超高压输电公司柳州局 一种基于辅助粒子滤波的igbt剩余寿命估算方法
CN112782558A (zh) * 2020-12-29 2021-05-11 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) 集成电路失效率获取方法

Similar Documents

Publication Publication Date Title
CN105373660A (zh) 基于等效电路的晶体管可靠性表征方法
CN102542077B (zh) AlGaN/GaN HEMT小信号模型的参数提取方法
CN106529102B (zh) AlGaN/GaN HEMT小信号模型及其参数的提取方法
CN103063995B (zh) 一种预测soi mosfet器件可靠性寿命的方法
CN102338846B (zh) 一种GaN基HEMT器件的可靠性评估方法
CN110376500A (zh) 一种功率mos器件开启过程中瞬态温升在线测量方法
CN102393501A (zh) 一种mosfet可靠性测试分析系统及方法
CN104316855B (zh) 一种hemt器件结温的测试方法
CN102103184B (zh) 一种提取晶体管非线性热阻的方法
CN105046066A (zh) AlGaN/GaN HETM小信号模型及其提参方法
CN103033728B (zh) 经时击穿矩阵测试电路及测试方法
CN104376161A (zh) 一种建立AlGaN/GaN HEMT器件直流模型的方法
CN106646180B (zh) 一种wat阈值电压测试方法及系统
CN105468828A (zh) 一种ⅲ-ⅴ族hemt表面势基集约型模型的建模方法
CN105787210A (zh) 氮化镓高电子迁移率晶体管小信号电路模型参数提取方法
CN106199366A (zh) 一种功率mos器件在线测温的方法
CN103954899A (zh) 一种实时测量二极管瞬态温升的方法
Gibiino et al. A GaN HEMT global large-signal model including charge trapping for multibias operation
CN111398764B (zh) 功率管电压测试方法、装置和测试电路
CN101459093B (zh) 验证非对称高压场效应管漂移区电阻的方法
CN107861042A (zh) 一种用于宽禁带半导体功率器件的测试方法
Hefyene et al. Bias-dependent drift resistance modeling for accurate DC and AC simulation of asymmetric HV-MOSFET
CN113505504B (zh) 一种提取GaN HEMT器件热源模型的方法
Luo et al. Reliable GaN HEMT modeling based on Chalmers model and pulsed S-parameter measurements
CN104716065A (zh) 金属氧化物半导体场效应晶体管电容-电压特性修正方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 610000 Sichuan, Shuangliu County, Southwest Airport Economic Development Zone, the Internet of things industry park

Applicant after: CHENGDU HIWAFER TECHNOLOGY CO., LTD.

Address before: 610000 Sichuan, Shuangliu County, Southwest Airport Economic Development Zone, the Internet of things industry park

Applicant before: CHENGDU GASTONE TECHNOLOGY CO., LTD.

COR Change of bibliographic data
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20160302