HMX不同晶型热膨胀特性及相变的ReaxFF分子动力学模拟

HMX（octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine）晶体不同晶型对温度的响应对于深入认识含能材料在外加载荷下的感度、稳定性、相变等具有重要意义．采用ReaxFF—lg势函数和等温等压分子动力学（NPT-MD）方法,研究了β-,δ-和α-HMX三种固态晶型在T=303—503K内的晶体结构和分子结构．结果表明,ReaxFF—lg势函数能够合理地描述HMX晶体的热膨胀行为,计算得到的晶体结构和热膨胀系数与实验数据比较符合．三种晶型的线膨胀系数表明,β-HMX晶体具有明显的热膨胀各向异性;δ-HMX晶体C轴的热膨胀特性与a,b轴略有不同;α-HMX晶体三个方向的热膨胀特性基本相同．三种晶型的体膨胀系数为8-HMX〉α-HMX〉β-HMX,表明δ-HMX晶体对温度最敏感,这可能是δ-HMX晶体具有更高感度的原因．在T=443K时,β-HMX晶体的晶胞参数突变,且部分分子构型由“椅式”转变为“船式-椅式”．采用两相热力学模型（2PT）得到的亥姆赫兹自由能表明,β—HMX晶体在T=423—443K时发生相变．β-,δ-和α-HMX晶体分别在T=303—423K,T=443—503K和T=363-423K内稳定存在．     

高压下Fe从bcc到hcp结构相变机理的第一性原理计算

采用基于密度泛函理论的第一性原理方法, 分别研究了压力作用下Fe从体心立方 (bcc,α 相) 结构到六角密排(hcp, ε相) 结构相变的两种不同的相变机理: 相变过程中出现亚稳定的面心立方(fcc) 结构(bcc-fcc-hcp) , 以及相变过程中没有出现亚稳定的fcc结构(bcc-hcp) . 计算结果表明: 静水压力条件下, 相变过程中并不会产生亚稳定的fcc结构, 这与最近的原位XRD实验测量结果相一致. 随着压力的增加, fcc-hcp的相变势垒逐渐增加, 压力趋向于阻止Fe从fcc结构到hcp结构的相变. 计算得到了相变过程中原子磁性和结构的详细信息, 分析表明相变过程中涉及复杂的磁性转变, 并且讨论了原子磁性对结构转变影响的物理机理. 此外, 对分子动力学模拟中产生亚稳定的fcc结构的原因也进行了讨论. 关键词： 相变机理 静水压力 第一性原理 铁     

Structural,elastic and thermodynamic properties under pressure and temperature effects of MgIn2S4 and CdIn2S4

A density functional-based method is used to investigate the structural, elastic and thermodynamic properties of the cubic spinel semiconductors MgIn₂S₄ and CdIn₂S₄ at different pressures and temperatures. Computed ground structural parameters are in good agreement with the available experimental data. Single-crystal elastic parameters are calculated for pressure up to 10 GPa and temperature up to 1200 K. The obtained elastic constants values satisfy the requirement of mechanical stability, indicating that MgIn₂S₄ and CdIn₂S₄ compounds could be stable in the investigated pressure range. Isotropic elastic parameters for ideal polycrystalline MgIn₂S₄ and CdIn₂S₄ aggregates are computed in the framework of the Voigt–Reuss–Hill approximation. Pressure and thermal effects on some macroscopic properties such as lattice constant, volume expansion coefficient and heat capacities are predicted using the quasi-harmonic Debye model in which the lattice vibrations are taken into account.     

Half-metallic magnetism of Co2CrX (X=As, Sb) Heusler compounds: An ab initio study

In this study, we present the electronic, magnetic, and structural properties of two novel half-metallic full-Heusler compounds, Co₂CrAs and Co₂CrSb, in cubic L2₁ geometry. The calculations are based on the density functional theory within plane-wave pseudopotential method and spin-polarized generalized gradient approximation of the exchange-correlation functional. The electronic band structures and density of states of the systems indicate half-metallic behavior with vanishing electronic density of states of minority spins at Fermi level, which yields perfect spin polarization. The calculated magnetic moments of both systems in L2₁ structure are 5.00 μ_B, which are largely localized on the chromium site. The energy gaps in minority spin states are restricted by the 3d-states of cobalt atoms on two different sublattices. The formation enthalpies for both structures are negative indicating stability of these systems against decomposition into stable solid compounds.     

