In order to simulate changes in the water quality of the Miyun Reservoir dueto continuous descent of surface water level, a 3-D ecological hydrodynamic model was developedthrough coupling the water quality analysis si...In order to simulate changes in the water quality of the Miyun Reservoir dueto continuous descent of surface water level, a 3-D ecological hydrodynamic model was developedthrough coupling the water quality analysis simulation program (WASP) with the environmental fluiddynamics code (EFDC). The model was then calibrated and verified. Four scenarios (S1, S21, S22 andS23) were simulated using the model. Results show that the water quality of the Miyun Reservoirunder conditions of low surface water level is apparently affected by different amounts of inflowand different total phosphorus (TP) loadings. The chlorophyll-a concentration might exceed 10 μg/Lin many areas of the Miyun Reservoir (This limitative value is seen as a critical value ofeutrophication) when large loadings of TP enter due to the amount of inflow increasing. Results ofscenario S23 indicate that control of TP loadings can decrease chlorophyll-a concentrationeffectively, and the water quality of the Miyun Reservoir will improve or retain its status quo.展开更多
Kinetics of dissociative O2 adsorption, OHad desorption, and oxygen reduction reaction (ORR) at Pt(111) electrode in 0.1 mol/L HClO4 has been investigated. Reversible OHad adsorption/desorption occurs at potential...Kinetics of dissociative O2 adsorption, OHad desorption, and oxygen reduction reaction (ORR) at Pt(111) electrode in 0.1 mol/L HClO4 has been investigated. Reversible OHad adsorption/desorption occurs at potentials from 0.6 V to 1.0 V (vs. RHE) with the exchange current density of ca. 50 mA/cm^2 at 0.8 V, the fast kinetics of OHad desorption indicates that it should not be the rate determining step for ORR. In the kineticor kinetic-mass transport mix controlled potential region, ORR current at constant potential displays slight decrease with reaction time. ORR current in the positive-going potential scan is slightly larger than that in the subsequent negative-going scan with electrode rotation speed (〉800 r/min) and slow potential scan rate (〈100 mV/s). The open circuit potential of Pt/0.1 mol/L HClO4 interface increases promptly from 0.9 V to 1.0 V after switch from O2 free- to O2-saturated solution. The increase of open circuit potential as well as ORR current decays under potential control due to the accumulation of OHad from dissociative adsorption of O2. It indicates that at Pt(111) the net rate for O2 decomposition to OHad is slightly faster than that for OHad removal, one cannot simply use the assumption of rate determining step to discuss ORR kinetics. Instead, the ORR kinetics is determined by both the kinetics for O2 decomposition to OHad as well as the thermo-equilibrium of OHad+H^++e→←H2O.展开更多
Exploring the production and application of clean energy has always been the core of sustainable development.As a clean and sustainable technology,electrocatalysis has been receiving widespread attention.It is crucial...Exploring the production and application of clean energy has always been the core of sustainable development.As a clean and sustainable technology,electrocatalysis has been receiving widespread attention.It is crucial to achieve efficient,stable and cheap electrocatalysts.However,the traditional“trial and error”method is time-consuming,laborious and costly.In recent years,with the significant increase in computing power,computations have played an important role in electrocatalyst design.Nevertheless,it is still difficult to search for advanced electrocatalysts in the vast chemical space through traditional density functional theory(DFT)computations.Fortunately,the development of machine learning and interdisciplinary integration will inject new impetus into targeted design of electrocatalysts.Machine learning is able to predict electrochemical performances with an accuracy close to DFT.Here we provide an overview of the application of machine learning in electrocatalyst design,including the prediction of structure,thermodynamic properties and kinetic barriers.We also discuss the potential of explicit solvent model combined with machine learning molecular dynamics in this field.Finally,the favorable circumstances and challenges are outlined for the future development of machine learning in electrocatalysis.The studies on electrochemical processes by machine learning will further realize targeted design of high-efficiency electrocatalysts.展开更多
We further study the validity of the Monte Carlo Hamiltonian method. The advantage of the method, in comparison with the standard Monte Carlo Lagrangian approach, is its capability to study the excited states. We cons...We further study the validity of the Monte Carlo Hamiltonian method. The advantage of the method, in comparison with the standard Monte Carlo Lagrangian approach, is its capability to study the excited states. We consider two quantum mechanical models: a symmetric one ; and an asymmetric one , for and , for . The results for the spectrum, wave functions and thermodynamical observables are in agreement with the analytical or Runge–Kutta calculations.展开更多
The influence of the confining potential strength and temperature on the structures and dynamics of a two-dimensional (2D) dusty plasma system is investigated through molecular dynamic (MD) simulation. The circular sy...The influence of the confining potential strength and temperature on the structures and dynamics of a two-dimensional (2D) dusty plasma system is investigated through molecular dynamic (MD) simulation. The circular symmetric confining potential leads to the nonuniform packing of particles, that is, an inner core with a hexagon lattice surrounded by a few outer circular shells. Under the appropriate confining potential and temperature, the particle trajectories on middle shells form a series of concentric and nested hexagons due to tangential movements of particles.Mean square displacement, self-diffusion constant, pair correlation function, and the nearest bond are used to characterize the structural and dynamical properties of the system. With the increase of the confining potential, the radial and tangential movements of particles have different behaviors. With the increase of temperature, the radial and tangential motions strengthen, particle trajectories gradually become disordered, and the system gradually changes from a crystal or liquid state to a gas state.展开更多
Photodissociation of HOBr is an important step in the reaction network of the depletion of ozone in stratosphere.Here,we report the first three-dimensional potential energy surfaces for the lowest three singlet states...Photodissociation of HOBr is an important step in the reaction network of the depletion of ozone in stratosphere.Here,we report the first three-dimensional potential energy surfaces for the lowest three singlet states for HOBr,based on high level multi reference configuration interaction calculations.Quantum dynamics calculations are performed with a real wavepacket method,yielding not only absorption spectra but also internal state and angular distributions of the photodissociation fragments.Our results agree quantitatively with the measured total absorption cross sections of HOBr in the ultraviolet region and reproduce well the observed vibrationally cold and rotationally hot OH/OD fragments via photodissociation of HOBr/DOBr at 266 nm.In addition,we predict that the recoil anisotropy parameters for OH/OD are close to the limiting value of a parallel transition,suggesting a rapid dissociation process at 266 nm following an in-plane transition from the ground state(1^1A')to the 21A'state.This is consistent with the experimental conclusion derived from the measured rotational alignment.However,spin and electronic angular momenta need to be taken into account in the future to achieve a more quantitative agreement with experiment.Our work is expected to motivate further experimental investigations for this benchmark system.展开更多
The reaction H+SO_(2)→OH+SO is important in the combustion and atmospheric chemistry,as well as the interstellar medium.It also represents a typical complex-forming reaction with deep complexes,serving as an ideal ca...The reaction H+SO_(2)→OH+SO is important in the combustion and atmospheric chemistry,as well as the interstellar medium.It also represents a typical complex-forming reaction with deep complexes,serving as an ideal candidate for testing various kinetics theories and providing interesting reaction dynamical phenomena.In this work,we reported a quasiclassical trajectory study of this reaction on our previously developed accurate full-dimensional potential energy surface.The experimental thermal rate coefficients over the temperature range 1400 K≤T≤2200 K were well reproduced.For the reactant SO_(2)being sampled at the ground ro-vibrational state,the calculated integral cross sections increased slightly along the collision energy ranging from 31.0 kcal/mol to 40.0 kcal/mol,and then became essentially flat at the collision energy within 40.0−55.0 kcal/mol.The product angular distributions are almost symmetric with nearly identical backward-forward double peak structure.The products OH and SO vibrational state distributions were also analyzed.展开更多
The quark potential model is extended to include the sea quark excitation using the random phase approximation. The effective quark interaction preserves the important QCD properties — chiral symmetry and confinement...The quark potential model is extended to include the sea quark excitation using the random phase approximation. The effective quark interaction preserves the important QCD properties — chiral symmetry and confinement simultaneously. A primary qualitative analysis shows that the π meson as a well-known typical Goldstone boson and the other mesons made up of valence quark pair such as the ρ meson can also be described in this extended quark potential model.展开更多
Unidirectional transport of a particle in a spatially periodic and symmetric potential under a periodic force with broken temporal symmetry is studied. With a collaboration of the potential field and the asymmetric ac...Unidirectional transport of a particle in a spatially periodic and symmetric potential under a periodic force with broken temporal symmetry is studied. With a collaboration of the potential field and the asymmetric ac force, a dc current can be observed. Resonant current steps are found for a finite period of the ac force. A phase diagram of these resonant steps is given. Stochastic-resonance-like directional transport induced by thermal noises is revealed.展开更多
Using coarse-grained molecular dynamics simulations based on Gay-Berne potential model, we have simulated the cooling process of liquid n-butanol. A new set of GB parameters are obtained by fitting the results of dens...Using coarse-grained molecular dynamics simulations based on Gay-Berne potential model, we have simulated the cooling process of liquid n-butanol. A new set of GB parameters are obtained by fitting the results of density functional theory calculations. The simulations are carried out in the range of 290-50 K with temperature decrements of 10 K. The cooling characteristics are determined on the basis of the variations of the density, the potential energy and orientational order parameter with temperature, whose slopes all show discontinuity. Both the radial distribution function curves and the second-rank orientational correlation function curves exhibit splitting in the second peak. Using the discontinuous change of these thermodynamic and structure properties, we obtain the glass transition at an estimate of temperature Tg=1204.10 K, which is in good agreement with experimental results 1104-1 K.展开更多
In this paper, we give a direct method for calculating the partition function, and hence the equation of state (EOS) of QCD at finite chemical potential and zero temperature. In the EOS derived in this paper the pre...In this paper, we give a direct method for calculating the partition function, and hence the equation of state (EOS) of QCD at finite chemical potential and zero temperature. In the EOS derived in this paper the pressure density is the sum of two terms: the first term P(μ)|μ=0 (the pressure density at μ = 0) is a μ-independent constant; the second term, which is totally determined by G[μ] (p) (the dressed quark propagator at finite μ), contains all the nontrivial μ-dependence. By applying a general result in the rainbow-ladder approximation of the Dyson-Schwinger approach obtained in our previous study [Phys. Rev. C 71 (2005) 015205], G[μ](p) is calculated from the meromorphic quark propagator proposed in [Phys. Rev. D 67 (2003) 054019]. From this the full analytic expression of the EOS of QCD at finite μ and zero T is obtained (apart from the constant term P(μ)|μ=0, which can in principle be caJculated from the CJT effective action). A comparison between our EOS and the cold, perturbative EOS of QCD of Fraga, Pisarski and Schaffner-Bielich is made. It is expected that our EOS can provide a possible new approach for the study of neutron stars.展开更多
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.展开更多
Complex-forming reactions widely exist in gas-phase chemical reactions. Various complexforming bimolecular reactions have been investigated and interesting phenomena have been discovered. The complex-forming reactions...Complex-forming reactions widely exist in gas-phase chemical reactions. Various complexforming bimolecular reactions have been investigated and interesting phenomena have been discovered. The complex-forming reactions usually have small or no barrier in the entrance channel, which leads to obvious differences in kinetic and dynamic characteristics compared with direct reactions. Theoretically, quantum state-resolved reaction dynamics can provide the most detailed microscopic dynamic mechanisms and is now feasible for a direct reaction with only one potential barrier. However, it is of great challenge to construct accurate potential energy surfaces and perform accurate quantum dynamics calculations for a complex polyatomic reaction involving deep potential wells and multi-channels. This paper reviews the most recent progress in two prototypical oxyhydrogen complex-forming reaction systems, HO2 and HO3, which are significant in combustion, atmospheric, and interstellar chemistry. We will present a brief survey of both computational and experimental work and emphasize on some unsolved problems existing in these systems.展开更多
文摘In order to simulate changes in the water quality of the Miyun Reservoir dueto continuous descent of surface water level, a 3-D ecological hydrodynamic model was developedthrough coupling the water quality analysis simulation program (WASP) with the environmental fluiddynamics code (EFDC). The model was then calibrated and verified. Four scenarios (S1, S21, S22 andS23) were simulated using the model. Results show that the water quality of the Miyun Reservoirunder conditions of low surface water level is apparently affected by different amounts of inflowand different total phosphorus (TP) loadings. The chlorophyll-a concentration might exceed 10 μg/Lin many areas of the Miyun Reservoir (This limitative value is seen as a critical value ofeutrophication) when large loadings of TP enter due to the amount of inflow increasing. Results ofscenario S23 indicate that control of TP loadings can decrease chlorophyll-a concentrationeffectively, and the water quality of the Miyun Reservoir will improve or retain its status quo.
基金V. ACKNOWLEDGMENTS This work was supported by the National Natural Science Foundation of China (No.20773116), the National Instrumentation Program (No.2011YQ03012416), and 973 Program from the Ministry of Science and Technology of China (No.2010CB923302).
文摘Kinetics of dissociative O2 adsorption, OHad desorption, and oxygen reduction reaction (ORR) at Pt(111) electrode in 0.1 mol/L HClO4 has been investigated. Reversible OHad adsorption/desorption occurs at potentials from 0.6 V to 1.0 V (vs. RHE) with the exchange current density of ca. 50 mA/cm^2 at 0.8 V, the fast kinetics of OHad desorption indicates that it should not be the rate determining step for ORR. In the kineticor kinetic-mass transport mix controlled potential region, ORR current at constant potential displays slight decrease with reaction time. ORR current in the positive-going potential scan is slightly larger than that in the subsequent negative-going scan with electrode rotation speed (〉800 r/min) and slow potential scan rate (〈100 mV/s). The open circuit potential of Pt/0.1 mol/L HClO4 interface increases promptly from 0.9 V to 1.0 V after switch from O2 free- to O2-saturated solution. The increase of open circuit potential as well as ORR current decays under potential control due to the accumulation of OHad from dissociative adsorption of O2. It indicates that at Pt(111) the net rate for O2 decomposition to OHad is slightly faster than that for OHad removal, one cannot simply use the assumption of rate determining step to discuss ORR kinetics. Instead, the ORR kinetics is determined by both the kinetics for O2 decomposition to OHad as well as the thermo-equilibrium of OHad+H^++e→←H2O.
文摘Exploring the production and application of clean energy has always been the core of sustainable development.As a clean and sustainable technology,electrocatalysis has been receiving widespread attention.It is crucial to achieve efficient,stable and cheap electrocatalysts.However,the traditional“trial and error”method is time-consuming,laborious and costly.In recent years,with the significant increase in computing power,computations have played an important role in electrocatalyst design.Nevertheless,it is still difficult to search for advanced electrocatalysts in the vast chemical space through traditional density functional theory(DFT)computations.Fortunately,the development of machine learning and interdisciplinary integration will inject new impetus into targeted design of electrocatalysts.Machine learning is able to predict electrochemical performances with an accuracy close to DFT.Here we provide an overview of the application of machine learning in electrocatalyst design,including the prediction of structure,thermodynamic properties and kinetic barriers.We also discuss the potential of explicit solvent model combined with machine learning molecular dynamics in this field.Finally,the favorable circumstances and challenges are outlined for the future development of machine learning in electrocatalysis.The studies on electrochemical processes by machine learning will further realize targeted design of high-efficiency electrocatalysts.
文摘We further study the validity of the Monte Carlo Hamiltonian method. The advantage of the method, in comparison with the standard Monte Carlo Lagrangian approach, is its capability to study the excited states. We consider two quantum mechanical models: a symmetric one ; and an asymmetric one , for and , for . The results for the spectrum, wave functions and thermodynamical observables are in agreement with the analytical or Runge–Kutta calculations.
基金the National Natural Science Foundation of China under
文摘The influence of the confining potential strength and temperature on the structures and dynamics of a two-dimensional (2D) dusty plasma system is investigated through molecular dynamic (MD) simulation. The circular symmetric confining potential leads to the nonuniform packing of particles, that is, an inner core with a hexagon lattice surrounded by a few outer circular shells. Under the appropriate confining potential and temperature, the particle trajectories on middle shells form a series of concentric and nested hexagons due to tangential movements of particles.Mean square displacement, self-diffusion constant, pair correlation function, and the nearest bond are used to characterize the structural and dynamical properties of the system. With the increase of the confining potential, the radial and tangential movements of particles have different behaviors. With the increase of temperature, the radial and tangential motions strengthen, particle trajectories gradually become disordered, and the system gradually changes from a crystal or liquid state to a gas state.
基金supported by the National Key R&D Program of China (2017YFA0303500)Anhui Initiative in Quantum Information Technologies(AHY090200).
文摘Photodissociation of HOBr is an important step in the reaction network of the depletion of ozone in stratosphere.Here,we report the first three-dimensional potential energy surfaces for the lowest three singlet states for HOBr,based on high level multi reference configuration interaction calculations.Quantum dynamics calculations are performed with a real wavepacket method,yielding not only absorption spectra but also internal state and angular distributions of the photodissociation fragments.Our results agree quantitatively with the measured total absorption cross sections of HOBr in the ultraviolet region and reproduce well the observed vibrationally cold and rotationally hot OH/OD fragments via photodissociation of HOBr/DOBr at 266 nm.In addition,we predict that the recoil anisotropy parameters for OH/OD are close to the limiting value of a parallel transition,suggesting a rapid dissociation process at 266 nm following an in-plane transition from the ground state(1^1A')to the 21A'state.This is consistent with the experimental conclusion derived from the measured rotational alignment.However,spin and electronic angular momenta need to be taken into account in the future to achieve a more quantitative agreement with experiment.Our work is expected to motivate further experimental investigations for this benchmark system.
基金supported by the National Natural Science Foundation of China(No.21973009)Chongqing Municipal Natural Science Foundation(No.cstc2019jcyj-msxm X0087)。
文摘The reaction H+SO_(2)→OH+SO is important in the combustion and atmospheric chemistry,as well as the interstellar medium.It also represents a typical complex-forming reaction with deep complexes,serving as an ideal candidate for testing various kinetics theories and providing interesting reaction dynamical phenomena.In this work,we reported a quasiclassical trajectory study of this reaction on our previously developed accurate full-dimensional potential energy surface.The experimental thermal rate coefficients over the temperature range 1400 K≤T≤2200 K were well reproduced.For the reactant SO_(2)being sampled at the ground ro-vibrational state,the calculated integral cross sections increased slightly along the collision energy ranging from 31.0 kcal/mol to 40.0 kcal/mol,and then became essentially flat at the collision energy within 40.0−55.0 kcal/mol.The product angular distributions are almost symmetric with nearly identical backward-forward double peak structure.The products OH and SO vibrational state distributions were also analyzed.
文摘The quark potential model is extended to include the sea quark excitation using the random phase approximation. The effective quark interaction preserves the important QCD properties — chiral symmetry and confinement simultaneously. A primary qualitative analysis shows that the π meson as a well-known typical Goldstone boson and the other mesons made up of valence quark pair such as the ρ meson can also be described in this extended quark potential model.
文摘Unidirectional transport of a particle in a spatially periodic and symmetric potential under a periodic force with broken temporal symmetry is studied. With a collaboration of the potential field and the asymmetric ac force, a dc current can be observed. Resonant current steps are found for a finite period of the ac force. A phase diagram of these resonant steps is given. Stochastic-resonance-like directional transport induced by thermal noises is revealed.
文摘Using coarse-grained molecular dynamics simulations based on Gay-Berne potential model, we have simulated the cooling process of liquid n-butanol. A new set of GB parameters are obtained by fitting the results of density functional theory calculations. The simulations are carried out in the range of 290-50 K with temperature decrements of 10 K. The cooling characteristics are determined on the basis of the variations of the density, the potential energy and orientational order parameter with temperature, whose slopes all show discontinuity. Both the radial distribution function curves and the second-rank orientational correlation function curves exhibit splitting in the second peak. Using the discontinuous change of these thermodynamic and structure properties, we obtain the glass transition at an estimate of temperature Tg=1204.10 K, which is in good agreement with experimental results 1104-1 K.
基金supported in part by the National Natural Science Foundation of China under Grant No.10575050the Research Fund for the Doctoral Program of Higher Education under Grant No.20060284020
文摘In this paper, we give a direct method for calculating the partition function, and hence the equation of state (EOS) of QCD at finite chemical potential and zero temperature. In the EOS derived in this paper the pressure density is the sum of two terms: the first term P(μ)|μ=0 (the pressure density at μ = 0) is a μ-independent constant; the second term, which is totally determined by G[μ] (p) (the dressed quark propagator at finite μ), contains all the nontrivial μ-dependence. By applying a general result in the rainbow-ladder approximation of the Dyson-Schwinger approach obtained in our previous study [Phys. Rev. C 71 (2005) 015205], G[μ](p) is calculated from the meromorphic quark propagator proposed in [Phys. Rev. D 67 (2003) 054019]. From this the full analytic expression of the EOS of QCD at finite μ and zero T is obtained (apart from the constant term P(μ)|μ=0, which can in principle be caJculated from the CJT effective action). A comparison between our EOS and the cold, perturbative EOS of QCD of Fraga, Pisarski and Schaffner-Bielich is made. It is expected that our EOS can provide a possible new approach for the study of neutron stars.
基金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 (No.91641104, No.21733006, and No.21590802)
文摘Complex-forming reactions widely exist in gas-phase chemical reactions. Various complexforming bimolecular reactions have been investigated and interesting phenomena have been discovered. The complex-forming reactions usually have small or no barrier in the entrance channel, which leads to obvious differences in kinetic and dynamic characteristics compared with direct reactions. Theoretically, quantum state-resolved reaction dynamics can provide the most detailed microscopic dynamic mechanisms and is now feasible for a direct reaction with only one potential barrier. However, it is of great challenge to construct accurate potential energy surfaces and perform accurate quantum dynamics calculations for a complex polyatomic reaction involving deep potential wells and multi-channels. This paper reviews the most recent progress in two prototypical oxyhydrogen complex-forming reaction systems, HO2 and HO3, which are significant in combustion, atmospheric, and interstellar chemistry. We will present a brief survey of both computational and experimental work and emphasize on some unsolved problems existing in these systems.