A dynamical simulation method is presented to model the cornea deformation in surgery of laser thermokeratoplasty. The virtual cornea is constructed as a mass-spring system. The corneal surface tension is simulated by...A dynamical simulation method is presented to model the cornea deformation in surgery of laser thermokeratoplasty. The virtual cornea is constructed as a mass-spring system. The corneal surface tension is simulated by damping spring stretch between mass points on the cornea model. The aqueous humor in the eyeball is modeled as ideal gas, and the intraocular pressure is simulated by gas pressure. The coagulation force is exerted on each photocoagulation spot to demonstrate its collapse caused by the condensation of corneal soft tissue irradiated by laser. An extra viscous drag force is added to each mass point to weaken the mass point oscillation. The use of the effective time-corrected Verlet integral method brings about flowing and stable dynamic simulation procedures. The simulation results show that, comparing to the undeformed model, the curvature of the region between the optical center and photocoagulation spot increases obviously. Moreover, the shape of the deformed virtual cornea is much similar to that of the real cornea after surgery.展开更多
INS （Inelastic Neutron Scattering） spectrum of methane hydrate was measured on MARI （a direct-geometry chopper spectrometer） at Ruther-ford Appleton Laboratory, UK. Compared with ice Ih, it is found that the whole...INS （Inelastic Neutron Scattering） spectrum of methane hydrate was measured on MARI （a direct-geometry chopper spectrometer） at Ruther-ford Appleton Laboratory, UK. Compared with ice Ih, it is found that the whole spectrum of methane hydrate moves toward high-energy by about 1.5 meV. Using lattice dynamical （LD） technique, computer simulations of methane hydrate were carried out. In the simulations, four potential models （BE TIP3P, TIP4P, MCY） were employed to calculate the phonon density of states （PDOS）. Comparing the calculated PDOS spectrum with the experimental spectrum, it is found that BF, TIP4P, and TIP3P potential lattices give out well-separated translational and librational bands while MCY potential lattice is unstable to do so and this model is not suitable to describe hydrate system.展开更多
The pre-wetting of aggregate surface is a means to improve the interface performance of SBS modified asphalt and aggregate.The effect of pre-wetting technology on the interaction between SBS modified asphalt and aggre...The pre-wetting of aggregate surface is a means to improve the interface performance of SBS modified asphalt and aggregate.The effect of pre-wetting technology on the interaction between SBS modified asphalt and aggregate was analyzed by molecular dynamics simulation.The diffusion coefficient and concentration distribution of SBS modified asphalt on aggregate surface are included.The simulation results show that the diffusion coefficient of the aggregate surface of SBS modified asphalt is increased by 47.6%and 70.5%respectively after 110#asphalt and 130#asphalt are pre-wetted.The concentration distribution of SBS modified asphalt on the aggregate surface after pre-wetting is more uniform.According to the results of interface energy calculation,the interface energy of SBS modified bitumen and aggregate can be increased by about 5%after pre-wetting.According to the results of molecular dynamics simulation,the pre-wetting technology can effectively improve the interface workability of SBS modified bitumen and aggregate,so as to improve the interface performance.展开更多
We calculate the electrical and thermal conductivity of hydrogen for a wide range of densities and temperatures by using molecular dynamics simulations informed by density functional theory.On the basis of the corresp...We calculate the electrical and thermal conductivity of hydrogen for a wide range of densities and temperatures by using molecular dynamics simulations informed by density functional theory.On the basis of the corresponding extended ab initio data set,we construct interpolation formulas covering the range from low-density,high-temperature to high-density,low-temperature plasmas.Our conductivity model repro-duces the well-known limits of the Spitzer and Ziman theory.We compare with available experimental data andfind very good agreement.The new conductivity model can be applied,for example,in dynamo simulations for magneticfield generation in gas giant planets,brown dwarfs,and stellar envelopes.展开更多
Stimuli-responsive polymers capable of rapidly altering their chain conformation in response to external stimuli exhibit broad applica-tion prospects.Experiments have shown that pressure plays a pivotal role in regula...Stimuli-responsive polymers capable of rapidly altering their chain conformation in response to external stimuli exhibit broad applica-tion prospects.Experiments have shown that pressure plays a pivotal role in regulating the microscopic chain conformation of polymers in mixed solvents,and one notable finding is that increasing the pressure can lead to the vanishing of the co-nonsolvency effect.However,the mecha-nisms underlying this phenomenon remain unclear.In this study,we systematically investigated the influence of pressure on the co-nonsolvency effect of single-chain and multi-chain homopolymers in binary mixed good-solvent systems using molecular dynamics simulations.Our results show that the co-nonsolvency-induced chain conformation transition and aggregation behavior significantly depend on pressure in allsingle-chain and multi-chain systems.In single-chain systems,at low pressures,the polymer chain maintains a collapsed state over a wide range of co-solvent fractions(x-range)owing to the co-nonsolvency effect.As the pressure increases,the x-range of the collapsed state gradually narrows,ac-companied by a progressive expansion of the chain.In multichain systems,polymer chains assemble into approximately spherical aggregates over a broad x-range at low pressures owing to the co-nonsolvency effect.Increasing the pressure reduces the x-range for forming aggregates and leads to the formation of loose aggregates or even to a state of dispersed chains at some x-range.These findings indicate that increasing the pressure can weaken or even offset the co-nonsolvency effect in some x-range,which is in good agreement with the experimental observations.Quantitative analysis of the radial density distributions and radial distribution functions reveals that,with increasing pressure,(1)the densities of both polymers and co-solvent molecules within aggregates decrease,while that of the solvent molecule increases;and(2)the effective interac-tions between the polymer and the co-solvent weaken,whereas those between the polymer and solvent strengthen.This enhances the incorpo-ration of solvent molecules within the chains,thereby weakening or even suppressing the chain aggregation.Our study not only elucidates the regulatory mechanism of pressure on the microscopic chain conformations and aggregation behaviors of polymers,but also may provide theo-retical guidance for designing smart polymericmaterials based on mixed solvents.展开更多
Carbon nanotubes(CNTs),black phosphorus nanotubes(BPNTs),and graphene derivatives exhibit significant promise for applications in nano-electromechanical systems(NEMS),energy storage,and sensing technologies due to the...Carbon nanotubes(CNTs),black phosphorus nanotubes(BPNTs),and graphene derivatives exhibit significant promise for applications in nano-electromechanical systems(NEMS),energy storage,and sensing technologies due to their exceptional mechanical,electrical,and thermal properties.This review summarizes recent advances in understanding the dynamic behaviors of these nanomaterials,with a particular focus on insights gained from molecular dynamics(MD)simulations.Key areas discussed include the oscillatory and rotational dynamics of double-walled CNTs,fabrication and stability challenges associated with BPNTs,and the emerging potential of graphyne nanotubes(GNTs).The review also outlines design strategies for enhancing nanodevice performance and underscores the importance of future efforts in experimental validation,multi-scale coupling analyses,and the development of novel nanocomposites to accelerate practical deployment.展开更多
Eu^(2+)doped fluorosilicate glass-ceramics containing BaF_(2) nanocrystals have high potential as spectral conversion materials for organic solar cells.However,it is difficult to realize the efficient design of BaF_(2...Eu^(2+)doped fluorosilicate glass-ceramics containing BaF_(2) nanocrystals have high potential as spectral conversion materials for organic solar cells.However,it is difficult to realize the efficient design of BaF_(2):Eu^(2+)doped fluorosilicate glass and to vividly observe the glass microstructure in experiment through traditional trial-and-error glass preparation method.BaF_(2):Eu^(2+)doped fluorosilicate glassceramics with high transparency,and high photoluminescence(PL)performance were predicted,designed and prepared via molecular dynamics(MD)simulation method.By MD simulation prediction,self-organized nanocrystallization was realized to inhibit the abnormal growth of nanocrystals due to[AlO_(4)]tetrahedra formed in the fluoride-oxide interface.The introduction of NaF reduces the effective phonon energy of the glass because Na+will prompt Al^(3+)to migrate from the fluoride phase to the silicate phase and interface.The local environment of Eu^(2+)is optimized by predicting the doping concentration of EuF_(3) and 2 mol%EuF3 is the best concentration in this work.Glass-ceramics sample GC2Eu as spectral conversion layer was successfully applied on organic solar cells to obtain more available visible phonons with a high photoelectric conversion efficiency(PCE).This work confirms the guidance of molecular dynamics simulation methods for fluorosilicate glasses design.展开更多
In the production process of silicone sealant,mineral oil is used to replace methyl silicone oil plasticizer in silicone sealant to reduce costs and increase efficiency.However,the silicone sealant content in mineral ...In the production process of silicone sealant,mineral oil is used to replace methyl silicone oil plasticizer in silicone sealant to reduce costs and increase efficiency.However,the silicone sealant content in mineral oil is prone to premature aging,which significantly reduces the mechanical properties of the silicone sealant and severely affects its service life.At the same time,there are few reports on the simulation research of the performance of silicone sealant.In this study,three mixed system models of crosslinking silicone sealant/plasticizer are constructed by the molecular dynamics simulationmethod,and the effect of three influencing factors,namely,crosslinking degree of silicone sealant,plasticizer content and external temperature on the mechanical properties of silicone sealant system is analyzed.The results show that at room temperature,the mechanical properties of the silicone sealant system are enhanced with the increase of its crosslinking degree;At a high crosslinking degree,with the increase of plasticizer content,themechanical properties of the silicone sealant system show an overall decreasing trend.When the methyl silicone oil in the range of 20%,themechanical properties of the silicone sealant appeared tobe a small degree of enhancement;As the temperature increases,the doped mineral oil mechanical properties of silicone sealant declined significantly,while doped with methyl silicone oil silicone sealant and doped with double-ended vinyl silicone oil silicone sealant mechanical properties have better heat resistance.It will provide scientific theoretical guidance for improving and predicting the mechanical properties of silicone sealant.展开更多
Recent advancements in nanotechnology have spotlighted the catalytic potential of nanozymes, particularly single-atom nanozymes(SANs), which are pivotal for innovations in biosensing and medical diagnostics. Among oth...Recent advancements in nanotechnology have spotlighted the catalytic potential of nanozymes, particularly single-atom nanozymes(SANs), which are pivotal for innovations in biosensing and medical diagnostics. Among others, DNA stands out as an ideal biological regulator. Its inherent programmability and interaction capabilities allow it to significantly modulate nanozyme activity. This study delves into the dynamic interplay between DNA and molybdenum-zinc single-atom nanozymes(Mo-Zn SANs). Using molecular dynamics simulations, we uncover how DNA influences the peroxidase-like activities of Mo-Zn SANs, providing a foundational understanding that broadens the application scope of SANs in biosensing.With these insights as a foundation, we developed and demonstrated a model aptasensor for point-ofcare testing(POCT), utilizing a label-free colorimetric approach that leverages DNA-nanozyme interactions to achieve high-sensitivity detection of lysozyme. Our work elucidates the nuanced control DNA exerts over nanozyme functionality and illustrates the application of this molecular mechanism through a smartphone-assisted biosensing platform. This study not only underscores the practical implications of DNA-regulated Mo-Zn SANs in enhancing biosensing platforms, but also highlights the potential of single-atom nanozyme technology to revolutionize diagnostic tools through its inherent versatility and sensitivity.展开更多
A super-radiant terahertz free-electron laser(THz-FEL)light source was developed for the first time in Thailand and Southeast Asia at the PBP-CMU Electron Linac Laboratory(PCELL)of Chiang Mai University.This radiation...A super-radiant terahertz free-electron laser(THz-FEL)light source was developed for the first time in Thailand and Southeast Asia at the PBP-CMU Electron Linac Laboratory(PCELL)of Chiang Mai University.This radiation source requires relatively ultrashort electron bunches to produce intense coherent THz pulses.Three electron bunch compression processes are utilized in the PCELL accelerator system comprising pre-bunch compression in an alpha magnet,velocity bunching in a radio-frequency(RF)linear accelerator(linac),and magnetic bunch compression in a 180°acromat system.Electron bunch compression in the magnetic compressor system poses considerable challenges,which are addressed through the use of three quadrupole doublets.The strengths of the quadrupole fields significantly influence the rotation of the beam line longitudinal phase space distribution along the bunch compressor.Start-to-end beam dynamics simulations using the ASTRA code were performed to optimize the electron beam properties for generating super-radiant THz-FEL radiation.The operational parameters considered in the simulations comprise the alpha magnet gradient,linac RF phase,and quadrupole field strengths.The optimization results show that 10-16MeV femtosecond electron bunches with a low energy spread(~0.2%),small normalized emittance(~15πmm·mrad),and high peak current(165-247A)can be produced by the PCELL accelerator system at the optimal parameters.A THz-FEL with sub-microjoule pulse energies can thus be obtained at the optimized electron beam parameters.The physical and conceptual design of the THz-FEL beamline were completed based on the beam dynamics simulation results.The construction and installation of this beamline are currently underway and expected to be completed by mid-2024.The commissioning of the beamline will then commence.展开更多
Eutectic high entropy alloys are noted for their excellent castability and comprehensive mechanical properties.The excellent mechanical properties are closely related to the activation and evolution of deformation mec...Eutectic high entropy alloys are noted for their excellent castability and comprehensive mechanical properties.The excellent mechanical properties are closely related to the activation and evolution of deformation mechanisms at the atomic scale.In this work,AlCoCrFeNi2.1 alloy is taken as the research object.The mechanical behaviors and deformation mechanisms of the FCC and B2 single crystals with different orientations and the FCC/B2 composites with K-S orientation relationship during nanoindentation processes are systematically studied by molecular dynamics simulations.The results show that the mechanical behaviors of FCC single crystals are significantly orientation-dependent,meanwhile,the indentation force of[110]single crystal is the lowest at the elastic-plastic transition point,and that for[100]single crystal is the lowest in plastic deformation stage.Compared with FCC,the stress for B2 single crystals at the elastic-plastic transition point is higher.However,more deformation systems such as stacking faults,twins and dislocation loops are activated in FCC single crystal during the plastic deformation process,resulting in higher indentation force.For composites,the flow stress increases with the increase of B2 phase thickness during the initial stage of deformation.When indenter penetrates heterogeneous interface,the significantly increased deformation system in FCC phase leads to a significant increase in indentation force.The mechanical behaviors and deformation mechanisms depend on the component single crystal.When the thickness of the component layer is less than 15 nm,the heterogeneous interfaces fail to prevent the dislocation slip and improve the indentation force.The results will enrich the plastic deformation mechanisms of multi-principal eutectic alloys and provide guidance for the design of nanocrystalline metallic materials.展开更多
Objective To evaluate the antibacterial potential of bioactive compounds from Persicaria hydropiper(L.)(P.hydropiper)against bacterial virulence proteins through molecular docking(MD)and experimental validation.Method...Objective To evaluate the antibacterial potential of bioactive compounds from Persicaria hydropiper(L.)(P.hydropiper)against bacterial virulence proteins through molecular docking(MD)and experimental validation.Methods Six bioactive compounds from P.hydropiper were investigated:catechin(CAT1),hyperin(HYP1),ombuin(OMB1),pinosylvin(PSV1),quercetin 3-sulfate(QSF1),and scutellarein(SCR1).Their binding affinities and potential binding pockets were assessed through MD against four bacterial target proteins with Protein Data Bank identifiers(PDB IDs):topoisomerase IV from Escherichia coli(E.coli)(PDB ID:3FV5),Staphylococcus aureus(S.aureus)gyrase ATPase binding domain(PDB ID:3U2K),CviR from Chromobacterium violaceum(C.violaceum)(PDB ID:3QP1),and glycosyl hydrolase from Pseudomonas aeruginosa(P.aeruginosa)(PDB ID:5BX9).Molecular dynamics simulations(MDS)were performed on the most promising compound-protein complexes for 50 nanoseconds(ns).Drug-likeness was evaluated using Lipinski's Rule of Five(RO5),followed by absorption,distribution,metabolism,excretion,and toxicity(ADMET)analysis using SwissADME and pkCSM web servers.Antibacterial activity was evaluated through disc diffusion assays,testing both individual compounds and combinations with conventional antibiotics[cefotaxime(CTX1,30μg/disc),ceftazidime(CAZ1,30μg/disc),and piperacillin(PIP1,100μg/disc)].Results MD revealed strong binding affinity(ranging from-9.3 to-5.9 kcal/mol)for all compounds,with CAT1 showing exceptional binding to 3QP1(-9.3 kcal/mol)and 5BX9(-8.4 kcal/mol).MDS confirmed the stability of CAT1-protein complexes with binding free energies of-84.71 kJ/mol(5BX9-CAT1)and-95.59 kJ/mol(3QP1-CAT1).Five compounds(CAT1,SCR1,PSV1,OMB1,and QSF1)complied with Lipinski's RO5 and showed favorable ADMET profiles.All compounds were non-carcinogenic,with CAT1 classified in the lowest toxicity class(VI).In antibacterial assays,CAT1 demonstrated significant activity against both gram-positive bacteria[Streptococcus pneumoniae(S.pneumoniae),S.aureus,and Bacillus cereus(B.cereus)][zone diameter of inhibition(ZDI):10-22 mm]and gram-negative bacteria[Acinetobacter baumannii(A.baumannii),E.coli,and P.aeruginosa](ZDI:14-27 mm).Synergistic effects were observed when CAT1 was combined with antibiotics and the growth inhibitory indices(GII)was 0.69-1.00.Conclusion P.hydropiper bioactive compounds,particularly CAT1,show promising antibacterial potential through multiple mechanisms,including direct inhibition of bacterial virulence proteins and synergistic activity with conventional antibiotics.The favorable pharmacological properties and low toxicity profiles support their potential development as therapeutic agents against bacterial infections.展开更多
Uranium–molybdenum(U–Mo) alloys are critical for nuclear power generation and propulsion because of their superior thermal conductivity, irradiation stability, and anti-swelling properties. This study explores the p...Uranium–molybdenum(U–Mo) alloys are critical for nuclear power generation and propulsion because of their superior thermal conductivity, irradiation stability, and anti-swelling properties. This study explores the plastic deformation mechanisms of γ-phase U–Mo alloys using molecular dynamics(MD) simulations. In the slip model, the generalized stacking fault energy(GSFE) and the modified Peierls–Nabarro(P–N) model are used to determine the competitive relationships among different slip systems. In the twinning model, the generalized plane fault energy(GPFE) is assessed to evaluate the competition between slip and twinning. The findings reveal that among the three slip systems, the {110}<111>slip system is preferentially activated, while in the {112}<111> system, twinning is favored over slip, as confirmed by MD tensile simulations conducted in various directions. Additionally, the impact of Mo content on deformation behavior is emphasized. Insights are provided for optimizing process conditions to avoid γ → α′′ transitions, thereby maintaining a higher proportion of γ-phase U–Mo alloys for practical applications.展开更多
In nature,cavitation bubbles typically appear in clusters,engaging in interactions that create a variety of dynamicmotion patterns.To better understand the behavior ofmultiple bubble collapses and the mechanisms of in...In nature,cavitation bubbles typically appear in clusters,engaging in interactions that create a variety of dynamicmotion patterns.To better understand the behavior ofmultiple bubble collapses and the mechanisms of interbubble interaction,this study employs molecular dynamics simulation combined with a coarse-grained force field.By focusing on collapsemorphology,local density,and pressure,it elucidates how the number and arrangement of bubbles influence the collapse process.The mechanisms behind inter-bubble interactions are also considered.The findings indicate that the collapse speed of unbounded bubbles located in lateral regions is greater than that of the bubbles in the center.Moreover,it is shown that asymmetrical bubble distributions lead to a shorter collapse time overall.展开更多
The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critica...The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critical speed points and maintain operational stability across the entire working speed range.However,the support stiffness and damping of rubber O-rings exhibit significant nonlinear frequency dependence.Conventional experimental methods for deriving equivalent stiffness and damping,based on the principle of the forced non-resonance method,require fabricating custom setups for each O-ring specification and conducting vibration tests at varying frequencies,resulting in low efficiency and high costs.This study proposes a hybrid simulation-experimental method for dynamic parameter identification.Firstly,the frequency-dependent dynamic parameters of a specific O-ring support system are experimentally obtained.Subsequently,a corresponding parametric finite element model is established to simulate and solve the equivalent elastic modulus and equivalent stiffness-damping coefficient of this O-ring support system.Ultimately,after iterative simulation,the simulated and experimental results achieve a 99.7%agreement.The parametric finite element model developed herein can directly simulate and inversely estimate frequency-dependent dynamic parameters for O-rings of different specifications but identical elastic modulus.展开更多
In the organic Rankine cycle,the refrigerant inevitably interacts with the lubricating oil.This study investigates the interaction mechanism between the fourth-generation refrigerant R1336mzz(Z)and the polyol ester(PO...In the organic Rankine cycle,the refrigerant inevitably interacts with the lubricating oil.This study investigates the interaction mechanism between the fourth-generation refrigerant R1336mzz(Z)and the polyol ester(POE)which is a representative component of the lubricating oil,using molecular dynamics simulations.The research focuses on pentaerythritol ester(PEC)with medium to long chain lengths,specifically PEC9.Relevant parameters such as solubility parameters,diffusion coefficients,binding energies,and radial distribution functions were calculated to elucidate the interaction dynamics.The variation in solubility parameters suggests that the miscibility of PEC9 and R1336mzz(Z)diminishes as the number of PEC9 chains increases.Additionally,the compatibility between these two components deteriorates with rising temperature,which is accompanied by a reduction in their binding energy.The simulation results presented in this study offer theoretical insights into the behavior of refrigerant R1336mzz(Z)upon contact with lubricating oil during actual operation,as well as implications for the operational efficiency of the equipment.展开更多
The strategic dispersion of carbon nanotubes(CNTs)within triblock copolymer matrix is key to fabricating nanocomposites with the desired electrical properties.This study investigated the self-assembly and electrical b...The strategic dispersion of carbon nanotubes(CNTs)within triblock copolymer matrix is key to fabricating nanocomposites with the desired electrical properties.This study investigated the self-assembly and electrical behavior of a polystyrene-polybutadiene-polystyrene(SBS)matrix with CNTs of different aspect ratios using hybrid particle-field molecular dynamics simulations.Structural factor analysis of the nanocomposites indicated that CNTs with higher aspect ratios promoted the transition of the SBS matrix from a bicontinuous to a lamellar phase.The resistor network algorithm method showed that the electrical conductivity of SBS and CNTs nanocomposites was influenced by the interplay between the CNTs aspect ratios,concentrations,and domain sizes of the triblock copolymer SBS.Our research sheds light on the relationship between CNTs dispersion and the electrical behavior of SBS/CNTs nanocomposites,guiding the engineering of materials to achieve desired electrical properties through the modulation of CNTs aspect ratios and tailored sizing of triblock copolymer domains.展开更多
Using molecular dynamics methods,simulations of collision cascades in polycrystalline tungsten(W)have been conducted in this study,including different primary-knock-on atom(PKA)directions,grain sizes,and PKA energies ...Using molecular dynamics methods,simulations of collision cascades in polycrystalline tungsten(W)have been conducted in this study,including different primary-knock-on atom(PKA)directions,grain sizes,and PKA energies between 1 keV and 150 keV.The results indicate that a smaller grain size leads to more defects forming in grain boundary regions during cascade processes.The impact of high-energy PKA may cause a certain degree of distortion of the grain boundaries,which has a higher probability in systems with smaller grain sizes and becomes more pronounced as the PKA energy increases.The direction of PKA can affect the formation and diffusion pathways of defects.When the PKA direction is perpendicular to the grain boundary,defects preferentially form near the grain boundary regions;by contrast,defects are more inclined to form in the interior of the grains.These results are of great significance for comprehending the changes in the performance of polycrystalline W under the high-energy fusion environments and can provide theoretical guidance for further optimization and application of W-based plasma materials.展开更多
The experiment explored the Fe_(2)O_(3) reduction process with H_(2)/CO mixed gas and confirmed a promoting effect from CO when its volume proportion in mixed gas is 20% at 850℃.The ReaxFF molecular dynamics(MD)simul...The experiment explored the Fe_(2)O_(3) reduction process with H_(2)/CO mixed gas and confirmed a promoting effect from CO when its volume proportion in mixed gas is 20% at 850℃.The ReaxFF molecular dynamics(MD)simulation method was used to observe the reduction process and provide an atomic-level explanation.The accuracy of the parameters used in the simulation was verified by the density functional theory(DFT)calculation.The simulation shows that the initial reduction rate of H_(2) is much faster than that of CO(from 800 to 950℃).As the reduction proceeds,cementite,obtained after CO participates in the reduction at 850℃,will appear on the iron surface.Due to the active properties of C atoms in cementite,they are easy to further react with the O atoms in Fe_(2)O_(3).The generation of internal CO may destroy the dense structure of the surface layer,thereby affecting the overall reduction swelling of Fe_(2)O_(3).However,excess CO is detrimental to the reaction rate,mainly because of the poor thermodynamic conditions of CO in the temperature range and the molecular diffusion capacity is not as good as that of H_(2).Furthermore,the surface structures obtained after H_(2) and CO reduction have been compared,and it was found that the structure obtained by CO reduction has a larger surface area,thus promoting the sub sequent reaction of H_(2).展开更多
基金Supported by the Specialized Research Fund for the Doctoral Program of Higher Education of China (200901431-10005)
文摘A dynamical simulation method is presented to model the cornea deformation in surgery of laser thermokeratoplasty. The virtual cornea is constructed as a mass-spring system. The corneal surface tension is simulated by damping spring stretch between mass points on the cornea model. The aqueous humor in the eyeball is modeled as ideal gas, and the intraocular pressure is simulated by gas pressure. The coagulation force is exerted on each photocoagulation spot to demonstrate its collapse caused by the condensation of corneal soft tissue irradiated by laser. An extra viscous drag force is added to each mass point to weaken the mass point oscillation. The use of the effective time-corrected Verlet integral method brings about flowing and stable dynamic simulation procedures. The simulation results show that, comparing to the undeformed model, the curvature of the region between the optical center and photocoagulation spot increases obviously. Moreover, the shape of the deformed virtual cornea is much similar to that of the real cornea after surgery.
基金supported by the National Natural Science Foundation of China (Grant No. 10474085)
文摘INS （Inelastic Neutron Scattering） spectrum of methane hydrate was measured on MARI （a direct-geometry chopper spectrometer） at Ruther-ford Appleton Laboratory, UK. Compared with ice Ih, it is found that the whole spectrum of methane hydrate moves toward high-energy by about 1.5 meV. Using lattice dynamical （LD） technique, computer simulations of methane hydrate were carried out. In the simulations, four potential models （BE TIP3P, TIP4P, MCY） were employed to calculate the phonon density of states （PDOS）. Comparing the calculated PDOS spectrum with the experimental spectrum, it is found that BF, TIP4P, and TIP3P potential lattices give out well-separated translational and librational bands while MCY potential lattice is unstable to do so and this model is not suitable to describe hydrate system.
基金Funded by the Research Funds of China University of Mining and Technology(No.102523215)。
文摘The pre-wetting of aggregate surface is a means to improve the interface performance of SBS modified asphalt and aggregate.The effect of pre-wetting technology on the interaction between SBS modified asphalt and aggregate was analyzed by molecular dynamics simulation.The diffusion coefficient and concentration distribution of SBS modified asphalt on aggregate surface are included.The simulation results show that the diffusion coefficient of the aggregate surface of SBS modified asphalt is increased by 47.6%and 70.5%respectively after 110#asphalt and 130#asphalt are pre-wetted.The concentration distribution of SBS modified asphalt on the aggregate surface after pre-wetting is more uniform.According to the results of interface energy calculation,the interface energy of SBS modified bitumen and aggregate can be increased by about 5%after pre-wetting.According to the results of molecular dynamics simulation,the pre-wetting technology can effectively improve the interface workability of SBS modified bitumen and aggregate,so as to improve the interface performance.
基金supported by the Priority Program SPP 1992 of the German Science Foundation(DFG)The Diversity of Exoplanets under project number 362460292.
文摘We calculate the electrical and thermal conductivity of hydrogen for a wide range of densities and temperatures by using molecular dynamics simulations informed by density functional theory.On the basis of the corresponding extended ab initio data set,we construct interpolation formulas covering the range from low-density,high-temperature to high-density,low-temperature plasmas.Our conductivity model repro-duces the well-known limits of the Spitzer and Ziman theory.We compare with available experimental data andfind very good agreement.The new conductivity model can be applied,for example,in dynamo simulations for magneticfield generation in gas giant planets,brown dwarfs,and stellar envelopes.
基金support provided by the National Natural Science Foundation of China(Nos.22173051,21829301,21774066),PCSIRT(IRT1257)the College Discipline Innovation and Intelligence Introduction Program(111 Project(B16027)+2 种基金the International Cooperation Base(No.2016D01025)Tianjin International Joint Research and Development Center)P.Zhang acknowledges the financial support provided by NSFC(No.22473024).
文摘Stimuli-responsive polymers capable of rapidly altering their chain conformation in response to external stimuli exhibit broad applica-tion prospects.Experiments have shown that pressure plays a pivotal role in regulating the microscopic chain conformation of polymers in mixed solvents,and one notable finding is that increasing the pressure can lead to the vanishing of the co-nonsolvency effect.However,the mecha-nisms underlying this phenomenon remain unclear.In this study,we systematically investigated the influence of pressure on the co-nonsolvency effect of single-chain and multi-chain homopolymers in binary mixed good-solvent systems using molecular dynamics simulations.Our results show that the co-nonsolvency-induced chain conformation transition and aggregation behavior significantly depend on pressure in allsingle-chain and multi-chain systems.In single-chain systems,at low pressures,the polymer chain maintains a collapsed state over a wide range of co-solvent fractions(x-range)owing to the co-nonsolvency effect.As the pressure increases,the x-range of the collapsed state gradually narrows,ac-companied by a progressive expansion of the chain.In multichain systems,polymer chains assemble into approximately spherical aggregates over a broad x-range at low pressures owing to the co-nonsolvency effect.Increasing the pressure reduces the x-range for forming aggregates and leads to the formation of loose aggregates or even to a state of dispersed chains at some x-range.These findings indicate that increasing the pressure can weaken or even offset the co-nonsolvency effect in some x-range,which is in good agreement with the experimental observations.Quantitative analysis of the radial density distributions and radial distribution functions reveals that,with increasing pressure,(1)the densities of both polymers and co-solvent molecules within aggregates decrease,while that of the solvent molecule increases;and(2)the effective interac-tions between the polymer and the co-solvent weaken,whereas those between the polymer and solvent strengthen.This enhances the incorpo-ration of solvent molecules within the chains,thereby weakening or even suppressing the chain aggregation.Our study not only elucidates the regulatory mechanism of pressure on the microscopic chain conformations and aggregation behaviors of polymers,but also may provide theo-retical guidance for designing smart polymericmaterials based on mixed solvents.
文摘Carbon nanotubes(CNTs),black phosphorus nanotubes(BPNTs),and graphene derivatives exhibit significant promise for applications in nano-electromechanical systems(NEMS),energy storage,and sensing technologies due to their exceptional mechanical,electrical,and thermal properties.This review summarizes recent advances in understanding the dynamic behaviors of these nanomaterials,with a particular focus on insights gained from molecular dynamics(MD)simulations.Key areas discussed include the oscillatory and rotational dynamics of double-walled CNTs,fabrication and stability challenges associated with BPNTs,and the emerging potential of graphyne nanotubes(GNTs).The review also outlines design strategies for enhancing nanodevice performance and underscores the importance of future efforts in experimental validation,multi-scale coupling analyses,and the development of novel nanocomposites to accelerate practical deployment.
基金Project supported by the National Natural Science Foundation of China(52172008,51872255)the Key Research and Development Project of Zhejiang Province(2021C01174)。
文摘Eu^(2+)doped fluorosilicate glass-ceramics containing BaF_(2) nanocrystals have high potential as spectral conversion materials for organic solar cells.However,it is difficult to realize the efficient design of BaF_(2):Eu^(2+)doped fluorosilicate glass and to vividly observe the glass microstructure in experiment through traditional trial-and-error glass preparation method.BaF_(2):Eu^(2+)doped fluorosilicate glassceramics with high transparency,and high photoluminescence(PL)performance were predicted,designed and prepared via molecular dynamics(MD)simulation method.By MD simulation prediction,self-organized nanocrystallization was realized to inhibit the abnormal growth of nanocrystals due to[AlO_(4)]tetrahedra formed in the fluoride-oxide interface.The introduction of NaF reduces the effective phonon energy of the glass because Na+will prompt Al^(3+)to migrate from the fluoride phase to the silicate phase and interface.The local environment of Eu^(2+)is optimized by predicting the doping concentration of EuF_(3) and 2 mol%EuF3 is the best concentration in this work.Glass-ceramics sample GC2Eu as spectral conversion layer was successfully applied on organic solar cells to obtain more available visible phonons with a high photoelectric conversion efficiency(PCE).This work confirms the guidance of molecular dynamics simulation methods for fluorosilicate glasses design.
基金supported by The Guangxi Scholarship Fund of Guangxi Education Department(GED),Guangxi Key Research and Development Project(Grant No.Guike AB24010217)the Major Special Project of Guangxi Science and Technology(GrantNo.Guike AA23062020)+1 种基金the Guangxi Science and Technology Base and Talent Project(Grant No.Guike AD20297016)the Guangxi Minzu University Startup Project for Talent Introduction in 2019(Grant No.2019KJQD11).
文摘In the production process of silicone sealant,mineral oil is used to replace methyl silicone oil plasticizer in silicone sealant to reduce costs and increase efficiency.However,the silicone sealant content in mineral oil is prone to premature aging,which significantly reduces the mechanical properties of the silicone sealant and severely affects its service life.At the same time,there are few reports on the simulation research of the performance of silicone sealant.In this study,three mixed system models of crosslinking silicone sealant/plasticizer are constructed by the molecular dynamics simulationmethod,and the effect of three influencing factors,namely,crosslinking degree of silicone sealant,plasticizer content and external temperature on the mechanical properties of silicone sealant system is analyzed.The results show that at room temperature,the mechanical properties of the silicone sealant system are enhanced with the increase of its crosslinking degree;At a high crosslinking degree,with the increase of plasticizer content,themechanical properties of the silicone sealant system show an overall decreasing trend.When the methyl silicone oil in the range of 20%,themechanical properties of the silicone sealant appeared tobe a small degree of enhancement;As the temperature increases,the doped mineral oil mechanical properties of silicone sealant declined significantly,while doped with methyl silicone oil silicone sealant and doped with double-ended vinyl silicone oil silicone sealant mechanical properties have better heat resistance.It will provide scientific theoretical guidance for improving and predicting the mechanical properties of silicone sealant.
基金supported by the Science and Technology Research Project from Education Department of Jilin Province (No. JJKH20231296KJ)the Natural Science Foundation of Science and Technology Department of Jilin Province (Joint Fund Project) (No. YDZJ202201ZYTS340)+9 种基金the Fundamental Research Funds for the Central Universities (No. 2412022ZD013)the Science and Technology Development Plan Project of Jilin Province (Nos. SKL202302030, SKL202402017, 20210204126YY, 20230204113YY, 20240602003RC, 20210402059GH)the National Natural Science Foundation of China (Nos. 22174137, 22322410, 92372102 and 22073094)the Cooperation Funding of Changchun with Chinese Academy of Sciences (No. 22SH13)the Capital Construction Fund Projects within the Budget of Jilin Province (No. 2023C042–5)the University Level Scientific Research Projects of Ordinary Universities in Xinjiang Uygur Autonomous Region (No. 2022YQSN002)the State Key Laboratory of Molecular Engineering of Polymers (Fudan University) (No. K2024–11)the Program for Young Scholars in Regional Development of CASthe essential support of the Network and Computing Center, CIAC, CASthe Computing Center of Jilin Province。
文摘Recent advancements in nanotechnology have spotlighted the catalytic potential of nanozymes, particularly single-atom nanozymes(SANs), which are pivotal for innovations in biosensing and medical diagnostics. Among others, DNA stands out as an ideal biological regulator. Its inherent programmability and interaction capabilities allow it to significantly modulate nanozyme activity. This study delves into the dynamic interplay between DNA and molybdenum-zinc single-atom nanozymes(Mo-Zn SANs). Using molecular dynamics simulations, we uncover how DNA influences the peroxidase-like activities of Mo-Zn SANs, providing a foundational understanding that broadens the application scope of SANs in biosensing.With these insights as a foundation, we developed and demonstrated a model aptasensor for point-ofcare testing(POCT), utilizing a label-free colorimetric approach that leverages DNA-nanozyme interactions to achieve high-sensitivity detection of lysozyme. Our work elucidates the nuanced control DNA exerts over nanozyme functionality and illustrates the application of this molecular mechanism through a smartphone-assisted biosensing platform. This study not only underscores the practical implications of DNA-regulated Mo-Zn SANs in enhancing biosensing platforms, but also highlights the potential of single-atom nanozyme technology to revolutionize diagnostic tools through its inherent versatility and sensitivity.
基金support from the NSRF via the Program Management Unit for Human Resources&Institutional Development,Research,and Innovation(No.B05F650022),as well as from Chiang Mai University.
文摘A super-radiant terahertz free-electron laser(THz-FEL)light source was developed for the first time in Thailand and Southeast Asia at the PBP-CMU Electron Linac Laboratory(PCELL)of Chiang Mai University.This radiation source requires relatively ultrashort electron bunches to produce intense coherent THz pulses.Three electron bunch compression processes are utilized in the PCELL accelerator system comprising pre-bunch compression in an alpha magnet,velocity bunching in a radio-frequency(RF)linear accelerator(linac),and magnetic bunch compression in a 180°acromat system.Electron bunch compression in the magnetic compressor system poses considerable challenges,which are addressed through the use of three quadrupole doublets.The strengths of the quadrupole fields significantly influence the rotation of the beam line longitudinal phase space distribution along the bunch compressor.Start-to-end beam dynamics simulations using the ASTRA code were performed to optimize the electron beam properties for generating super-radiant THz-FEL radiation.The operational parameters considered in the simulations comprise the alpha magnet gradient,linac RF phase,and quadrupole field strengths.The optimization results show that 10-16MeV femtosecond electron bunches with a low energy spread(~0.2%),small normalized emittance(~15πmm·mrad),and high peak current(165-247A)can be produced by the PCELL accelerator system at the optimal parameters.A THz-FEL with sub-microjoule pulse energies can thus be obtained at the optimized electron beam parameters.The physical and conceptual design of the THz-FEL beamline were completed based on the beam dynamics simulation results.The construction and installation of this beamline are currently underway and expected to be completed by mid-2024.The commissioning of the beamline will then commence.
基金supported by the Natural Science Foundation of Hebei Province(E2024209052)the Youth Scholars Promotion Plan of North China University of Science and Technology(QNTJ202307).
文摘Eutectic high entropy alloys are noted for their excellent castability and comprehensive mechanical properties.The excellent mechanical properties are closely related to the activation and evolution of deformation mechanisms at the atomic scale.In this work,AlCoCrFeNi2.1 alloy is taken as the research object.The mechanical behaviors and deformation mechanisms of the FCC and B2 single crystals with different orientations and the FCC/B2 composites with K-S orientation relationship during nanoindentation processes are systematically studied by molecular dynamics simulations.The results show that the mechanical behaviors of FCC single crystals are significantly orientation-dependent,meanwhile,the indentation force of[110]single crystal is the lowest at the elastic-plastic transition point,and that for[100]single crystal is the lowest in plastic deformation stage.Compared with FCC,the stress for B2 single crystals at the elastic-plastic transition point is higher.However,more deformation systems such as stacking faults,twins and dislocation loops are activated in FCC single crystal during the plastic deformation process,resulting in higher indentation force.For composites,the flow stress increases with the increase of B2 phase thickness during the initial stage of deformation.When indenter penetrates heterogeneous interface,the significantly increased deformation system in FCC phase leads to a significant increase in indentation force.The mechanical behaviors and deformation mechanisms depend on the component single crystal.When the thickness of the component layer is less than 15 nm,the heterogeneous interfaces fail to prevent the dislocation slip and improve the indentation force.The results will enrich the plastic deformation mechanisms of multi-principal eutectic alloys and provide guidance for the design of nanocrystalline metallic materials.
基金Research Grants of Senior Research Fellowship in favor of first author(Gloak Majumdar)from Council of Scientific and Industrial Research(CSIR,New Delhi,Government of India)(CSIR-SRF)with Award No.09/1151/(0008)2020-EMR-I.
文摘Objective To evaluate the antibacterial potential of bioactive compounds from Persicaria hydropiper(L.)(P.hydropiper)against bacterial virulence proteins through molecular docking(MD)and experimental validation.Methods Six bioactive compounds from P.hydropiper were investigated:catechin(CAT1),hyperin(HYP1),ombuin(OMB1),pinosylvin(PSV1),quercetin 3-sulfate(QSF1),and scutellarein(SCR1).Their binding affinities and potential binding pockets were assessed through MD against four bacterial target proteins with Protein Data Bank identifiers(PDB IDs):topoisomerase IV from Escherichia coli(E.coli)(PDB ID:3FV5),Staphylococcus aureus(S.aureus)gyrase ATPase binding domain(PDB ID:3U2K),CviR from Chromobacterium violaceum(C.violaceum)(PDB ID:3QP1),and glycosyl hydrolase from Pseudomonas aeruginosa(P.aeruginosa)(PDB ID:5BX9).Molecular dynamics simulations(MDS)were performed on the most promising compound-protein complexes for 50 nanoseconds(ns).Drug-likeness was evaluated using Lipinski's Rule of Five(RO5),followed by absorption,distribution,metabolism,excretion,and toxicity(ADMET)analysis using SwissADME and pkCSM web servers.Antibacterial activity was evaluated through disc diffusion assays,testing both individual compounds and combinations with conventional antibiotics[cefotaxime(CTX1,30μg/disc),ceftazidime(CAZ1,30μg/disc),and piperacillin(PIP1,100μg/disc)].Results MD revealed strong binding affinity(ranging from-9.3 to-5.9 kcal/mol)for all compounds,with CAT1 showing exceptional binding to 3QP1(-9.3 kcal/mol)and 5BX9(-8.4 kcal/mol).MDS confirmed the stability of CAT1-protein complexes with binding free energies of-84.71 kJ/mol(5BX9-CAT1)and-95.59 kJ/mol(3QP1-CAT1).Five compounds(CAT1,SCR1,PSV1,OMB1,and QSF1)complied with Lipinski's RO5 and showed favorable ADMET profiles.All compounds were non-carcinogenic,with CAT1 classified in the lowest toxicity class(VI).In antibacterial assays,CAT1 demonstrated significant activity against both gram-positive bacteria[Streptococcus pneumoniae(S.pneumoniae),S.aureus,and Bacillus cereus(B.cereus)][zone diameter of inhibition(ZDI):10-22 mm]and gram-negative bacteria[Acinetobacter baumannii(A.baumannii),E.coli,and P.aeruginosa](ZDI:14-27 mm).Synergistic effects were observed when CAT1 was combined with antibiotics and the growth inhibitory indices(GII)was 0.69-1.00.Conclusion P.hydropiper bioactive compounds,particularly CAT1,show promising antibacterial potential through multiple mechanisms,including direct inhibition of bacterial virulence proteins and synergistic activity with conventional antibiotics.The favorable pharmacological properties and low toxicity profiles support their potential development as therapeutic agents against bacterial infections.
基金Project supported by the National Natural Science Foundation of China (Grant No. 52271105)。
文摘Uranium–molybdenum(U–Mo) alloys are critical for nuclear power generation and propulsion because of their superior thermal conductivity, irradiation stability, and anti-swelling properties. This study explores the plastic deformation mechanisms of γ-phase U–Mo alloys using molecular dynamics(MD) simulations. In the slip model, the generalized stacking fault energy(GSFE) and the modified Peierls–Nabarro(P–N) model are used to determine the competitive relationships among different slip systems. In the twinning model, the generalized plane fault energy(GPFE) is assessed to evaluate the competition between slip and twinning. The findings reveal that among the three slip systems, the {110}<111>slip system is preferentially activated, while in the {112}<111> system, twinning is favored over slip, as confirmed by MD tensile simulations conducted in various directions. Additionally, the impact of Mo content on deformation behavior is emphasized. Insights are provided for optimizing process conditions to avoid γ → α′′ transitions, thereby maintaining a higher proportion of γ-phase U–Mo alloys for practical applications.
基金funded by the Natural Science Foundation of China[U20A20292]Shandong Province Science andTechnology SMES InnovationAbility Improvement Project[2023TSGC0005]China Postdoctoral Science Foundation[2024M752697].
文摘In nature,cavitation bubbles typically appear in clusters,engaging in interactions that create a variety of dynamicmotion patterns.To better understand the behavior ofmultiple bubble collapses and the mechanisms of interbubble interaction,this study employs molecular dynamics simulation combined with a coarse-grained force field.By focusing on collapsemorphology,local density,and pressure,it elucidates how the number and arrangement of bubbles influence the collapse process.The mechanisms behind inter-bubble interactions are also considered.The findings indicate that the collapse speed of unbounded bubbles located in lateral regions is greater than that of the bubbles in the center.Moreover,it is shown that asymmetrical bubble distributions lead to a shorter collapse time overall.
基金National Key R&D Program of China(No.2017YFB1304000)Fundamental Research Funds for the Central Universities,China(No.2232023G-05-1)。
文摘The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critical speed points and maintain operational stability across the entire working speed range.However,the support stiffness and damping of rubber O-rings exhibit significant nonlinear frequency dependence.Conventional experimental methods for deriving equivalent stiffness and damping,based on the principle of the forced non-resonance method,require fabricating custom setups for each O-ring specification and conducting vibration tests at varying frequencies,resulting in low efficiency and high costs.This study proposes a hybrid simulation-experimental method for dynamic parameter identification.Firstly,the frequency-dependent dynamic parameters of a specific O-ring support system are experimentally obtained.Subsequently,a corresponding parametric finite element model is established to simulate and solve the equivalent elastic modulus and equivalent stiffness-damping coefficient of this O-ring support system.Ultimately,after iterative simulation,the simulated and experimental results achieve a 99.7%agreement.The parametric finite element model developed herein can directly simulate and inversely estimate frequency-dependent dynamic parameters for O-rings of different specifications but identical elastic modulus.
基金supported by Hainan Provincial Natural Science Foundation of China(No.422CXTD509).
文摘In the organic Rankine cycle,the refrigerant inevitably interacts with the lubricating oil.This study investigates the interaction mechanism between the fourth-generation refrigerant R1336mzz(Z)and the polyol ester(POE)which is a representative component of the lubricating oil,using molecular dynamics simulations.The research focuses on pentaerythritol ester(PEC)with medium to long chain lengths,specifically PEC9.Relevant parameters such as solubility parameters,diffusion coefficients,binding energies,and radial distribution functions were calculated to elucidate the interaction dynamics.The variation in solubility parameters suggests that the miscibility of PEC9 and R1336mzz(Z)diminishes as the number of PEC9 chains increases.Additionally,the compatibility between these two components deteriorates with rising temperature,which is accompanied by a reduction in their binding energy.The simulation results presented in this study offer theoretical insights into the behavior of refrigerant R1336mzz(Z)upon contact with lubricating oil during actual operation,as well as implications for the operational efficiency of the equipment.
基金financially supported by the National Natural Science Foundation of China(Nos.52273019,62173065,22133002,22273031,and 12274056)Fundamental Research Funds for the Central Universities(No.04442024074)+2 种基金NationalKey R&D Program of China(No.2022YFB3707300)Beijing Natural Science Foundation(No.4242040)Scientific Research Funds Project of Liaoning Provincial Department of Education(No.LJKZ0034)。
文摘The strategic dispersion of carbon nanotubes(CNTs)within triblock copolymer matrix is key to fabricating nanocomposites with the desired electrical properties.This study investigated the self-assembly and electrical behavior of a polystyrene-polybutadiene-polystyrene(SBS)matrix with CNTs of different aspect ratios using hybrid particle-field molecular dynamics simulations.Structural factor analysis of the nanocomposites indicated that CNTs with higher aspect ratios promoted the transition of the SBS matrix from a bicontinuous to a lamellar phase.The resistor network algorithm method showed that the electrical conductivity of SBS and CNTs nanocomposites was influenced by the interplay between the CNTs aspect ratios,concentrations,and domain sizes of the triblock copolymer SBS.Our research sheds light on the relationship between CNTs dispersion and the electrical behavior of SBS/CNTs nanocomposites,guiding the engineering of materials to achieve desired electrical properties through the modulation of CNTs aspect ratios and tailored sizing of triblock copolymer domains.
基金Project supported by the National MCF Energy Research and Development Program of China(Grant No.2018YFE0308101)the National Key Research and Development Program of China(Grant No.2018YFB0704000)+1 种基金the Suqian Science and Technology Program(Grant No.K202337)the Natural Science Foundation of the Jiangsu Higher Education Institutions of China(Grant No.23KJD490001).
文摘Using molecular dynamics methods,simulations of collision cascades in polycrystalline tungsten(W)have been conducted in this study,including different primary-knock-on atom(PKA)directions,grain sizes,and PKA energies between 1 keV and 150 keV.The results indicate that a smaller grain size leads to more defects forming in grain boundary regions during cascade processes.The impact of high-energy PKA may cause a certain degree of distortion of the grain boundaries,which has a higher probability in systems with smaller grain sizes and becomes more pronounced as the PKA energy increases.The direction of PKA can affect the formation and diffusion pathways of defects.When the PKA direction is perpendicular to the grain boundary,defects preferentially form near the grain boundary regions;by contrast,defects are more inclined to form in the interior of the grains.These results are of great significance for comprehending the changes in the performance of polycrystalline W under the high-energy fusion environments and can provide theoretical guidance for further optimization and application of W-based plasma materials.
基金financial support from the National Natural Science Foundation of China(Nos.52204335 and 52374319)the National Nature Science Foundation of China(No.52174291)the Central Universities Foundation of China(No.06500170)。
文摘The experiment explored the Fe_(2)O_(3) reduction process with H_(2)/CO mixed gas and confirmed a promoting effect from CO when its volume proportion in mixed gas is 20% at 850℃.The ReaxFF molecular dynamics(MD)simulation method was used to observe the reduction process and provide an atomic-level explanation.The accuracy of the parameters used in the simulation was verified by the density functional theory(DFT)calculation.The simulation shows that the initial reduction rate of H_(2) is much faster than that of CO(from 800 to 950℃).As the reduction proceeds,cementite,obtained after CO participates in the reduction at 850℃,will appear on the iron surface.Due to the active properties of C atoms in cementite,they are easy to further react with the O atoms in Fe_(2)O_(3).The generation of internal CO may destroy the dense structure of the surface layer,thereby affecting the overall reduction swelling of Fe_(2)O_(3).However,excess CO is detrimental to the reaction rate,mainly because of the poor thermodynamic conditions of CO in the temperature range and the molecular diffusion capacity is not as good as that of H_(2).Furthermore,the surface structures obtained after H_(2) and CO reduction have been compared,and it was found that the structure obtained by CO reduction has a larger surface area,thus promoting the sub sequent reaction of H_(2).