Mechanism of reversible photo-induced magnetization in prussian blue analogues

Photo-induced magnetic phase transitions observed in Co—Fe Prussian blue analogues are discussed from a theoretical point of view. In this material a magnetization is induced by visible light, and the induced magnetization can be completely diminished not only by thermal annealing but also by irradiation of infra red light. After a brief survey of relevant experiments, we present an overview of our theoretical approaches to clarify the underlying mechanism of this interesting phenomenon. The ligand field at the Co site is shown to be the most relevant variable controlling spin states. The presence of Fe vacancies and the resultant ligand substitution by water molecules turns to be an essential prerequisite for the photoinduced magnetization.     

Ab initio calculation of elastic properties of rock-salt and zinc-blend MgS under pressure

The elastic properties of zinc-blend (ZB) and rock-salt (RS) MgS are calculated by ab initio method. The calculation shows that the enthalpy for RS structure and for ZB structure essentially is the same at ambient pressure. The ZB structure becomes unstable above 5 GPa. For these two structures, the pressure dependences of typical elastic properties, i.e. the bulk modulus, the shear modulus, the Young's modulus, the Poisson's ratio, and the anisotropy factor, are presented. The Debye temperature and sound velocity under high pressure have also been calculated. Debye approximation is used to estimate the zero-point vibrational energy.     

An ab initio calculation of 17O and 29Si NMR parameters for SiO2 polymorphs

Ab initio molecular orbital calculations (Hartree–Fock, HF and density functional theories, DFTs) have been carried out for SiO₂ polymorphs coesite, low cristobalite, and -quartz, in order to investigate the reliability of this method for predicting ²⁹Si and ¹⁷O nuclear magnetic resonance (NMR) properties of silicates. Oxygen- and silicon-centered clusters consisting of one (1T) to three tetrahedral (3T) shells (one to four atomic shells), taken from real crystal structure, have been investigated. It is found that for reasonable predication of both the ²⁹Si and ¹⁷O chemical shifts (δ_i^Si and δ_i^O), the minimum cluster is one that gives the correct second neighbors to the nucleus of interest. Both the δ_i^Si and δ_i^O have reached convergence with respect to cluster size at the OH-terminated two tetrahedral (2T) shell (three atomic shells around Si and four atomic shells around O) model. At convergence, the calculated δ_i^Si values agree well (within ±1 ppm) with experimental data. The calculated ¹⁷O electric field gradient (EFG)-related parameters also agree with experimental data within experimental uncertainties. The calculation also reproduces small differences in δ_i^O for O sites with similar tetrahedral connectivities, but shows deviations up to about 10 ppm in relative difference for O sites with different tetrahedral connectivities. The poor performance for the latter is mainly due to the approximations of the HF method. Our study thus suggests that the ab initio calculation method is a reliable mean for predicting ²⁹Si and ¹⁷O NMR parameters for silicates. Such an approach should find application not only to well-ordered crystalline phases, but also to disordered materials, by combining with other techniques, such as the molecular dynamics simulation method.     

FLiNaK-ZrF_4体系的Raman光谱及量子化学计算研究

本文采用量子化学从头计算方法对FLiNaK-ZrF_4体系中锆的不同配位结构及不同阳离子的影响进行了模拟计算,锆的配位数与其特征峰的峰位呈线性关系。对不同配比的FLiNaK-ZrF_4体系样品进行拉曼光谱实验,结果表明,随着ZrF_4含量的增加,体系中自由F-浓度减小,导致锆的配位数降低并出现桥氟结构。高温熔融状态下,FLiNaK-0.37ZrF_4样品中六配位结构分解形成五配位和七配位结构。     

Superconducting properties of barium in three phases under high pressure from first principles

