改进的本征态展开方法求解含时薛定谔方程

本文改进了原有的本征态展开方法.通过从能量E到q=( )2E的表象变换,不仅准确地计算了在此方法中起重要作用的低能电子布居,而且大幅度地减少了计算时间.利用这个高效的方法,我们计算了在强激光作用下一个模型原子的高次谐波发射谱.     

量子点中Stark效应及微观电极化率的理论研究

以半导体材料InSb为研究对象,利用有效质量近似方法,分别研究了在强弱电场下含有两个电子的量子点中的能级移动.首先,利用微扰法计算了库仑相互作用所引起的能量变化ΔE,计算由于弱电场而引起的Stark能级移动ΔE1,并研究了ΔE1与特征长度L的关系;其次,在强电场作用下,用解薛定谔方程的方法求出了能量变化ΔE2并将其与弱电场下的情况进行对比;最后,在此基础上利用非线性光学知识研究了强电场作用下微观电极化率χ是随着L的增加而增大的.     

高剥离类钠离子3s2S—3d2电四极矩E2光谱跃迁的理论计算

 利用全相对论性多组态Dirac-Fock平均能级方法系统地计算了高剥离类钠离子3s²S—3d²D(Z=14～103）电四极矩E2光谱跃迁的能级间隔，跃迁几率和振子强度，计算中考虑了核的有限体积效应，Breit修正和QED修正，所得结果和最近的实验数据及理论计算值进行了比较，计算结果表明：高原子序数的高荷电离子的电四极矩E2光谱跃迁的跃迁几率和中性原子的电偶极E1的相当，在ICF和MCF高温高密度激光等离子体中，电四极矩E2光谱跃迁过程不容被忽视。     

100～5000eV下电子被若干分子散射的总截面

利用可加性规则，使用Hartree-Fock波函数，在100～5000eV内采用由束缚原子概念修正过的复光学势，对电子被约40个分子散射的总截面进行了计算，并将计算结果与实验结果及其它理论计算结果进行了比较。得出如下结论：①对大多数散射体系而言，仅当入射能量E≥100eV时，计算结果即与实验结果符合得较好，少数较复杂的体系，当E≥200～300eV时，也能与实验结果符合得较好。但采用未修正的复光学势进行计算，仅当E≥500～800eV时才与实验结果符合得较好，不少体系只有在E≥1000eV时才能得到较好的结果。这说明本修正是成功的；②在较高的能量下，对于同类型分子构成的散射体系，在同一能量下其总截面与分子内的总电子数成正比，且电子数越多，其线性度也就越好；③在较高能量下总截面与能量间存在QT≈aE^-b的关系。②与③相结合，可推算出一些复杂分子的总截面。     

路径积分的动量表象计算和库仑散射

在合理的近似下，通过在过量表象中计算路径积分，较为方便地求出了库仑散射的相移表达式，同时，利用本方法可容易求出散射振幅和截面。     

高剥离类镁等电子序列的磁偶极M1和电四极E2光谱跃迁理论研究

利用全相对论性多组态Dirac-Fock广义平均能级方法，系统地计算了类镁离子3s3p磁偶极Ml^3P2--^3P1和电四极E2 ^2P2--^3P0(Z=20-103)光谱跃迁的能级间隔、跃迁几率和振子强度。计算中考虑了原子核的有限体积效应，进行了高阶Breit修正和QED修正，所得到的能级间隔和最近的实验数据及理论计算值进行了比较。计算结果表明：高原子序数的高荷电离子的磁偶极矩M1和电四极矩E2跃迁几率和中性原子的电偶极E1的相当，在ICF和MCF高温激光等离子体中，磁偶极矩M1和电四极矩E2跃迁过程不容被忽视。     

石墨Raman活性模E2g的压缩行为及C—C键力常数的压力效应

 本文在0~13 GPa压力范围内测量了石墨Raman活性模E_2g的频率随压力的变化，结果表明在4 GPa附近dν/dp出现了不连续点，该点前的dν/dp大于该点后的dν/dp，这是由于在压力的作用下石墨六角网格发生畸变所致。本工作还导出了高压状态下C—C键力常数的表达式，对E. Fitzer在常压下用Raman频率计算C—C键长的计算公式进行了修正。     

N2介质阻挡放电中有效电场的光谱测定

采用蒙特卡罗方法，模拟了N2非平衡态直流放电等离子体过程。计算了不同E／N条件下电子在氮气中的漂移速度，结果与实验数据符合得很好。模拟了电子激发态粒子数随E／N的变化关系。提出在介质阻挡放电研究和应用中，通过监测两条紫外谱线145和201nm强度之比来监测有效电场的方法。     

LiNbO_3在反常色散区光学常数的测定

我们用170SX傅里叶变换红外光谱仪,测量了LiNbO_3在30cm~(-1)～2800cm~(-1)波段内的反射光谱。应用Kramers-Kronig关系的积分形式,快速而准确地计算出LiNbO_3在同一波段内的各种光学常数的谱图。由此,可直接得到A_1模和E模的横频率ω_T、纵频率ω_L和高频介电常数ε_∞。计算了振子强度。讨论了振子间的相互耦合作用。并认为670cm~(-1)模就是E模的一个基频。     

球形粒子复合介质非线性有效电导交流响应

利用摄动方法，研究了在外加交流电场作用下，含球形杂质的非线性复合介质有效电导的交流响应问题。在一类形如E1sinωt +E3sin3ωt的正弦电场作用下，分别导出了非线性复合介质中杂质粒子及基质所在区域的局部电势。在杂质浓度较低时，进一步得到了复合介质非线性有效电导在基频和三阶谐波下的响应，并讨论了二者之间的关系。     

电子入射双电离He原子的理论研究

完成了高入射能(E05.6keV),低敲出能(E1=E2=10eV)情形下,电子入射双电离He原子五重微分截面(FDCS)的理论计算.计算中使用了三个电子在He+2离子场中的四体末态波函数.计算结果与其他理论结果进行比较发现:所得结果与最新的(e,3e)实验数据较好地一致.     

具有复折射率的多层光波导的改善的解析方法

从理论上研究具有复折射率的多层光波导的关键是求解复数超越本征方程，计算工作量是庞大的，本文提出了改善计算这种光波导中导模折射率的普遍公式。将它与叠代法结合起来，即使对于高损耗或高增益的结构，也能迅速地求得很精确的结果。文中以四层和五层的金属覆盖层波导作为实例说明这一计算方法，并将计算结果与其它方法所得的结果进行了比较。     

分层海洋介质中本征声线的快速计算

本文提供一种在恒定水深的水平分层海洋中快速计算本征声线的方法，该方法引入声线跨度的概念，将本征声线分解为若干声线跨度的组合，并以声线常数为自变量建立声线组合方程。利用牛顿下山法迭代求解该方程得到本征声线，极大地提高了计算效率，给定海深，声剖面，声源距离、深度及接收器深度等参数时，应用本方法可计算出若干条本征声线的出射角、到达接收位置时的入射角和传播时间。若附加声源的辐射频率和海底底质等信息，还可     

Monte Carlo方法模拟低能电子在多元介质中散射——基于平均散射截面方法

应用基于平均散射截面低能电子在多元介质中散射Monte Carlo方法,模拟E≤5eV低能电子在多种多元介质中散射。计算了电子背散射系数,背散射电子能谱、角分布,入射电子、背散射电子在介质中的作用范围、沉积能分布,并与确定散射中心方法的结果比较。两种方法计算结果广泛一致,进一步证明基于平均散射截面方法的有效性和可靠性。入射电子能量较低,介质平均原子序数较大时,计算的背散射电子角分析不服从余弦分布律。     

BaVS3晶格动力学研究

采用基于密度泛函理论的平面波赝势方法,在局域密度近似下采用线性响应的密度泛函微扰理论研究了具有六角结构的BaVS3化合物的品格动力学性质,得到了整个声子谱.计算得到的г点拉曼频率和实验数据进行了比较,其中E22g模、A1g模和实验测量值符合得比较好.对于E32g模,采用线性响应计算的结果与实验值差别较大.对该模应用冻结声子方法研究后认为差异主要是由于E32g模的较强的非谐性引起的.此外V原子在平面内的振动模E12g出现了虚频.虚频的出现预示着六角相的BaVS3结构的不稳定性,从而很好地解释了该材料由六角相到正交相的结构相变.     

非线性系统的任意项精细积分外插多步法及其在混沌数值分析中的应用

对非线性系统提出了高精度的精细积分任意项外插多步法的计算公式.本方法只需增加插值 项数即可提高计算精度，同时不会增加过大的计算量，发展完善了精细积分法.将本方法应 用于混沌方程中，取得了较好的效果.数值计算结果表明，该方法是一种高精度、高效率的 方法，在求解混沌系统上比传统方法有很大的优势. 关键词： 非线性系统 任意项精细积分多步法 混沌系统     

电子阻止本领S_e(E)的z振荡及其计算

本文应用Hartree-Fock-Slater势能代替Thomas-Fermi势能,得出电子阻止本领S_e(E)的归一化关系式;同时,结合有效原子序数概念,得出能量对S_e(E)影响的关系式,利用前一关系式,定量计算了S_e(E)随离子和靶原子的原子序数z变化的振荡关系,计算结果与实验值相当符合,利用后一关系式,定量计算了多种离子-靶体系的S_e(E),与文献[1]的关系式比较,本文结果更好地符合实验值,而且有更广的能量适用范围。     

钨合金的高压本构研究

通过静高压实验、动高压实验及理论计算相结合的方法，确定了钨合金的一种高压本构方程 （Steinberg模型方程）中的各参数.用数值模拟计算方法，采用几种不同的本构模型进行了 计算，并与实测结果进行了比较.结果显示，确定的高压本构方程计算值与实测值符合最好. 关键词： 高压本构方程 Steinberg 模型 数值模拟计算 粒子速度波形     

等效磁荷方法及其在三维磁场计算中的应用

本从普通物理教学角度讨论了等效磁荷概念；以及等效磁荷方法在包含有非线性高磁导介质的三维磁场计算中的应用，给出了用此方法计算电磁除铁器磁场的实例，最后讨论此方法的优缺点。     

冲击载荷下混凝土本构模型构建研究

 在对混凝土动态力学性能和现有本构模型综合分析的基础上,构建了一个新的适用于冲击响应问题数值分析的混凝土本构模型。该本构模型全面考虑了压力、应力第三不变量、变形的硬化和软化、应变率强化以及拉伸损伤等各个影响因素。将其加入LTZ-2D程序,确定了本构模型参数,对混凝土靶板的穿透问题进行了数值验证分析。计算得到的弹体剩余速度同实验结果基本一致,同时得到了混凝土靶板破裂的计算图像。计算结果及其分析表明,所构建的本构模型能够较好地反映冲击载荷作用下混凝土动态响应的主要特性。