The density functional theory （DFT） is the most popular method for evaluating bond dis- sociation enthalpies （BDEs） of most molecules. Thus, we are committed to looking for alternative methods that can balance the...The density functional theory （DFT） is the most popular method for evaluating bond dis- sociation enthalpies （BDEs） of most molecules. Thus, we are committed to looking for alternative methods that can balance the computational cost and higher precision to the best for large systems. The performance of DFT, double-hybrid DFT, and high-level com- posite methods are examined. The tested sets contain monocyclic and polycyclic aromatic molecules, branched hydrocarbons, small inorganic molecules, etc. The results show that the mPW2PLYP and G4MP2 methods achieve reasonable agreement with the benchmark val- ues for most tested molecules, and the mean absolute deviations are 2.43 and 1.96 kcal/mol after excluding the BDEs of branched hydrocarbons. We recommend the G4MP2 is the most appropriate method for small systems （atoms number≤20）; the double-hybrid DFT methods are advised for large aromatic molecules in medium size （20≤atoms number≤50）, and the double-hybrid DFT methods with empirical dispersion correction are recommended for long-chain and branched hydrocarbons in the same size scope; the DFT methods are ad- vised to apply for large systems （atoms number〉50）, and the M06-2X and B3P86 methods are also favorable. Moreover, the differences of optimized geometry of different methods are discussed and the effects of basis sets for various methods are investigated.展开更多
A large database is desired for machine learning(ML) technology to make accurate predictions of materials physicochemical properties based on their molecular structure.When a large database is not available,the develo...A large database is desired for machine learning(ML) technology to make accurate predictions of materials physicochemical properties based on their molecular structure.When a large database is not available,the development of proper featurization method based on physicochemical nature of target proprieties can improve the predictive power of ML models with a smaller database.In this work,we show that two new featurization methods,volume occupation spatial matrix and heat contribution spatial matrix,can improve the accuracy in predicting energetic materials' crystal density(ρ_(crystal)) and solid phase enthalpy of formation(H_(f,solid)) using a database containing 451 energetic molecules.Their mean absolute errors are reduced from 0.048 g/cm~3 and 24.67 kcal/mol to 0.035 g/cm~3 and 9.66 kcal/mol,respectively.By leave-one-out-cross-validation,the newly developed ML models can be used to determine the performance of most kinds of energetic materials except cubanes.Our ML models are applied to predict ρ_(crystal) and H_(f,solid) of CHON-based molecules of the 150 million sized PubChem database,and screened out 56 candidates with competitive detonation performance and reasonable chemical structures.With further improvement in future,spatial matrices have the potential of becoming multifunctional ML simulation tools that could provide even better predictions in wider fields of materials science.展开更多
The parameters of embedded atom method for elements Ce, Th and Yb were determined by fitting the lattice constants, the cohesive energy, the monovacancy formation energy and the bulk modulus of elements. The alloy pot...The parameters of embedded atom method for elements Ce, Th and Yb were determined by fitting the lattice constants, the cohesive energy, the monovacancy formation energy and the bulk modulus of elements. The alloy potential was taken as the form of Johnson′s. The formation enthalpies of Th-Ce, Th-Yb and Ce-Yb binary alloys systems and Ce-Th-Yb ternary alloy were calculated with the present embedded atom potentials. The calculations for binary alloys are in good agreement with the results calculated with Miedema′s theory. As for the ternary alloy, the calculated formation enthalpies are in good agreement with those extrapolated from the formation enthalpies of constitutive binary alloys by Toop′s model.展开更多
The modified embedded atom method proposed by authors has been applied to calculating the enthalpies of formation of random alloys and the ordered intermetallic compounds for noble metal binary systems bearing Rh or I...The modified embedded atom method proposed by authors has been applied to calculating the enthalpies of formation of random alloys and the ordered intermetallic compounds for noble metal binary systems bearing Rh or Ir. The present results are in good agreement with those of Miedema theory, available experiments and the first-principles quantum mechanics calculations. The present results indicate that Cu-Rh, Cu-lr, Ag-Rh, Ag-lr, Au-Rh, Au-lr, Pd-Rh and Pd-lr systems are repulsive, however, IMi-Rh, Ni-lr, Pt-lr, Pt-Rh and Rh-lr systems form solid solutions and Ni-Rh, Ni-lr and Pt-Rh show ordering tendency.展开更多
The purpose of this paper is to present a novel way to building quantitative structure-property relationship(QSPR) models for predicting the gas-to-benzene solvation enthalpy(ΔHSolv) of 158 organic compounds based on...The purpose of this paper is to present a novel way to building quantitative structure-property relationship(QSPR) models for predicting the gas-to-benzene solvation enthalpy(ΔHSolv) of 158 organic compounds based on molecular descriptors calculated from the structure alone. Different kinds of descriptors were calculated for each compounds using dragon package. The variable selection technique of enhanced replacement method(ERM) was employed to select optimal subset of descriptors. Our investigation reveals that the dependence of physico-chemical properties on solvation enthalpy is a nonlinear observable fact and that ERM method is unable to model the solvation enthalpy accurately. The standard error value of prediction set for support vector machine(SVM) is 1.681 kJ ? mol^(-1) while it is 4.624 kJ ? mol^(-1) for ERM. The results established that the calculated ΔHSolvvalues by SVM were in good agreement with the experimental ones, and the performances of the SVM models were superior to those obtained by ERM one. This indicates that SVM can be used as an alternative modeling tool for QSPR studies.展开更多
A numerical analysis was carried out to study the moving boundary problem in the physical process of pulsed Nd-YAG laser surface melting prior to vaporization.The enthalpy method was applied to solve this two-phase ax...A numerical analysis was carried out to study the moving boundary problem in the physical process of pulsed Nd-YAG laser surface melting prior to vaporization.The enthalpy method was applied to solve this two-phase axisymmetrical mehing problem.Computational results of tempera ture fields were obtained,which provide useful information to practical lair treatment processing. The validity of enthalpy method in solving such problems is presented.展开更多
In this work, the group contribution method of Chickos et al. was applied to estimate the fusion enthalpy of isonicotinic acid, and the obtained result(29.2 k J·mol^(-1)) showed a large difference with the value(...In this work, the group contribution method of Chickos et al. was applied to estimate the fusion enthalpy of isonicotinic acid, and the obtained result(29.2 k J·mol^(-1)) showed a large difference with the value(135 k J·mol^(-1)) as reported from literatures and as determined by differential scanning calorimetry(DSC). The results of DSC/TG measurement showed that the phase transition of isonicotinic acid from 187.27 °C to277.47 °C underwent a sublimation process, with a sublimation enthalpy of 128.03 k J·mol^(-1). An efficient analytical technique combining pyrolysis and gas chromatography/mass spectrometry(Py-GC/MS) was used to prove this conclusion.展开更多
文摘The density functional theory （DFT） is the most popular method for evaluating bond dis- sociation enthalpies （BDEs） of most molecules. Thus, we are committed to looking for alternative methods that can balance the computational cost and higher precision to the best for large systems. The performance of DFT, double-hybrid DFT, and high-level com- posite methods are examined. The tested sets contain monocyclic and polycyclic aromatic molecules, branched hydrocarbons, small inorganic molecules, etc. The results show that the mPW2PLYP and G4MP2 methods achieve reasonable agreement with the benchmark val- ues for most tested molecules, and the mean absolute deviations are 2.43 and 1.96 kcal/mol after excluding the BDEs of branched hydrocarbons. We recommend the G4MP2 is the most appropriate method for small systems （atoms number≤20）; the double-hybrid DFT methods are advised for large aromatic molecules in medium size （20≤atoms number≤50）, and the double-hybrid DFT methods with empirical dispersion correction are recommended for long-chain and branched hydrocarbons in the same size scope; the DFT methods are ad- vised to apply for large systems （atoms number〉50）, and the M06-2X and B3P86 methods are also favorable. Moreover, the differences of optimized geometry of different methods are discussed and the effects of basis sets for various methods are investigated.
基金support from the Ministry of Education(MOE) Singapore Tier 1 (RG8/20)。
文摘A large database is desired for machine learning(ML) technology to make accurate predictions of materials physicochemical properties based on their molecular structure.When a large database is not available,the development of proper featurization method based on physicochemical nature of target proprieties can improve the predictive power of ML models with a smaller database.In this work,we show that two new featurization methods,volume occupation spatial matrix and heat contribution spatial matrix,can improve the accuracy in predicting energetic materials' crystal density(ρ_(crystal)) and solid phase enthalpy of formation(H_(f,solid)) using a database containing 451 energetic molecules.Their mean absolute errors are reduced from 0.048 g/cm~3 and 24.67 kcal/mol to 0.035 g/cm~3 and 9.66 kcal/mol,respectively.By leave-one-out-cross-validation,the newly developed ML models can be used to determine the performance of most kinds of energetic materials except cubanes.Our ML models are applied to predict ρ_(crystal) and H_(f,solid) of CHON-based molecules of the 150 million sized PubChem database,and screened out 56 candidates with competitive detonation performance and reasonable chemical structures.With further improvement in future,spatial matrices have the potential of becoming multifunctional ML simulation tools that could provide even better predictions in wider fields of materials science.
文摘The parameters of embedded atom method for elements Ce, Th and Yb were determined by fitting the lattice constants, the cohesive energy, the monovacancy formation energy and the bulk modulus of elements. The alloy potential was taken as the form of Johnson′s. The formation enthalpies of Th-Ce, Th-Yb and Ce-Yb binary alloys systems and Ce-Th-Yb ternary alloy were calculated with the present embedded atom potentials. The calculations for binary alloys are in good agreement with the results calculated with Miedema′s theory. As for the ternary alloy, the calculated formation enthalpies are in good agreement with those extrapolated from the formation enthalpies of constitutive binary alloys by Toop′s model.
基金This work is partly supported by The National Natural Science Foundation of China(59762001)Guangxi Natural Science Foundation(9811031)the Foundation of Ten-Hundred-Thousand Project of Personnel Department of Guangxi(2000209).
文摘The modified embedded atom method proposed by authors has been applied to calculating the enthalpies of formation of random alloys and the ordered intermetallic compounds for noble metal binary systems bearing Rh or Ir. The present results are in good agreement with those of Miedema theory, available experiments and the first-principles quantum mechanics calculations. The present results indicate that Cu-Rh, Cu-lr, Ag-Rh, Ag-lr, Au-Rh, Au-lr, Pd-Rh and Pd-lr systems are repulsive, however, IMi-Rh, Ni-lr, Pt-lr, Pt-Rh and Rh-lr systems form solid solutions and Ni-Rh, Ni-lr and Pt-Rh show ordering tendency.
文摘The purpose of this paper is to present a novel way to building quantitative structure-property relationship(QSPR) models for predicting the gas-to-benzene solvation enthalpy(ΔHSolv) of 158 organic compounds based on molecular descriptors calculated from the structure alone. Different kinds of descriptors were calculated for each compounds using dragon package. The variable selection technique of enhanced replacement method(ERM) was employed to select optimal subset of descriptors. Our investigation reveals that the dependence of physico-chemical properties on solvation enthalpy is a nonlinear observable fact and that ERM method is unable to model the solvation enthalpy accurately. The standard error value of prediction set for support vector machine(SVM) is 1.681 kJ ? mol^(-1) while it is 4.624 kJ ? mol^(-1) for ERM. The results established that the calculated ΔHSolvvalues by SVM were in good agreement with the experimental ones, and the performances of the SVM models were superior to those obtained by ERM one. This indicates that SVM can be used as an alternative modeling tool for QSPR studies.
基金the National Natural Science Foundation of China and the Chinese Academy of Sciences
文摘A numerical analysis was carried out to study the moving boundary problem in the physical process of pulsed Nd-YAG laser surface melting prior to vaporization.The enthalpy method was applied to solve this two-phase axisymmetrical mehing problem.Computational results of tempera ture fields were obtained,which provide useful information to practical lair treatment processing. The validity of enthalpy method in solving such problems is presented.
文摘In this work, the group contribution method of Chickos et al. was applied to estimate the fusion enthalpy of isonicotinic acid, and the obtained result(29.2 k J·mol^(-1)) showed a large difference with the value(135 k J·mol^(-1)) as reported from literatures and as determined by differential scanning calorimetry(DSC). The results of DSC/TG measurement showed that the phase transition of isonicotinic acid from 187.27 °C to277.47 °C underwent a sublimation process, with a sublimation enthalpy of 128.03 k J·mol^(-1). An efficient analytical technique combining pyrolysis and gas chromatography/mass spectrometry(Py-GC/MS) was used to prove this conclusion.
