Electron-phonon induced superconductivity in highT c superconductors

It is proposed that phonon induced mechanism is capable of explaining the observed range ofT _c for the newly discovered high temperature superconductor within the Allen and Dynes strong coupling theory.     

Comparisons between statistics of low energy collision cascades generated by full molecular dynamics and by its binary collision approximation

Abstract

Collision cascades in Cu, Au and Cu₃Au are generated by full molecular dynamics (MD) and by its binary collision approximation (BCA) with the Marlowe program. Cu and Au primaries have 1 keV initial energy.

The same Molière repulsive potential is used in both models for close encounters. In the MD model, this potential is carefully splined to the pair component of the N-body potential developed by Ackland and Vitek. In the BCA, this N-body interaction is roughly modeled by a constant isotropic 4 eV binding energy of the target atoms to their rest positions.

Time distributions of the number of atoms moving with a total energy higher than a threshold value E _d are compared and discussed. Recoil range distributions during the cascade development are discussed as well. The agreement between MD and BCA is fairly good in all cases for E _d larger than about 3 eV. In the case of smaller E _d-values, the BCA may result in an overestimate of the number of moving atoms in the late development of the cascades. This discrepancy is suggested to originate in the lack of attractive forces between the moving particles and the surrounding atoms in the BCA.     

Bi4Se3薄膜的声子结构及电声子相互作用的第一性原理研究

基于密度泛函微扰理论,运用第一性原理研究两种终结面的Bi₄Se₃薄膜的声子结构和电-声子相互作用.结果表明两种终结面的Bi₄Se₃薄膜体系都是动力学稳定的. Bi₄Se₃薄膜中插入的Bi₂双原子层与Bi₂Se₃五原子层的声子投影态密度并不完全匹配,这会阻碍部分声子的输运,降低热导,从而有利于提高材料的热电性能.另外,两种终结面的Bi₄Se₃薄膜的电-声子耦合系数都不太大（约0.278）,有利于制备基于室温工作的电子学器件.     

Electron-phonon coupling in cuprate and iron-based superconductors revealed by Raman scattering

Electron-phonon coupling （EPC） in cuprate and iron-based superconducting systems, as revealed by Raman scat- tering, is briefly reviewed. We introduce how to extract the coupling information through phonon lineshape. Then we discuss the strength of EPC in different high-temperature superconductor （HTSC） systems and possible factors affecting the strength. A comparative study between Raman phonon theories and experiments allows us to gain insight into some crucial electronic properties, especially superconductivity. Finally, we summarize and compare EPC in the two existing HTSC systems, and discuss what role it may play in the HTSC.     

Study of the antimicrobial peptide indolicidin and mutants in eukaryotic modelled membrane by molecular dynamics simulations

In this work the interaction of the antimicrobial peptide indolicidin (IND) and its mutants CP10A and CP11 with a eukaryotic membrane model was examined by molecular dynamics simulations. The aim was to analyse the behaviour of these antimicrobial peptides when they interact with a eukaryotic modelled membrane, thereby obtaining atomic detailed observations that are not experimentally available. In the simulations, the widely studied dipalmitoylphosphatidylcholine hydrated bilayer was used as a eukaryotic membrane model. In agreement with experimental observations, the peptides IND, CP10A, and CP11 insert into the bilayer differently; the peptides that insert more deeply present the major hemolytic activities. The hydrophobic residues are responsible for the insertion, but some Trp residues of the peptides remain at the bilayer/water interface because they interact with the bilayer choline groups by cation-π interactions that should be important for recognition of eukaryotic membrane by the three studied peptides.     

Molecular dynamics simulations of displacement cascades in Fe-10%Cr systems

