Both current and pressure coupling schemes have been adopted in the hybrid kinetic–magnetohydrodynamic code CLT-K recently.Numerical equivalences between these two coupling schemes are strictly verified under differe...Both current and pressure coupling schemes have been adopted in the hybrid kinetic–magnetohydrodynamic code CLT-K recently.Numerical equivalences between these two coupling schemes are strictly verified under different approximations.First,when considering only the perturbed distribution function of energetic particles(EPs),the equivalence can be proved analytically.Second,when both the variations of the magnetic field and the EP distribution function are included,the current and pressure coupling schemes numerically produce the same result in the nonlinear simulations.On this basis,the influences of co-/counter-passing and trapped EPs on the linear stabilities of the m/n=2/1 tearing mode(TM)have been investigated(where m and n represent the poloidal and toroidal mode numbers,respectively).The results of scanningβh of EPs show that the co-passing and trapped EPs are found to stabilize the TM,while the counter-passing EPs tend to destabilize the TM.The behind(de)stabilization mechanisms of the TM by EPs are carefully analyzed.Furthermore,after exceeding critical EP betas,the same branch of the high-frequency mode is excited by co-/counterpassing and trapped EPs,which is identified as the m/n=2/1 energetic particle mode.展开更多
A coupled system simulating both firebox and reactor is established to study the naphtha pyrolysis in an industrial tubular furnace.The firebox model is based on zone method including combustion,radiation,and convecti...A coupled system simulating both firebox and reactor is established to study the naphtha pyrolysis in an industrial tubular furnace.The firebox model is based on zone method including combustion,radiation,and convection to simulate heat transfer in the furnace.A two-dimensional recirculation model is proposed to estimate the flow field in furnace.The reactor model integrates the feedstock reconstruction model,an auto-generator of detail kinetic schemes,and the reactor simulation model to simulate the reaction process in the tubular coil.The coupled simulation result is compared with industrial process and shows agreement within short computation time.展开更多
We present an investigation of the one-dimensional ferromagnetism in Au–Co nanowires deposited on the Cu（110）surface. By using the density functional theory, the influence of the nonmagnetic copper substrate Cu（11...We present an investigation of the one-dimensional ferromagnetism in Au–Co nanowires deposited on the Cu（110）surface. By using the density functional theory, the influence of the nonmagnetic copper substrate Cu（110） on the magnetic properties of the bimetallic Au–Co nanowires is studied. The results show the emergence of magnetic anisotropy in the supported Au–Co nanowires. The magnetic anisotropy energy has the same order of magnitude as the exchange interaction energy between Co atoms in the wire. Our electronic structure calculation reveals the emergence of new hybridized bands between Au and Co atoms and surface Cu atoms. The Curie temperature of the Au–Co wires is calculated by means of kinetic Monte Carlo simulation. The strong size effect of the Curie temperature is demonstrated.展开更多
The effects of nanosecond discharge on ignition characteristics of a stoichiometric methane–air mixture without inert diluent gas were studied by numerical simulation at 0.1 MPa and an initial temperature of 1300 K. ...The effects of nanosecond discharge on ignition characteristics of a stoichiometric methane–air mixture without inert diluent gas were studied by numerical simulation at 0.1 MPa and an initial temperature of 1300 K. A modified non-equilibrium plasma kinetic model was developed to simulate the temporal evolution of particles produced during nanosecond discharge and its afterglow. As important roles in ignition, path fluxes of O and H radicals were analyzed in detail. Different strength of E/N and different discharge duration were applied to the discharge process in this study. And the results presented that a deposited energy of 1–30 m J·cm^(-3) could dramatically reduce the ignition delay time. Furthermore, temperature and radicals analysis was conducted to investigate the effect of non-equilibrium plasma on production of intermediate radicals. Finally, sensitivity analysis was employed to have further understanding on ignition chemistries of the mixture under nanosecond discharge.展开更多
Metal additive manufacturing (AM) is a disruptive manufacturing technology that takes into account the needs of complex structural forming and high-performance component forming. At present, the understanding of metal...Metal additive manufacturing (AM) is a disruptive manufacturing technology that takes into account the needs of complex structural forming and high-performance component forming. At present, the understanding of metal additive manufacturing simulation methods is not thorough enough, which restricts the development of metal additive manufacturing. Present work discusses the evolution of KMC method simulation results for simulating metal additive manufacturing at different length ratios and different scanning speeds. The results reveal that as the scanning speed increases, the main grains in the simulation results are transformed from coarse columnar grains to crescent-shaped grains, which are in good agreement with the existing experimental results. Besides, as the ratio of unit physical length to unit simulation length increases, the ratio of unit physical time to unit simulation time gradually decreases.展开更多
We develop a coarse grained （CG） approach for efficiently simulating calcium dynamics in the endoplasmic reticulum membrane based on a fine stochastic lattice gas model. By grouping neighboring microscopic sites tog...We develop a coarse grained （CG） approach for efficiently simulating calcium dynamics in the endoplasmic reticulum membrane based on a fine stochastic lattice gas model. By grouping neighboring microscopic sites together into CG cells and deriving CG reaction rates using local mean field approximation, we perform CG kinetic Monte Carlo （kMC） simulations and find the results of CG-kMC simulations are in excellent agreement with that of the microscopic ones. Strikingly, there is an appropriate range of coarse proportion rn, corresponding to the minimal deviation of the phase transition point compared to the microscopic one. For fixed m, the critical point increases monotonously as the system size increases, especially, there exists scaling law between the deviations of the phase transition point and the system size. Moreover, the CG approach provides significantly faster Monte Carlo simulations which are easy to implement and are directly related to the microscopics, so that one can study the system size effects at the cost of reasonable computational time.展开更多
The processes of multilayer thin Cu films grown on Cu （100） surfaces at elevated temperature （250-400K） are simulated by mean of kinetic Monte Carlo （KMC） method, where the realistic growth model and physical pa...The processes of multilayer thin Cu films grown on Cu （100） surfaces at elevated temperature （250-400K） are simulated by mean of kinetic Monte Carlo （KMC） method, where the realistic growth model and physical parameters are used. The effects of small island （dimer and trimer） diffusion, edge diffusion along the islands, exchange of the adatom with an atom in the existing island, as well as mass transport between interlayers are included in the simulation model. Emphasis is placed on revealing the influence of the Ehrlic-Schwoebel （ES） barrier on growth mode and morphology during multilayer thin film growth. We present numerical evidence that the ES barrier does exist for the Cu/Cu（100） system and an ES barrier EB 〉 0.125eV is estimated from a comparison of the KMC simulation with the realistic experimental images. The transitions of growth modes with growth conditions and the influence of exchange barrier on growth mode are also investigated.展开更多
This paper presents the experimental progress of laser-focused Cr atomic deposition and the experimental condition. The result is an accurate array of lines with a periodicity of 212.8±0.2 nm and mean full-width ...This paper presents the experimental progress of laser-focused Cr atomic deposition and the experimental condition. The result is an accurate array of lines with a periodicity of 212.8±0.2 nm and mean full-width at half maximum as approximately 95 nm. Surface growth in laser-focused Cr atomic deposition is modeled and studied by kinetic Monte Carlo simulation including two events： the one is that atom trajectories in laser standing wave are simulated with the semiclassical equations of motion to obtain the deposition position; the other is that adatom diffuses by considering two major diffusion processes, namely, terrace diffusion and step-edge descending. Comparing with experimental results （Anderson W R, Bradley C C, McClelland J J and Celotta R J 1999 Phys. Rev. A 59 2476）, it finds that the simulated trend of dependence on feature width is in agreement with the power of standing wave, the other two simulated trends are the same in the initial stage. These results demonstrate that some surface diffusion processes play important role in feature width broadening. Numerical result also shows that high incoming beam flux of atoms deposited redounds to decrease the distance between adatoms which can diffuse to minimize the feature width and enhance the contrast.展开更多
LiFePO_(4),one of the mainstream cathode materials of current EV batteries,exhibits experimental diffusion coefficients(D_(c))of Li^(+)which are not only several orders of magnitude lower than those predicted by the i...LiFePO_(4),one of the mainstream cathode materials of current EV batteries,exhibits experimental diffusion coefficients(D_(c))of Li^(+)which are not only several orders of magnitude lower than those predicted by the ionic hopping barriers obtained from theoretical calculations and spectroscopic measurements,but also span several orders from 10^(-14)to 10^(-18)cm^(2)s^(-1)under different states of charge(SOC)and the charging rates(C-rates).Atomic level understanding of such sluggishness and diversity of Li^(+)transport kinetics would be of significance in improving the rate performance of LiFePO_(4)through material and operation optimization but remain challenging.Herein,we show that the high sensitivity of Li^(+)hopping barriers on the local Li–Li coordination environments(numbers and configurations)plays a key role in the ion transport kinetics.This is due a neural network-based deep potential(DP)which allows accurate and efficient calculation of hopping barriers of Li^(+)in LiFePO_(4)with various Li–Li coordination environments,with which the kinetic Monte-Carlo(KMC)method was employed to determine the D_(c)values at various C-rates and SOC across a broad spectrum.Especially,an accelerated KMC simulation strategy is proposed to obtain the D_(c)values under a wide range of SOC at low C-rates,which agree well with that obtained from the galvanostatic intermittent titration technique(GITT).The present study provides accurate descriptions of Li^(+)transport kinetics at both very high and low C-rates,which remains challenging to experiments and first-principles calculations,respectively.Finally,it is revealed that the gradient distributions of Li^(+)density along the diffusion path result in great asymmetry in the barriers of the forward and backward hopping,causing very slow diffusion of Li^(+)and the diverse variation of D_(c).展开更多
In supersonic flowing plasmas,the auto-resonant behavior of ion acoustic waves driven by stimulated Brillouin backscattering is self-consistently investigated.A nature of absolute instability appears in the evolution ...In supersonic flowing plasmas,the auto-resonant behavior of ion acoustic waves driven by stimulated Brillouin backscattering is self-consistently investigated.A nature of absolute instability appears in the evolution of the stimulated Brillouin backscattering.By adopting certain form of incident lights combined by two perpendicular linear polarization lasers or polarization rotation lasers,the absolute instability is suppressed significantly.The suppression of auto-resonant stimulated Brillouin scattering is verified with the fully kinetic Vlasov code.展开更多
A comprehensive simulation model -- deposition, diffusion, rotation, reaction and aggregation model is presented to simulate the formation processes of ramified clusters on liquid surfaces, where clusters can disuse a...A comprehensive simulation model -- deposition, diffusion, rotation, reaction and aggregation model is presented to simulate the formation processes of ramified clusters on liquid surfaces, where clusters can disuse and rotate easily. The mobility （including diffusion and rotation） of clusters is related to its mass, which is given by Dm = Dos^-γD and θm = θos^-γθ, respectively. The influence of the reaction probability on the kinetics and structure formation is included in the simulation model. We concentrate on revealing dynamic scaling during ramified cluster formation. For this purpose, the time evolution of the cluster density and the weight-average cluster size as well as the cluster-size distribution scaling function at different time are determined for various conditions. The dependence of the cluster density on the deposition flux and time-dependence of fractal dimension are also investigated. The obtained results are helpful in understanding the formation of clusters or thin film growth on liquid surfaces.展开更多
A realistic kinetic Monte Carlo （KMC） simulation model with physical parameters is developed, which well reproduces the heteroepitaxial growth of multilayered Ni thin film on Cu（100） surfaces at room temperature. ...A realistic kinetic Monte Carlo （KMC） simulation model with physical parameters is developed, which well reproduces the heteroepitaxial growth of multilayered Ni thin film on Cu（100） surfaces at room temperature. The effects of mass transport between interlayers and edge diffusion of atoms along the islands are included in the simulation model, and the surface roughness and the layer distribution versus total coverage are calculated. Speeially, the simulation model reveals the transition of growth mode with coverage and the difference between the Ni heteroepitaxy on Cu（100） and the Ni homoepitaxy on Ni（100）. Through comparison of KMC simulation with the real scanning tunneling microscopy （STM） experiments, the Ehrlich-Schwoebel （ES） barrier Ees is estimated to be 0.18±0.02 eV for Ni/Cu（100） system while 0.28 eV for Ni/Ni（100）. The simulation also shows that the growth mode depends strongly on the thickness of thin film and the surface temperature, and the critical thickness of growth mode transition is dependent on the growth condition such as surface temperature and deposition flux as well.展开更多
Low-dimensional materials have excellent properties which are closely related to their dimensionality.However,the growth mechanism underlying tunable dimensionality from 2D triangles to 1D ribbons of such materials is...Low-dimensional materials have excellent properties which are closely related to their dimensionality.However,the growth mechanism underlying tunable dimensionality from 2D triangles to 1D ribbons of such materials is still unrevealed.Here,we establish a general kinetic Monte Carlo model for transition metal dichalcogenides(TMDs) growth to address such an issue.Our model is able to reproduce several key findings in experiments,and reveals that the dimensionality is determined by the lattice mismatch and the interaction strength between TMDs and the substrate.We predict that the dimensionality can be well tuned by the interaction strength and the geometry of the substrate.Our work deepens the understanding of tunable dimensionality of low-dimensional materials and may inspire new concepts for the design of such materials with expected dimensionality.展开更多
The heteroepitaxial growth of multilayer Cu/Pd（100） thin film via pulse laser deposition （PLD） at room temperature is simulated by using kinetic Monte Carlo （KMC） method with realistic physical parameters. The e...The heteroepitaxial growth of multilayer Cu/Pd（100） thin film via pulse laser deposition （PLD） at room temperature is simulated by using kinetic Monte Carlo （KMC） method with realistic physical parameters. The effects of mass transport between interlayers, edge diffusion of adatoms along the islands and instantaneous deposition are considered in the simulation model, Emphasis is placed on revealing the details of multilayer Cu/Pd（100） thin film growth and estimating the Ehrlich-Schwoebel （ES） barrier. It is shown that the instantaneous deposition in the PLD growth gives rise to the layer-by-layer growth mode, persisting up to about 9 monolayers （ML） of Cu/Pd（100）. The ES barriers of 0.08 ± 0.01 eV is estimated by comparing the KMC simulation results with the real scanning tunnelling microscopy （STM） measurements,展开更多
The turbulent flow over a channel bed roughened by three layers of closely packed spheres with a Reynolds number of Re= 15 000 is investigated using the large eddy simulation（LES） and the double-averaging（DA） meth...The turbulent flow over a channel bed roughened by three layers of closely packed spheres with a Reynolds number of Re= 15 000 is investigated using the large eddy simulation（LES） and the double-averaging（DA） method. The DA velocity is compared with the results of the corresponding laboratory experiments to validate the LES results. The existence of the types of vortex structures is demonstrated by the Q-criterion above the permeable bed. The turbulent kinetic energy（TKE） fluxes and budget are quantified and discussed. The results show that the TKE fluxes are directed downward and downstream near the virtual bed level. In the TKE budget, the form-induced diffusion rate is significant in the vicinity of the crest bed level, and the TKE production rate and the dissipation rate attain their peaks at the crest bed level and decrease sharply below it.展开更多
We report kinetic Monte-Karlo(KMC)simulation of self-assembled synthesis of nanocrystals by physical vapor deposition(PVD),which is one of most flexible,efficient,and clean techniques to fabricate nanopatterns.In part...We report kinetic Monte-Karlo(KMC)simulation of self-assembled synthesis of nanocrystals by physical vapor deposition(PVD),which is one of most flexible,efficient,and clean techniques to fabricate nanopatterns.In particular,self-assembled arrays of nanocrystals can be synthesized by PVD.However size,shape and density of self-assembled nanocrystals are highly sensitive to the process conditions such as duration of deposition,temperature,substrate material,etc.To efficiently synthesize nanocrystalline arrays by PVD,the process control factors should be understood in detail.KMC simulations of film deposition are an important tool for understanding the mechanisms of film deposition.In this paper,we report a KMC modeling that explicitly represents PVD synthesis of self-assembled nanocrystals.We study how varying critical process parameters such as deposition rate,duration,temperature,and substrate type affect the lateral 2D morphologies of self-assembled metallic islands on substrates,and compare our results with experimentally observed surface morphologies generated by PVD.Our simulations align well with experimental results reported in the literature.展开更多
The conventional low-temperature method of solution polycondensation was developed to realize the reaction of p-phenylenediamin and terephthaloyl chloride for the preparation of poly（p-phenylene terephthalamide）（P...The conventional low-temperature method of solution polycondensation was developed to realize the reaction of p-phenylenediamin and terephthaloyl chloride for the preparation of poly（p-phenylene terephthalamide）（PPTA）. Some main factors influencing this process were investigated to determine the optimum condition for high molecular weight. Experiment showed significant slowing of the reaction and gradual deviation of second-order reaction kinetics due to diffusion control. These phenomena were studied theoretically via dynamic Monte Carlo simulation. A concise expression,n～c0^-0.88·t^0.37, was proposed to describe the diffusioncontrolled polycondensation process as a function of the monomer concentration and reaction time. The theoretical results provided a good description of diffusion-effected kinetics for the polycondensation process of PPTA, and demonstrated good agreement with the experimental data. Some differences of scaling relations between model and experiment results were also discussed.展开更多
We present simulations of the mechanism of secondary nucleation of polymer crystallization,based on a new model accounting for the microscopic kinetics of attaching and detaching.As the key feature of the model,we int...We present simulations of the mechanism of secondary nucleation of polymer crystallization,based on a new model accounting for the microscopic kinetics of attaching and detaching.As the key feature of the model,we introduced multibody-interaction parameters that establish correlations between the attaching and detaching rate constants and the resulting thickness and width of the crystalline lamella.Using MATLAB and Monte Carlo method,we followed the evolution of the secondary nuclei as a function of various multibody-interaction parameters.We identified three different growth progressions of the crystal：（i） Widening,（ii） thickening and（iii） simultaneously thickening and widening of lamellar crystals,controlled by the corresponding kinetic parameters.展开更多
4H-silicon carbides deposited by diamond films have wide applications in many fields such as semiconductor heterojunction,heat sink and mechanical sealing.Nucleation plays a critical role in the deposition of the diam...4H-silicon carbides deposited by diamond films have wide applications in many fields such as semiconductor heterojunction,heat sink and mechanical sealing.Nucleation plays a critical role in the deposition of the diamond film on 4H-silicon carbides.Nevertheless,as a typical polar material,the fundamental mechanism of diamond nucleation on different faces of 4H-silicon carbides has not been fully understood yet.In this contribution,nucleation of diamond was performed on the carbon-and silicon-faces of 4H-silicon carbides in a direct current chemical vapor deposition device.The nucleation density on the carbon-face is higher by 2-3 orders of magnitude compared to the silicon-face.Transmission electron microscopy verifies that there are high density diamond nuclei on the interface between the carbon-face and the diamond film,which is different from columnar diamond growth structure on the silicon-face.Transition state theory calculation reveals that the unprecedented distinction of the nucleation density between the carbon-face and the silicon-face is attributed to different desorption rates of the absorbed hydrocarbon radicals.In addition,kinetic model simulations demonstrate that it is more difficult to form CH2(s)-CH2(s)dimers on silicon-faces than carbon-faces,resulting in much lower nucleation densities on silicon-faces.展开更多
基金supported by National Natural Science Foundation of China(No.11835010)by the National Key Research and Development Program of China(Nos.2019YFE03030004 and 2022YFE031000001)。
文摘Both current and pressure coupling schemes have been adopted in the hybrid kinetic–magnetohydrodynamic code CLT-K recently.Numerical equivalences between these two coupling schemes are strictly verified under different approximations.First,when considering only the perturbed distribution function of energetic particles(EPs),the equivalence can be proved analytically.Second,when both the variations of the magnetic field and the EP distribution function are included,the current and pressure coupling schemes numerically produce the same result in the nonlinear simulations.On this basis,the influences of co-/counter-passing and trapped EPs on the linear stabilities of the m/n=2/1 tearing mode(TM)have been investigated(where m and n represent the poloidal and toroidal mode numbers,respectively).The results of scanningβh of EPs show that the co-passing and trapped EPs are found to stabilize the TM,while the counter-passing EPs tend to destabilize the TM.The behind(de)stabilization mechanisms of the TM by EPs are carefully analyzed.Furthermore,after exceeding critical EP betas,the same branch of the high-frequency mode is excited by co-/counterpassing and trapped EPs,which is identified as the m/n=2/1 energetic particle mode.
基金Supported by the National Natural Science Foundation of China(U1462206)
文摘A coupled system simulating both firebox and reactor is established to study the naphtha pyrolysis in an industrial tubular furnace.The firebox model is based on zone method including combustion,radiation,and convection to simulate heat transfer in the furnace.A two-dimensional recirculation model is proposed to estimate the flow field in furnace.The reactor model integrates the feedstock reconstruction model,an auto-generator of detail kinetic schemes,and the reactor simulation model to simulate the reaction process in the tubular coil.The coupled simulation result is compared with industrial process and shows agreement within short computation time.
基金Project supported by the Russian Foundation of Basic Researches
文摘We present an investigation of the one-dimensional ferromagnetism in Au–Co nanowires deposited on the Cu（110）surface. By using the density functional theory, the influence of the nonmagnetic copper substrate Cu（110） on the magnetic properties of the bimetallic Au–Co nanowires is studied. The results show the emergence of magnetic anisotropy in the supported Au–Co nanowires. The magnetic anisotropy energy has the same order of magnitude as the exchange interaction energy between Co atoms in the wire. Our electronic structure calculation reveals the emergence of new hybridized bands between Au and Co atoms and surface Cu atoms. The Curie temperature of the Au–Co wires is calculated by means of kinetic Monte Carlo simulation. The strong size effect of the Curie temperature is demonstrated.
基金Supported by the National Natural Science Foundation of China(No.51376021)the Fundamental Research Funds for the Central Universities(No.2015YJS146)
文摘The effects of nanosecond discharge on ignition characteristics of a stoichiometric methane–air mixture without inert diluent gas were studied by numerical simulation at 0.1 MPa and an initial temperature of 1300 K. A modified non-equilibrium plasma kinetic model was developed to simulate the temporal evolution of particles produced during nanosecond discharge and its afterglow. As important roles in ignition, path fluxes of O and H radicals were analyzed in detail. Different strength of E/N and different discharge duration were applied to the discharge process in this study. And the results presented that a deposited energy of 1–30 m J·cm^(-3) could dramatically reduce the ignition delay time. Furthermore, temperature and radicals analysis was conducted to investigate the effect of non-equilibrium plasma on production of intermediate radicals. Finally, sensitivity analysis was employed to have further understanding on ignition chemistries of the mixture under nanosecond discharge.
文摘Metal additive manufacturing (AM) is a disruptive manufacturing technology that takes into account the needs of complex structural forming and high-performance component forming. At present, the understanding of metal additive manufacturing simulation methods is not thorough enough, which restricts the development of metal additive manufacturing. Present work discusses the evolution of KMC method simulation results for simulating metal additive manufacturing at different length ratios and different scanning speeds. The results reveal that as the scanning speed increases, the main grains in the simulation results are transformed from coarse columnar grains to crescent-shaped grains, which are in good agreement with the existing experimental results. Besides, as the ratio of unit physical length to unit simulation length increases, the ratio of unit physical time to unit simulation time gradually decreases.
基金This work was supported by the National Natural Science Foundation of China （No.11205002）. Chuansheng Shen was also supported by the Key Scientific Research Fund of Anhui Provincial Education Department （No.KJ2012A189）.
文摘We develop a coarse grained （CG） approach for efficiently simulating calcium dynamics in the endoplasmic reticulum membrane based on a fine stochastic lattice gas model. By grouping neighboring microscopic sites together into CG cells and deriving CG reaction rates using local mean field approximation, we perform CG kinetic Monte Carlo （kMC） simulations and find the results of CG-kMC simulations are in excellent agreement with that of the microscopic ones. Strikingly, there is an appropriate range of coarse proportion rn, corresponding to the minimal deviation of the phase transition point compared to the microscopic one. For fixed m, the critical point increases monotonously as the system size increases, especially, there exists scaling law between the deviations of the phase transition point and the system size. Moreover, the CG approach provides significantly faster Monte Carlo simulations which are easy to implement and are directly related to the microscopics, so that one can study the system size effects at the cost of reasonable computational time.
基金Project supported by the Natural Science Foundation for Young Scientists of Zhejiang Province, China （Grant No RC02069）.
文摘The processes of multilayer thin Cu films grown on Cu （100） surfaces at elevated temperature （250-400K） are simulated by mean of kinetic Monte Carlo （KMC） method, where the realistic growth model and physical parameters are used. The effects of small island （dimer and trimer） diffusion, edge diffusion along the islands, exchange of the adatom with an atom in the existing island, as well as mass transport between interlayers are included in the simulation model. Emphasis is placed on revealing the influence of the Ehrlic-Schwoebel （ES） barrier on growth mode and morphology during multilayer thin film growth. We present numerical evidence that the ES barrier does exist for the Cu/Cu（100） system and an ES barrier EB 〉 0.125eV is estimated from a comparison of the KMC simulation with the realistic experimental images. The transitions of growth modes with growth conditions and the influence of exchange barrier on growth mode are also investigated.
基金Project supported by the Nanoscience Foundation of Shanghai,China (Grant Nos. 0852nm07000 and 0952nm07000)the National Natural Science Foundation of China (Grant No. 10804084)the National Science & Technology Support Project (Grant No. 2006BAF06B08)
文摘This paper presents the experimental progress of laser-focused Cr atomic deposition and the experimental condition. The result is an accurate array of lines with a periodicity of 212.8±0.2 nm and mean full-width at half maximum as approximately 95 nm. Surface growth in laser-focused Cr atomic deposition is modeled and studied by kinetic Monte Carlo simulation including two events： the one is that atom trajectories in laser standing wave are simulated with the semiclassical equations of motion to obtain the deposition position; the other is that adatom diffuses by considering two major diffusion processes, namely, terrace diffusion and step-edge descending. Comparing with experimental results （Anderson W R, Bradley C C, McClelland J J and Celotta R J 1999 Phys. Rev. A 59 2476）, it finds that the simulated trend of dependence on feature width is in agreement with the power of standing wave, the other two simulated trends are the same in the initial stage. These results demonstrate that some surface diffusion processes play important role in feature width broadening. Numerical result also shows that high incoming beam flux of atoms deposited redounds to decrease the distance between adatoms which can diffuse to minimize the feature width and enhance the contrast.
基金financially supported by the National Natural Science Foundation of China(22272122,21832004 and 21673163)。
文摘LiFePO_(4),one of the mainstream cathode materials of current EV batteries,exhibits experimental diffusion coefficients(D_(c))of Li^(+)which are not only several orders of magnitude lower than those predicted by the ionic hopping barriers obtained from theoretical calculations and spectroscopic measurements,but also span several orders from 10^(-14)to 10^(-18)cm^(2)s^(-1)under different states of charge(SOC)and the charging rates(C-rates).Atomic level understanding of such sluggishness and diversity of Li^(+)transport kinetics would be of significance in improving the rate performance of LiFePO_(4)through material and operation optimization but remain challenging.Herein,we show that the high sensitivity of Li^(+)hopping barriers on the local Li–Li coordination environments(numbers and configurations)plays a key role in the ion transport kinetics.This is due a neural network-based deep potential(DP)which allows accurate and efficient calculation of hopping barriers of Li^(+)in LiFePO_(4)with various Li–Li coordination environments,with which the kinetic Monte-Carlo(KMC)method was employed to determine the D_(c)values at various C-rates and SOC across a broad spectrum.Especially,an accelerated KMC simulation strategy is proposed to obtain the D_(c)values under a wide range of SOC at low C-rates,which agree well with that obtained from the galvanostatic intermittent titration technique(GITT).The present study provides accurate descriptions of Li^(+)transport kinetics at both very high and low C-rates,which remains challenging to experiments and first-principles calculations,respectively.Finally,it is revealed that the gradient distributions of Li^(+)density along the diffusion path result in great asymmetry in the barriers of the forward and backward hopping,causing very slow diffusion of Li^(+)and the diverse variation of D_(c).
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11875091 and 11975059)the Science Challenge Project,China(Grant No.TZ2016005).
文摘In supersonic flowing plasmas,the auto-resonant behavior of ion acoustic waves driven by stimulated Brillouin backscattering is self-consistently investigated.A nature of absolute instability appears in the evolution of the stimulated Brillouin backscattering.By adopting certain form of incident lights combined by two perpendicular linear polarization lasers or polarization rotation lasers,the absolute instability is suppressed significantly.The suppression of auto-resonant stimulated Brillouin scattering is verified with the fully kinetic Vlasov code.
基金The project supported by the Natural Science Foundation for Young Scientists of Zhejiang Province of China under Grant No. RC02069 .
文摘A comprehensive simulation model -- deposition, diffusion, rotation, reaction and aggregation model is presented to simulate the formation processes of ramified clusters on liquid surfaces, where clusters can disuse and rotate easily. The mobility （including diffusion and rotation） of clusters is related to its mass, which is given by Dm = Dos^-γD and θm = θos^-γθ, respectively. The influence of the reaction probability on the kinetics and structure formation is included in the simulation model. We concentrate on revealing dynamic scaling during ramified cluster formation. For this purpose, the time evolution of the cluster density and the weight-average cluster size as well as the cluster-size distribution scaling function at different time are determined for various conditions. The dependence of the cluster density on the deposition flux and time-dependence of fractal dimension are also investigated. The obtained results are helpful in understanding the formation of clusters or thin film growth on liquid surfaces.
基金the State Key Basic Research Development Project of China under Grant No.2006CB708612the Key Project of Science and Technology of Zhejiang Province under Grant No.2007C21120
文摘A realistic kinetic Monte Carlo （KMC） simulation model with physical parameters is developed, which well reproduces the heteroepitaxial growth of multilayered Ni thin film on Cu（100） surfaces at room temperature. The effects of mass transport between interlayers and edge diffusion of atoms along the islands are included in the simulation model, and the surface roughness and the layer distribution versus total coverage are calculated. Speeially, the simulation model reveals the transition of growth mode with coverage and the difference between the Ni heteroepitaxy on Cu（100） and the Ni homoepitaxy on Ni（100）. Through comparison of KMC simulation with the real scanning tunneling microscopy （STM） experiments, the Ehrlich-Schwoebel （ES） barrier Ees is estimated to be 0.18±0.02 eV for Ni/Cu（100） system while 0.28 eV for Ni/Ni（100）. The simulation also shows that the growth mode depends strongly on the thickness of thin film and the surface temperature, and the critical thickness of growth mode transition is dependent on the growth condition such as surface temperature and deposition flux as well.
基金supported by the Ministry of Science and Technology (No.2018YFA0208702)the National Natural Science Foundation of China (No.32090044,No. 21973085,No.21833007,No.21790350)+1 种基金Anhui Initiative in Quantum Information Technologies (AHY 090200)the Fundamental Research Funds for the Central Universities (WK2340000104)。
文摘Low-dimensional materials have excellent properties which are closely related to their dimensionality.However,the growth mechanism underlying tunable dimensionality from 2D triangles to 1D ribbons of such materials is still unrevealed.Here,we establish a general kinetic Monte Carlo model for transition metal dichalcogenides(TMDs) growth to address such an issue.Our model is able to reproduce several key findings in experiments,and reveals that the dimensionality is determined by the lattice mismatch and the interaction strength between TMDs and the substrate.We predict that the dimensionality can be well tuned by the interaction strength and the geometry of the substrate.Our work deepens the understanding of tunable dimensionality of low-dimensional materials and may inspire new concepts for the design of such materials with expected dimensionality.
基金Project supported by the State Key Development Program for Basic Research of China （Grant No 2006CB708612） and Natural Science Foundation for Young Scientists of Zhejiang Province, China （Grant No RC02069）.Acknowledgment We would like to thank Dr Ling-wang Wang of the Computational Research Division at the Lawrence Berkeley National Lab and Dr Xiangrong Ye of Department of Material Science and Chemical Engineering, University of California at San Diego for helpful discussion.
文摘The heteroepitaxial growth of multilayer Cu/Pd（100） thin film via pulse laser deposition （PLD） at room temperature is simulated by using kinetic Monte Carlo （KMC） method with realistic physical parameters. The effects of mass transport between interlayers, edge diffusion of adatoms along the islands and instantaneous deposition are considered in the simulation model, Emphasis is placed on revealing the details of multilayer Cu/Pd（100） thin film growth and estimating the Ehrlich-Schwoebel （ES） barrier. It is shown that the instantaneous deposition in the PLD growth gives rise to the layer-by-layer growth mode, persisting up to about 9 monolayers （ML） of Cu/Pd（100）. The ES barriers of 0.08 ± 0.01 eV is estimated by comparing the KMC simulation results with the real scanning tunnelling microscopy （STM） measurements,
基金Project supported by the National Natural Science Foun-dation of China(Grant No.11372161,51209230)
文摘The turbulent flow over a channel bed roughened by three layers of closely packed spheres with a Reynolds number of Re= 15 000 is investigated using the large eddy simulation（LES） and the double-averaging（DA） method. The DA velocity is compared with the results of the corresponding laboratory experiments to validate the LES results. The existence of the types of vortex structures is demonstrated by the Q-criterion above the permeable bed. The turbulent kinetic energy（TKE） fluxes and budget are quantified and discussed. The results show that the TKE fluxes are directed downward and downstream near the virtual bed level. In the TKE budget, the form-induced diffusion rate is significant in the vicinity of the crest bed level, and the TKE production rate and the dissipation rate attain their peaks at the crest bed level and decrease sharply below it.
文摘We report kinetic Monte-Karlo(KMC)simulation of self-assembled synthesis of nanocrystals by physical vapor deposition(PVD),which is one of most flexible,efficient,and clean techniques to fabricate nanopatterns.In particular,self-assembled arrays of nanocrystals can be synthesized by PVD.However size,shape and density of self-assembled nanocrystals are highly sensitive to the process conditions such as duration of deposition,temperature,substrate material,etc.To efficiently synthesize nanocrystalline arrays by PVD,the process control factors should be understood in detail.KMC simulations of film deposition are an important tool for understanding the mechanisms of film deposition.In this paper,we report a KMC modeling that explicitly represents PVD synthesis of self-assembled nanocrystals.We study how varying critical process parameters such as deposition rate,duration,temperature,and substrate type affect the lateral 2D morphologies of self-assembled metallic islands on substrates,and compare our results with experimentally observed surface morphologies generated by PVD.Our simulations align well with experimental results reported in the literature.
基金financially supported by the National Natural Science Foundation of China (Nos.21204011,51603120 and 21404023)the National Basic Research Program of the China 973 Program (No.2011CB606101)+1 种基金Fundamental Research Funds for the Central Universities,Innovation Program of Shanghai Science and Technology Commission (No.14521100605)the Innovation Program of Shanghai Municipal Education Commission
文摘The conventional low-temperature method of solution polycondensation was developed to realize the reaction of p-phenylenediamin and terephthaloyl chloride for the preparation of poly（p-phenylene terephthalamide）（PPTA）. Some main factors influencing this process were investigated to determine the optimum condition for high molecular weight. Experiment showed significant slowing of the reaction and gradual deviation of second-order reaction kinetics due to diffusion control. These phenomena were studied theoretically via dynamic Monte Carlo simulation. A concise expression,n～c0^-0.88·t^0.37, was proposed to describe the diffusioncontrolled polycondensation process as a function of the monomer concentration and reaction time. The theoretical results provided a good description of diffusion-effected kinetics for the polycondensation process of PPTA, and demonstrated good agreement with the experimental data. Some differences of scaling relations between model and experiment results were also discussed.
基金financially supported by the National Natural Science Foundation of China(No.21374054)the Sino-German Center for Research Promotion
文摘We present simulations of the mechanism of secondary nucleation of polymer crystallization,based on a new model accounting for the microscopic kinetics of attaching and detaching.As the key feature of the model,we introduced multibody-interaction parameters that establish correlations between the attaching and detaching rate constants and the resulting thickness and width of the crystalline lamella.Using MATLAB and Monte Carlo method,we followed the evolution of the secondary nuclei as a function of various multibody-interaction parameters.We identified three different growth progressions of the crystal：（i） Widening,（ii） thickening and（iii） simultaneously thickening and widening of lamellar crystals,controlled by the corresponding kinetic parameters.
基金the National Key Research and Development Project(No.2017YFE0128600)Ningbo 3315 Innovation Team(No.2019A-18-C)+9 种基金Science and Technology Innovation 2025 Major Project of Ningbo(No.2018023)National Defense Key Laboratory Fund(No.6142807180511)Innovation Funding of State Oceanic Administration(No.NBHY-2017-Z3)Ningbo Industrial Technology Innovation Project(No.2016B10038)‘13th Five-Year’Equipment Pre-research Sharing Project(No.E1710161)‘Key Talents’Senior Engineer Project of Ningbo Institute of Materials Technology and Engineeringthe financial support of CONICYT in the project Fondecyt 11180121the financial support given the VID in the framework of U-IniciaUI013/2018 and the academic direction of the University of Chilefinancial support from the Swedish Government Strategic Research Area in Materials Science on Functional Materials at Linkoping University(Faculty Grant SFO Mat LiU No.200900971)the Swedish Research Council(VR)。
文摘4H-silicon carbides deposited by diamond films have wide applications in many fields such as semiconductor heterojunction,heat sink and mechanical sealing.Nucleation plays a critical role in the deposition of the diamond film on 4H-silicon carbides.Nevertheless,as a typical polar material,the fundamental mechanism of diamond nucleation on different faces of 4H-silicon carbides has not been fully understood yet.In this contribution,nucleation of diamond was performed on the carbon-and silicon-faces of 4H-silicon carbides in a direct current chemical vapor deposition device.The nucleation density on the carbon-face is higher by 2-3 orders of magnitude compared to the silicon-face.Transmission electron microscopy verifies that there are high density diamond nuclei on the interface between the carbon-face and the diamond film,which is different from columnar diamond growth structure on the silicon-face.Transition state theory calculation reveals that the unprecedented distinction of the nucleation density between the carbon-face and the silicon-face is attributed to different desorption rates of the absorbed hydrocarbon radicals.In addition,kinetic model simulations demonstrate that it is more difficult to form CH2(s)-CH2(s)dimers on silicon-faces than carbon-faces,resulting in much lower nucleation densities on silicon-faces.
