期刊文献+
共找到18篇文章
< 1 >
每页显示 20 50 100
钙钛矿超晶格SrTiO_(3)/BaTiO_(3)的挠曲电效应
1
作者 陈许敏 叶盼 +2 位作者 王继光 霍德璇 曹东兴 《物理学报》 SCIE EI CAS CSCD 北大核心 2022年第20期207-215,共9页
挠曲电效应是应变梯度与电极化之间的机电耦合,存在尺寸效应,在许多微纳米尺寸结构中起着非常重要的作用.本文采用密度泛函理论,对交错层SrTiO_(3)/BaTiO_(3)超晶格进行了系统的挠曲电效应研究,通过探究超晶格在给定的应变梯度下的力电... 挠曲电效应是应变梯度与电极化之间的机电耦合,存在尺寸效应,在许多微纳米尺寸结构中起着非常重要的作用.本文采用密度泛函理论,对交错层SrTiO_(3)/BaTiO_(3)超晶格进行了系统的挠曲电效应研究,通过探究超晶格在给定的应变梯度下的力电响应,独立得到其纵向挠曲电系数、横向挠曲电系数和剪切挠曲电系数.结果表明:超晶格的横向、剪切挠曲电系数和纵向挠曲电系数较其组分材料有不同的变化,其中超晶格的横向和剪切挠曲电系数分量分别较块体BaTiO_(3)提升约6倍,较块体SrTiO_(3)提升4.2倍和1.3倍;纵向挠曲电系数较其组成材料基本不变;这种挠曲电系数分量不同程度提升的综合效果能够使超晶格SrTiO_(3)/BaTiO_(3)较其单一组分材料的挠曲电效应产生数倍提升.本文对寻找性能优异的复合挠曲电材料具有一定理论指导意义. 展开更多
关键词 挠曲电效应 第一性原理 钙钛矿超晶胞
下载PDF
第一性原理研究ZnS掺V的光学性质和电子结构 被引量:20
2
作者 何开华 余飞 +2 位作者 姬广富 颜其礼 郑澍奎 《高压物理学报》 EI CAS CSCD 北大核心 2006年第1期56-60,共5页
运用密度泛函平面波赝势方法(PWP)和广义梯度近似(GGA),对替代式掺杂钒(V)的闪锌矿(ZnS)的超晶胞电子结构进行了计算。研究了ZnS掺杂的光学性质及其电子结构,通过分析发现,光吸收的计算结果与运用配位场理论得到的计算结果以及实验数据... 运用密度泛函平面波赝势方法(PWP)和广义梯度近似(GGA),对替代式掺杂钒(V)的闪锌矿(ZnS)的超晶胞电子结构进行了计算。研究了ZnS掺杂的光学性质及其电子结构,通过分析发现,光吸收的计算结果与运用配位场理论得到的计算结果以及实验数据符合较好;同时还对引入杂质前后电子结构的异同以及价键的一些性质进行了对比分析:因为杂质V的引入,态密度出现了几个新的峰,并且本体能带向低能方向偏移了大约2.5 eV,V所带正电荷为0.28,比任何一类Zn原子都要小,S—V键的共价性最强、键长最短。 展开更多
关键词 超晶胞 杂质 光学性质
下载PDF
氢氧化铝晶体生长习性的Morphology分析 被引量:4
3
作者 吴争平 尹周澜 +1 位作者 陈启元 李洁 《中南大学学报(自然科学版)》 EI CAS CSCD 北大核心 2009年第4期897-903,共7页
对氢氧化铝晶体生长习性进行Morphology理论研究,明确氢氧化铝晶粒长大后的晶粒形态特征及氢氧化铝的理论晶体生长习性。首先,用BFDH法说明了生长过程中最重要的面;其次,采用晶体平衡形态(Equilibrium morphology)法,用Morphology程序... 对氢氧化铝晶体生长习性进行Morphology理论研究,明确氢氧化铝晶粒长大后的晶粒形态特征及氢氧化铝的理论晶体生长习性。首先,用BFDH法说明了生长过程中最重要的面;其次,采用晶体平衡形态(Equilibrium morphology)法,用Morphology程序对氢氧化铝3×4×2超晶胞及其(001),(100),(010),(011),(110),(101)和(112)面真空slab模型进行平衡形态计算,计算结果与实验结果吻合,很好地预测了各显露面族(001),(100),(010),(011),(110),(101)和(112)面族。理论生长习性计算结果表明,氢氧化铝3×4×2超晶胞及其各表面真空slab模型的生长习性差异显著;氢氧化铝(gibbsite)超晶胞趋向于生长为准六棱柱或厚的准六角板状晶体;(001),(100)和(010)面slab基本上是生长为楔板状晶体;(110),(101)和(112)slab的晶体形态则趋向于生长为较细长的棒状晶体;(011)slab的生长形态介于楔板状和棒状之间,很可能生长为较薄的片状晶体。 展开更多
关键词 氢氧化铝 超晶胞 生长习性 晶体形态
下载PDF
过渡金属掺杂金红石相TiO_2能带结构的第一性原理计算 被引量:16
4
作者 陈琦丽 唐超群 《材料科学与工程学报》 CAS CSCD 北大核心 2006年第4期514-516,共3页
本文采用第一性原理能带计算方法和超晶胞模型计算金红石相TiO2掺杂过渡金属元素的电子结构。计算结果表明,Zn掺杂对TiO2的带隙宽度影响不明显,V、Cr、Mn、Fe、Co、Ni、Cu的掺杂都有可能使TiO2吸收带出现红移现象或产生在可见光区的吸收... 本文采用第一性原理能带计算方法和超晶胞模型计算金红石相TiO2掺杂过渡金属元素的电子结构。计算结果表明,Zn掺杂对TiO2的带隙宽度影响不明显,V、Cr、Mn、Fe、Co、Ni、Cu的掺杂都有可能使TiO2吸收带出现红移现象或产生在可见光区的吸收,其中杂质原子的t2g态起了重要作用。 展开更多
关键词 TIO2 密度泛函理论 超晶胞 过渡金属掺杂
下载PDF
ZnS掺Mn^(2+)电子结构和光学性质的第一性原理计算 被引量:10
5
作者 万淼 郑广 +2 位作者 何开华 洪汉烈 陈琦丽 《原子与分子物理学报》 CAS CSCD 北大核心 2009年第1期157-162,共6页
在广义梯度近似下,利用平面波赝势对ZnS(闪锌矿,F-43m)和Mn^(2+)掺杂的ZnS超晶胞的电子态密度、原子间电子云重叠布局数和光学性质等进行了自恰计算.自旋极化的计算结果显示,掺入Mn^(2+)离子后态密度整体向低能方向移动,在禁带中出现了... 在广义梯度近似下,利用平面波赝势对ZnS(闪锌矿,F-43m)和Mn^(2+)掺杂的ZnS超晶胞的电子态密度、原子间电子云重叠布局数和光学性质等进行了自恰计算.自旋极化的计算结果显示,掺入Mn^(2+)离子后态密度整体向低能方向移动,在禁带中出现了由Mn 3d、Zn3d与S 3p组成的新态.电子云重叠布局分析了掺杂前后价键性质的变化,解释了磁性离子Mn^(2+)导致Zn和S出现不对称的自旋向上和自旋向下态的机理.复介电系数谱图向高能方向迁移,并且在0到2.7 eV范围内出现一个新的尖峰,利用晶体场理论和态密度,对改介电峰进行了指认,为研究此类材料光学性质提供了一定理论依据. 展开更多
关键词 超晶胞 掺杂 电子结构 光学性质
下载PDF
Sn_(1-x)N_xO_2材料光电性质的研究
6
作者 张国莲 逯瑶 +2 位作者 冯现徉 张昌文 王培吉 《人工晶体学报》 EI CAS CSCD 北大核心 2011年第6期1623-1627,共5页
基于密度泛函理论的第一性原理,采用广义梯度近似(GGA)下的线性缀加平面波(FP-LAPW)方法,应用WIEN2K软件计算了超晶胞结构Sn1-xNxO2材料的总态密度、能带结构和光学折射率及介电函数虚部。计算结果表明掺杂后费米能级向低能方向移动,随... 基于密度泛函理论的第一性原理,采用广义梯度近似(GGA)下的线性缀加平面波(FP-LAPW)方法,应用WIEN2K软件计算了超晶胞结构Sn1-xNxO2材料的总态密度、能带结构和光学折射率及介电函数虚部。计算结果表明掺杂后费米能级向低能方向移动,随着掺杂量的增加,Sn1-xNxO2材料的价带和导带的分裂程度增强,禁带宽度逐渐减小,并且在1.35~2.50 eV的能量范围上形成了杂质带,其主要来源是N的2p态上的电子。分析Sn1-xNxO2材料的能带结构可知掺杂前后均是直接跃迁半导体,掺杂后其介电函数谱和折射率也与带隙相对应地发生红移,介电谱的跃迁峰与电子从价带到导带的跃迁有关,从理论上指出光学性质与电子结构之间的内在关系。 展开更多
关键词 超晶胞SnO2 光电性质 第一性原理
下载PDF
Sn_(1-x)N_xO_2材料的电子结构和光学性质
7
作者 张国莲 王培吉 张昌文 《济南大学学报(自然科学版)》 CAS 北大核心 2011年第4期396-399,共4页
采用第一性原理,以W IEN2K软件为平台对Sn1-xNxO2超晶胞体系的态密度(DOS)、能带结构、介电函数和吸收系数进行模拟计算,从理论上指出光学特性与电子结构之间的内在联系。分析结果表明:掺入杂质后体系带隙减小了0.59eV,费米能级向低能... 采用第一性原理,以W IEN2K软件为平台对Sn1-xNxO2超晶胞体系的态密度(DOS)、能带结构、介电函数和吸收系数进行模拟计算,从理论上指出光学特性与电子结构之间的内在联系。分析结果表明:掺入杂质后体系带隙减小了0.59eV,费米能级向低能方向移动进入价带,并且由半导体变为半金属材料,N的掺入有助于p型SnO2的实现;掺杂前后体系为直接跃迁半导体,其介电函数谱和吸收谱与带隙相对应均发生了红移,并且光学吸收边变宽,增大了光学响应,Sn1-xNxO2材料可广泛应用于红外发光器件。 展开更多
关键词 超晶胞SnO2 光电性质 第一性原理
下载PDF
磷掺杂金红石TiO2电子结构的第一性原理研究 被引量:4
8
作者 宋桥台 杨苹 +2 位作者 王锦标 鲁雄 黄楠 《材料科学与工程学报》 CAS CSCD 北大核心 2008年第3期435-437,440,共4页
TiO2作为一种光催化材料,具有广阔的应用领域。近几年,利用掺杂改性来改善TiO2光催化性已成为当今研究的热点。本文采用第一性原理方法和超晶胞模型计算了磷元素掺杂金红石相TiO2的电子结构,并分析了三种不同掺杂位置的电子结构对TiO2... TiO2作为一种光催化材料,具有广阔的应用领域。近几年,利用掺杂改性来改善TiO2光催化性已成为当今研究的热点。本文采用第一性原理方法和超晶胞模型计算了磷元素掺杂金红石相TiO2的电子结构,并分析了三种不同掺杂位置的电子结构对TiO2光催化性能的影响。计算结果表明,P置换O和P在晶格间隙位置的掺杂对TiO2的带隙影响比较明显,在带隙中均形成了中间掺杂能级,都有可能使TiO2吸收带出现红移现象或产生在可见光区的吸收,但P在晶格间隙位置的掺杂所形成的中间掺杂能级比较平稳,有可能改善TiO2光催化性能,其中磷原子的3p态起了重要作用。 展开更多
关键词 密度泛函理论 金红石相TiO2 超晶胞 光催化
下载PDF
二维复式晶格完全带隙光子晶体缺陷频率的分析 被引量:2
9
作者 年秀芝 陈鹤鸣 《光学与光电技术》 2009年第4期23-26,共4页
提出了一种具有完全带隙的二维复式晶格光子晶体,该晶体是在二维正方形格子中,旋转截面为正方形的柱子,同时在每个原胞中心引入圆形截面的柱子而形成的,并在其中引入点缺陷。运用平面波展开法并结合超晶胞理论分析此缺陷态光子晶体的频... 提出了一种具有完全带隙的二维复式晶格光子晶体,该晶体是在二维正方形格子中,旋转截面为正方形的柱子,同时在每个原胞中心引入圆形截面的柱子而形成的,并在其中引入点缺陷。运用平面波展开法并结合超晶胞理论分析此缺陷态光子晶体的频率特性。仿真结果表明,通过调节缺陷的尺寸、角度等结构参量,可以改变缺陷态频率的位置,使TE和TM模缺陷态频率一致,TE模和TM模同时谐振,处于缺陷态频率的入射光就能够完全耦合进点缺陷,具有较高的耦合效率。这种结构的复式晶格完全带隙光子晶体为制作完全带隙光子晶体谐振腔提供了理论基础。 展开更多
关键词 光子晶体 复式晶格结构 完全带隙 点缺陷 平面波展开法 超晶胞
下载PDF
钽掺杂金红石TiO2电子结构的第一性原理研究
10
作者 宋桥台 杨苹 +2 位作者 王锦标 鲁雄 黄楠 《功能材料》 EI CAS CSCD 北大核心 2007年第A10期4116-4118,共3页
采用第一性原理中的CASTEP模块计算了金红石相TiO2掺杂过渡金属钽元素的能带结构以及态密度。计算结果表明,Ta原子替换Ti原子后,Ta的5d轨道对TiO2的导带影响较明显,且使TiO2的带隙宽度减小,这样有可能使金红石TiO2吸收带出现红移现... 采用第一性原理中的CASTEP模块计算了金红石相TiO2掺杂过渡金属钽元素的能带结构以及态密度。计算结果表明,Ta原子替换Ti原子后,Ta的5d轨道对TiO2的导带影响较明显,且使TiO2的带隙宽度减小,这样有可能使金红石TiO2吸收带出现红移现象或产生在可见光区的吸收,其中Ta原子的t2g态起了重要作用。 展开更多
关键词 金红石相TiO2 密度泛函理论 超晶胞 钽掺杂
下载PDF
空位缺陷对单层硅烯薄膜力学性能的影响 被引量:3
11
作者 郑超 徐建刚 +1 位作者 张云光 李桂霞 《原子与分子物理学报》 CAS 北大核心 2018年第2期321-327,共7页
采用基于分子动力学理论的Forcite模拟软件包对含不同浓度的单、双空位缺陷硅烯薄膜的超晶胞体系进行优化,并对其力学性能进行了计算和分析.结果表明:随着空位缺陷浓度的增加,硅烯薄膜的拉梅常数、泊松比、体弹模量和剪切模量呈线性递... 采用基于分子动力学理论的Forcite模拟软件包对含不同浓度的单、双空位缺陷硅烯薄膜的超晶胞体系进行优化,并对其力学性能进行了计算和分析.结果表明:随着空位缺陷浓度的增加,硅烯薄膜的拉梅常数、泊松比、体弹模量和剪切模量呈线性递减趋势,而由于空位缺陷附近键长的缩减导致硅烯薄膜"硬化"与空位缺陷浓度的增加导致硅烯晶格中硅原子密度降低,两种体制的竞争使得硅烯杨氏模量表现出先升高在降低的趋势. 展开更多
关键词 空位缺陷 超晶胞体系 分子动力学 力学性质
下载PDF
Fe掺杂锐钛矿型TiO_2光催化剂第一性原理研究
12
作者 赵栩莹 张先兵 +2 位作者 吴江 赵镇 彭小龙 《上海电力学院学报》 CAS 2016年第6期561-564,共4页
采用了基于密度泛函理论的分子模拟软件Materials Studio的Dmol^3计算模块,对铁(Fe)掺杂锐钛矿型TiO_2光催化剂的超晶胞体系进行了几何结构优化.模拟计算过程中,使用广义梯度近似的方法处理交换关联能,得到了TiO_2在掺杂Fe前后的能带结... 采用了基于密度泛函理论的分子模拟软件Materials Studio的Dmol^3计算模块,对铁(Fe)掺杂锐钛矿型TiO_2光催化剂的超晶胞体系进行了几何结构优化.模拟计算过程中,使用广义梯度近似的方法处理交换关联能,得到了TiO_2在掺杂Fe前后的能带结构和电子态密度的变化状况,为光催化剂的选择与制备提供了理论依据. 展开更多
关键词 光催化剂 第一性原理 超晶胞
下载PDF
线缺陷对硅烯力学性能的影响
13
作者 李丽萍 李朔 +2 位作者 张聪颖 张云光 李桂霞 《原子与分子物理学报》 CAS 北大核心 2021年第2期61-66,共6页
硅烯是一种类似石墨烯的二维硅,由于其与石墨烯类似的电子特性以及与现有半导体技术的兼容性而引起了人们的广泛关注.采用基于分子动力学理论的Forcite模拟软件包对含不同条数、不同位置的线缺陷硅烯的超晶胞体系进行优化,并计算和分析... 硅烯是一种类似石墨烯的二维硅,由于其与石墨烯类似的电子特性以及与现有半导体技术的兼容性而引起了人们的广泛关注.采用基于分子动力学理论的Forcite模拟软件包对含不同条数、不同位置的线缺陷硅烯的超晶胞体系进行优化,并计算和分析了其力学性能.结果表明:随着线缺陷位置的变化,硅烯的拉梅常数、泊松比、体积模量、杨氏模量和剪切模量变化趋势基本吻合,说明硅烯的力学性能受线缺陷位置影响不敏感;伴随着线缺陷条数的增加,五环和八环密度升高,硅烯力学性能拉梅常数、泊松比、体积模量和剪切模量变化幅度小,影响不大,而扶手椅方向的杨氏模量线性增加,锯齿型方向的杨氏模量线性减小. 展开更多
关键词 线缺陷 超晶胞体系 分子动力学理论 力学性质
下载PDF
α-Al_2O_3(0001)基片表面结构与能量研究 被引量:17
14
作者 杨春 李言荣 +4 位作者 薛卫东 陶佰万 刘兴钊 张鹰 黄玮 《物理学报》 SCIE EI CAS CSCD 北大核心 2003年第9期2268-2273,共6页
对α Al2 O3( 0 0 0 1)晶体表层三种不同终止原子结构的计算模型 ,在三维周期边界条件下的κ空间中 ,采用超软赝势平面波函数描述多电子体系 .应用基于密度泛函理论的局域密度近似 ,计算了不同表层结构的体系能量 ,表明最表层终止原子... 对α Al2 O3( 0 0 0 1)晶体表层三种不同终止原子结构的计算模型 ,在三维周期边界条件下的κ空间中 ,采用超软赝势平面波函数描述多电子体系 .应用基于密度泛函理论的局域密度近似 ,计算了不同表层结构的体系能量 ,表明最表层终止原子为单层Al的表面结构最稳定 .对由 10个原子组成的菱形原胞进行了结构优化 ,得到晶胞参数值(a0 =0 4 8178nm)与实验报道值误差小于 1 3% .进一步计算了超晶胞 ( 2× 2 )表面弛豫 ,弛豫后原第 2层O原子层成为最表层 ;对不同表层O ,Al原子最外层电子进行了布居分析 ,表面电子有更大的概率被定域在O原子的周围 。 展开更多
关键词 α-Al2O3(0001)基片 表面结构 能量 电子表面态 电子功能薄膜材料 超晶胞 表面弛豫 软赝势
原文传递
氧气在Al(001)面吸附的第一原理研究 被引量:11
15
作者 张芳英 朱圣龙 滕英元 《科学通报》 EI CAS CSCD 北大核心 2004年第16期1687-1690,共4页
采用基于广义梯度近似(GGA)交换关联近似的超软(ultrasoft)赝势和具有三维周期性边界条件的超晶胞模型, 用第一原理计算方法, 计算并分析了氧分子在Al(001)面吸附的价键结构和局域电子结构. 超晶胞表面为2×2原胞, 共有14层原子, ... 采用基于广义梯度近似(GGA)交换关联近似的超软(ultrasoft)赝势和具有三维周期性边界条件的超晶胞模型, 用第一原理计算方法, 计算并分析了氧分子在Al(001)面吸附的价键结构和局域电子结构. 超晶胞表面为2×2原胞, 共有14层原子, 其中铝原子9层, 真空5层. 氧分子层与表面铝原子层的距离为一倍面间距, 选用了几种不同的初始吸附形态. 结果表明, 氧气在Al(001)面时, 分子键平行于铝表面时容易被吸附, 分子键垂直于铝表面时不容易被吸附. 吸附过程中价键的分析表明, 氧分子在O2/Al(001)界面的吸附过程有两种形式:2222O(O)2O2O---和222O(O)OO--?, 吸附过程与氧气在Al(001)面初始吸附形态有关. 展开更多
关键词 第一原理 超晶胞 氧分子 吸附 广义梯度近似 Al(001)面 分子键 腐蚀
原文传递
合金元素对β-Ti25Nb合金的结构稳定性与弹性性质的影响 被引量:3
16
作者 万小军 吴昌义 林建国 《金属热处理》 CAS CSCD 北大核心 2011年第4期16-20,共5页
采用密度泛函理论的缀加平面波加局域轨道的第一性原理和超晶胞方法对β结构的Ti75Nb25以及Ti68.75Nb25X6.25(X为Mo,Sn,Ta和Zr)合金的能量、电子结构以及弹性性质进行了理论计算。通过对比这几种合金的结构稳定性和弹性性质,发现Mo和Ta... 采用密度泛函理论的缀加平面波加局域轨道的第一性原理和超晶胞方法对β结构的Ti75Nb25以及Ti68.75Nb25X6.25(X为Mo,Sn,Ta和Zr)合金的能量、电子结构以及弹性性质进行了理论计算。通过对比这几种合金的结构稳定性和弹性性质,发现Mo和Ta元素既能起到稳定β-Ti25Nb相作用,又能使其弹性常量增加。而Sn对β结构的稳定性影响较小,Zr却降低了β结构的稳定性,且两者对弹性常数的影响均较小。 展开更多
关键词 Ti-25Nb合金 结构稳定性 弹性性质 第一性原理 超晶胞
原文传递
六角蜂窝结构光子晶体异质结带隙特性研究 被引量:3
17
作者 关春颖 苑立波 《物理学报》 SCIE EI CAS CSCD 北大核心 2006年第3期1244-1247,共4页
利用蜂窝结构光子晶体具有两个范围较大绝对带隙的特性,设计新型六角蜂窝结构光子晶体异质结,采用平面波展开和超晶胞相结合的方法来研究异质结的能带结构特性.给出异质结的结构和相应的能带图,分析异质结界面传导模,研究两侧结构作横... 利用蜂窝结构光子晶体具有两个范围较大绝对带隙的特性,设计新型六角蜂窝结构光子晶体异质结,采用平面波展开和超晶胞相结合的方法来研究异质结的能带结构特性.给出异质结的结构和相应的能带图,分析异质结界面传导模,研究两侧结构作横向拉开和侧向滑移时对传导模的影响,讨论这些结构的实际可行性.计算结果表明,没有任何晶格移动,此种结构异质结在绝对带隙中就有导模存在,两边晶格横向拉开对导模影响较大,而侧向滑移的影响则较小. 展开更多
关键词 光子晶体 异质结 光子带隙 超晶胞
原文传递
Spatial and temporal charac-teristics of VHF radiation source produced by lightning in supercell thunderstorms 被引量:14
18
作者 ZHANGYijun MENGQing +2 位作者 P.R.Krehbiel LIUXinsheng ZHOUXiuji 《Chinese Science Bulletin》 SCIE EI CAS 2004年第6期624-631,共8页
The three-dimensional temporal and spatial characteristics of VHF radiation events produced by light-ning discharges in three supercell thunderstorms have been analyzed based on the data measured by the lightning map-... The three-dimensional temporal and spatial characteristics of VHF radiation events produced by light-ning discharges in three supercell thunderstorms have been analyzed based on the data measured by the lightning map-ping array system with high time and space resolution. The results indicate that lightning hole (lighting free region) with about 5—6 km in diameter or lighting ring (annular lighting free region) is associated with the strong updraft in thunder-storm. The lasting time of lightning holes is either short or long, being about 20 min in a tornado-producing thunder-storm. The lightning holes appear before the occurrence of tornado. The lightning hole is the most obvious during the occurrence of tornado and some self-existent lighting radia-tion events appear at a height of 15—16 km. The lightning channels of inter-cloud (IC) lightning discharge exhibit clockwise rotary structures and do not have clear bilevel structures in the vicinity of the tornado. The lightning holes are corresponding to the strong updraft region. The tempo-ral and spatial distribution of lightning radiation events re-veals the structure of strong updraft in supercell thunder-storms. Positive cloud-to-ground (CG) lightning discharges dominate in these thunderstorms and the peak of positive CG lightning flash rate appears, with the maximum of 6 per minute, after or before the occurrence of tornado. 展开更多
关键词 时间特性 空间特性 VHF辐射源 超晶胞雷暴 龙卷风 闪电
原文传递
上一页 1 下一页 到第
使用帮助 返回顶部