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.展开更多
Explorations into new electrolytes have highlighted the critical impact of solvation structure on battery performance,Classical molecular dynamics(CMD)using semi-empirical force fields has become an essential tool for...Explorations into new electrolytes have highlighted the critical impact of solvation structure on battery performance,Classical molecular dynamics(CMD)using semi-empirical force fields has become an essential tool for simulating solvation structures.However,mainstream force fields often lack accuracy in describing strong ion-solvent interactions,causing disparities between CMD simulations and experimental observations.Although some empirical methods have been employed in some of the studies to address this issue,their effectiveness has been limited.Our CMD research,supported by quantum chemical calculations and experimental data,reveals that the solvation structure is influenced not only by the charge model but also by the polarization description.Previous empirical approaches that focused solely on adjusting ion-solvent interaction strengths overlooked the importance of polarization effects.Building on this insight,we propose integrating the Drude polarization model into mainstream force fields and verify its feasibility in carbonate,ether,and nitrile electrolytes.Our experimental results demonstrate that this approach significantly enhances the accuracy of CMD-simulated solvation structures.This work is expected to provide a more reliable CMD method for electrolyte design,shielding researchers from the pitfalls of erroneous simulation outcomes.展开更多
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.展开更多
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.展开更多
The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope ...The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope and image processing technology were employed to achieve a quantitative grain size distribution of NiTi alloys so as to provide experimental data for molecular dynamics modeling at the atomic scale.Considering the size effect of molecular dynamics model on material properties,a reasonable modeling size was provided by taking into account three characteristic dimensions from the perspective of macro,meso,and micro scales according to the Buckinghamπtheorem.Then,the corresponding MD simulation on deformation and fracture behavior was investigated to derive a parameterized traction-separation(T-S)law,and then it was embedded into cohesive elements of finite element software.Thus,the crack propagation behavior in NiTi alloys was reproduced by the finite element method(FEM).The experimental results show that the predicted initiation fracture toughness is in good agreement with experimental data.In addition,it is found that the dynamics initiation fracture toughness increases with decreasing grain size and increasing loading velocity.展开更多
The effects of different Al_(2)O_(3)/SiO_(2)(Al/Si)ratios on the structure and tensile strength of Na_(2)O-CaO-MgO-Al_(2)O_(3)-SiO_(2)glass fiber were investigated by Raman,tensile strength tests and molecular dynamic...The effects of different Al_(2)O_(3)/SiO_(2)(Al/Si)ratios on the structure and tensile strength of Na_(2)O-CaO-MgO-Al_(2)O_(3)-SiO_(2)glass fiber were investigated by Raman,tensile strength tests and molecular dynamics simulation.The results showed that Al^(3+)mainly existed in the form of[AlO_(4)]within the glass network.With the increase of Al/Si ratio,the Si-O-Al linkage gradually became the main connection mode of glass network.The increase of bridging oxygen content and variation of Q^(n) indicated that a higher degree of network polymerization was formed.The tensile strength of the glass fibers obtained through experiments increased from 2653.56 to 2856.83 MPa,which was confirmed by the corresponding molecular dynamics simulation.During the stretching process,the Si-O bonds in the Si-O-Al linkage tended to break regardless of the compositional changes,and the increase of fractured Si-O-Al and Al-O-Al linkage absorbed more energy to resist the destroy.展开更多
We present a general numerical simulation method to solve non-Markovian dynamics of an open quantum system influenced by quantum Brownian motion.Based on the determined memory kernel function,this method enables the r...We present a general numerical simulation method to solve non-Markovian dynamics of an open quantum system influenced by quantum Brownian motion.Based on the determined memory kernel function,this method enables the resolution of non-Markovian dynamics for a wide range of system Hamiltonians and spectral densities.The system dynamics are described by exact integro-differential operator equations without any common approximations and they are simulated in this work by definite-number equations with stochastic initial conditions.This approach ensures the applicability of mature numerical methods and maintains computational complexity that remains largely invariant,even when dealing with more complex models.The high accuracy of our simulation is evident from a comparison with the results obtained from corresponding exact master equations,underscoring the reliability and precision of our method.展开更多
The formation of donut-shaped penetration pore upon membrane fusion in a closed lipid membrane system is of biological significance,since such the structures extensively exist in living body with various functions.How...The formation of donut-shaped penetration pore upon membrane fusion in a closed lipid membrane system is of biological significance,since such the structures extensively exist in living body with various functions.However,the related formation dynamics is unclear because of the limitation of experimental techniques.This work developed a new model of intra-vesicular fusion to elaborate the formation and stabilization of penetration pores by employing molecular dynamics simulations,based on simplified spherical lipid vesicle system,and investigated the regulation of membrane lipid composition.Results showed that penetration pore could be successfully formed based on the strategy of membrane fusion.The ease of intra-vesicular fusion and penetration pore formation was closely correlated with the lipid curvature properties,where negative spontaneous curvature of lipids seemed to be unfavorable for intra-vesicle fusion.Furthermore,the inner membrane tension around the pore was much larger than other regions,which governed the penetration pore size and stability.This work provided basic understanding for vesicle penetration pore formation and stabilization mechanisms.展开更多
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.展开更多
The effects of temperature and Re content on the mechanical properties,dislocation morphology,and deformation mechanism of γ-γ′phases nickel-based single crystal superalloys are investigated by using the molecular ...The effects of temperature and Re content on the mechanical properties,dislocation morphology,and deformation mechanism of γ-γ′phases nickel-based single crystal superalloys are investigated by using the molecular dynamics method through the model of γ-γ′phases containing hole defect.The addition of Re makes the dislocation distribution tend towards the γ phase.The higher the Re content,the earlier theγphase yields,while the γ′phase yields later.Dislocation bends under the combined action of the applied force and the resistance of the Re atoms to form a bend point.The Re atoms are located at the bend points and strengthen the alloy by fixing the dislocation and preventing it from cutting the γ′phase.Dislocations nucleate first in the γ phase,causing theγphase to deform plastically before the γ′phase.As the strain increases,the dislocation length first remains unchanged,then increases rapidly,and finally fluctuates and changes.The dislocation lengths in the γ phase are larger than those in the γ′phase at different temperatures.The dislocation length shows a decreasing tendency with the increase of the temperature.Temperature can affect movement of the dislocation,and superalloys have different plastic deformation mechanisms at low,medium and high temperatures.展开更多
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.展开更多
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.展开更多
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.展开更多
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 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 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).展开更多
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.展开更多
The hybrid CO_(2) thermal technique has achieved considerable success globally in extracting residual heavy oil from reserves following a long-term steam stimulation process.Using microscopic visualization experiments...The hybrid CO_(2) thermal technique has achieved considerable success globally in extracting residual heavy oil from reserves following a long-term steam stimulation process.Using microscopic visualization experiments and molecular dynamics(MD)simulations,this study investigates the microscopic enhanced oil recovery(EOR)mechanisms underlying residual oil removal using hybrid CO_(2) thermal systems.Based on the experimental models for the occurrence of heavy oil,this study evaluates the performance of hybrid CO_(2) thermal systems under various conditions using MD simulations.The results demonstrate that introducing CO_(2) molecules into heavy oil can effectively penetrate and decompose dense aggregates that are originally formed on hydrophobic surfaces.A stable miscible hybrid CO_(2) thermal system,with a high effective distribution ratio of CO_(2),proficiently reduces the interaction energies between heavy oil and rock surfaces,as well as within heavy oil.A visualization analysis of the interactions reveals that strong van der Waals(vdW)attractions occur between CO_(2) and heavy oil molecules,effectively promoting the decomposition and swelling of heavy oil.This unlocks the residual oil on the hydrophobic surfaces.Considering the impacts of temperature and CO_(2) concentration,an optimal gas-to-steam injection ratio(here,the CO_(2):steam ratio)ranging between 1:6 and 1:9 is recommended.This study examines the microscopic mechanisms underlying the hybrid CO_(2) thermal technique at a molecular scale,providing a significant theoretical guide for its expanded application in EOR.展开更多
CrCoNi medium entropy alloy(MEA)fabricated by laser powder bed fusion(LPBF)benefits from its distinctive hierarchical microstructure and has great potential as a structural material.However,while the intriguing chemic...CrCoNi medium entropy alloy(MEA)fabricated by laser powder bed fusion(LPBF)benefits from its distinctive hierarchical microstructure and has great potential as a structural material.However,while the intriguing chemical short-range order(CSRO)widely exists in high/medium entropy alloys,its formation in the LPBF-built samples still lacks enough understanding.In this study,we verified its existence by fine transmission electron microscopy characterizations and utilized hybrid Monte Carlo/molecular dynamics simulations to investigate the features and effects of CSRO in LPBF-built CrCoNi MEA(AM model).Results showed that the CSRO fraction and the stacking fault energy of the AM model lie between those of the well-annealed and random solid solution counterparts.Among these models,the AM model exhibited the best strain hardening ability due to its highest capability to generate and store sessile dislocations.The results agreed well with existing data and provide guidance to the future development of LPBF-built CrCoNi MEA.展开更多
基金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.
基金supported by the Science and Technology Project of State Grid Corporation of China(5419-202199552A-0-5-ZN).
文摘Explorations into new electrolytes have highlighted the critical impact of solvation structure on battery performance,Classical molecular dynamics(CMD)using semi-empirical force fields has become an essential tool for simulating solvation structures.However,mainstream force fields often lack accuracy in describing strong ion-solvent interactions,causing disparities between CMD simulations and experimental observations.Although some empirical methods have been employed in some of the studies to address this issue,their effectiveness has been limited.Our CMD research,supported by quantum chemical calculations and experimental data,reveals that the solvation structure is influenced not only by the charge model but also by the polarization description.Previous empirical approaches that focused solely on adjusting ion-solvent interaction strengths overlooked the importance of polarization effects.Building on this insight,we propose integrating the Drude polarization model into mainstream force fields and verify its feasibility in carbonate,ether,and nitrile electrolytes.Our experimental results demonstrate that this approach significantly enhances the accuracy of CMD-simulated solvation structures.This work is expected to provide a more reliable CMD method for electrolyte design,shielding researchers from the pitfalls of erroneous simulation outcomes.
基金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.
基金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.
基金Funded by the National Natural Science Foundation of China Academy of Engineering Physics and Jointly Setup"NSAF"Joint Fund(No.U1430119)。
文摘The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope and image processing technology were employed to achieve a quantitative grain size distribution of NiTi alloys so as to provide experimental data for molecular dynamics modeling at the atomic scale.Considering the size effect of molecular dynamics model on material properties,a reasonable modeling size was provided by taking into account three characteristic dimensions from the perspective of macro,meso,and micro scales according to the Buckinghamπtheorem.Then,the corresponding MD simulation on deformation and fracture behavior was investigated to derive a parameterized traction-separation(T-S)law,and then it was embedded into cohesive elements of finite element software.Thus,the crack propagation behavior in NiTi alloys was reproduced by the finite element method(FEM).The experimental results show that the predicted initiation fracture toughness is in good agreement with experimental data.In addition,it is found that the dynamics initiation fracture toughness increases with decreasing grain size and increasing loading velocity.
基金Funded by National Natural Science Foundation of China(No.52172019)Shandong Provincial Youth Innovation Team Development Plan of Colleges and Universities(No.2022KJ100)。
文摘The effects of different Al_(2)O_(3)/SiO_(2)(Al/Si)ratios on the structure and tensile strength of Na_(2)O-CaO-MgO-Al_(2)O_(3)-SiO_(2)glass fiber were investigated by Raman,tensile strength tests and molecular dynamics simulation.The results showed that Al^(3+)mainly existed in the form of[AlO_(4)]within the glass network.With the increase of Al/Si ratio,the Si-O-Al linkage gradually became the main connection mode of glass network.The increase of bridging oxygen content and variation of Q^(n) indicated that a higher degree of network polymerization was formed.The tensile strength of the glass fibers obtained through experiments increased from 2653.56 to 2856.83 MPa,which was confirmed by the corresponding molecular dynamics simulation.During the stretching process,the Si-O bonds in the Si-O-Al linkage tended to break regardless of the compositional changes,and the increase of fractured Si-O-Al and Al-O-Al linkage absorbed more energy to resist the destroy.
基金supported by the National Natural Science Foundation of China(Grant No.12304389,12274053)the Science and Technology Research Project of Xiamen University of Technology(Grant No.YKJ19025R)the Scientific Research Foundation of NEU(Grant No.01270021920501*115)。
文摘We present a general numerical simulation method to solve non-Markovian dynamics of an open quantum system influenced by quantum Brownian motion.Based on the determined memory kernel function,this method enables the resolution of non-Markovian dynamics for a wide range of system Hamiltonians and spectral densities.The system dynamics are described by exact integro-differential operator equations without any common approximations and they are simulated in this work by definite-number equations with stochastic initial conditions.This approach ensures the applicability of mature numerical methods and maintains computational complexity that remains largely invariant,even when dealing with more complex models.The high accuracy of our simulation is evident from a comparison with the results obtained from corresponding exact master equations,underscoring the reliability and precision of our method.
基金supported by the National Natural Science Foundation of China(Grants Nos.T2394512,32130061,and 12172366)the Scientific Instrument Developing Project of the Chinese Academy of Sciences(Grant No.GJJSTD20220002).
文摘The formation of donut-shaped penetration pore upon membrane fusion in a closed lipid membrane system is of biological significance,since such the structures extensively exist in living body with various functions.However,the related formation dynamics is unclear because of the limitation of experimental techniques.This work developed a new model of intra-vesicular fusion to elaborate the formation and stabilization of penetration pores by employing molecular dynamics simulations,based on simplified spherical lipid vesicle system,and investigated the regulation of membrane lipid composition.Results showed that penetration pore could be successfully formed based on the strategy of membrane fusion.The ease of intra-vesicular fusion and penetration pore formation was closely correlated with the lipid curvature properties,where negative spontaneous curvature of lipids seemed to be unfavorable for intra-vesicle fusion.Furthermore,the inner membrane tension around the pore was much larger than other regions,which governed the penetration pore size and stability.This work provided basic understanding for vesicle penetration pore formation and stabilization mechanisms.
基金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.
基金Project supported by the Xi’an Science and Technology Plan Project of Shaanxi Province of China(Grant No.23GXFW0086).
文摘The effects of temperature and Re content on the mechanical properties,dislocation morphology,and deformation mechanism of γ-γ′phases nickel-based single crystal superalloys are investigated by using the molecular dynamics method through the model of γ-γ′phases containing hole defect.The addition of Re makes the dislocation distribution tend towards the γ phase.The higher the Re content,the earlier theγphase yields,while the γ′phase yields later.Dislocation bends under the combined action of the applied force and the resistance of the Re atoms to form a bend point.The Re atoms are located at the bend points and strengthen the alloy by fixing the dislocation and preventing it from cutting the γ′phase.Dislocations nucleate first in the γ phase,causing theγphase to deform plastically before the γ′phase.As the strain increases,the dislocation length first remains unchanged,then increases rapidly,and finally fluctuates and changes.The dislocation lengths in the γ phase are larger than those in the γ′phase at different temperatures.The dislocation length shows a decreasing tendency with the increase of the temperature.Temperature can affect movement of the dislocation,and superalloys have different plastic deformation mechanisms at low,medium and high temperatures.
基金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.
基金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.
基金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.
基金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.
基金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.
基金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).
基金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.
基金financially supported by the National Natural Science Foundation of China(No.U20B6003)the China Scholarship Council(No.202306440015)a project of the China Petroleum&Chemical Corporation(No.P22174)。
文摘The hybrid CO_(2) thermal technique has achieved considerable success globally in extracting residual heavy oil from reserves following a long-term steam stimulation process.Using microscopic visualization experiments and molecular dynamics(MD)simulations,this study investigates the microscopic enhanced oil recovery(EOR)mechanisms underlying residual oil removal using hybrid CO_(2) thermal systems.Based on the experimental models for the occurrence of heavy oil,this study evaluates the performance of hybrid CO_(2) thermal systems under various conditions using MD simulations.The results demonstrate that introducing CO_(2) molecules into heavy oil can effectively penetrate and decompose dense aggregates that are originally formed on hydrophobic surfaces.A stable miscible hybrid CO_(2) thermal system,with a high effective distribution ratio of CO_(2),proficiently reduces the interaction energies between heavy oil and rock surfaces,as well as within heavy oil.A visualization analysis of the interactions reveals that strong van der Waals(vdW)attractions occur between CO_(2) and heavy oil molecules,effectively promoting the decomposition and swelling of heavy oil.This unlocks the residual oil on the hydrophobic surfaces.Considering the impacts of temperature and CO_(2) concentration,an optimal gas-to-steam injection ratio(here,the CO_(2):steam ratio)ranging between 1:6 and 1:9 is recommended.This study examines the microscopic mechanisms underlying the hybrid CO_(2) thermal technique at a molecular scale,providing a significant theoretical guide for its expanded application in EOR.
基金financially supported by the National Key R&D Program of China(No.2022YFB4602102)the National Natural Science Foundation of China(Grant No.51971144)the Natural Science Foundation of Shanghai(Grant No.19ZR1425200)。
文摘CrCoNi medium entropy alloy(MEA)fabricated by laser powder bed fusion(LPBF)benefits from its distinctive hierarchical microstructure and has great potential as a structural material.However,while the intriguing chemical short-range order(CSRO)widely exists in high/medium entropy alloys,its formation in the LPBF-built samples still lacks enough understanding.In this study,we verified its existence by fine transmission electron microscopy characterizations and utilized hybrid Monte Carlo/molecular dynamics simulations to investigate the features and effects of CSRO in LPBF-built CrCoNi MEA(AM model).Results showed that the CSRO fraction and the stacking fault energy of the AM model lie between those of the well-annealed and random solid solution counterparts.Among these models,the AM model exhibited the best strain hardening ability due to its highest capability to generate and store sessile dislocations.The results agreed well with existing data and provide guidance to the future development of LPBF-built CrCoNi MEA.