Molecular dynamics simulations of the displacement cascades in Fe-10%Cr systems are used to simulate the primary knocked-on atom events of the irradiation damage at temperatures 300,600,and 750 K with primary knockedon atom energies between 1 and 15 keV.The results indicate that the vacancies produced by the cascade are all in the central region of the displacement cascade.During the cascade,all recoil Fe and Cr atoms combine with each other to form Fe-Cr or Fe-Fe interstitial dumbbells as well as interstitial clusters.The number and the size of interstitial clusters increase with the energy of the primary knocked-on atom and the temperature.A few large clusters consist of a large number of Fe interstitials with a few Cr atoms,the rest are Fe-Cr clusters with small and medium sizes.The interstitial dumbbells of Fe-Fe and Fe-Cr are in the 111 and 110 series directions,respectively.     

Molecular dynamics simulations of displacement cascades in Feben 10%Cr systems

Molecular dynamics simulations of the displacement cascades in Fe-10%Cr systems are used to simulate the primary knocked-on atom events of the irradiation damage at temperatures 300, 600, and 750 K with primary knocked-on atom energies between 1 and 15 keV. The results indicate that the vacancies produced by the cascade are all in the central region of the displacement cascade. During the cascade, all recoil Fe and Cr atoms combine with each other to form Fe-Cr or Fe-Fe interstitial dumbbells as well as interstitial clusters. The number and the size of interstitial clusters increase with the energy of the primary knocked-on atom and the temperature. A few large clusters consist of a large number of Fe interstitials with a few Cr atoms, the rest are Fe-Cr clusters with small and medium sizes. The interstitial dumbbells of Fe-Fe and Fe-Cr are in the lan111ran and lan110ran series directions, respectively.     

Energy flow in thermostatted non-equilibrium molecular dynamics simulations

Energy transfers between internal kinetic and potential energy reservoirs in a simple liquid are studied by setting the temperature of one energy reservoir to a different value from that of the others and computing the resultant energy flows. In the first set of simulations, the x-directional kinetic temperature was artificially raised above the other five, and in the second, the x-directional configurational temperature was artificially raised above the other five. In both cases, external energy flows balanced, but unexpected energy flows between different directional components of the potential energy were observed. Additional simulations showed that these energy flows occurred regardless of the arrangement of thermostats imposed on the six degrees of freedom and the addition of shear. Heat flow between degrees of freedom that were ostensibly at the same temperature was anomalously observed. It was concluded that a different breakdown of the contributions to the configurational energy that is consistent with the definition of the directional configurational temperatures is required.     

Predicting XAFS scattering path cumulants and XAFS spectra for metals (Cu,Ni, Fe,Ti, Au) using molecular dynamics simulations

The ability of molecular dynamics (MD) simulations to support the analysis of X‐ray absorption fine‐structure (XAFS) data for metals is evaluated. The low‐order cumulants (ΔR, σ², C₃) for XAFS scattering paths are calculated for the metals Cu, Ni, Fe, Ti and Au at 300 K using 28 interatomic potentials of the embedded‐atom method type. The MD cumulant predictions were evaluated within a cumulant expansion XAFS fitting model, using global (path‐independent) scaling factors. Direct simulations of the corresponding XAFS spectra, χ(R), are also performed using MD configurational data in combination with the FEFFab initio code. The cumulant scaling parameters compensate for differences between the real and effective scattering path distributions, and for any errors that might exist in the MD predictions and in the experimental data. The fitted value of ΔR is susceptible to experimental errors and inadvertent lattice thermal expansion in the simulation crystallites. The unadjusted predictions of σ² vary in accuracy, but do not show a consistent bias for any metal except Au, for which all potentials overestimate σ². The unadjusted C₃ predictions produced by different potentials display only order‐of‐magnitude consistency. The accuracy of direct simulations of χ(R) for a given metal varies among the different potentials. For each of the metals Cu, Ni, Fe and Ti, one or more of the tested potentials was found to provide a reasonable simulation of χ(R). However, none of the potentials tested for Au was sufficiently accurate for this purpose.     

正常金属-分子量子点-超导耦合系统的介观输运研究

