原子与双原子分子的碰撞反应

一个原子与宁一个同核双原子分子反应碰撞的过渡态构型是线型的, 计算了平衡临界键长的统计分布函数及其平均值、反应活化能和反应速率常数, 理论结果与实验值相符.     

Molecular Dynamical Simulation of Ice Phase Transition： Ice Ih to High-Density Amorphous

We put 5kbar and 12kbar on perfect ice Ih lattice at 77K and 180K. After 30000 simulation steps （in units of 10^-15 s）, high-density amorphous ice is formed. Four-site simple-pair potential TIP4P is used for molecular interactions and the rigid molecular model is employed. Phase transition processes are fitted by an exponential function, and different phase transition times τ are obtained from O-O radial distribution functions （366 and 359fs for 77K and 180K） and O-O-O angle distribution functions （126 and 116fs for 77K and 180K）.     

Molecular dynamics study of ice structural evolution

Molecular dynamics simulation is employed to study the structural evolution of low density amorphous ice during its compression from one atmosphere to 2.5 GPa. Calculated results show that high density amorphous ice is formed at an intermediate pressure of -1.0 GPa; the O-O-O bond angle ranges from 83° to 113°, and the O-H… O bond is bent from 112° to 160°. Very high density amorphous ice is obtained by quenching to 80 K and decompressing the ice to ambient pressure from 160 K/1.3 GPa or 160 K/1.7 GPa; and the next-nearest O-O length is found to be 0.310 nm, just 0.035 nm beyond the nearest O-O distance of 0.275 nm.     

有机自旋器件磁性渗透层中双极化子对自旋极化输运的影响

根据实验发现的有机器件如Co/有机半导体/La_0.7Sr_0.3MnO₃中磁性原子渗透现象, 利用自旋漂移-扩散方程, 理论研究了磁性渗透层中极化子-双极化子的转化对自旋极化输运的影响. 研究发现: 磁性渗透层具有不同于纯净有机层的迁移率和自旋反转时间, 都将影响极化子-双极化子的转化, 进而影响自旋极化的输运; 在磁性渗透层中极化子自旋反转时间的劈裂是引起自旋弛豫的主要因素, 而极化子和双极化子之间的转化是重要因素.     

Lattice Dynamical Simulation of Guest-Host Interaction in N2 Clathrate Hydrate

We perform the lattice dynamical simulation studies of hydrate host lattice interacting with Xe, Ar, and N2 atoms/molecules. The calculated results show that the well-defined peaks （2.0meV and 3.8meV） and another peak （6.2 me V） are assigned to the vibrations of N2 molecules in large and small cages, respectively. It is confirmed that the double N2 molecule occupancies of large cage lead to filling of the mode gap between the small cage and the large cage.     

Structural Evolution of Compressing Amorphous Ice

Molecular dynamics simulation is employed to study structural evolution during compressing low density amorphous ice from one atmosphere to 2.5 GPa. The calculated results show that high density amorphous ice is formed under intermediate pressure of about 1.0GPa and O-O-O angle ranges from about 83°to 113° and O-H…O is bent from 112° to 160°The very high density amorphous ice is also formed under the pressure larger than 1.4 GPa and interstitial molecules are found in 0.3-0.4 A just beyond the nearest O-O distance. Low angle O-H… O disappears and it is believed that these hydrogen bonds are broken or re-bonded under high pressures.     

Lattice dynamics study of low energy guest-host coupling in clathrate hydrate

Our lattice dynamics simulation of Xe-hydrate with four-site TIP4P oxygen-shell model can accurately reproduce each peak position in the inelastic incoherent neutron scattering spectrum at the acoustic band （below 15 meV） and yield correct relative intensity. Based on the results, the uncertain profile at ～6 meV is assigned to anharmonic guest modes coupled strongly to small cages. Blue shift is proposed in phonon dispersion sheet in the case of anticrossing and found to be an evident signal for guest-host coupling that explains the anomalous thermal conductivity of clathrate hydrate.     

NEUTRON SCATTERING AND LATTICE DYNAMICAL STUDIES OF THE HIGH-PRESSURE PHASE ICE (II)

Lattice dynamical calculations have been carried out for ice II based on the force field constructed for ice Ih. In order to fully understand ice II inelastic neutron scattering spectra, the decomposed phonon density of states was shown mode by mode. Calculated results have shown that the hydrogen bond force constant between the six-molecule rings is significantly weaker, 75eV/nm², compared with the force constant, 220eV/nm², within the rings. Inelastic neutron scattering spectra of clathrate hydrate H₂O+He are almost the same as ice II. This means that the absorption of He atoms cannot affect the bond strengths of the ice II host lattice. Based on the force field model for ice II, the van der Waals interactions between water molecules and helium atoms are considered. The results obtained are consistent with experimental data. Lattice dynamical calculations have been carried out for ice II using seven rigid pairwise potentials. It was found that MCY makes the stretching and bending interactions in ice II too weak and makes the O-O bond length too long (～5%), thus its lattice densities are obviously lower than other potential lattices or experimental values.