J. Cent. South Univ. (2017) 24: 2542-2549
DOI: https://doi.org/10.1007/s11771-017-3667-6
![](/web/fileinfo/upload/magazine/12630/317131/image002.jpg)
Volatilization behaviors of molybdenum and sulfur in vacuum decomposition of molybdenite concentrate
ZHOU Yue-zhen(周岳珍)1, 2, LU Yong(卢勇)1, 2, LIU Da-chun(刘大春)1, 2, 3,
CHEN Xiu-min(陈秀敏)1, 2, 3, LI Hui(李慧)1, 2, LI Wei(李玮)1, 2
1. National Engineering Laboratory of Vacuum Metallurgy, Kunming University of Science and Technology,Kunming 650093, China;
2. Key Laboratory of Vacuum Metallurgy of Non-ferrous Metals of Yunnan Province, Kunming 650093, China;
3. State Key Laboratory of Complex Non-ferrous Metal Resources Clear Utilization in Yunnan Province,Kunming 650093, China
Central South University Press and Springer-Verlag GmbH Germany, part of Springer Nature 2017
Abstract: Thermodynamic calculation, ab initio molecular dynamics (AIMD) and vacuum decomposition experiments were performed to study the volatilization behaviors of Mo and S from molybdenite concentrate by vacuum decomposition. In thermodynamic calculation, starting decomposition temperatures of reactions were calculated, and saturated vapor pressures of Mo, S and MoS2 were also analyzed. In AIMD, geometries of the Sn (n≤8), Mom (m≤8) and MomSn (m+n≤8) clusters have been optimized using density functional theory (DFT) with generalized gradient approximation (GGA). And these clusters were simulated in DFT with Cambridge Sequential Total Energy Package (CASTEP) code of Material Studio software. Structures and stabilities of these clusters before and after molecular dynamics simulations were discussed, and diffusion coefficients were also calculated. In vacuum decomposition experiments, relationship between heat preservation time and volatilization rate of Mo and S was obtained, while the constant temperature and chamber pressure were 1823 K and 5–35 Pa, respectively. Above all, both the theoretical and experimental results showed that volatilization behaviors of Mo and S during vacuum decomposition process of molybdenite concentrate were as follows: Mo could partly evaporate into the condensate in the form of clusters, and S could easily evaporate into the condensate.
Key words: volatilization behavior; vacuum decomposition; thermodynamics; ab initio molecular dynamics; clusters
1 Introduction
Molybdenite is the major industrial mineral for molybdenum [1]. At present, molybdenite concentrate is converted to molybdenum oxide by oxidizing roasting firstly, and then the resultant oxidized concentrate is purified by conversion of the molybdenum oxide to ammonium molybdate, finally commercial molybdenum metal can be obtained by reduction of molybdenum oxide with H2 in two or more stages [2]. However, there are several serious problems in processing molybdenite by this method: valuable elements in the concentrate such as rhenium and selenium cannot be recovered easily to high degrees and air pollution may occur due to the emission of SO2 gas [3].
Comparing to traditional metallurgical process, vacuum metallurgy owns many advantages, such as high metal recovery rate, less pollution, less energy consuming. Vacuum distillation and decomposition are the important parts of vacuum metallurgy, the former is regarded as one of the most effective and environment- friendly methods for metal separation, preparation of high purity metal and recycling of secondary metal resources,and the latter is mainly applied to thermal decomposition of compounds [4–7].
Many scholars made several meaningful researches in vacuum decomposition of molybdenite concentrate. On the view of experiments, CHEN [8] investigated vacuum decomposition process of analytic grade molybdenum disulfide and molybdenite concentrate respectively, and useful experimental parameters were obtained. WANG et al [9] studied the key steps of thermal decomposition process of molybdenum concentrate in vacuum and verified it by vacuum decomposition experiment. On the view of molecular dynamics simulation, structural and other properties of Mon(n=2–55) [10], Mon(n=2–8) [11], Sn(n=2–8) [12], Mo5Sn(n=5–15) [13], MomSn(n=1–6, m=n–3n) [14] and MoS6 clusters [15] were calculated based on DFT. And LIU et al [16] optimized MoS2 crystal structure and simulated the thermal decomposition of MoS2 by DFT.
The main purpose of this research is to study the volatilization behaviors of Mo and S during vacuum decomposition process of molybdenite concentrate, from the view of both thermodynamics and molecular dynamics simulation, and further verify it through vacuum decomposition experiments.
2 Theoretical analysis of volatilization behaviors
2.1 Thermodynamic calculation
The possible mechanism of thermal decomposition of MoS2 under vacuum was as follows [17], and the reaction’s Gibbs free energy at different pressures were shown in Table 1 [9].
4MoS2=2Mo2S3+S2 (1)
Mo2S3=2Mo+1.5S2 (2)
Table 1 Starting decomposition temperature of MoS2 and Mo2S3 at different pressures
![](/web/fileinfo/upload/magazine/12630/317131/image004.jpg)
As seen in Table 1, the starting decomposition temperature kept decreasing as pressure fell. And comparing to reactions happened under condition of atmospheric pressure, whose starting decomposition temperatures were at least above 2000 K, the temperatures were 1491 K and 1575 K while the pressure fell to 10 Pa. Thus, it is feasible to separate Mo and S in molybdenite concentrate by vacuum decomposition.
In order to investigate the volatilization behaviors of Mo and S during vacuum decomposition process, saturated vapor pressure of Mo, S and MoS2 were calculated with the van LAAR equation expressed as follows [18]:
lgPΘ=AT–1+BlgT+CT+D (3)
lgPΘ=–AT–1+B (4)
where PΘ was the saturated vapor pressure, Pa; A, B, C and D were evaporation constants [19]; T was the thermodynamic temperature, K. The melting points and saturated vapor pressure of Mo, S and its compound were shown in Table 2 and Fig. 1.
Table 2 Melting points of Mo, S and its compound
![](/web/fileinfo/upload/magazine/12630/317131/image005.jpg)
![](/web/fileinfo/upload/magazine/12630/317131/image007.jpg)
Fig. 1 Relationship between temperature and saturated vapor pressure
As seen in Table 2 and Fig. 1, saturated vapor pressure of S was much larger than 10 Pa as temperature rose, which indicated that S can easily evaporate into the condensate during the whole vacuum decomposition process. While the temperature was below starting decomposition temperature (1491 K or 1575 K), the saturated vapor pressure of MoS2 was smaller than the pressure in furnace. Namely, Mo could not evaporate into the condensate in the form of MoS2 before the vacuum decomposition reactions started. While the temperature was above starting decomposition temperature, Mo was still in solid-state. Namely, Mo could not evaporate into the condensate in the form of pure substance after reactions started.
According to thermodynamic calculation results, volatilization behaviors of Mo and S during vacuum decomposition process of molybdenite concentrate were as follows: Mo nearly enriched in the residual, and S could easily evaporate into the condensate.
2.2 AIMD
2.2.1 Computational methods
Based on the authors’ previous theoretical calculation, ground-state structures of Sn(n≤8) clusters, Mom(m≤8) clusters and MomSn(m+n≤8) clusters were obtained [20]. The computational details were as follows: calculations were performed with DFT implemented in the Dmol3 package in the Materials Studio software [21]. In the geometric optimizations, the exchange correlation interaction was treated within the GGA using the PW91 functional [22]. The convergent threshold was set to 10–5Ha for energy change, 0.002 Ha/
for forces and 0.005
for displacement. In the electronic structure calculations, the double numerical basis set including d-polarization function (DND) was utilized to describe the valence electrons and the core electrons were described with effective core potential. Self-consistent field (SCF) calculations were done with a convergence criterion of 10–6 Ha on the total energy. The direct inversion in an iterative subspace (DIIS) approach was used to speed up SCF convergence. A 0.005 Ha of smearing was applied to the orbital occupation.
The dynamics simulation calculations of the optimized Sn clusters, Mom clusters and MomSn clusters were carried out with DFT by using the program package CASTEP [23]. Exchange and correlation effects were treated with GGA implemented in PW91 as exchange- correlation functional. It was simulated by NVT system, with the simulation time 10.0 ps and time step size 1.0 fs, total 10000 step. And the simulation temperature was 1823 K. The Nosé thermostat [24] and Andersen barostat [25] were used in the dynamics simulation. Brillouin- zone sampling of 1×1×1 k-points and cutoff energy of 180.0 eV in plane wave basis sets were verified to be sufficient to obtain good numerical accuracy with reasonable computational costs.
2.2.2 Structures and stabilities of clusters
Ground-state structures and simulated structures of Mom clusters are shown in Fig. 2. Structural changes were as follows: 1) Comparing to the ground-state structures of Mom (m=2, 5, 6) clusters, structures after dynamics simulation were basically unchanged. Though bond length changed to a certain degree, Mo–Mo bonds did not break in these clusters; 2) Due to the relatively great changing bond length, several Mo–Mo bonds broke or formed, and structures of Mom (m=3, 4, 7, 8) clusters after dynamics simulation changed obviously.
Thus, the structural stability of Mom (m=2, 5, 6) clusters was larger than it of Mom (m=3, 4, 7, 8) clusters. Namely, Mom (m=2, 5, 6) clusters were more likely to steadily exist under the condition of vacuum and high temperature.
Ground-state structures and simulated structures of Sn clusters were shown in Fig. 3, structural changes were as follows: 1) Comparing to the ground-state structures of Sn (n=2, 3, 4, 6) clusters, structures after dynamics simulation were basically unchanged. 2) Comparing to the ground-state structures of Sn (n=5, 7, 8) clusters, structures after dynamics simulation changed obviously. The spatial structure of S5 cluster transformed from cyclic structure to chain structure because several S—S bonds in this cluster broke. Sn (n=7, 8) clusters’ structures changed greatly, and they broke into several parts which consisted with less atoms.
Thus, the structural stability of Sn (n=2, 3, 4, 6) clusters was larger than it of Sn (n=7, 8) clusters. Namely, Sn (n=2, 3, 4, 6) clusters were more likely to steadily exist under the condition of vacuum and high temperature.
Part ground-state structures and simulated structures of MomSn clusters were shown in Fig. 4, structural changes were as follows: 1) To a certain degree, all kinds of the MomSn (m+n≤8) clusters’ structures changed, and only 10 kinds of these clusters were basically unchanged after dynamics simulation. These MomSn clusters were MoS, MoS2, MoS3, MoS4, Mo2S2, Mo2S4, Mo3S2, Mo4S4, Mo5S2 and Mo6S; 2) Due to the conspicuous changing bond of Mo—Mo, S—S and Mo—S, the other 18 kinds of MomSn clusters’ structures changed obviously, such as Mo2S, Mo2S5, Mo3S4 and Mo7S.
![](/web/fileinfo/upload/magazine/12630/317131/image009.jpg)
Fig. 2 Structures of Mom clusters before (ground-state) and after (1823 K) simulations
![](/web/fileinfo/upload/magazine/12630/317131/image011.jpg)
Fig. 3 Structures of Sn clusters before (ground-state) and after (1823 K) simulations
Thus, the structural stability of the former 10 kinds of MomSn clusters was larger than that of the latter 18 kinds of MomSn clusters. Namely, the former 10 kinds of MomSn clusters were more likely to steadily exist under the condition of vacuum and high temperature.
In conclusion, Mom (m=2, 5, 6) clusters, Sn (n=2, 3, 4, 6) clusters and the 10 kinds of MomSn clusters mentioned before were more likely to steadily exist under the condition of vacuum and high temperature.
2.2.3 Diffusion property of clusters
In order to investigate the volatilization behaviors of the 17 kinds of clusters which could steadily exist under the condition of vacuum and high temperature, diffusion property was calculated based on the Einstein law as follows [26, 27]:
(5)
(6)
where dMS is the mean square value of particles’ displacement; ri(t) is location of particle i while time equaled to t; N is the sum of particles; < > was the ensemble average of all particles during the dynamics simulation process; D is diffusion coefficient. The diffusion coefficients of clusters under the condition of vacuum and high temperature are shown in Tables 3 and 4.
As seen in Tables 3 and 4, diffusion coefficients of these clusters kept decreasing with the increasing atomic number of the clusters. And while the ratio between atomic number of molybdenum and sulfur increasing, diffusion coefficients kept decreasing too.
According to the diffusion coefficients of the 17 kinds of clusters, which were more likely to steadily exist under the condition of vacuum and high temperature, several rules were found as follows: 1) Diffusion coefficient of Mo2 cluster was relatively larger than that of other clusters in Mom clusters; 2) Diffusion coefficient of S2, S3 and S4 clusters were relatively larger than that of other clusters in Sn clusters, and their values were about the double of Mo2 cluster’s. 3) Diffusion coefficients of MoS, MoS2, MoS4 and Mo2S4 clusters were relatively larger than those of other clusters in MomSn clusters, and their values were approximately between the values of Mom clusters and Sn clusters.
![](/web/fileinfo/upload/magazine/12630/317131/image017.jpg)
Fig. 4 Structures of MomSn cluster before (ground-state) and after (1823 K) simulations
Table 3 Diffusion coefficients (D/(
2·s–1)) of Mom clusters and Sn clusters after simulations
![](/web/fileinfo/upload/magazine/12630/317131/image019.jpg)
Table 4 Diffusion coefficient (D/(
2·s–1)) of MomSn clusters after simulations
![](/web/fileinfo/upload/magazine/12630/317131/image021.jpg)
Thus, among all kinds of Mom, Sn and MomSn clusters simulated here, 8 kinds of clusters’ structures were basically unchanged after dynamics simulation, and the diffusion coefficients of these clusters were relatively large. Namely, these clusters could steadily exist under the condition of vacuum and high temperature, and could easily evaporate into the condensate relatively easily. These clusters were Mo2, S2, S3, S4, MoS, MoS2, MoS4 and Mo2S4 clusters.
According to AIMD results, volatilization behaviors of Mo and S during vacuum decomposition process of molybdenite concentrate were as follows: Mo and S could evaporate into the condensate in the form of clusters, and the volatilization ability of Mo was much weaker than that of S.
In conclusion, theoretical analysis of volatilization behaviors of Mo and S during vacuum decomposition process of molybdenite concentrate was as follows: 1) Thermodynamic calculation results demonstrated that Mo nearly enriched in the residual, but AIMD results indicated that Mo could evaporate into the condensate in the form of clusters, the conclusions draw here were contradictory; 2) Both thermodynamic calculation and AIMD results proved that S could easily evaporate into the condensate.
Thus, it was necessary to perform experimental study to further investigate the volatilization behaviors of Mo and S during vacuum decomposition process.
3 Vacuum decomposition experiments of molybdenite concentrate
3.1 Experimental conditions and materials
3.1.1 Experimental conditions
In order to make experimental results as a reliable evidence to further prove the volatilization behaviors of Mo and S during vacuum decomposition process, appropriate experimental conditions were very important.
We selected 1823 K as the constant temperature while studying the influence of heat preservation time on the volatilization behaviors of Mo and S, and the chamber pressure was 5–35 Pa. Firstly, the starting decomposition temperatures of vacuum decomposition reactions were 1491 K and 1575 K in thermodynamic calculation, and the simulation temperature was set 1823 K in AIMD. Thus, it was appropriate to select 1823 K as the constant temperature. Secondly, the pressure obtained in our experimental equipment was 5–35 Pa, and this narrow range of pressure was basically consistent with the theoretical analysis.
3.1.2 Experimental materials
The chemical composition of molybdenite concentrate used in this experiment was shown in Table 5.
Table 5 Chemical composition of molybdenite concentrate (mass fraction, %)
![](/web/fileinfo/upload/magazine/12630/317131/image022.jpg)
As seen in Table 5, the main chemical components were Mo and S. And several impurities, such as O, Al, Si, Fe and Cu, were found in molybdenite concentrate. Thus, the experimental materials used here was high grade molybdenite concentrate, and its main compound was MoS2.
Here, we selected molybdenite concentrate as experimental materials for the following reasons: Firstly, the experimental results could be a reliable evidence to further prove the theoretical analysis results because the molybdenite concentrate mainly contained Mo and S; Secondly, both experimental and theoretical results could further provide heat preservation time selection in vacuum decomposition experiments of molybdenite concentrate.
3.2 Experimental equipment and methods
Vacuum decomposition experiments were carried out in a laboratory-scale vertical vacuum distillation furnace, as shown in Fig. 5. It mainly consists of a chamber, condensing zone, heating zone and vacuum pump systems [28].
A cylindrical material, with diameter and thickness of 23 mm and 18 mm, respectively, was obtained through tableting. Then the material was removed into crucible and put into the vacuum furnace. The feeding material was heated to different temperatures at a heating rate of 15 K/min and each temperature was kept for different time when the chamber pressure was (5–35) Pa. Then, the heater power was switched off after the end of heating period and vacuum was unnecessary to maintain while the temperature in furnace below 323 K. In the end, residual and condensate were gathered, weighed and chemically analyzed, respectively.
![](/web/fileinfo/upload/magazine/12630/317131/image024.jpg)
Fig. 5 Schematic diagram of vertical vacuum distillation furnace
3.3 Experimental results and discussion
Vacuum decomposition experiments were carried under the pressure of (5–35) Pa at 1823 K for different heat preservation time, from 15 to 120 min. Volatilization rate of Mo or S was calculated as follow:
(7)
where r is the volatilization rate, %; m0 is the mass of Mo or S in molybdenite concentrate, g; mr is the mass of Mo or S in residual, g. The effect of heat preservation time on volatilization rate of Mo and S is shown in Fig. 6.
![](/web/fileinfo/upload/magazine/12630/317131/image028.jpg)
Fig. 6 Relationship between heat preservation time and volatilization rate
As seen in Fig. 6, volatilization rate of S increased from 58.75% to 99.82% while the heat preservation time extended from 15 to 120 min. Namely, S could easily evaporate into the condensate and its volatilization rate nearly reached to 100%. Thus, experimental results were consistent with the thermodynamic calculation and AIMD results about the volatilization behavior of S during vacuum decomposition process. While extending the heat preservation time, volatilization rate of Mo increased from 0.93% to 6.42%. Namely, Mo could partly evaporate into the condensate and its volatilization rate was much less than that of S. Although thermodynamic calculation results indicated that Mo nearly enriched in the residual, it still could safely draw the conclusion that Mo could partly evaporate into the condensate, according to both AIMD and experimental results.
According to theoretical and experimental results, volatilization behaviors of Mo and S during vacuum decomposition process of molybdenite concentrate were as follows: Mo could partly evaporate into the condensate in the form of clusters, and S could easily evaporate into the condensate.
Volatilization rate of S increased from 58.75% to 99.30% while the heat preservation time extended from 15 to 60 min, and it slightly increased to 99.82% while doubled the heat preservation time. In the meantime, volatilization rate of Mo still increased from 5.96% to 6.42%. Thus, on the premise of desulphurization, it was reasonable to select 60 min as the heat preservation time while the vacuum decomposition experiments were carried at 1823 K.
4 Conclusions
Thermodynamic calculation results showed that Mo nearly enriched in the residual, and S could easily evaporate into the condensate. However, AIMD results indicated that Mo and S could evaporate into the condensate in the form of clusters, and the volatilization ability of Mo was much weaker than that of S. In order to further investigate the volatilization behaviors of Mo and S, relationship between heat preservation time and volatilization rate of Mo and S was obtained through vacuum decomposition experiments of molybdenite concentrate. Taking both theoretical and experimental results into consideration, volatilization behaviors of Mo and S were as follows: Mo could partly evaporate into the condensate in the form of clusters, and S could easily evaporate into the condensate. And heat preservation time selection in vacuum decomposition experiments was also investigated as follows: On the premise of desulphurization, it was reasonable to select 60 min as the heat preservation time while the vacuum decomposition experiments were carried at 1823 K.
References
[1] ZHANG Qi-xiu, ZHAO Qin-sheng. Metallurgy of molybdenum and tungsten [M]. Beijing: Metallurgical Industry Press, 2005: 59–61. (in Chinese)
[2] WANG Lei, GUO Pei-min, PANG Jian-ming, LUO Lin-geng, ZHAO Pei. Phase change and kinetics of vacuum decomposition of molybdenite concentrate [J]. Vacuum, 2015, 116: 77–81.
[3] BLANCO E, SOHN H Y, HAN G, HAKOBYAN K Y. The kinetics of oxidation of molybdenite concentrate by water vapor [J]. Metall Mater Trans B, 2007, 38(4): 689–693.
[4] LI Dong-sheng, DAI Yong-nian, YANG Bin, LIU Da-chun, DENG Yong. Purification of indium by vacuum distillation and its analysis [J]. Journal of Central South University, 2013, 20(2): 337–341.
[5] JIANG Wen-long, DENG Yong, YANG Bin, LIU Da-chun, DAI Yong-nian, XU Bao-qiang. Application of vacuum distillation in refining crude indium [J]. Rare Met, 2013, 32(6): 627–631.
[6] KONG Xiang-feng, YANG Bin, XIONG Heng, LIU Da-chun, XU Bao-qiang. Removal of impurities from crude lead with high impurities by vacuum distillation and its analysis [J]. Vacuum, 2014, 105:17–20.
[7] YUAN Hai-bin, YANG Bin, XU Bao-qiang, YU Qing-chun, FENG Yue-bin, DAI Yong-nian. Aluminum production by carbothermo- chlorination reduction of alumina in vacuum [J]. Trans Nonferrous Met Soc China, 2010, 20(8): 1505–1510.
[8] CHEN Jie. Studies on vacuum decomposition of molybdenite concentrate [D]. Kunming: Kunming University of Science and Technology, 2009: 16–49. (in Chinese)
[9] WANG Lei, GUO Pei-min, PANG Jian-ming, ZHAO Pei. Thermodynamic analysis of vacuum decomposition process of molybdenum concentrate [J]. The Chinese Journal of Nonferrous Metals, 2015, 25(1): 190–196. (in Chinese)
[10] ZHANG Wen-qin, RAN Xiao-rong, ZHAO Hai-tao, WANG Li-chang. The nonmetallicity of molybdenum clusters [J]. J Chem Phys, 2004, 121(16): 7717–7724.
[11] LEI Xue-ling. Theoretical study of small Mo clusters and molecular nitrogen adsorption on Mo clusters [J]. Chin Phys B, 2010, 19(10): 17103.
[12] BAI Yu-lin, CHEN Xiang-rong, YANG Xiang-dong, LU Peng-fei. Structures of small sulfur clusters Sn (n=2–8) from Langevin molecular dynamics methods [J]. Acta Phys Chim Sin, 2003, 19(12): 1102–1107.
[13] MURUGAN P, KUMAR V, KAWAZOE Y, OTA N. Understanding the structural stability of compound Mo-S clusters at sub-nanometer level Mater [J]. Trans, 2007, 48(4): 658–661.
[14] MURUGAN P, KUMAR V, KAWAZOE Y, OTA N. Ab initio study of structural stability of Mo-S clusters and size specific stoichiometries of magic clusters [J]. J Phys Chem A, 2007, 111(14): 2778–2782.
[15] WU Ni, ZHANG Chang-fu, ZHOU Qi, HUANG Xin, ZHANG Yong-fan, DING Kai-ning, WANG Bin. DFT study on the electronic and structural properties of MoS6-/0 clusters [J]. Chinese J Struc Chem, 2013, 32(7): 1046–1054.
[16] LIU Da-chun, CHEN Xiu-min, LI Dian-jun, WANG Fei, LUO Xiao-gang, YANG Bin. Simulation of MoS2 crystal structure and the experimental study of thermal decomposition [J]. J Mol Struct, 2010, 980(1–3): 66–71.
[17] BUKER D O. Process for thermal dissociation of molybdenum U.S. Patent No. 3,966,459. [P] 1976.
[18] KONG Xiang-feng, YANG Bin, XIONG Heng, KONG Ling-xin, LIU Da-chun, XU Bao-qiang. Thermodynamics of removing impurities from crude lead by vacuum distillation refining [J]. Trans Nonferrous Met Soc China, 2014, 24(6): 1946–1950.
[19] DAI Yong-nian, YANG Bin. Vacuum metallurgy of nonferrous metal [M]. Beijing: Metallurgical Industry Press, 2009: 38–40. (in Chinese)
[20] ZHOU Yue-zhen, LU Yong, LIU Da-chun, CHEN Xiu-min, LI Wei, XU Jun-jie. Geometries, stabilities, and electronic properties of MomSn (m+n≤8) clusters [J]. The Chinese Journal of Nonferrous Metals, 2015, 25(8): 2251–2258. (in Chinese)
[21] DELLEY B. An all-electron numerical method for solving the local density functional for polyatomic molecules [J]. J Chem Phys, 1990; 92(1): 508–517.
[22] PERDEW J P, WANG Yue. Pair-distribution function and its coupling-constant average for the spin-polarized electron gas [J]. Phys Rev B, 1992, 46(20): 12947.
[23] SEGALL M D, LINDA P J D, PROBERT M J, PICKARD C J, HASNIP P J, CLARK S J, PAYNE M C. First-principles simulation: Ideas, illustrations and the CASTEP code [J]. J Phys: Condens Matter, 2002, 14(11): 2717–2744.
[24] NOSE S. A unified formulation of the constant temperature molecular dynamics methods [J]. J Chem Phys, 1984, 81(1): 511–519.
[25] ANDERSEN H C. Molecular dynamics simulations at constant pressure and/or temperature [J]. J Chem Phys, 1980, 72(4): 2384–2393.
[26] AHMED A, WU Jiang-tao. Critical anomaly and finite size scaling of the self-diffusion coeffcient for Lennard Jones fluids by non- equilibrium molecular dynamic simulation [J]. Chin Phys B, 2011, 20(10): 106601.
[27] GAO Song, WANG Fei, HE Dong-qi, ZHANG Shou-yu, BAO Shang-lian. An efficient non-iterative algorithm to estimate the apparent diffusion coefficient of water molecules in the magnetic resonance diffusion imaging experiments [J]. Chin Sin Bull, 2014, 59(8): 734–739.
[28] LI Liang, TIAN Yang, LIU Da-chun, ZHOU Hou-jun, DAI Yong-nian, YANG Bin. Pretreatment of lead anode slime with low silver by vacuum distillation for concentrating silver [J]. Journal of Central South University, 2013, 20(3): 615–621.
(Edited by HE Yun-bin)
Cite this article as: ZHOU Yue-zhen, LU Yong, LIU Da-chun, CHEN Xiu-min, LI Hui, LI Wei. Volatilization behaviors of molybdenum and sulfur in vacuum decomposition of molybdenite concentrate [J]. Journal of Central South University, 2017, 24(11): 2542–2549. DOI:https://doi.org/10.1007/s11771-017-3667-6.
Foundation item: Projects(1202271, 51104078) supported by the National Natural Science Foundation of China; Project(IRT1250) supported by the Program for Innovative Research Team in University of Ministry of Education of China
Received date: 2015-10-10; Accepted date: 2017-01-20
Corresponding author: LIU Da-chun, PhD; Tel: +86–871–65114017; E-mail: lcd_2002@sina.com