运用非平衡格林函数理论、正则变换以及BCS平均场理论研究了正常金属-分子量子点-超导耦合系统（N-MQD-S）的介观输运,得到了系统的电流公式并选择适当的参数进行了数值计算。数值计算结果表明：电声子耦合强度 与线宽函数 对系统的输运行为有较大影响。当 增大时,由声子辅助隧穿所产生的共振峰将高于分子量子点自身能级产生的共振峰;当线宽函数 增大时,在Andreev反射共振峰的两侧将出现新的边峰。     

正常金属-分子量子点-超导耦合系统的介观输运研究

运用非平衡格林函数理论、正则变换以及BCS平均场理论研究了正常金属-分子量子点-超导耦合系统（N-MQD-S）的介观输运,得到了系统的电流公式并选择适当的参数进行了数值计算。数值计算结果表明：电声子耦合强度 与线宽函数 对系统的输运行为有较大影响。当 增大时,由声子辅助隧穿所产生的共振峰将高于分子量子点自身能级产生的共振峰;当线宽函数 增大时,在Andreev反射共振峰的两侧将出现新的边峰。     

刚性多原子分子的正则系综分子动力学算法

通过将Nose所定义的扩展系统哈密顿量推广到刚性多原子分子系统中，严格地推导了刚性多原子分子在正则系综下的运动方程.在此基础上，证明了依靠所给出的运动方程能得到平衡态的正则分布，并给出了该方程所对应的可逆而守恒的积分格式. 关键词： 刚性多原子分子 正则系综 分子动力学算法     

双壁碳纳米管电浸润现象的分子动力学模拟

在单壁碳纳米管电浸润现象原子模拟的基础上,对双壁碳纳米管的电浸润现象进行了计算机模拟.运用经典分子动力学方法结合一个宏观的电毛细管模型模拟了双壁碳纳米管在水银中的电浸润过程,对不同内管尺寸情况下的浸润现象作了研究和比较.计算结果表明双壁碳管和单壁碳管的电浸润过程存在很大的不同,双壁碳管的内管在电浸润过程中起到重要的作用：当改变双壁碳管中内管的尺寸时,浸润现象会产生很大的改变. 关键词： 双壁碳纳米管 电浸润 分子动力学     

Atomic collision cascades on void evolution in vanadium

Molecular dynamics simulations were performed to study void evolution subject to unidirectional self-bombardment and radiation-induced variation of mechanical properties in single crystalline vanadium. 3D simulation cells of perfect body-centered cubic (BCC) vanadium, as well as those with one, two, four, and six voids, were investigated. For the no void case, the maximum number of defects, maximum volumetric swelling, and the number of defects left in bulk after a sufficiently long recovery period increased with higher primary recoil energy. For the cases containing voids, a primary recoil energy was carefully assigned to an atom so as to initiate a dense collision spike in the voids center, where some self-interstitial atoms gained kinetic energy by secondary replacement collision sequence traveling along the ? 111? direction. It is found that the larger or the greater the number of voids contained initially in the box, the larger the normalized void volume, and the smaller the volumetric swelling after sufficient recovery of systems. In the single void case, the void became elongated along the bombarding direction; in the multiple void cases, the voids coalesced only when the intervoid ligament distance was short. After sufficient relaxation of the irradiated specimen, a hydrostatic tension was exerted on the box, where the voids were treated as dislocation sources. It is shown that with higher primary recoil energy, the yield stress dropped in cases with smaller or fewer voids but rose in those with larger or greater number of voids. This radiation-induced softening to hardening transition with increasing dislocation density can be attributed to the combined effects of the defect-induced dislocation nucleation and the resistance of defects to dislocation motion. Moreover, as the primary recoil energy increased, the ductility of vanadium in the no void case decreased, but was only slightly changed in the cases containing void.     

Hard-sphere heat conductivity via nonequilibrium molecular dynamics

