The micro-capsules used for drug delivery are fabricated using polylactic acid（PLA）,which is a biomedical material approved by the FDA.A coarse-grained model of long-chain PLA was built,and molecular dynamics（MD）s...The micro-capsules used for drug delivery are fabricated using polylactic acid（PLA）,which is a biomedical material approved by the FDA.A coarse-grained model of long-chain PLA was built,and molecular dynamics（MD）simulations of the model were performed using a MARTINI force field.Based on the nonlocal theory,the formula for the initial elastic modulus of polymers considering the nonlocal effect was derived,and the scaling law of internal characteristic length of polymers was proposed,which was used to adjust the cut-off radius in the MD simulations of PLA.The results show that the elastic modulus should be computed using nonlinear regression.The nonlocal effect has a certain influence on the simulation results of PLA.According to the scaling law,the cut-off radius was determined and applied to the MD simulations,the results of which reflect the influence of the molecular weight change on the elastic moduli of PLA,and are in agreement with the experimental outcome.展开更多
Compared to the well-studied two-dimensional(2D)ferroelectricity,the appearance of 2D antiferroelectricity is much rarer,where local dipoles from the nonequivalent sublattices within 2D monolayers are oppositely orien...Compared to the well-studied two-dimensional(2D)ferroelectricity,the appearance of 2D antiferroelectricity is much rarer,where local dipoles from the nonequivalent sublattices within 2D monolayers are oppositely oriented.Using NbOCl_(2) monolayer with competing ferroelectric(FE)and antiferroelectric(AFE)phases as a 2D material platform,we demonstrate the emergence of intrinsic antiferroelectricity in NbOCl_(2) monolayer under experimentally accessible shear strain,along with new functionality associated with electric field-induced AFE-to-FE phase transition.Specifically,the complex configuration space accommodating FE and AFE phases,polarization switching kinetics,and finite temperature thermodynamic properties of 2D NbOCl_(2) are all accurately predicted by large-scale molecular dynamics simulations based on deep learning interatomic potential model.Moreover,room temperature stable antiferroelectricity with low polarization switching barrier and one-dimensional collinear polarization arrangement is predicted in shear-deformed NbOCl_(2) monolayer.The transition from AFE to FE phase in 2D NbOCl_(2) can be triggered by a low critical electric field,leading to a double polarization–electric(P–E)loop with small hysteresis.A new type of optoelectronic device composed of AFE-NbOCl_(2) is proposed,enabling electric“writing”and nonlinear optical“reading”logical operation with fast operation speed and low power consumption.展开更多
The secondary structure of different Iβcellulose was analyzed by a molecular dynamics sim-ulation with MARTINI coarse-grained force field,where each chain of the cellulose includes 40 D-glucoses units.Calculation giv...The secondary structure of different Iβcellulose was analyzed by a molecular dynamics sim-ulation with MARTINI coarse-grained force field,where each chain of the cellulose includes 40 D-glucoses units.Calculation gives a satisfied description about the secondary structure of the cellulose.As the chain number increasing,the cellulose becomes the form of a helix,with the diameter of screw growing and spiral rising.Interestingly,the celluloses with chain number N of 4,6,24 and 36 do show right-hand twisting.On the contrast,the celluloses with N of 8,12,16 chains are left-hand twisting.These simulations indicate that the cellulose with chain number larger than 36 will break down to two parts.Besides,the result indicates that 36-chains cellulose model is the most stable among all models.Furthermore,the Lennard-Jones potential determines the secondary structure.In addition,an equation was set up to analyze the twisting structure.展开更多
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.展开更多
To overcome the limitations of microscale experimental techniques and molecular dynamics(MD)simulations,a coarse-grained molecular dynamics(CGMD)method was used to simulate the wetting processes of clay aggregates.Bas...To overcome the limitations of microscale experimental techniques and molecular dynamics(MD)simulations,a coarse-grained molecular dynamics(CGMD)method was used to simulate the wetting processes of clay aggregates.Based on the evolution of swelling stress,final dry density,water distribution,and clay arrangements under different target water contents and dry densities,a relationship between the swelling behaviors and microstructures was established.The simulated results showed that when the clay-water well depth was 300 kcal/mol,the basal spacing from CGMD was consistent with the X-ray diffraction(XRD)data.The effect of initial dry density on swelling stress was more pronounced than that of water content.The anisotropic swelling characteristics of the aggregates are related to the proportion of horizontally oriented clay mineral layers.The swelling stress was found to depend on the distribution of tactoids at the microscopic level.At lower initial dry density,the distribution of tactoids was mainly controlled by water distribution.With increase in the bound water content,the basal spacing expanded,and the swelling stresses increased.Free water dominated at higher water contents,and the particles were easily rotated,leading to a decrease in the number of large tactoids.At higher dry densities,the distances between the clay mineral layers decreased,and the movement was limited.When bound water enters the interlayers,there is a significant increase in interparticle repulsive forces,resulting in a greater number of small-sized tactoids.Eventually,a well-defined logarithmic relationship was observed between the swelling stress and the total number of tactoids.These findings contribute to a better understanding of coupled macro-micro swelling behaviors of montmorillonite-based materials,filling a study gap in clay-water interactions on a micro scale.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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 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 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.展开更多
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 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.展开更多
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 global rapid transition towards sustainable energy systems has heightened the demand for highperformance lithium metal batteries(LMBs),where understanding interfacial phenomena is paramount.In this contribution,we...The global rapid transition towards sustainable energy systems has heightened the demand for highperformance lithium metal batteries(LMBs),where understanding interfacial phenomena is paramount.In this contribution,we present an on-the-fly machine learning molecular dynamics(OTF-MLMD)approach to probe the complex side reactions at lithium metal anode–electrolyte interfaces with exceptional accuracy and computational efficiency.The machine learning force field(MLFF)was firstly validated in a bulk-phase system comprising twenty 1,2-dimethoxyethane(DME)molecules,demonstrating energy fluctuations and structural parameters in close agreement with ab initio molecular dynamics(AIMD)benchmarks.Subsequent simulations of lithium–DME and lithium–electrolyte interfaces revealed minimal discrepancies in energy,bond lengths,and net charge variations(notably in FSI-species),underscoring the method's DFT-level precision of the approach.A further small-scale interfacial model enabled on-the-fly training over a mere of 340 fs,which was then successfully transferred to a large-scale simulation encompassing nearly 300,000 atoms,representing the largest interfacial model in LMB research up to date.The hierarchical validation strategy not only establishes the robustness of the MLFF in capturing both interfacial and bulk-phase chemistry but also paves the way for statistically meaningful simulations of battery interfaces.The fruitful findings highlight the transformative potential of OTF-MLMD in bridging the gap between atomistic accuracy and macroscopic modeling,affording a universal approach to understand interfacial reactions in LMBs.展开更多
基金Project supported by the National Natural Science Foundation of China(no.11272360)the Natural Science Foundation of Guangdong Province(no.2014A030313793)+1 种基金the Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund(the second phase)National Supercomputer Center in Guangzhou
文摘The micro-capsules used for drug delivery are fabricated using polylactic acid（PLA）,which is a biomedical material approved by the FDA.A coarse-grained model of long-chain PLA was built,and molecular dynamics（MD）simulations of the model were performed using a MARTINI force field.Based on the nonlocal theory,the formula for the initial elastic modulus of polymers considering the nonlocal effect was derived,and the scaling law of internal characteristic length of polymers was proposed,which was used to adjust the cut-off radius in the MD simulations of PLA.The results show that the elastic modulus should be computed using nonlinear regression.The nonlocal effect has a certain influence on the simulation results of PLA.According to the scaling law,the cut-off radius was determined and applied to the MD simulations,the results of which reflect the influence of the molecular weight change on the elastic moduli of PLA,and are in agreement with the experimental outcome.
基金supported by the National Natural Science Foundation of China (Grant No.11574244 for G.Y.G.)the XJTU Research Fund for AI Science (Grant No.2025YXYC011 for G.Y.G.)the Hong Kong Global STEM Professorship Scheme (for X.C.Z.)。
文摘Compared to the well-studied two-dimensional(2D)ferroelectricity,the appearance of 2D antiferroelectricity is much rarer,where local dipoles from the nonequivalent sublattices within 2D monolayers are oppositely oriented.Using NbOCl_(2) monolayer with competing ferroelectric(FE)and antiferroelectric(AFE)phases as a 2D material platform,we demonstrate the emergence of intrinsic antiferroelectricity in NbOCl_(2) monolayer under experimentally accessible shear strain,along with new functionality associated with electric field-induced AFE-to-FE phase transition.Specifically,the complex configuration space accommodating FE and AFE phases,polarization switching kinetics,and finite temperature thermodynamic properties of 2D NbOCl_(2) are all accurately predicted by large-scale molecular dynamics simulations based on deep learning interatomic potential model.Moreover,room temperature stable antiferroelectricity with low polarization switching barrier and one-dimensional collinear polarization arrangement is predicted in shear-deformed NbOCl_(2) monolayer.The transition from AFE to FE phase in 2D NbOCl_(2) can be triggered by a low critical electric field,leading to a double polarization–electric(P–E)loop with small hysteresis.A new type of optoelectronic device composed of AFE-NbOCl_(2) is proposed,enabling electric“writing”and nonlinear optical“reading”logical operation with fast operation speed and low power consumption.
文摘The secondary structure of different Iβcellulose was analyzed by a molecular dynamics sim-ulation with MARTINI coarse-grained force field,where each chain of the cellulose includes 40 D-glucoses units.Calculation gives a satisfied description about the secondary structure of the cellulose.As the chain number increasing,the cellulose becomes the form of a helix,with the diameter of screw growing and spiral rising.Interestingly,the celluloses with chain number N of 4,6,24 and 36 do show right-hand twisting.On the contrast,the celluloses with N of 8,12,16 chains are left-hand twisting.These simulations indicate that the cellulose with chain number larger than 36 will break down to two parts.Besides,the result indicates that 36-chains cellulose model is the most stable among all models.Furthermore,the Lennard-Jones potential determines the secondary structure.In addition,an equation was set up to analyze the twisting structure.
基金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 National Natural Science Foundation of China(Grant No.42172308)the Youth Innovation Promotion Association CAS(Grant No.2022331)the Key Research and Development Program of Hubei Province(Grant No.2022BAA036).
文摘To overcome the limitations of microscale experimental techniques and molecular dynamics(MD)simulations,a coarse-grained molecular dynamics(CGMD)method was used to simulate the wetting processes of clay aggregates.Based on the evolution of swelling stress,final dry density,water distribution,and clay arrangements under different target water contents and dry densities,a relationship between the swelling behaviors and microstructures was established.The simulated results showed that when the clay-water well depth was 300 kcal/mol,the basal spacing from CGMD was consistent with the X-ray diffraction(XRD)data.The effect of initial dry density on swelling stress was more pronounced than that of water content.The anisotropic swelling characteristics of the aggregates are related to the proportion of horizontally oriented clay mineral layers.The swelling stress was found to depend on the distribution of tactoids at the microscopic level.At lower initial dry density,the distribution of tactoids was mainly controlled by water distribution.With increase in the bound water content,the basal spacing expanded,and the swelling stresses increased.Free water dominated at higher water contents,and the particles were easily rotated,leading to a decrease in the number of large tactoids.At higher dry densities,the distances between the clay mineral layers decreased,and the movement was limited.When bound water enters the interlayers,there is a significant increase in interparticle repulsive forces,resulting in a greater number of small-sized tactoids.Eventually,a well-defined logarithmic relationship was observed between the swelling stress and the total number of tactoids.These findings contribute to a better understanding of coupled macro-micro swelling behaviors of montmorillonite-based materials,filling a study gap in clay-water interactions on a micro scale.
基金Research Grants of Senior Research Fellowship in favor of first author(Gloak Majumdar)from Council of Scientific and Industrial Research(CSIR,New Delhi,Government of India)(CSIR-SRF)with Award No.09/1151/(0008)2020-EMR-I.
文摘Objective To evaluate the antibacterial potential of bioactive compounds from Persicaria hydropiper(L.)(P.hydropiper)against bacterial virulence proteins through molecular docking(MD)and experimental validation.Methods Six bioactive compounds from P.hydropiper were investigated:catechin(CAT1),hyperin(HYP1),ombuin(OMB1),pinosylvin(PSV1),quercetin 3-sulfate(QSF1),and scutellarein(SCR1).Their binding affinities and potential binding pockets were assessed through MD against four bacterial target proteins with Protein Data Bank identifiers(PDB IDs):topoisomerase IV from Escherichia coli(E.coli)(PDB ID:3FV5),Staphylococcus aureus(S.aureus)gyrase ATPase binding domain(PDB ID:3U2K),CviR from Chromobacterium violaceum(C.violaceum)(PDB ID:3QP1),and glycosyl hydrolase from Pseudomonas aeruginosa(P.aeruginosa)(PDB ID:5BX9).Molecular dynamics simulations(MDS)were performed on the most promising compound-protein complexes for 50 nanoseconds(ns).Drug-likeness was evaluated using Lipinski's Rule of Five(RO5),followed by absorption,distribution,metabolism,excretion,and toxicity(ADMET)analysis using SwissADME and pkCSM web servers.Antibacterial activity was evaluated through disc diffusion assays,testing both individual compounds and combinations with conventional antibiotics[cefotaxime(CTX1,30μg/disc),ceftazidime(CAZ1,30μg/disc),and piperacillin(PIP1,100μg/disc)].Results MD revealed strong binding affinity(ranging from-9.3 to-5.9 kcal/mol)for all compounds,with CAT1 showing exceptional binding to 3QP1(-9.3 kcal/mol)and 5BX9(-8.4 kcal/mol).MDS confirmed the stability of CAT1-protein complexes with binding free energies of-84.71 kJ/mol(5BX9-CAT1)and-95.59 kJ/mol(3QP1-CAT1).Five compounds(CAT1,SCR1,PSV1,OMB1,and QSF1)complied with Lipinski's RO5 and showed favorable ADMET profiles.All compounds were non-carcinogenic,with CAT1 classified in the lowest toxicity class(VI).In antibacterial assays,CAT1 demonstrated significant activity against both gram-positive bacteria[Streptococcus pneumoniae(S.pneumoniae),S.aureus,and Bacillus cereus(B.cereus)][zone diameter of inhibition(ZDI):10-22 mm]and gram-negative bacteria[Acinetobacter baumannii(A.baumannii),E.coli,and P.aeruginosa](ZDI:14-27 mm).Synergistic effects were observed when CAT1 was combined with antibiotics and the growth inhibitory indices(GII)was 0.69-1.00.Conclusion P.hydropiper bioactive compounds,particularly CAT1,show promising antibacterial potential through multiple mechanisms,including direct inhibition of bacterial virulence proteins and synergistic activity with conventional antibiotics.The favorable pharmacological properties and low toxicity profiles support their potential development as therapeutic agents against bacterial infections.
基金Project supported by the National 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.
基金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 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.
基金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.
基金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.
基金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.
基金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 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.
基金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.
基金supported by the National Key Research and Development Program(2021YFB2500300)the National Natural Science Foundation of China(T2322015,92472101,22393903,22393900,and 52394170)+1 种基金the Beijing Municipal Natural Science Foundation(L247015 and L233004)Tsinghua University Initiative Scientific Research Program。
文摘The global rapid transition towards sustainable energy systems has heightened the demand for highperformance lithium metal batteries(LMBs),where understanding interfacial phenomena is paramount.In this contribution,we present an on-the-fly machine learning molecular dynamics(OTF-MLMD)approach to probe the complex side reactions at lithium metal anode–electrolyte interfaces with exceptional accuracy and computational efficiency.The machine learning force field(MLFF)was firstly validated in a bulk-phase system comprising twenty 1,2-dimethoxyethane(DME)molecules,demonstrating energy fluctuations and structural parameters in close agreement with ab initio molecular dynamics(AIMD)benchmarks.Subsequent simulations of lithium–DME and lithium–electrolyte interfaces revealed minimal discrepancies in energy,bond lengths,and net charge variations(notably in FSI-species),underscoring the method's DFT-level precision of the approach.A further small-scale interfacial model enabled on-the-fly training over a mere of 340 fs,which was then successfully transferred to a large-scale simulation encompassing nearly 300,000 atoms,representing the largest interfacial model in LMB research up to date.The hierarchical validation strategy not only establishes the robustness of the MLFF in capturing both interfacial and bulk-phase chemistry but also paves the way for statistically meaningful simulations of battery interfaces.The fruitful findings highlight the transformative potential of OTF-MLMD in bridging the gap between atomistic accuracy and macroscopic modeling,affording a universal approach to understand interfacial reactions in LMBs.
