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.展开更多
Using molecular dynamics methods,simulations of collision cascades in polycrystalline tungsten(W)have been conducted in this study,including different primary-knock-on atom(PKA)directions,grain sizes,and PKA energies ...Using molecular dynamics methods,simulations of collision cascades in polycrystalline tungsten(W)have been conducted in this study,including different primary-knock-on atom(PKA)directions,grain sizes,and PKA energies between 1 keV and 150 keV.The results indicate that a smaller grain size leads to more defects forming in grain boundary regions during cascade processes.The impact of high-energy PKA may cause a certain degree of distortion of the grain boundaries,which has a higher probability in systems with smaller grain sizes and becomes more pronounced as the PKA energy increases.The direction of PKA can affect the formation and diffusion pathways of defects.When the PKA direction is perpendicular to the grain boundary,defects preferentially form near the grain boundary regions;by contrast,defects are more inclined to form in the interior of the grains.These results are of great significance for comprehending the changes in the performance of polycrystalline W under the high-energy fusion environments and can provide theoretical guidance for further optimization and application of W-based plasma materials.展开更多
The 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.展开更多
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.展开更多
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.展开更多
According to the complex differential accumulation history of deep marine oil and gas in superimposed basins,the Lower Paleozoic petroleum system in Tahe Oilfield of Tarim Basin is selected as a typical case,and the p...According to the complex differential accumulation history of deep marine oil and gas in superimposed basins,the Lower Paleozoic petroleum system in Tahe Oilfield of Tarim Basin is selected as a typical case,and the process of hydrocarbon generation and expulsion,migration and accumulation,adjustment and transformation of deep oil and gas is restored by means of reservoine-forming dynamics simulation.The thermal evolution history of the Lower Cambrian source rocks in Tahe Oilfield reflects the obvious differences in hydrocarbon generation and expulsion process and intensity in different tectonic zones,which is the main reason controlling the differences in deep oil and gas phases.The complex transport system composed of strike-slip fault and unconformity,etc.controlled early migration and accumulation and late adjustment of deep oil and gas,while the Middle Cambrian gypsum-salt rock in inner carbonate platform prevented vertical migration and accumulation of deep oil and gas,resulting in an obvious"fault-controlled"feature of deep oil and gas,in which the low potential area superimposed by the NE-strike-slip fault zone and deep oil and gas migration was conducive to accumulation,and it is mainly beaded along the strike-slip fault zone in the northeast direction.The dynamic simulation of reservoir formation reveals that the spatio-temporal configuration of"source-fault-fracture-gypsum-preservation"controls the differential accumulation of deep oil and gas in Tahe Oilfield.The Ordovician has experienced the accumulation history of multiple periods of charging,vertical migration and accumulation,and lateral adjustment and transformation,and deep oil and gas have always been in the dynamic equilibrium of migration,accumulation and escape.The statistics of residual oil and gas show that the deep stratum of Tahe Oilfield still has exploration and development potential in the Ordovician Yingshan Formation and Penglaiba Formation,and the Middle and Upper Cambrian ultra-deep stratum has a certain oil and gas resource prospect.This study provides a reference for the dynamic quantitative evaluation of deep oil and gas in the Tarim Basin,and also provides a reference for the study of reservoir formation and evolution in carbonate reservoir of paleo-craton basin.展开更多
In order to ensure the penetrability of double-cased perforation in offshore oil and gas fields and to maximize the capacity of perforation completion, This study establishes a dynamic model of double-cased perforatio...In order to ensure the penetrability of double-cased perforation in offshore oil and gas fields and to maximize the capacity of perforation completion, This study establishes a dynamic model of double-cased perforation using ANSYS/LS-DYNA simulation technology. The combination of critical perforation parameters for double casing is obtained by studying the influencing factors of the jet-forming process,perforation depth, diameter, and stress changes of the inner and outer casing. The single-target perforation experiments under high-temperature and high-pressure(HTHP) conditions and ground full-scale ring target perforation tests are designed to verify the accuracy of numerical simulation results. The reduced factor is adopted as the quantitative measure of perforation depth and diameter for different types of perforation charge under different conditions. The results show that the perforation depth reduction increases with temperature and pressure, and the reduced factor is between 0.67 and 0.87 under HTHP conditions of 130℃/44 MPa and 137℃/60 MPa. Comparing the results of the numerical simulation and the full-scale test correction, the maximum error is less than 8.91%, and this numerical simulation has strong reliability. This research provides a basis for a reasonable range of double-cased perforation parameters and their optimal selection.展开更多
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.展开更多
Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced...Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced critical load and hardness compared to nc-Ni,where perfect,stair-rod and Shockley dislocations are activated at (111),(111) and (111) slip planes in nt-Ni compared to only SSockley dislocation nucleation at (111) and (111) slip planes of nc-Ni.The nt-Ni exhibits a less significant indentation size effect in comparison with nc-Ni due to the dislocation slips hindrance of the twin boundary.The atomic deformation associated with the indentation size effect is investigated during dislocation transmission.Different from the decreasing partial slips parallel to the indenter surface in nc-Ni with increasing temperature,the temperaturedependent atomic deformation of nt-Ni is closely related to the twin boundary:from the partial slips parallel to the twin boundary (~10 K),to increased confined layer slips and decreased twin migration(300 K–600 K),to decreased confined layer slips and increased dislocation interaction of dislocation pinning and dissociation (900 K–1200 K).Dislocation density and atomic structure types through quantitative analysis are implemented to further reveal the above-mentioned dislocation motion and atomic structure alteration.Our study is helpful for understanding the temperature-dependent plasticity of twin boundary in nanotwinned materials.展开更多
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.展开更多
We employ the Green–Kubo(G-K)and Einstein relations to estimate the self-diffusion coefficients(denoted as D_(G)and D_(E),respectively)in two-dimensional(2D)strongly coupled dusty plasmas(SC-DPs)via equilibrium molec...We employ the Green–Kubo(G-K)and Einstein relations to estimate the self-diffusion coefficients(denoted as D_(G)and D_(E),respectively)in two-dimensional(2D)strongly coupled dusty plasmas(SC-DPs)via equilibrium molecular dynamics(EMD)simulations.D_(G)and D_(E)are computed for a broad domain of screening length(κ)and coupling parameters(Г)along with different system sizes.It is observed that both D_(G)and D_(E)decrease linearly with increasing G in warm liquid states and increase with increasingκ.In cold liquid states,the Einstein relation accurately predicts D_(E)in 2D SC-DPs because diffusion motion is close to normal diffusion,but the G-K relation provides overestimations of D_(G),because VACF indicates anomalous diffusion;thus,D_(G)is not accurate.Our new simulation outcomes reveal that D_(G)and D_(E)remain independent of system sizes.Furthermore,our investigations demonstrate that at higher temperatures,D_(G)and D_(E)converge,suggesting diffusion motion close to normal diffusion,while at lower temperatures,these two values diverge.We find reasonable agreement by comparing current and existing numerical,theoretical and experimental data.Moreover,when normalizing diffusion coefficients by the Einstein frequency and testing against the universal temperature scaling law,D_(G)deviates from theoretical curves at low temperatures and k,whereas D_(E)only disagrees with theory at very smallκ(■0.10).These findings provide valuable insight into diagnosing dust component parameters within 2D DP systems and contribute to the broader understanding of diffusion processes in DP environments.展开更多
The track model used in the dynamic analysis and design system software is investigated. A home made tank is taken as an example to illustrate the method for modeling an integral tracked vehicle and perform the dynam...The track model used in the dynamic analysis and design system software is investigated. A home made tank is taken as an example to illustrate the method for modeling an integral tracked vehicle and perform the dynamic simulation. The obtained results have demonstrated that the simulation method has the advantage of high efficiency, more convenience and more insight into the dynamical behavior of the system.展开更多
In this paper molecular dynamics simulations are performed to study the accumulation behaviour of N2 and H2 at water/graphite interface under ambient temperature and pressure. It finds that both N2 and H2 molecules ca...In this paper molecular dynamics simulations are performed to study the accumulation behaviour of N2 and H2 at water/graphite interface under ambient temperature and pressure. It finds that both N2 and H2 molecules can accumulate at the interface and form one of two states according to the ratio of gas molecules number to square of graphite surface from our simulation results： gas films （pancake-like） for a larger ratio and nanobubbles for a smaller ratio. In addition, we discuss the stabilities of nanobubbles at different environment temperatures. Surprisingly, it is found that the density of both kinds of gas states can be greatly increased, even comparable with that of the liquid N2 and liquid H2. The present results are expected to be helpful for the understanding of the stable existence of gas film （pancake-like） and nanobubbles.展开更多
Selectivity control is a difficult scientific and industrial challenge in methanol-to-olefins(MTO)conversion.It has been experimentally established that the topology of zeolite catalysts influenced the distribution of...Selectivity control is a difficult scientific and industrial challenge in methanol-to-olefins(MTO)conversion.It has been experimentally established that the topology of zeolite catalysts influenced the distribution of products.Besides the topology effect on reaction kinetics,the topology influences the diffusion of reactants and products in catalysts as well.In this work,by using COMPASS force-field molecular dynamics method,we investigated the intracrystalline diffusion of ethene and propene in four different zeolites,CHA,MFI,BEA and FAU,at different temperatures.The self-diffusion coefficients and diffusion activation barriers were calculated.A strong restriction on the diffusion of propene in CHA was observed because the self-diffusion coefficient ratio of ethene to propene is larger than 18 and the diffusion activation barrier of propene is more than 20 kJ/mol in CHA.This ratio decreases with the increase of temperature in the four investigated zeolites.The shape selectivity on products from diffusion perspective can provide some implications on the understanding of the selectivity difference between HSAPO-34 and HZSM-5 catalysts for the MTO conversion.展开更多
This paper performs molecular dynamics simulations to investigate the role of the monovalent cations K, Na and the divalent cation Ca on the stability and swelling of montmorillonite. The recently developed CLAYFF for...This paper performs molecular dynamics simulations to investigate the role of the monovalent cations K, Na and the divalent cation Ca on the stability and swelling of montmorillonite. The recently developed CLAYFF force field is used to predict the basal spacing as a function of the water content in the interlayer. The simulations reproduced the swelling pattern of these montmorillonites, suggesting a mechanism of their hydration different （K＋ 〈 Na＋ 〈 Ca2＋） from that of K＋-, Na＋-, and Ca2＋-montmorillonites. In particular, these results indicate that the valence of the cations has the larger impact on the behaviour of clay water systems. It also finds that the differences in size and hydration energy of K＋, Na＋ and Ca2＋ ions have strong implications for the structure of interlayer. This leads to the differences in the layer spacings of the simulated K＋-, Na＋-, and Ca2＋-montmorillonites. Furthermore, these simulations show that the K cations interact strongly with the clay sheets for the dehydrated clay sheets, but for the hydrated clays the Ca cations interact clearly strongly with the clay sheets.展开更多
Molecular dynamics simulations are performed to study the growth mechanism of CH4-CO2 mixed hydrate in xco2 = 75%, xco2 = 50%, and zco2 = 25% systems at T = 250 K, 255 K and 260 K, respectively. Our simulation results...Molecular dynamics simulations are performed to study the growth mechanism of CH4-CO2 mixed hydrate in xco2 = 75%, xco2 = 50%, and zco2 = 25% systems at T = 250 K, 255 K and 260 K, respectively. Our simulation results show that the growth rate of CH4-CO2 mixed hydrate increases as the CO2 concentration in the initial solution phase increases and the temperature decreases. Via hydrate formation, the composition of CO2 in hydrate phase is higher than that in initial solution phase and the encaging capacity of CO2 in hydrates increases with the decrease in temperature. By analysis of the cage occupancy ratio of CH4 molecules and CO2 molecules in large cages to small cages, we find that CO2 molecules are preferably encaged into the large cages of the hydrate crystal as compared with CH4 molecules. Interestingly, CH4 molecules and CO2 molecules frequently replace with each other in some particular cage sites adjacent to hydrate/solution interface during the crystal growth process. These two species of guest molecules eventually act to stabilize the newly formed hydrates, with CO2 molecules occupying large cages and CH4 molecules occupying small cages in hydrate.展开更多
We use non-equilibrium molecular dynamics simulations to calculate the self-diffusion coefficient, D, of a Lennard Jones fluid over a wide density and temperature range. The change in self-diffusion coefficient with t...We use non-equilibrium molecular dynamics simulations to calculate the self-diffusion coefficient, D, of a Lennard Jones fluid over a wide density and temperature range. The change in self-diffusion coefficient with temperature decreases by increasing density. For density ρ＊ = ρσ3 = 0.84 we observe a peak at the value of the self-diffusion coefficient and the critical temperature T＊ = kT/ε = 1.25. The value of the self-diffusion coefficient strongly depends on system size. The data of the self-diffusion coefficient are fitted to a simple analytic relation based on hydrodynamic arguments. This correction scales as N-α, where α is an adjustable parameter and N is the number of particles. It is observed that the values of a 〈 1 provide quite a good correction to the simulation data. The system size dependence is very strong for lower densities, but it is not as strong for higher densities. The self-diffusion coefficient calculated with non-equilibrium molecular dynamic simulations at different temperatures and densities is in good agreement with other calculations fronl the literature.展开更多
Currently,wire bonding is the most popular first-level interconnection technology used between the die and package terminals,but even with its long-term and excessive usage,the mechanism of wire bonding has not been c...Currently,wire bonding is the most popular first-level interconnection technology used between the die and package terminals,but even with its long-term and excessive usage,the mechanism of wire bonding has not been completely evaluated.Therefore,fundamental research is still needed.In this study,the mechanism of microweld formation and breakage during Cu-Cu wire bonding was investigated by using molecular dynamics simulation.The contact model for the nanoindentation process between the wire and substrate was developed to simulate the contact process of the Cu wire and Cu substrate.Elastic contact and plastic instability were investigated through the loading and unloading processes.Moreover,the evolution of the indentation morphology and distributions of the atomic stress were also investigated.It was shown that the loading and unloading curves do not coincide,and the unloading curve exhibited hysteresis.For the substrate,in the loading process,the main force changed from attractive to repulsive.The maximum von Mises stress increased and shifted from the center toward the edge of the contact area.During the unloading process,the main force changed from repulsive to attractive.The Mises stress reduced first and then increased.Stress concentration occurs around dislocations in the middle area of the Cu wire.展开更多
A deep understanding of explosive sensitivities and their factors is important for safe and reliable applications.However,quantitative prediction of the sensitivities is difficult.Here,reactive molecular dynamics simu...A deep understanding of explosive sensitivities and their factors is important for safe and reliable applications.However,quantitative prediction of the sensitivities is difficult.Here,reactive molecular dynamics simulation models for high-speed piston impacts on explosive supercells were established.Simulations were also performed to investigate shock-induced reactions of various high-energy explosives.The fraction of reacted explosive molecules in an initial supercell changed linearly with the propagation distance of the shock-wave front.The corresponding slope could be used as a reaction rate for a specific shock-loading velocity.Reaction rates that varied with the shock-loading pressure exhibited two-stage linearities with different slopes.The two inflection points corresponded to the initial and accelerated reactions,which respectively correlated to the thresholds of shock-induced ignition and detonation.Therefore,the ignition and detonation critical pressures could be determined.The sensitivity could then be a quantitative prediction of the critical pressure.The accuracies of the quantitative shock sensitivity predictions were verified by comparing the impact and shock sensitivities of common explosives and the characteristics of anisotropic shock-induced reactions.Molecular dynamics simulations quantitatively predict and rank shock sensitivities by using only crystal structures of the explosives.Overall,this method will enable the design and safe use of explosives.展开更多
基金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.
基金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 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 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.
文摘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 Sichuan Province Regional Innovation Cooperation Project(21QYCX0048)Sinopec Science and Technology Department Project(P21048-3)。
文摘According to the complex differential accumulation history of deep marine oil and gas in superimposed basins,the Lower Paleozoic petroleum system in Tahe Oilfield of Tarim Basin is selected as a typical case,and the process of hydrocarbon generation and expulsion,migration and accumulation,adjustment and transformation of deep oil and gas is restored by means of reservoine-forming dynamics simulation.The thermal evolution history of the Lower Cambrian source rocks in Tahe Oilfield reflects the obvious differences in hydrocarbon generation and expulsion process and intensity in different tectonic zones,which is the main reason controlling the differences in deep oil and gas phases.The complex transport system composed of strike-slip fault and unconformity,etc.controlled early migration and accumulation and late adjustment of deep oil and gas,while the Middle Cambrian gypsum-salt rock in inner carbonate platform prevented vertical migration and accumulation of deep oil and gas,resulting in an obvious"fault-controlled"feature of deep oil and gas,in which the low potential area superimposed by the NE-strike-slip fault zone and deep oil and gas migration was conducive to accumulation,and it is mainly beaded along the strike-slip fault zone in the northeast direction.The dynamic simulation of reservoir formation reveals that the spatio-temporal configuration of"source-fault-fracture-gypsum-preservation"controls the differential accumulation of deep oil and gas in Tahe Oilfield.The Ordovician has experienced the accumulation history of multiple periods of charging,vertical migration and accumulation,and lateral adjustment and transformation,and deep oil and gas have always been in the dynamic equilibrium of migration,accumulation and escape.The statistics of residual oil and gas show that the deep stratum of Tahe Oilfield still has exploration and development potential in the Ordovician Yingshan Formation and Penglaiba Formation,and the Middle and Upper Cambrian ultra-deep stratum has a certain oil and gas resource prospect.This study provides a reference for the dynamic quantitative evaluation of deep oil and gas in the Tarim Basin,and also provides a reference for the study of reservoir formation and evolution in carbonate reservoir of paleo-craton basin.
基金the support of the Foundation of Natural Science Foundation of Shaanxi Province, Grant/ Award nos. 2023-JC-YB-361National Natural Science Foundation (Number 52104033)。
文摘In order to ensure the penetrability of double-cased perforation in offshore oil and gas fields and to maximize the capacity of perforation completion, This study establishes a dynamic model of double-cased perforation using ANSYS/LS-DYNA simulation technology. The combination of critical perforation parameters for double casing is obtained by studying the influencing factors of the jet-forming process,perforation depth, diameter, and stress changes of the inner and outer casing. The single-target perforation experiments under high-temperature and high-pressure(HTHP) conditions and ground full-scale ring target perforation tests are designed to verify the accuracy of numerical simulation results. The reduced factor is adopted as the quantitative measure of perforation depth and diameter for different types of perforation charge under different conditions. The results show that the perforation depth reduction increases with temperature and pressure, and the reduced factor is between 0.67 and 0.87 under HTHP conditions of 130℃/44 MPa and 137℃/60 MPa. Comparing the results of the numerical simulation and the full-scale test correction, the maximum error is less than 8.91%, and this numerical simulation has strong reliability. This research provides a basis for a reasonable range of double-cased perforation parameters and their optimal selection.
基金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 (Grant No.11572090)。
文摘Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced critical load and hardness compared to nc-Ni,where perfect,stair-rod and Shockley dislocations are activated at (111),(111) and (111) slip planes in nt-Ni compared to only SSockley dislocation nucleation at (111) and (111) slip planes of nc-Ni.The nt-Ni exhibits a less significant indentation size effect in comparison with nc-Ni due to the dislocation slips hindrance of the twin boundary.The atomic deformation associated with the indentation size effect is investigated during dislocation transmission.Different from the decreasing partial slips parallel to the indenter surface in nc-Ni with increasing temperature,the temperaturedependent atomic deformation of nt-Ni is closely related to the twin boundary:from the partial slips parallel to the twin boundary (~10 K),to increased confined layer slips and decreased twin migration(300 K–600 K),to decreased confined layer slips and increased dislocation interaction of dislocation pinning and dissociation (900 K–1200 K).Dislocation density and atomic structure types through quantitative analysis are implemented to further reveal the above-mentioned dislocation motion and atomic structure alteration.Our study is helpful for understanding the temperature-dependent plasticity of twin boundary in nanotwinned materials.
文摘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.
基金support of the Fundamental Research Funds for the Central Universities of China(Grant No.2019ZDPY16).
文摘We employ the Green–Kubo(G-K)and Einstein relations to estimate the self-diffusion coefficients(denoted as D_(G)and D_(E),respectively)in two-dimensional(2D)strongly coupled dusty plasmas(SC-DPs)via equilibrium molecular dynamics(EMD)simulations.D_(G)and D_(E)are computed for a broad domain of screening length(κ)and coupling parameters(Г)along with different system sizes.It is observed that both D_(G)and D_(E)decrease linearly with increasing G in warm liquid states and increase with increasingκ.In cold liquid states,the Einstein relation accurately predicts D_(E)in 2D SC-DPs because diffusion motion is close to normal diffusion,but the G-K relation provides overestimations of D_(G),because VACF indicates anomalous diffusion;thus,D_(G)is not accurate.Our new simulation outcomes reveal that D_(G)and D_(E)remain independent of system sizes.Furthermore,our investigations demonstrate that at higher temperatures,D_(G)and D_(E)converge,suggesting diffusion motion close to normal diffusion,while at lower temperatures,these two values diverge.We find reasonable agreement by comparing current and existing numerical,theoretical and experimental data.Moreover,when normalizing diffusion coefficients by the Einstein frequency and testing against the universal temperature scaling law,D_(G)deviates from theoretical curves at low temperatures and k,whereas D_(E)only disagrees with theory at very smallκ(■0.10).These findings provide valuable insight into diagnosing dust component parameters within 2D DP systems and contribute to the broader understanding of diffusion processes in DP environments.
文摘The track model used in the dynamic analysis and design system software is investigated. A home made tank is taken as an example to illustrate the method for modeling an integral tracked vehicle and perform the dynamic simulation. The obtained results have demonstrated that the simulation method has the advantage of high efficiency, more convenience and more insight into the dynamical behavior of the system.
基金supported in part by National Natural Science Foundation of China (Grant Nos 10474109 and 10674146)supported is part by the Shanghai Supercomputer Center of China
文摘In this paper molecular dynamics simulations are performed to study the accumulation behaviour of N2 and H2 at water/graphite interface under ambient temperature and pressure. It finds that both N2 and H2 molecules can accumulate at the interface and form one of two states according to the ratio of gas molecules number to square of graphite surface from our simulation results： gas films （pancake-like） for a larger ratio and nanobubbles for a smaller ratio. In addition, we discuss the stabilities of nanobubbles at different environment temperatures. Surprisingly, it is found that the density of both kinds of gas states can be greatly increased, even comparable with that of the liquid N2 and liquid H2. The present results are expected to be helpful for the understanding of the stable existence of gas film （pancake-like） and nanobubbles.
基金supported by the National Basic Research Program of China (2009CB623504)the National Science Foundation of China (21103231)Shanghai Science Foundation (11ZR1449700)
文摘Selectivity control is a difficult scientific and industrial challenge in methanol-to-olefins(MTO)conversion.It has been experimentally established that the topology of zeolite catalysts influenced the distribution of products.Besides the topology effect on reaction kinetics,the topology influences the diffusion of reactants and products in catalysts as well.In this work,by using COMPASS force-field molecular dynamics method,we investigated the intracrystalline diffusion of ethene and propene in four different zeolites,CHA,MFI,BEA and FAU,at different temperatures.The self-diffusion coefficients and diffusion activation barriers were calculated.A strong restriction on the diffusion of propene in CHA was observed because the self-diffusion coefficient ratio of ethene to propene is larger than 18 and the diffusion activation barrier of propene is more than 20 kJ/mol in CHA.This ratio decreases with the increase of temperature in the four investigated zeolites.The shape selectivity on products from diffusion perspective can provide some implications on the understanding of the selectivity difference between HSAPO-34 and HZSM-5 catalysts for the MTO conversion.
基金Project supported by the Key Laboratory of Mountain Hazards and Earth Surface Processes, the Chinese Academy of Sciences
文摘This paper performs molecular dynamics simulations to investigate the role of the monovalent cations K, Na and the divalent cation Ca on the stability and swelling of montmorillonite. The recently developed CLAYFF force field is used to predict the basal spacing as a function of the water content in the interlayer. The simulations reproduced the swelling pattern of these montmorillonites, suggesting a mechanism of their hydration different （K＋ 〈 Na＋ 〈 Ca2＋） from that of K＋-, Na＋-, and Ca2＋-montmorillonites. In particular, these results indicate that the valence of the cations has the larger impact on the behaviour of clay water systems. It also finds that the differences in size and hydration energy of K＋, Na＋ and Ca2＋ ions have strong implications for the structure of interlayer. This leads to the differences in the layer spacings of the simulated K＋-, Na＋-, and Ca2＋-montmorillonites. Furthermore, these simulations show that the K cations interact strongly with the clay sheets for the dehydrated clay sheets, but for the hydrated clays the Ca cations interact clearly strongly with the clay sheets.
基金supported by the National Natural Science Foundation of China(No.51176192)CAS Program(KGZD-EW-301)NOG Program(GHZ2012006003)
文摘Molecular dynamics simulations are performed to study the growth mechanism of CH4-CO2 mixed hydrate in xco2 = 75%, xco2 = 50%, and zco2 = 25% systems at T = 250 K, 255 K and 260 K, respectively. Our simulation results show that the growth rate of CH4-CO2 mixed hydrate increases as the CO2 concentration in the initial solution phase increases and the temperature decreases. Via hydrate formation, the composition of CO2 in hydrate phase is higher than that in initial solution phase and the encaging capacity of CO2 in hydrates increases with the decrease in temperature. By analysis of the cage occupancy ratio of CH4 molecules and CO2 molecules in large cages to small cages, we find that CO2 molecules are preferably encaged into the large cages of the hydrate crystal as compared with CH4 molecules. Interestingly, CH4 molecules and CO2 molecules frequently replace with each other in some particular cage sites adjacent to hydrate/solution interface during the crystal growth process. These two species of guest molecules eventually act to stabilize the newly formed hydrates, with CO2 molecules occupying large cages and CH4 molecules occupying small cages in hydrate.
基金supported by the National Natural Science Foundation of China (Grant No. 51076128)the National High Technology Research and Development Program of China (Grant No. 2009AA05Z107)
文摘We use non-equilibrium molecular dynamics simulations to calculate the self-diffusion coefficient, D, of a Lennard Jones fluid over a wide density and temperature range. The change in self-diffusion coefficient with temperature decreases by increasing density. For density ρ＊ = ρσ3 = 0.84 we observe a peak at the value of the self-diffusion coefficient and the critical temperature T＊ = kT/ε = 1.25. The value of the self-diffusion coefficient strongly depends on system size. The data of the self-diffusion coefficient are fitted to a simple analytic relation based on hydrodynamic arguments. This correction scales as N-α, where α is an adjustable parameter and N is the number of particles. It is observed that the values of a 〈 1 provide quite a good correction to the simulation data. The system size dependence is very strong for lower densities, but it is not as strong for higher densities. The self-diffusion coefficient calculated with non-equilibrium molecular dynamic simulations at different temperatures and densities is in good agreement with other calculations fronl the literature.
基金the National Key R&D Program of China(Grant No.2019YFB1704600)the Hubei Provincial Natural Science Foundation of China(Grant No.2020CFA032).
文摘Currently,wire bonding is the most popular first-level interconnection technology used between the die and package terminals,but even with its long-term and excessive usage,the mechanism of wire bonding has not been completely evaluated.Therefore,fundamental research is still needed.In this study,the mechanism of microweld formation and breakage during Cu-Cu wire bonding was investigated by using molecular dynamics simulation.The contact model for the nanoindentation process between the wire and substrate was developed to simulate the contact process of the Cu wire and Cu substrate.Elastic contact and plastic instability were investigated through the loading and unloading processes.Moreover,the evolution of the indentation morphology and distributions of the atomic stress were also investigated.It was shown that the loading and unloading curves do not coincide,and the unloading curve exhibited hysteresis.For the substrate,in the loading process,the main force changed from attractive to repulsive.The maximum von Mises stress increased and shifted from the center toward the edge of the contact area.During the unloading process,the main force changed from repulsive to attractive.The Mises stress reduced first and then increased.Stress concentration occurs around dislocations in the middle area of the Cu wire.
基金supported by the National Natural Science Foundation of China(Grant No.11832006).
文摘A deep understanding of explosive sensitivities and their factors is important for safe and reliable applications.However,quantitative prediction of the sensitivities is difficult.Here,reactive molecular dynamics simulation models for high-speed piston impacts on explosive supercells were established.Simulations were also performed to investigate shock-induced reactions of various high-energy explosives.The fraction of reacted explosive molecules in an initial supercell changed linearly with the propagation distance of the shock-wave front.The corresponding slope could be used as a reaction rate for a specific shock-loading velocity.Reaction rates that varied with the shock-loading pressure exhibited two-stage linearities with different slopes.The two inflection points corresponded to the initial and accelerated reactions,which respectively correlated to the thresholds of shock-induced ignition and detonation.Therefore,the ignition and detonation critical pressures could be determined.The sensitivity could then be a quantitative prediction of the critical pressure.The accuracies of the quantitative shock sensitivity predictions were verified by comparing the impact and shock sensitivities of common explosives and the characteristics of anisotropic shock-induced reactions.Molecular dynamics simulations quantitatively predict and rank shock sensitivities by using only crystal structures of the explosives.Overall,this method will enable the design and safe use of explosives.
