The commonly used oxide-supported metal catalysts are usually prepared in aqueous phase,which then often need to undergo calcination before usage.Therefore,the surface hydration and dehydration of oxide supports are c...The commonly used oxide-supported metal catalysts are usually prepared in aqueous phase,which then often need to undergo calcination before usage.Therefore,the surface hydration and dehydration of oxide supports are critical for the realistic modeling of supported metal catalysts.In this work,by ab initio molecular dynamics(AIMD)simulations,the initial anhydrous monoclinic ZrO_(2)(111)surfaces are evaluated within explicit solvents in aqueous phase at mild temperatures.During the simulations,all the two-fold-coordinated O sites will soon be protonated to form the acidic hydroxyls(HO_(L)),remaining the basic hydroxyls(HO^(∗))on Zr.The basic hydroxyls(HO^(∗))can easily diffuse on surfaces via the active proton exchange with the undissociated adsorption water(H_(2)O^(∗)).Within the temperatures ranging from 273 K to 373 K,in aqueous phase a certain representative equilibrium hydrated m-ZrO_(2)(¯111)surface is obtained with the coverage(θ)of 0.75 on surface Zr atoms.Later,free energies on the stepwise surface water desorption are calculated by density functional theory to mimic the surface dehydration under the mild calcination temperatures lower than 800 K.By obtaining the phase diagrams of surface dehydration,the representative partially hydrated m-ZrO_(2)(111)surfaces(0.25≤θ<0.75)at various calcination temperatures are illustrated.These hydrated m-ZrO_(2)(111)surfaces can be crucial and readily applied for more realistic modeling of ZrO_(2) catalysts and ZrO_(2)-supported metal catalysts.展开更多
The mixing enthalpies and structural order in liquid Mg−Si system were investigated via ab-initio molecular dynamics at 1773 K.By calculating the transferred charges and electron density differences,the dominance of S...The mixing enthalpies and structural order in liquid Mg−Si system were investigated via ab-initio molecular dynamics at 1773 K.By calculating the transferred charges and electron density differences,the dominance of Si−Si interactions in the chemical environments around Si was demonstrated,which determined that the mixing enthalpy reached the minimum on Mg-rich side.In terms of Honeycutt and Anderson(HA)bond pairs based on the partial pair correlation functions,the attraction between Si−Si pairs and Mg atoms was revealed,and the evolution of structural order with Si content was characterized as a process of constituting frame structures by Si−Si pairs that dispersed Mg atoms.Focusing on tetrahedral order of local Si-configurations,a correlation between the mixing enthalpy and structural order was uncovered ultimately,which provided a new perspective combining the energetics with geometry to understand the liquid Mg−Si binary system.展开更多
We study the binding of molecular oxygen to a (5, 0) single walled SiC nanotube, by means of density functional calculations. The center of a hexagon of silicon and carbon atoms in sites on SiCNT surfaces is the mos...We study the binding of molecular oxygen to a (5, 0) single walled SiC nanotube, by means of density functional calculations. The center of a hexagon of silicon and carbon atoms in sites on SiCNT surfaces is the most stable adsorption site for 02 molecule, with a binding energy of -38.22 eV and an average Si-O binding distance of 1.698 A. We have also tested the stability of the 02-adsorbed SiCNT/CNT with ab initio molecular dynamics simulation which have been carried out at room temperature. Furthermore, the adsorption of 02 on the single walled carbon nanotubes has been investigated. Our first-principles calculations predict that the 02 adsorptive capability of silicon carbide nanotubes is much better than that of carbon nanotubes. This might have potential for gas detection and energy storage.展开更多
The irradiation-induced sputtering and the structural damage at tungsten surface are investigated by using molecular dynamics simulations at the level of quantum mechanics. Our simulations indicate that the sputtered ...The irradiation-induced sputtering and the structural damage at tungsten surface are investigated by using molecular dynamics simulations at the level of quantum mechanics. Our simulations indicate that the sputtered atoms appear when the energy of incident primary knock-on atom (PKA) is more than 200 eV and the incident angle of the PKA is larger than 65°. Meanwhile, the irradiation-induced vacancies are less when the incident angle of PKA is in the range of 45°-65°. So, the optimum incident angles of PKA are suggested to reduce the irradiation-induced damage of the W surface. Furthermore, we find that the interstitials contained in the systems accelerate the sputtering whereas the intrinsic vacancies suppress the sputtering when the PKA is near the defects.展开更多
We investigate the structure,trapping,and diffusion behaviors of helium(He) at vacancy in a Fe single crystal using first-principles simulations.Vacancy with more space can provide the lower electron density region fo...We investigate the structure,trapping,and diffusion behaviors of helium(He) at vacancy in a Fe single crystal using first-principles simulations.Vacancy with more space can provide the lower electron density region for He binding in comparison with intrinsic Fe,causing He to diffuse into the vacancy inner easily.We provide the quantitative microscopic studies related to the atomic-level thermo-kinetic trapping processes.Moreover,such physical viewpoint can be applied to other vacancy-like defects such as vacancy clusters,void and grain boundaries which can open a space with reduced electron density region to increase He binding in metals and metal alloys.展开更多
Despite their rich information content,electronic structure data amassed at high volumes in ab initio molecular dynamics simulations are generally under-utilized.We introduce a transferable high-fidelity neural networ...Despite their rich information content,electronic structure data amassed at high volumes in ab initio molecular dynamics simulations are generally under-utilized.We introduce a transferable high-fidelity neural network representation of such data in the form of tight-binding Hamiltonians for crystalline materials.This predictive representation of ab initio electronic structure,combined with machinelearning boosted molecular dynamics,enables efficient and accurate electronic evolution and sampling.When it is applied to a one-dimension charge-density wave material,carbyne,we are able to compute the spectral function and optical conductivity in the canonical ensemble.The spectral functions evaluated during soliton-antisoliton pair annihilation process reveal significant renormalization of low-energy edge modes due to retarded electron-lattice coupling beyond the Born-Oppenheimer limit.The availability of an efficient and reusable surrogate model for the electronic structure dynamical system will enable calculating many interesting physical properties,paving the way to previously inaccessible or challenging avenues in materials modeling.展开更多
Surface eigenstress and eigendisplacement models were used to investigate the surface stress, surface relaxation and surface elasticity of thin films with different surface orientations. Molecular dynamics simulations...Surface eigenstress and eigendisplacement models were used to investigate the surface stress, surface relaxation and surface elasticity of thin films with different surface orientations. Molecular dynamics simulations and first-principles calculations were conducted on face-centered cubic Au films with the focus on relaxation induced nonlinear initial deformation. The simu- lation results verify the theoretical predictions of the size dependency of surface energy density and surface stress, and the non- linear scaling law of the size-dependent Young's modulus of thin films. The mechanism of the size-dependent behaviors was further explored at the atomic bonding level with the charge density field. The Au atomic bonding at surfaces is enhanced compared to its interior counterpart and therefore the nominal Young's modulus of the Au thin films is larger when the film thickness is smaller.surface elasticity,展开更多
A laboratory-based aerobic incubation was conducted to investigate nitrogen (N) isotopic fractionation related to nitrification in five agricultural soils after application of ammonium sulfate ((NH4)2804). The s...A laboratory-based aerobic incubation was conducted to investigate nitrogen (N) isotopic fractionation related to nitrification in five agricultural soils after application of ammonium sulfate ((NH4)2804). The soil samples were collected from a subtropical barren land soil derived from granite (RGB), three subtropical upland soils derived from granite (RQU), Quaternary red earth (RGU), Quaternary Xiashu loess (YQU) and a temperate upland soil generated from alluvial deposit (FAU). The five soils varied in nitrification potential, being in the order of FAU 〉 YQU 〉 RGU 〉 RQU 〉 RGB. Significant N isotopic fractionation accompanied nitrification of NH4+. 615N values of NH4+ increased with enhanced nitrification over time in the four upland soils with NH4+ addition, while those of NO3 decreased consistently to the minimum and thereafter increased. 515N values of NH4+ showed a significantly negative linear relationship with NH4+-N concentration, but a positive linear relationship with NO3-N concentration. The apparent isotopic fractionation factor calculated based on the loss of NH4+ was 1.036 for RQU, 1.022 for RGU, 1.016 for YQU, and 1.020 for FAU, respectively. Zero- and first-order reaction kinetics seemed to have their limitations in describing the nitrification process affected by NH4+ input in the studied soils. In contrast, N kinetic isotope fractionation was closely related to the nitrifying activity, and might serve as an alternative tool for estimating the nitrification capacity of agricultural soils.展开更多
First-principles calculations,which are based on the plane-wave pseudopotential approach to density functional perturbation theory within the local density approximation,have been performed to investigate the structur...First-principles calculations,which are based on the plane-wave pseudopotential approach to density functional perturbation theory within the local density approximation,have been performed to investigate the structural,lattice dynamical and thermodynamic properties of zinc blende(B3) structure magnesium chalcogenides:MgS,MgSe and MgTe.The results of ground state parameters and phonon dispersion are compared and agree well with the experimental data available and other calculations.We obtain the change of Born effective charge and LO-TO splitting under hydrostatic pressure.Finally,by the calculations of phonon frequencies,some thermodynamic properties such as the entropy,heat capacity,internal energy,and free energy are also successfully obtained.展开更多
基金This work was supported by the National Natural Science Foundation of China(No.22022504,No.22003022)of ChinaNatural Science Foundation of Guangdong,China(No.2021A1515010213,No.2021A1515110406)+2 种基金Guangdong“Pearl River”Talent Plan(No.2019QN01L353)Higher Education Innovation Strong School Project of Guangdong Province of China(No.2020KTSCX122)Guangdong Provincial Key Laboratory of Catalysis(No.2020B121201002).Most calculations are performed on the CHEM Highperformance Computing Cluster(CHEM-HPC)located at the Department of Chemistry,Southern University of Science and Technology(SUSTech).The computational resources are also supported by the Center for Computational Science and Engineering at SUSTech.
文摘The commonly used oxide-supported metal catalysts are usually prepared in aqueous phase,which then often need to undergo calcination before usage.Therefore,the surface hydration and dehydration of oxide supports are critical for the realistic modeling of supported metal catalysts.In this work,by ab initio molecular dynamics(AIMD)simulations,the initial anhydrous monoclinic ZrO_(2)(111)surfaces are evaluated within explicit solvents in aqueous phase at mild temperatures.During the simulations,all the two-fold-coordinated O sites will soon be protonated to form the acidic hydroxyls(HO_(L)),remaining the basic hydroxyls(HO^(∗))on Zr.The basic hydroxyls(HO^(∗))can easily diffuse on surfaces via the active proton exchange with the undissociated adsorption water(H_(2)O^(∗)).Within the temperatures ranging from 273 K to 373 K,in aqueous phase a certain representative equilibrium hydrated m-ZrO_(2)(¯111)surface is obtained with the coverage(θ)of 0.75 on surface Zr atoms.Later,free energies on the stepwise surface water desorption are calculated by density functional theory to mimic the surface dehydration under the mild calcination temperatures lower than 800 K.By obtaining the phase diagrams of surface dehydration,the representative partially hydrated m-ZrO_(2)(111)surfaces(0.25≤θ<0.75)at various calcination temperatures are illustrated.These hydrated m-ZrO_(2)(111)surfaces can be crucial and readily applied for more realistic modeling of ZrO_(2) catalysts and ZrO_(2)-supported metal catalysts.
基金The authors are grateful for the financial supports from the National Key Research and Development Program of China(2016YFB0701202)the National Natural Science Foundation of China(51901117,51801116,51804190,and 11804179)the Shandong Provincial Key Research and Development Plan,China(2019GGX102047).
文摘The mixing enthalpies and structural order in liquid Mg−Si system were investigated via ab-initio molecular dynamics at 1773 K.By calculating the transferred charges and electron density differences,the dominance of Si−Si interactions in the chemical environments around Si was demonstrated,which determined that the mixing enthalpy reached the minimum on Mg-rich side.In terms of Honeycutt and Anderson(HA)bond pairs based on the partial pair correlation functions,the attraction between Si−Si pairs and Mg atoms was revealed,and the evolution of structural order with Si content was characterized as a process of constituting frame structures by Si−Si pairs that dispersed Mg atoms.Focusing on tetrahedral order of local Si-configurations,a correlation between the mixing enthalpy and structural order was uncovered ultimately,which provided a new perspective combining the energetics with geometry to understand the liquid Mg−Si binary system.
文摘We study the binding of molecular oxygen to a (5, 0) single walled SiC nanotube, by means of density functional calculations. The center of a hexagon of silicon and carbon atoms in sites on SiCNT surfaces is the most stable adsorption site for 02 molecule, with a binding energy of -38.22 eV and an average Si-O binding distance of 1.698 A. We have also tested the stability of the 02-adsorbed SiCNT/CNT with ab initio molecular dynamics simulation which have been carried out at room temperature. Furthermore, the adsorption of 02 on the single walled carbon nanotubes has been investigated. Our first-principles calculations predict that the 02 adsorptive capability of silicon carbide nanotubes is much better than that of carbon nanotubes. This might have potential for gas detection and energy storage.
基金This work is supported by the National Magnetic Confinement Fusion Program (No.2013GB107004), the National Natural Science Foundation of China (No.11275191). The Computational Center of USTC is acknowledged for computational support.
文摘The irradiation-induced sputtering and the structural damage at tungsten surface are investigated by using molecular dynamics simulations at the level of quantum mechanics. Our simulations indicate that the sputtered atoms appear when the energy of incident primary knock-on atom (PKA) is more than 200 eV and the incident angle of the PKA is larger than 65°. Meanwhile, the irradiation-induced vacancies are less when the incident angle of PKA is in the range of 45°-65°. So, the optimum incident angles of PKA are suggested to reduce the irradiation-induced damage of the W surface. Furthermore, we find that the interstitials contained in the systems accelerate the sputtering whereas the intrinsic vacancies suppress the sputtering when the PKA is near the defects.
基金supported by the National Natural Science Foundation of China (Grant No. 51101135)the Natural Science Foundation of Shandong (Grant No. ZR2010AM002)
文摘We investigate the structure,trapping,and diffusion behaviors of helium(He) at vacancy in a Fe single crystal using first-principles simulations.Vacancy with more space can provide the lower electron density region for He binding in comparison with intrinsic Fe,causing He to diffuse into the vacancy inner easily.We provide the quantitative microscopic studies related to the atomic-level thermo-kinetic trapping processes.Moreover,such physical viewpoint can be applied to other vacancy-like defects such as vacancy clusters,void and grain boundaries which can open a space with reduced electron density region to increase He binding in metals and metal alloys.
基金supported by the National Natural Science Foundation of China(11725415 and 11934001)the Ministry of Science and Technology of China(2018YFA0305601 and2016YFA0301004)+1 种基金by the Strategic Priority Research Program of the Chinese Academy of Sciences(XDB28000000)supported in part by the Center for Chemistry in Solution and at Interfaces(CSI)at Princeton University,funded by the DOE Award DE-SC0019394。
文摘Despite their rich information content,electronic structure data amassed at high volumes in ab initio molecular dynamics simulations are generally under-utilized.We introduce a transferable high-fidelity neural network representation of such data in the form of tight-binding Hamiltonians for crystalline materials.This predictive representation of ab initio electronic structure,combined with machinelearning boosted molecular dynamics,enables efficient and accurate electronic evolution and sampling.When it is applied to a one-dimension charge-density wave material,carbyne,we are able to compute the spectral function and optical conductivity in the canonical ensemble.The spectral functions evaluated during soliton-antisoliton pair annihilation process reveal significant renormalization of low-energy edge modes due to retarded electron-lattice coupling beyond the Born-Oppenheimer limit.The availability of an efficient and reusable surrogate model for the electronic structure dynamical system will enable calculating many interesting physical properties,paving the way to previously inaccessible or challenging avenues in materials modeling.
基金supported by the Hong Kong Research Grants Council(Grant No.622312)
文摘Surface eigenstress and eigendisplacement models were used to investigate the surface stress, surface relaxation and surface elasticity of thin films with different surface orientations. Molecular dynamics simulations and first-principles calculations were conducted on face-centered cubic Au films with the focus on relaxation induced nonlinear initial deformation. The simu- lation results verify the theoretical predictions of the size dependency of surface energy density and surface stress, and the non- linear scaling law of the size-dependent Young's modulus of thin films. The mechanism of the size-dependent behaviors was further explored at the atomic bonding level with the charge density field. The Au atomic bonding at surfaces is enhanced compared to its interior counterpart and therefore the nominal Young's modulus of the Au thin films is larger when the film thickness is smaller.surface elasticity,
基金Supported by the Natural Science Foundation of Jiangsu Province,China(No.BK2010612)the Foundation of State Key Laboratory of Soil and Sustainable Agriculture(No.Y052010034)the Knowledge Innovation Program of the Institute of Soil Science,Chinese Academy of Sciences(No.ISSASIP0723)
文摘A laboratory-based aerobic incubation was conducted to investigate nitrogen (N) isotopic fractionation related to nitrification in five agricultural soils after application of ammonium sulfate ((NH4)2804). The soil samples were collected from a subtropical barren land soil derived from granite (RGB), three subtropical upland soils derived from granite (RQU), Quaternary red earth (RGU), Quaternary Xiashu loess (YQU) and a temperate upland soil generated from alluvial deposit (FAU). The five soils varied in nitrification potential, being in the order of FAU 〉 YQU 〉 RGU 〉 RQU 〉 RGB. Significant N isotopic fractionation accompanied nitrification of NH4+. 615N values of NH4+ increased with enhanced nitrification over time in the four upland soils with NH4+ addition, while those of NO3 decreased consistently to the minimum and thereafter increased. 515N values of NH4+ showed a significantly negative linear relationship with NH4+-N concentration, but a positive linear relationship with NO3-N concentration. The apparent isotopic fractionation factor calculated based on the loss of NH4+ was 1.036 for RQU, 1.022 for RGU, 1.016 for YQU, and 1.020 for FAU, respectively. Zero- and first-order reaction kinetics seemed to have their limitations in describing the nitrification process affected by NH4+ input in the studied soils. In contrast, N kinetic isotope fractionation was closely related to the nitrifying activity, and might serve as an alternative tool for estimating the nitrification capacity of agricultural soils.
基金Supported by Education Department Foundation of Liaoning Province of China under Grant Nos. 201064145,2010397Education Science Foundation of Liaoning of China under Grant No. 201102166
文摘First-principles calculations,which are based on the plane-wave pseudopotential approach to density functional perturbation theory within the local density approximation,have been performed to investigate the structural,lattice dynamical and thermodynamic properties of zinc blende(B3) structure magnesium chalcogenides:MgS,MgSe and MgTe.The results of ground state parameters and phonon dispersion are compared and agree well with the experimental data available and other calculations.We obtain the change of Born effective charge and LO-TO splitting under hydrostatic pressure.Finally,by the calculations of phonon frequencies,some thermodynamic properties such as the entropy,heat capacity,internal energy,and free energy are also successfully obtained.