We use an Evans-Gillan driving forceF ^d, together with isokinetic and isoenergetic constraint forcesF ^c, to drive steady heat currents in periodic systems of 4 and 32 hard spheres. The additional driving and constraint forces produce curved trajectories as well as additional streaming and collisional contributions to the momentum and energy fluxes. Here we develop an analytic treatment of the collisions so that the simulation becomes approximately ten times faster than our previous numerical treatment. At low field strengths, for less than 0.4, where is the hard-sphere diameter, the 32-sphere conductivity is consistent with Alder, Gass, and Wainwright's 108-sphere value. At higher field strengths the conductivity varies roughly as ^1/2, in parallel with the logarithmic dependence found previously for three hard disks.     

氦热力学性质的分子动力学研究

 利用经典分子动力学和第一性原理分子动力学,研究了氦在高压下的熔化曲线、状态方程和非金属-金属转变。得到了氦在温度小于4.5 eV、 密度0.3~5.0 g/cm³范围内的状态方程,并把氦的熔化曲线的压强范围拓展到了50 GPa。氦的能隙宽度曲线表明,温度大大降低了氦的金属化密度。     

非绝热分子动力学的量子路径模拟

基于近期发展的经典-量子混合模拟非绝热分子动力学的量子路径方案,本文对5个典型势能面模型进行了模拟,包括单交叉模型、双交叉模型、拓展耦合模型、哑铃模型以及双弓模型.由于难以在严格意义上得到退相干速率,数值模拟中,我们比较了三个不同的退相干速率公式,包括冻结高斯波包近似退相干速率、能量分辨速率以及力分辨速率.在模拟过程中,我们恰当地处理了势能面跳跃时的能量守恒和力的反向问题.通过与全量子动力学模拟的精确结果进行对比发现,对于结构较简单的势能面模型,三种退相干速率都能得到较好的结果;然而对于较复杂的势能面模型,由于复杂量子干涉的原因,与其他混合经典-量子动力学方案类似,量子路径方案仍然难以得到较准确的结果.如何发展更加有效的混合经典-量子模拟方案,是未来研究的重要课题.     

非晶态石英的变电荷分子动力学模拟

基于迭代变电荷方法,用分子动力学模拟了非晶态石英的结构与振动特征.首先利用熔化-猝火方法得到了非晶石英的平衡结构.在此基础上获得了体系不同原子对之间的对关联函数、键角分布函数和振动频谱等,结果与实验数据均符合较好.变电荷方法的计算结果表明,非晶石英体系内粒子的电荷与石英晶体内粒子电荷显著不同,并且出现了较大的涨落. 关键词： 分子动力学 变电荷 非晶石英     

Isospin- and momentum-dependent interaction in isospin-dependent quantum molecular dynamics

An isospin-dependent quantum molecular dynamical model (IQMD) is developed, with the isospin degree of freedom in the momentum-dependent interaction(MDI) included in IQMD, to obtain an isospin- and momentum-dependent interaction (IMDI) in IQMD. We investigate the effect of IMDI on the isospin fractionation ratio and its dynamical mechanism in the intermediate energy heavy ion collisions. It is found that the IMDI induces the significant reductions in the isospin fractionation ratio for all of beam energies, impact parameters, neutron--proton ratios and mass number of colliding systems. However, the strong dependence of isospin fractionation ratio on the symmetrical potential is preserved, with the isospin degree of freedom included in the MDI, i.e. the isospin fractionation ratio is still a good probe for extracting the information about the equation of state of isospin asymmetrical nuclear matter.     

Molecular dynamics simulation of self-diffusion coefficients for liquid metals

The temperature-dependent coefficients of self-diffusion for liquid metals are simulated by molecular dynamics meth ods based on the embedded-atom-method （EAM） potential function. The simulated results show that a good inverse linear relation exists between the natural logarithm of self-diffusion coefficients and temperature, though the results in the litera ture vary somewhat, due to the employment of different potential functions. The estimated activation energy of liquid metals obtained by fitting the Arrhenius formula is close to the experimental data. The temperature-dependent shear-viscosities obtained from the Stokes-Einstein relation in conjunction with the results of molecular dynamics simulation are generally consistent with other values in the literature.