Electron-phonon coupling (EPC) in the three high-pressure phases of Ba is investigated using a pseudopotential plane-wave method based on density functional perturbation theory. The calculated values of superconducting critical temperature T c of Ba-I and Ba-II under pressure are consistent well with the trends observed experimentally. Moreover, Ba-V is found to be superconducting with a maximum T c exceeding 7.8 K at 45 GPa. With the increase of pressure, the values of T c increase in Ba I and Ba-II but the value of T c decreases in Ba-V. For Ba-I at pressures below 2 GPa, the increases of logarithmic average frequency ω log and electron-phonon coupling parameters λ both contribute to the enhancement of T c . For all the three phases at pressures above 2 GPa, T c is found to be primarily determined by λ . Further investigation reveals that for all the three phases, the change in λ with pressure can be explained mainly by change in the phonon frequency. Thus for Ba-II and Ba-V, although they exhibit completely different superconducting behaviors, their superconductivities have the same origin; the pressure dependence of T c is determined finally by the pressure dependence of phonon frequency.     

Phase transition and critical behavior of spin–orbital coupled spinel ZnV_2O_4

We present the temperature-dependent susceptibility and specific heat measurement of spinel ZnV_2O_4.The structural transition with orbital ordering and the antiferromagnetic transition with spin ordering were observed at 50 K and 37 K,respectively.By analysis of the hysteresis behavior between the specific heat curves obtained in warming and cooling processes,the structural transition was confirmed to be the first-order transition,while the antiferromagnetic transition was found to be of the second-order type.At the structural transition,the latent heat and entropy change were calculated from the excess specific heat,and the derivative of pressure with respect to temperature was obtained using the Clausius-Clapayron equation.At the magnetic transition,the width of the critical fluctuation region was obtained to be about 0.5 K by comparing with Gaussian fluctuations.In the critical region,the critical behavior was analyzed by using renormalization-group theory.The critical amplitude ratio A~+/A~- = 1.46,which deviates from the 3D Heisenburg model;while the critical exponent α is-0.011,which is close to the 3D XY model.We proposed that these abnormal critical behaviors can be attributed to strong spin-orbital coupling accompanied with the antiferromagnetic transition.Moreover,in the low temperature range(2-5 K),the Fermi energy,the density of states near the Fermi surface,and the low limit of Debye temperature were estimated to be2.42 eV,2.48 eV~(-1),and 240 K,respectively.     

Electron spin resonance study of a CuIr2S4 single crystal

A spinel sulphide CuIr₂S₄ single crystal, which exhibits an orbitally induced Peierls phase transition at ～230?K, is investigated by electron spin resonance (ESR) spectroscopy. The phase transition is clearly manifested on the ESR spectra. It is suggested that the ESR signals are produced by a few non-dimerized Ir⁴⁺ ions. Moreover, an extra ESR spectrum appears at low temperature in addition to the paramagnetic ESR signals of Ir⁴⁺ ions, which is suggested to be caused by the Jahn–Teller effect of the non-dimerized Ir⁴⁺ ions. From the ESR results, it is found that the Jahn–Teller splitting energy ΔE _JT is much smaller than the spin-dimerization gap.     

Structural phase transition in Bi2Te3 under high pressure

Structural change in Bi₂Te₃ under high pressure up to 16.6 GPa has been studied by powder x-ray diffraction. We observed two times of phase transitions at room temperature at the pressures of 8 and 14 GPa, respectively. According to our preliminary result on electrical resistance, it is reasonable to suppose that superconducting transition with T _c =2.8 K at the pressures of 10.2 GPa is observed in phase II. On the other hand, we found anomalies of the pressure dependences of lattice parameters and volume at around 2 GPa, which probably means the change in electrical structure on the Fermi surface.     

Structure instability of A2Al2B2O7 (A=K, Na) crystals

A theoretical study on the stabilities of the crystal structures of K₂Al₂B₂O₇ (KABO) and Na₂Al₂B₂O₇ (NABO) has been carried out using density functional theory with generalized-gradient corrections. All structures have been optimized by minimizing the total energies with respect to lattice constants and to the atomic coordinates within the unit cell. In the case of KABO, the structure with space group P321 always takes energy advantage over the whole volume range studied, whereas on both volume expansion and contraction, the original NABO structure becomes unstable, taking a structure transition from its space group to the P321 space group of KABO structure. The theoretical predictions were well realized in the experimental results on the structure changes of the (K_1−xNa_x)₂Al₂B₂O₇ solid solution system.     

First-principles study on electronic structure of Sr2Bi2O5 crystal

