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.展开更多
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 this study,circular dichroism(CD)and molecular dynamics(MD)simulation were used to investigate the thermal unfolding pathway of staphylococcal enterotoxin B(SEB)at temperatures of 298–371 and 298–500 K,and the re...In this study,circular dichroism(CD)and molecular dynamics(MD)simulation were used to investigate the thermal unfolding pathway of staphylococcal enterotoxin B(SEB)at temperatures of 298–371 and 298–500 K,and the relationship between the experimental and simulation results were explored.Our computational findings on the secondary structure of SEB showed that at room temperature,the CD spectroscopic results were highly consistent with the MD results.Moreover,under heating conditions,the changing trends of helix,sheet and random coil obtained by CD spectral fitting were highly consistent with those obtained by MD.In order to gain a deeper understanding of the thermal stability mechanism of SEB,the MD trajectories were analyzed in terms of root mean square deviation(RMSD),secondary structure assignment(SSA),radius of gyration(R_(g)),free energy surfaces(FES),solvent-accessible surface area(SASA),hydrogen bonds and salt bridges.The results showed that at low heating temperature,domain Ⅰ without loops(omitting the mobile loop region)mainly relied on hydrophobic interaction to maintain its thermal stability,whereas the thermal stability of domain Ⅱ was mainly controlled by salt bridges and hydrogen bonds.Under high heating temperature conditions,the hydrophobic interactions in domain Ⅰ without loops were destroyed and the secondary structure was almost completely lost,while domain Ⅱ could still rely on salt bridges as molecular staples to barely maintain the stability of the secondary structure.These results help us to understand the thermodynamic and kinetic mechanisms that maintain the thermal stability of SEB at the molecular level,and provide a direction for establishing safer and more effective food sterilization processes.展开更多
Nitrogen vacancy(NV)color centers in diamond have useful applications in quantum sensing andfluorescent marking.They can be gen-erated experimentally by ion implantation,femtosecond lasers,and chemical vapor deposition...Nitrogen vacancy(NV)color centers in diamond have useful applications in quantum sensing andfluorescent marking.They can be gen-erated experimentally by ion implantation,femtosecond lasers,and chemical vapor deposition.However,there is a lack of studies of the yield of NV color centers at the atomic scale.In the molecular dynamics simulations described in this paper,NV color centers are pre-pared by ion implantation in diamond with pre-doped nitrogen and subsequent annealing.The differences between the yields of NV color centers produced by implantation of carbon(C)and nitrogen(N)ions,respectively,are investigated.It is found that C-ion implantation gives a greater yield of NV color centers and superior location accuracy.The effects of different pre-doping concentrations(400–1500 ppm)and implantation energies(1.0–3.0 keV)on the NV color center yield are analyzed,and it is shown that a pre-doping concentra-tion of 1000 ppm with 2 keV C-ion implantation can produce a 13%yield of NV color centers after 1600 K annealing for 7.4 ns.Finally,a brief comparison of the NV color center identification methods is presented,and it is found that the error rate of an analysis utiliz-ing the identify diamond structure coordination analysis method is reduced by about 7%compared with conventional identification+methods.展开更多
Grain boundaries(GBs)play a crucial role on the structural stability and mechanical properties of Cu and its alloys.In this work,molecular dynamics(MD)simulations are employed to study the effects of Fe solutes on the...Grain boundaries(GBs)play a crucial role on the structural stability and mechanical properties of Cu and its alloys.In this work,molecular dynamics(MD)simulations are employed to study the effects of Fe solutes on the formation energy,excess volume,dislocations and melting behaviors of GBs in CuFe alloys.It is illustrated that Fe solute affects the structural stability of Cu GBs substantially,the formation energy of GBs is reduced,but the thickness and melting point of GBs are increased,that is,the structural stability of Cu GBs is significantly improved owing to the Fe solutes.A strong scaling law exists between the formation energy,excess volume,thickness and melting point of GBs.Therefore,Fe solid solute plays an important role in the characteristics of GBs in bi-crystal Cu.展开更多
Recently,there has been a growing prevalence in the utilization of graphdiyne(GDY)in the field of biomedicine,attributed to its distinctive physical structure and chemical properties.Additionally,its biocompatibility ...Recently,there has been a growing prevalence in the utilization of graphdiyne(GDY)in the field of biomedicine,attributed to its distinctive physical structure and chemical properties.Additionally,its biocompatibility has garnered increasing attention.However,there is a lack of research on the biological effects and physical mechanisms of GDYprotein interactions at the molecular scale.In this study,the villin headpiece subdomain(HP35)served as a representative protein model.Molecular dynamics simulations were employed to investigate the interaction process between the HP35 protein and GDY,as well as the structural evolution of the protein.The data presented in our study demonstrate that GDY can rapidly adsorb HP35 protein and induce denaturation to one of the a-helix structures of HP35 protein.This implies a potential cytotoxicity concern of GDY for biological systems.Compared to graphene,GDY induced less disruption to HP35 protein.This can be attributed to the presence of natural triangular vacancies in GDY,which prevents p–p stacking action and the limited interaction of GDY with HP35 protein is not conducive to the expansion of protein structures.These findings unveil the biological effects of GDY at the molecular level and provide valuable insights for the application of GDY in biomedicine.展开更多
Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear...Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear-thinning fluid in a microchannel.We validated the feasibility of our simulation method by evaluating the mean square displacement and Reynolds number of the solution layers.The results show that the change rule of the fluid system's velocity profile and interaction energy can reflect the shear-thinning characteristics of the fluids.The velocity profile resembles a top-hat shape,intensifying as the fluid's power law index decreases.The interaction energy between the wall and the fluid decreases gradually with increasing velocity,and a high concentration of non-Newtonian fluid reaches a plateau sooner.Moreover,the velocity profile of the fluid is related to the molecule number density distribution and their values are inversely proportional.By analyzing the radial distribution function,we found that the hydrogen bonds between solute and water molecules weaken with the increase in velocity.This observation offers an explanation for the shear-thinning phenomenon of the non-Newtonian flow from a micro perspective.展开更多
Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectil...Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectile impacts is thoroughly investigated using full-atomic molecular dynamics(MD)simulations.The study results show that the porous structure and density are key factors determining the mechanical response of GA under impact loading.Specifically,the impact-induced penetration of the projectile leads to the collapse of the pore structure,causing stretching and subsequent rupture of covalent bonds in graphene sheets.Moreover,the effects of temperature on the mechanical performance of GA have been proven to be minimal,thereby highlighting the mechanical stability of GA over a wide range of temperatures.Finally,the energy absorption density(EAD)and energy absorption efficiency(EAE)metrics are adopted to assess the energy absorption capacity of GA during projectile penetration.The research findings of this work demonstrate the significant potential of GA for energy absorption applications.展开更多
Studying the relationship between ionic interactions and salt solubility in seawater has implications for seawater desalination and mineral extraction.In this paper,a new method of expressing ion-to-ion interaction is...Studying the relationship between ionic interactions and salt solubility in seawater has implications for seawater desalination and mineral extraction.In this paper,a new method of expressing ion-to-ion interaction is proposed by using molecular dynamics simulation,and the relationship between ion-to-ion interaction and salt solubility in a simulated seawater water-salt system is investigated.By analyzing the variation of distance and contact time between ions in an electrolyte solution,from both spatial and temporal perspectives,new parameters were proposed to describe the interaction between ions:interaction distance(ID),and interaction time ratio(ITR).The best correlation between characteristic time ratio and solubility was found for a molar ratio of salt-to-water of 10:100 with a correlation coefficient of 0.96.For the same salt,a positive correlation was found between CTR and the molar ratio of salt and water.For type 1-1,type 2-1,type 1-2,and type 2-2 salts,the correlation coefficients between CTR and solubility were 0.93,0.96,0.92,and 0.98 for a salt-to-water molar ratio of 10:100,respectively.The solubility of multiple salts was predicted by simulations and compared with experimental values,yielding an average relative deviation of 12.4%.The new ion-interaction parameters offer significant advantages in describing strongly correlated and strongly hydrated electrolyte solutions.展开更多
To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitro...To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitrobenzene(DNB),solvents ligands dimethyl carbonate(DMC) and gamma-butyrolactone(GBL)with void were simulated,using molecular dynamics method and reactive force field.It is found that the CL-20 co-crystals with void defects will form hot spots when impacted,significantly affecting the decomposition of molecules around the void.The degree of molecular fragmentation is relatively low under the reflection velocity of 2 km/s,and the main reactions are the formation of dimer and the shedding of nitro groups.The existence of voids reduces the safety of CL-20 co-crystals,which induced the sensitivity of energetic co-crystals CL-20/TNT and CL-20/DNB to increase more significantly.Detonation has occurred under the reflection velocity of 4 km/s,energetic co-crystals are easier to polymerize than solvent co-crystals,and are not obviously affected by voids.The results show that the energy of the wave decreases after sweeping over the void,which reduces the chemical reaction frequency downstream of the void and affects the detonation performance,especially the solvent co-crystals.展开更多
The C–H bond activation in alkane dehydrogenation reactions is a key step in determining the reaction rate.To understand the impact of entropy,we performed ab initio static and molecular dynamics free energy simulati...The C–H bond activation in alkane dehydrogenation reactions is a key step in determining the reaction rate.To understand the impact of entropy,we performed ab initio static and molecular dynamics free energy simulations of ethane dehydrogenation over Co@BEA zeolite at different temperatures.AIMD simulations showed that a sharp decrease in free energy barrier as temperature increased.Our analysis of the temperature dependence of activation free energies uncovered an unusual entropic effect accompanying the reaction.The unique spatial structures around the Co active site at different temperatures influenced both the extent of charge transfer in the transition state and the arrangement of 3d orbital energy levels.We provided explanations consistent with the principles of thermodynamics and statistical physics.The insights gained at the atomic level have offered a fresh interpretation of the intricate long-range interplay between local chemical reactions and extensive chemical environments.展开更多
Microscopic pore structure in continental shale oil reservoirs is characterized by small pore throats and complex micro-structures.The adsorption behavior of hydrocarbons on the pore walls exhibits unique physical and...Microscopic pore structure in continental shale oil reservoirs is characterized by small pore throats and complex micro-structures.The adsorption behavior of hydrocarbons on the pore walls exhibits unique physical and chemical properties.Therefore,studying the adsorption morphology of hydrocarbon components in nanometer-sized pores and clarifying the exploitation limits of shale oil at the microscopic level are of great practical significance for the efficient development of continental shale oil.In this study,molecular dynamics simulations were employed to investigate the adsorption characteristics of various single-component shale oils in inorganic quartz fissures,and the influence of pore size and shale oil hydrocarbon composition on the adsorption properties in the pores was analyzed.The results show that different molecules have different adsorption capacities in shale oil pores,with lighter hydrocarbon components(C6H14)exhibiting stronger adsorption abilities.For the same adsorbed molecule,the adsorption amount linearly increases with the increase in pore diameter,but larger pores contribute more to shale oil adsorption.In shale pores,the thickness of the adsorption layer formed by shale oil molecules ranges from 0.4 to 0.5 nm,which is similar to the width of alkane molecules.Shale oil in the adsorbed state that is difficult to be exploited is mainly concentrated in the first adsorption layer.Among them,the volume fraction of adsorbed shale oil in 6 nm shale pores is 40.8%,while the volume fraction of shale oil that is difficult to be exploited is 16.2%.展开更多
Background:Based on network pharmacology and molecular docking,the present study investigated the mechanism of curcumin(CUR)in diabetic retinopathy treatment.Methods:Based on the DisGeNET,Swiss TargetPrediction,GeneCa...Background:Based on network pharmacology and molecular docking,the present study investigated the mechanism of curcumin(CUR)in diabetic retinopathy treatment.Methods:Based on the DisGeNET,Swiss TargetPrediction,GeneCards,Online Mendelian Inheritance in Man,Gene Expression Omnibus,and Comparative Toxicogenomics Database,the intersection core targets of CUR and diabetic retinopathy were identified.The intersection target was imported into the STRING database to obtain the protein-protein interaction map.According to the Database for Annotation,Visualization and Integrated Discovery database,the intersected targets were enriched in Gene Ontology(GO)and Kyoto Encyclopedia of Genes and Genomes pathways.Then Cytoscape 3.9.1 is used to make the drug-target-disease-pathway network.The mechanism of CUR and diabetic retinopathy was further verified by molecular docking and molecular dynamics simulation.Results:There were 203 intersecting targets of CUR and diabetic retinopathy identified.1320 GO entries were enriched for GO functions,which were primarily involved in the composition of cells such as identical protein binding,protein binding,enzyme binding,etc.It was found that 175 pathways were enriched using Kyoto Encyclopedia of Genes and Genomes pathway enrichment methods,which were mainly included in the lipid and atherosclerosis,AGE-RAGE signaling pathway in diabetic complications,pathways in cancer,etc.In the molecular docking analysis,CUR was found to have a good ability to bind to the core targets of albumin,IL-1B,and IL-6.The binding of albumin to CUR was further verified by molecular dynamics simulation.Conclusion:As a result of this study,CUR may exert a role in the treatment of diabetic retinopathy through multi-target and multi-pathway regulation,which indicates a possible direction of future research.展开更多
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.展开更多
This study investigates the mechanism by which baicalin inhibits cancer cell growth through estrogen receptor 1 (ESR1) using molecular dynamics simulations. The results show that baicalin primarily binds to the ligand...This study investigates the mechanism by which baicalin inhibits cancer cell growth through estrogen receptor 1 (ESR1) using molecular dynamics simulations. The results show that baicalin primarily binds to the ligand-binding domain (LBD) of ESR1, interacting through hydrogen bonds and hydrophobic interactions. After binding, the overall and local conformations of ESR1 change, affecting its interactions with other proteins and thus modulating the signaling pathways of cancer cells. Binding free energy analysis indicates that the binding of baicalin to ESR1 is spontaneous and relatively stable. Additionally, baicalin can inhibit the binding of ESR1 to estrogen, blocking the estrogen signaling pathway and thereby suppressing the growth and proliferation of cancer cells. This study provides theoretical and experimental foundations for the potential use of baicalin as an anticancer drug, offering new insights and methods for the development of novel anticancer drugs. However, the study has some limitations, such as limited simulation time and simplified systems. Future research can extend the simulation time and consider more physiological factors to more accurately simulate the interactions between baicalin and ESR1.展开更多
The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this wo...The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this work.The following properties were determined:density,self-diffusion coefficient,excess molar volume,and radial distribution function.The results show that with an increase in the mole fraction of[HPy][BF_(4)],the self-diffusion coefficient decreases.Additionally,the excess molar volume initially decreases,reaches a minimum,and then increases.The rules of radial distribution functions(RDFs)of characteristic atoms are different.With increasing the mole fraction of[HPy][BF_(4)],the first peak of the RDFs of HA1-F decreases,while that of CT6-CT6 rises at first and then decreases.This indicates that the solvent molecules affect the polar and non-polar regions of[HPy][BF_(4)]differently.展开更多
This paper develops a comprehensive computational modeling and simulation framework based on Complex Adaptive Systems(CAS)theory to unveil the underlying mechanisms of self-organization,nonlinear evolution,and emergen...This paper develops a comprehensive computational modeling and simulation framework based on Complex Adaptive Systems(CAS)theory to unveil the underlying mechanisms of self-organization,nonlinear evolution,and emergence in social systems.By integrating mathematical models,agent-based modeling,network dynamic analysis,and hybrid modeling approaches,the study applies CAS theory to case studies in economic markets,political decision-making,and social interactions.The experimental results demonstrate that local interactions among individual agents can give rise to complex global phenomena,such as market fluctuations,opinion polarization,and sudden outbreaks of social movements.This framework not only provides a more robust explanation for the nonlinear dynamics and abrupt transitions that traditional models often fail to capture,but also offers valuable decision-support tools for public policy formulation,social governance,and risk management.Emphasizing the importance of interdisciplinary approaches,this work outlines future research directions in high-performance computing,artificial intelligence,and real-time data integration to further advance the theoretical and practical applications of CAS in the social sciences.展开更多
The rapid solidification process of Mg7Zn3 alloy was simulated by the molecular dynamics method. The relationship between the local structure and the dynamics during the liquid-glass transition was deeply investigated...The rapid solidification process of Mg7Zn3 alloy was simulated by the molecular dynamics method. The relationship between the local structure and the dynamics during the liquid-glass transition was deeply investigated. It was found that the Mg-centered FK polyhedron and the Zn-centered icosahedron play a critical role in the formation of Mg7Zn3 metallic glass. The self-diffusion coefficients of Mg and Zn atoms deviate from the Arrhenius law near the melting temperature and then satisfy the power law. According to the time correlation functions of mean-square displacement, incoherent intermediate scattering function and non-Gaussian parameter, it was found that the β-relaxation in Mg7Zn3 supercooled liquid becomes more and more evident with decreasing temperature, and the α-relaxation time rapidly increases in the VFT law. Moreover, the smaller Zn atom has a faster relaxation behavior than the Mg atom. Some local atomic structures with short-range order have lower mobility, and they play a critical role in the appearance of cage effect in theβ-relaxation regime. The dynamics deviates from the Arrhenius law just at the temperature as the number of local atomic structures begins to rapidly increase. The dynamic glass transition temperature （Tc） is close to the glass transition point in structure （TgStr）.展开更多
The collecting power of tertiary amines（DRN,DEN and DPN） on kaolinite follows the order of DENDPNDRN.After reacting with DRN,DEN and DPN,the surface potentials of kaolinite increase remarkably,and the recruitments c...The collecting power of tertiary amines（DRN,DEN and DPN） on kaolinite follows the order of DENDPNDRN.After reacting with DRN,DEN and DPN,the surface potentials of kaolinite increase remarkably,and the recruitments caused by collectors also follow the order of DENDPNDRN.The results of dynamics simulation show that the geometries of substituent groups bonding to N are deflected and twisted,and some of bond angles are changed when tertiary amines cations adsorb on kaolinite（001） surface.Based on the results of dynamics simulations and quantum chemistry calculations,the electrostatic forces between three tertiary amines cations and 4×4×3（001） plane of kaolinite are 1.38×10？7 N（DRN12H＋）,1.44×10-6 N（DEN12H＋）,1.383×10-6 N（DPN12H＋）,respectively.展开更多
The distributions of local structural units of calcium silicate melts were quantified by means of classical molecular dynamics simulation and a newly constructed structural thermodynamic model. The distribution of fiv...The distributions of local structural units of calcium silicate melts were quantified by means of classical molecular dynamics simulation and a newly constructed structural thermodynamic model. The distribution of five kinds of Si-O tetrahedra Qi from these two methods was compared with each other and also with the experimental Raman spectra, an excellent agreement was achieved. These not only displayed the panorama distribution of microstructural units in the whole composition range, but also proved that the thermodynamic model is suitable for the utilization as the subsequent application model of spectral experiments for the thermodynamic calculation. Meanwhile, the five refined regions mastered by different disproportionating reactions were obtained. Finally, the distributions of two kinds of connections between Qi were obtained, denoted as Qi-Ca-Qj and Qi-[Ob]-Qj, from the thermodynamic model, and a theoretical verification was given that the dominant connections for any composition are equivalent connections.展开更多
基金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.
基金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.
文摘In this study,circular dichroism(CD)and molecular dynamics(MD)simulation were used to investigate the thermal unfolding pathway of staphylococcal enterotoxin B(SEB)at temperatures of 298–371 and 298–500 K,and the relationship between the experimental and simulation results were explored.Our computational findings on the secondary structure of SEB showed that at room temperature,the CD spectroscopic results were highly consistent with the MD results.Moreover,under heating conditions,the changing trends of helix,sheet and random coil obtained by CD spectral fitting were highly consistent with those obtained by MD.In order to gain a deeper understanding of the thermal stability mechanism of SEB,the MD trajectories were analyzed in terms of root mean square deviation(RMSD),secondary structure assignment(SSA),radius of gyration(R_(g)),free energy surfaces(FES),solvent-accessible surface area(SASA),hydrogen bonds and salt bridges.The results showed that at low heating temperature,domain Ⅰ without loops(omitting the mobile loop region)mainly relied on hydrophobic interaction to maintain its thermal stability,whereas the thermal stability of domain Ⅱ was mainly controlled by salt bridges and hydrogen bonds.Under high heating temperature conditions,the hydrophobic interactions in domain Ⅰ without loops were destroyed and the secondary structure was almost completely lost,while domain Ⅱ could still rely on salt bridges as molecular staples to barely maintain the stability of the secondary structure.These results help us to understand the thermodynamic and kinetic mechanisms that maintain the thermal stability of SEB at the molecular level,and provide a direction for establishing safer and more effective food sterilization processes.
基金supported by the National Natural Science Foundation of China(Grant Nos.52035009 and 51761135106)the State Key Laboratory of Precision Measuring Technology and Instruments(Pilt1705)+1 种基金the Henan Key Laboratory of Intelligent Manufacturing Equipment Integration for Superhard Materials(JDKJ2022-01)the“111”project by the State Administration of Foreign Experts Affairs and the Ministry of Education of China(Grant No.B07014).
文摘Nitrogen vacancy(NV)color centers in diamond have useful applications in quantum sensing andfluorescent marking.They can be gen-erated experimentally by ion implantation,femtosecond lasers,and chemical vapor deposition.However,there is a lack of studies of the yield of NV color centers at the atomic scale.In the molecular dynamics simulations described in this paper,NV color centers are pre-pared by ion implantation in diamond with pre-doped nitrogen and subsequent annealing.The differences between the yields of NV color centers produced by implantation of carbon(C)and nitrogen(N)ions,respectively,are investigated.It is found that C-ion implantation gives a greater yield of NV color centers and superior location accuracy.The effects of different pre-doping concentrations(400–1500 ppm)and implantation energies(1.0–3.0 keV)on the NV color center yield are analyzed,and it is shown that a pre-doping concentra-tion of 1000 ppm with 2 keV C-ion implantation can produce a 13%yield of NV color centers after 1600 K annealing for 7.4 ns.Finally,a brief comparison of the NV color center identification methods is presented,and it is found that the error rate of an analysis utiliz-ing the identify diamond structure coordination analysis method is reduced by about 7%compared with conventional identification+methods.
基金supported by National Key Research and Development Program of China(No.2021YFB3400800)National Natural Science Foundation of China(Grant No.52271136,51901177)Natural Science Foundation of Shaanxi Province(No.2021JC-06,2019TD-020).
文摘Grain boundaries(GBs)play a crucial role on the structural stability and mechanical properties of Cu and its alloys.In this work,molecular dynamics(MD)simulations are employed to study the effects of Fe solutes on the formation energy,excess volume,dislocations and melting behaviors of GBs in CuFe alloys.It is illustrated that Fe solute affects the structural stability of Cu GBs substantially,the formation energy of GBs is reduced,but the thickness and melting point of GBs are increased,that is,the structural stability of Cu GBs is significantly improved owing to the Fe solutes.A strong scaling law exists between the formation energy,excess volume,thickness and melting point of GBs.Therefore,Fe solid solute plays an important role in the characteristics of GBs in bi-crystal Cu.
基金Project supported by the National Natural Science Foundation of China(Grant No.52072132).
文摘Recently,there has been a growing prevalence in the utilization of graphdiyne(GDY)in the field of biomedicine,attributed to its distinctive physical structure and chemical properties.Additionally,its biocompatibility has garnered increasing attention.However,there is a lack of research on the biological effects and physical mechanisms of GDYprotein interactions at the molecular scale.In this study,the villin headpiece subdomain(HP35)served as a representative protein model.Molecular dynamics simulations were employed to investigate the interaction process between the HP35 protein and GDY,as well as the structural evolution of the protein.The data presented in our study demonstrate that GDY can rapidly adsorb HP35 protein and induce denaturation to one of the a-helix structures of HP35 protein.This implies a potential cytotoxicity concern of GDY for biological systems.Compared to graphene,GDY induced less disruption to HP35 protein.This can be attributed to the presence of natural triangular vacancies in GDY,which prevents p–p stacking action and the limited interaction of GDY with HP35 protein is not conducive to the expansion of protein structures.These findings unveil the biological effects of GDY at the molecular level and provide valuable insights for the application of GDY in biomedicine.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.51775077 and 51909023)。
文摘Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear-thinning fluid in a microchannel.We validated the feasibility of our simulation method by evaluating the mean square displacement and Reynolds number of the solution layers.The results show that the change rule of the fluid system's velocity profile and interaction energy can reflect the shear-thinning characteristics of the fluids.The velocity profile resembles a top-hat shape,intensifying as the fluid's power law index decreases.The interaction energy between the wall and the fluid decreases gradually with increasing velocity,and a high concentration of non-Newtonian fluid reaches a plateau sooner.Moreover,the velocity profile of the fluid is related to the molecule number density distribution and their values are inversely proportional.By analyzing the radial distribution function,we found that the hydrogen bonds between solute and water molecules weaken with the increase in velocity.This observation offers an explanation for the shear-thinning phenomenon of the non-Newtonian flow from a micro perspective.
基金supported by the National Natural Science Foundation of China(No.12102256).
文摘Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectile impacts is thoroughly investigated using full-atomic molecular dynamics(MD)simulations.The study results show that the porous structure and density are key factors determining the mechanical response of GA under impact loading.Specifically,the impact-induced penetration of the projectile leads to the collapse of the pore structure,causing stretching and subsequent rupture of covalent bonds in graphene sheets.Moreover,the effects of temperature on the mechanical performance of GA have been proven to be minimal,thereby highlighting the mechanical stability of GA over a wide range of temperatures.Finally,the energy absorption density(EAD)and energy absorption efficiency(EAE)metrics are adopted to assess the energy absorption capacity of GA during projectile penetration.The research findings of this work demonstrate the significant potential of GA for energy absorption applications.
基金supported by the National Natural Science Foundation of China(No.21776264).
文摘Studying the relationship between ionic interactions and salt solubility in seawater has implications for seawater desalination and mineral extraction.In this paper,a new method of expressing ion-to-ion interaction is proposed by using molecular dynamics simulation,and the relationship between ion-to-ion interaction and salt solubility in a simulated seawater water-salt system is investigated.By analyzing the variation of distance and contact time between ions in an electrolyte solution,from both spatial and temporal perspectives,new parameters were proposed to describe the interaction between ions:interaction distance(ID),and interaction time ratio(ITR).The best correlation between characteristic time ratio and solubility was found for a molar ratio of salt-to-water of 10:100 with a correlation coefficient of 0.96.For the same salt,a positive correlation was found between CTR and the molar ratio of salt and water.For type 1-1,type 2-1,type 1-2,and type 2-2 salts,the correlation coefficients between CTR and solubility were 0.93,0.96,0.92,and 0.98 for a salt-to-water molar ratio of 10:100,respectively.The solubility of multiple salts was predicted by simulations and compared with experimental values,yielding an average relative deviation of 12.4%.The new ion-interaction parameters offer significant advantages in describing strongly correlated and strongly hydrated electrolyte solutions.
基金supported by the National Natural Science Foundation of China (22275018)the Project of State Key Laboratory of Explosion Science and Technology (Beijing Institute of Technology)(Grant No.QNKT20-04)。
文摘To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitrobenzene(DNB),solvents ligands dimethyl carbonate(DMC) and gamma-butyrolactone(GBL)with void were simulated,using molecular dynamics method and reactive force field.It is found that the CL-20 co-crystals with void defects will form hot spots when impacted,significantly affecting the decomposition of molecules around the void.The degree of molecular fragmentation is relatively low under the reflection velocity of 2 km/s,and the main reactions are the formation of dimer and the shedding of nitro groups.The existence of voids reduces the safety of CL-20 co-crystals,which induced the sensitivity of energetic co-crystals CL-20/TNT and CL-20/DNB to increase more significantly.Detonation has occurred under the reflection velocity of 4 km/s,energetic co-crystals are easier to polymerize than solvent co-crystals,and are not obviously affected by voids.The results show that the energy of the wave decreases after sweeping over the void,which reduces the chemical reaction frequency downstream of the void and affects the detonation performance,especially the solvent co-crystals.
文摘The C–H bond activation in alkane dehydrogenation reactions is a key step in determining the reaction rate.To understand the impact of entropy,we performed ab initio static and molecular dynamics free energy simulations of ethane dehydrogenation over Co@BEA zeolite at different temperatures.AIMD simulations showed that a sharp decrease in free energy barrier as temperature increased.Our analysis of the temperature dependence of activation free energies uncovered an unusual entropic effect accompanying the reaction.The unique spatial structures around the Co active site at different temperatures influenced both the extent of charge transfer in the transition state and the arrangement of 3d orbital energy levels.We provided explanations consistent with the principles of thermodynamics and statistical physics.The insights gained at the atomic level have offered a fresh interpretation of the intricate long-range interplay between local chemical reactions and extensive chemical environments.
文摘Microscopic pore structure in continental shale oil reservoirs is characterized by small pore throats and complex micro-structures.The adsorption behavior of hydrocarbons on the pore walls exhibits unique physical and chemical properties.Therefore,studying the adsorption morphology of hydrocarbon components in nanometer-sized pores and clarifying the exploitation limits of shale oil at the microscopic level are of great practical significance for the efficient development of continental shale oil.In this study,molecular dynamics simulations were employed to investigate the adsorption characteristics of various single-component shale oils in inorganic quartz fissures,and the influence of pore size and shale oil hydrocarbon composition on the adsorption properties in the pores was analyzed.The results show that different molecules have different adsorption capacities in shale oil pores,with lighter hydrocarbon components(C6H14)exhibiting stronger adsorption abilities.For the same adsorbed molecule,the adsorption amount linearly increases with the increase in pore diameter,but larger pores contribute more to shale oil adsorption.In shale pores,the thickness of the adsorption layer formed by shale oil molecules ranges from 0.4 to 0.5 nm,which is similar to the width of alkane molecules.Shale oil in the adsorbed state that is difficult to be exploited is mainly concentrated in the first adsorption layer.Among them,the volume fraction of adsorbed shale oil in 6 nm shale pores is 40.8%,while the volume fraction of shale oil that is difficult to be exploited is 16.2%.
基金supported by the Hubei Province Research Innovation Team Project(T2021022)Scientific Research Projects of Hubei Health Commission(WJ2023M119).
文摘Background:Based on network pharmacology and molecular docking,the present study investigated the mechanism of curcumin(CUR)in diabetic retinopathy treatment.Methods:Based on the DisGeNET,Swiss TargetPrediction,GeneCards,Online Mendelian Inheritance in Man,Gene Expression Omnibus,and Comparative Toxicogenomics Database,the intersection core targets of CUR and diabetic retinopathy were identified.The intersection target was imported into the STRING database to obtain the protein-protein interaction map.According to the Database for Annotation,Visualization and Integrated Discovery database,the intersected targets were enriched in Gene Ontology(GO)and Kyoto Encyclopedia of Genes and Genomes pathways.Then Cytoscape 3.9.1 is used to make the drug-target-disease-pathway network.The mechanism of CUR and diabetic retinopathy was further verified by molecular docking and molecular dynamics simulation.Results:There were 203 intersecting targets of CUR and diabetic retinopathy identified.1320 GO entries were enriched for GO functions,which were primarily involved in the composition of cells such as identical protein binding,protein binding,enzyme binding,etc.It was found that 175 pathways were enriched using Kyoto Encyclopedia of Genes and Genomes pathway enrichment methods,which were mainly included in the lipid and atherosclerosis,AGE-RAGE signaling pathway in diabetic complications,pathways in cancer,etc.In the molecular docking analysis,CUR was found to have a good ability to bind to the core targets of albumin,IL-1B,and IL-6.The binding of albumin to CUR was further verified by molecular dynamics simulation.Conclusion:As a result of this study,CUR may exert a role in the treatment of diabetic retinopathy through multi-target and multi-pathway regulation,which indicates a possible direction of future research.
基金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.
文摘This study investigates the mechanism by which baicalin inhibits cancer cell growth through estrogen receptor 1 (ESR1) using molecular dynamics simulations. The results show that baicalin primarily binds to the ligand-binding domain (LBD) of ESR1, interacting through hydrogen bonds and hydrophobic interactions. After binding, the overall and local conformations of ESR1 change, affecting its interactions with other proteins and thus modulating the signaling pathways of cancer cells. Binding free energy analysis indicates that the binding of baicalin to ESR1 is spontaneous and relatively stable. Additionally, baicalin can inhibit the binding of ESR1 to estrogen, blocking the estrogen signaling pathway and thereby suppressing the growth and proliferation of cancer cells. This study provides theoretical and experimental foundations for the potential use of baicalin as an anticancer drug, offering new insights and methods for the development of novel anticancer drugs. However, the study has some limitations, such as limited simulation time and simplified systems. Future research can extend the simulation time and consider more physiological factors to more accurately simulate the interactions between baicalin and ESR1.
文摘The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this work.The following properties were determined:density,self-diffusion coefficient,excess molar volume,and radial distribution function.The results show that with an increase in the mole fraction of[HPy][BF_(4)],the self-diffusion coefficient decreases.Additionally,the excess molar volume initially decreases,reaches a minimum,and then increases.The rules of radial distribution functions(RDFs)of characteristic atoms are different.With increasing the mole fraction of[HPy][BF_(4)],the first peak of the RDFs of HA1-F decreases,while that of CT6-CT6 rises at first and then decreases.This indicates that the solvent molecules affect the polar and non-polar regions of[HPy][BF_(4)]differently.
文摘This paper develops a comprehensive computational modeling and simulation framework based on Complex Adaptive Systems(CAS)theory to unveil the underlying mechanisms of self-organization,nonlinear evolution,and emergence in social systems.By integrating mathematical models,agent-based modeling,network dynamic analysis,and hybrid modeling approaches,the study applies CAS theory to case studies in economic markets,political decision-making,and social interactions.The experimental results demonstrate that local interactions among individual agents can give rise to complex global phenomena,such as market fluctuations,opinion polarization,and sudden outbreaks of social movements.This framework not only provides a more robust explanation for the nonlinear dynamics and abrupt transitions that traditional models often fail to capture,but also offers valuable decision-support tools for public policy formulation,social governance,and risk management.Emphasizing the importance of interdisciplinary approaches,this work outlines future research directions in high-performance computing,artificial intelligence,and real-time data integration to further advance the theoretical and practical applications of CAS in the social sciences.
基金Project (51101022) supported by the National Natural Science Foundation of ChinaProject (CHD2012JC096) supported by the Fundamental Research Funds for the Central Universities,China
文摘The rapid solidification process of Mg7Zn3 alloy was simulated by the molecular dynamics method. The relationship between the local structure and the dynamics during the liquid-glass transition was deeply investigated. It was found that the Mg-centered FK polyhedron and the Zn-centered icosahedron play a critical role in the formation of Mg7Zn3 metallic glass. The self-diffusion coefficients of Mg and Zn atoms deviate from the Arrhenius law near the melting temperature and then satisfy the power law. According to the time correlation functions of mean-square displacement, incoherent intermediate scattering function and non-Gaussian parameter, it was found that the β-relaxation in Mg7Zn3 supercooled liquid becomes more and more evident with decreasing temperature, and the α-relaxation time rapidly increases in the VFT law. Moreover, the smaller Zn atom has a faster relaxation behavior than the Mg atom. Some local atomic structures with short-range order have lower mobility, and they play a critical role in the appearance of cage effect in theβ-relaxation regime. The dynamics deviates from the Arrhenius law just at the temperature as the number of local atomic structures begins to rapidly increase. The dynamic glass transition temperature （Tc） is close to the glass transition point in structure （TgStr）.
基金Project (2005CB623701) supported by the National Basic Research Program of China Project (201011031) supported by National Department Public Benefit Research Foundation from Ministry of Land and Resources of China+1 种基金 Project (2935) supported by the Foundation for the Author of Zhengzhou Institute of Multipurpose Utilization of Mineral Resources CAGS, China Project (1212011120304) supported by the Geological Surrey Program
文摘The collecting power of tertiary amines（DRN,DEN and DPN） on kaolinite follows the order of DENDPNDRN.After reacting with DRN,DEN and DPN,the surface potentials of kaolinite increase remarkably,and the recruitments caused by collectors also follow the order of DENDPNDRN.The results of dynamics simulation show that the geometries of substituent groups bonding to N are deflected and twisted,and some of bond angles are changed when tertiary amines cations adsorb on kaolinite（001） surface.Based on the results of dynamics simulations and quantum chemistry calculations,the electrostatic forces between three tertiary amines cations and 4×4×3（001） plane of kaolinite are 1.38×10？7 N（DRN12H＋）,1.44×10-6 N（DEN12H＋）,1.383×10-6 N（DPN12H＋）,respectively.
基金Project(2012CB722805)supported by the National Basic Research Program of ChinaProjects(50504010,50974083,51174131,51374141)supported by the National Natural Science Foundation of China+1 种基金Project(50774112)supported by the Joint Fund of NSFC and Baosteel,ChinaProject(07QA4021)supported by the Shanghai"Phosphor"Science Foundation,China
文摘The distributions of local structural units of calcium silicate melts were quantified by means of classical molecular dynamics simulation and a newly constructed structural thermodynamic model. The distribution of five kinds of Si-O tetrahedra Qi from these two methods was compared with each other and also with the experimental Raman spectra, an excellent agreement was achieved. These not only displayed the panorama distribution of microstructural units in the whole composition range, but also proved that the thermodynamic model is suitable for the utilization as the subsequent application model of spectral experiments for the thermodynamic calculation. Meanwhile, the five refined regions mastered by different disproportionating reactions were obtained. Finally, the distributions of two kinds of connections between Qi were obtained, denoted as Qi-Ca-Qj and Qi-[Ob]-Qj, from the thermodynamic model, and a theoretical verification was given that the dominant connections for any composition are equivalent connections.