The electronic structure of Sr₂Bi₂O₅ is calculated by the GGA approach. Both of the valence band maximum and the conduction band minimum are located at Γ-point. This means that Sr₂Bi₂O₅ is a direct band-gap material. The wide energy-band dispersions near the valence band maximum and the conduction band minimum predict that holes and electrons generated by band gap excitation have a high mobility. The conduction band is composed of Bi 6p, Sr 4d and O 2p energy states. On the other hand, the valence band can be divided into two energy regions ranging from −9.5 to −7.9 eV (lower valence band) and from −4.13 to 0 eV (upper valence band). The former mainly consists of Bi 6s states hybridizing with O 2s and O 2p states, and the latter is mainly constructed from O 2p states strongly interacting with Bi 6s and Bi 6p states.     

Ab initio predicted elastic and thermodynamic properties of Imm2-BN under high pressure

The structural parameters, elastic constants, thermodynamic properties of Imm2-BN under high pressure were calculated via the density functional theory in combination with quasi-harmonic Debye approach. The results showed that the pressure has the significant effect on the equilibrium lattice parameters, elastic and thermodynamic properties of Imm2-BN. The obtained ground state structural parameters are in good agreement with previous theoretical results. The elastic constants, elastic modulus, and elastic anisotropy were determined in the pressure range of 0–90?GPa. Furthermore, by analyzing the B/G ratio, the brittle/ductile behavior under high pressure is evaluated and the elastic anisotropy of the Imm2-BN up to 90?GPa is studied in detail. Moreover, the pressure and temperature dependence of thermal expansion coefficient, heat capacity, Debye temperature, and Grüneisen parameter are predicted in a wide pressure (0–90?GPa) and temperature (0–1600?K) ranges. The obtained results are expected to provide helpful guidance for the future synthesis and application of Imm2-BN.     

High pressure X-ray diffraction study of CaMnO3 perovskite

Using a diamond anvil cell device and synchrotron radiation,the in-situ high-pressure structure of CaMnO3 has been investigated.In the pressure up to 36.5 GPa,no pressure-induced phase transition is observed.The pressure dependence on the lattice parameters of CaMnO3 is reported,and the relationship of the axial compression coefficients is βa 〉 βc 〉 βb.The isothermal bulk modulus K298=224（25） GPa is also obtained by fitting the pressure-volume data using the Murnaghan equation of state.     

Electronic properties of MoS2 nanoribbons with disorder effects

A theoretical study of electronic properties on MoS₂ nanoribbon is made on focusing the calculation of zero bias transport in the presence of disorders. Disorders including intrinsic and extrinsic vacancies and also weak uniform scatter defects are considered. The calculations are based on the tight-binding Green's function formalism by including an iterative procedure. The Slater–Koster transformations are used to determine the parameters. This model reduces the numerical calculation time. The unsaturated atoms at the edge of armchair (zigzag) ribbon induce some mid-gap states with nearly high (low) localization, which act as scattering centers. The antiresonances of created quasi-localized states due to vacancy cause the conductance of the armchair nanoribbon to decrease. Finally, the zigzag ribbon provides the highest sensitivity as well as selectivity between the smaller energy range, in the presence of the single weak scatter with potential value of 2 eV at the edge of the ribbon.     

High pressure X-ray diffraction study of SrMnO3 perovskite

Using a diamond anvil cell device and synchrotron radiation, the in-situ high-pressure structure of SrMnO3 has been investigated. At pressure up to 28.6 GPa, no pressure-induced phase transition is observed. The lattice parameters as a function of pressure is reported, and the relationship of the axial compression coefficients is β<,a>> β<,c>. The isothermal bulk modulus K<,298>=266(4) GPa is also obtained by fitting the pressure- volume data using the Murnaghan equation of state.     

High pressure X-ray diffraction study of CaMnO3 perovskite

Using a diamond anvil cell device and synchrotron radiation,the in-situ high-pressure structure of CaMnO3 has been investigated.In the pressure up to 36.5 GPa,no pressure-induced phase transition is observed.The pressure dependence on the lattice parameters of CaMnO3 is reported,and the relationship of the axial compression coefficients is βa >βc > βb.The isothermal bulk modulus K298=224(25)GPa is also obtained by fitting the pressure-volume data using the Murnaghan equation of state.