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.展开更多
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.展开更多
The effects of tensile temperatures ranging from 100 K to 900 K on the phase transition of hexagonal close-packed(HCP)zirconium were investigated by molecular dynamics simulations,which were combined with experimental...The effects of tensile temperatures ranging from 100 K to 900 K on the phase transition of hexagonal close-packed(HCP)zirconium were investigated by molecular dynamics simulations,which were combined with experimental observation under high resolution transmission electron microscopy.The results show that externally applied loading first induced the HCP to body-centered cubic(BCC)phase transition in the Pitsch-Schrader(PS)orientation relationship(OR).Then,the face-centered cubic(FCC)structure transformed from the BCC phase in the Bain path.However,the HCP-to-BCC transition was incomplete at 100 K and 300 K,resulting in a prismatic-type OR between the FCC and original HCP phase.Additionally,at the temperature ranging from 100 K to 600 K,the inverse BCC-to-HCP transition occurred locally following other variants of the PS OR,resulting in a basal-type relation between the newly generated HCP and FCC phases.A higher tensile temperature promoted the amount of FCC phase transforming into the BCC phase when the strain exceeded 45%.Besides,the crystal stretched at lower temperatures exhibits relatively higher strength but by the compromise of plasticity.This study reveals the deformation mechanisms in HCP-Zr at different temperatures,which may provide a better understanding of the deformation mechanism of zirconium alloys under different application environments.展开更多
Molecular dynamics simulations are performed to observe the evolutions of 512 and 51262 cage-like water clusters filled with or without a methane molecule immersed in bulk liquid water at 250 K and 230 K. The lifetime...Molecular dynamics simulations are performed to observe the evolutions of 512 and 51262 cage-like water clusters filled with or without a methane molecule immersed in bulk liquid water at 250 K and 230 K. The lifetimes of these clusters are calculated according to their Lindemann index δ (t) using the criteria of δ≥0.07. For both the filled and empty clusters, we find the dynamics of bulk water determines the lifetimes of cage-like water clusters, and that the lifetime of 512 62 cage-like cluster is the same as that of 512 cage-like cluster. Although the methane molecule indeed makes the filled cage-like cluster more stable than the empty one, the empty cage-like cluster still has chance to be long-lived compared with the filled clusters. These observations support the labile cluster hypothesis on the formation mechanisms of gas hydrates.展开更多
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.展开更多
The characteristic of easy sintering of aluminum nanoparticle(ANP)limits its application in solid propellants.Coating ANP with fluoropolymer could effectively improve its combustion performance.To find out how the coa...The characteristic of easy sintering of aluminum nanoparticle(ANP)limits its application in solid propellants.Coating ANP with fluoropolymer could effectively improve its combustion performance.To find out how the coating layer inhibits sintering and promotes complete combustion of particles from an atomic view,a comparative study has been done for bare ANP and PTFE coated ANP by using reactive molecular dynamics simulations.The sintering process is quantified by shrinkage ratio and gyration radius.Our results show that,at the same heating rate and combustion temperatures,bare ANPs are sintered together after the temperature exceeds the melting point of aluminum but the decomposition of PTFE coating layer pushes particles away and increases reaction surface area by producing small Al-F clusters.The sintering of ANPs which are heated in PTFE is alleviated compared with particles heated in oxygen,but particles still sinter together due to the lack of intimate contact between PTFE and alumina surface.The effect of temperature on the combustion of PTFE coated ANPs is also studied from 1000 to3500 K.The number density analysis shows the particles will not be sintered at any temperature.Aluminum fluoride prefers diffusing to the external space and the remained particles are mainly composed of Al,C and O.Fast ignition simulations are performed by adopting micro canonical ensemble.With the expansion of aluminum core and the melting of alumina shell,bare ANPs are sintered into a liquid particle directly.For PTFE coated ANPs,the volatilization of gaseous aluminum fluoride products continually endows particles opposite momentum.展开更多
The interface properties between hydrated cement paste(hcp)and aggregates largely determine the various performances of concrete.In this work,molecular dynamics simulations were employed to explore the atomistic inter...The interface properties between hydrated cement paste(hcp)and aggregates largely determine the various performances of concrete.In this work,molecular dynamics simulations were employed to explore the atomistic interaction mechanisms between the commonly used aggregate phase calcite/silica and calcium silicate hydrates(C-S-H),as well as the effect of moisture.The results suggest that the C-S-H/calcite interface is relatively strong and stable under both dry and moist conditions,which is caused by the high-strength interfacial connections formed between calcium ions from calcite and high-polarity non-bridging oxygen atoms from the C-S-H surface.Silica can be also adsorbed on the dry C-S-H surface by the H-bonds;however,the presence of water molecules on the interface may substantially decrease the affinities.Furthermore,the dynamics interface separation tests of C-S-H/aggregates were also implemented by molecular dynamics.The shape of the calculated stress-separation distance curves obeys the quasi-static cohesive law obtained experimentally.The moisture conditions and strain rates were found to affect the separation process of C-S-H/silica.A wetter interface and smaller loading rate may lead to a lower adhesion strength.The mechanisms interpreted here may shed new lights on the understandings of hcp/aggregate interactions at a nano-length scale and creation of high performance cementitious materials.展开更多
Background Polygalacturonase inhibiting proteins(PGIPs)play a pivotal role in plant defense against plant patho-gens by inhibiting polygalacturonase(PG),an enzyme produced by pathogens to degrade plant cell wall pecti...Background Polygalacturonase inhibiting proteins(PGIPs)play a pivotal role in plant defense against plant patho-gens by inhibiting polygalacturonase(PG),an enzyme produced by pathogens to degrade plant cell wall pectin.PGIPs,also known as leucine-rich repeat pathogenesis-related(PR)proteins,activate the host’s defense response upon interaction with PG,thereby reinforcing the host defense against plant pathogens attacks.In Egyptian or extra-long staple cotton(Gossypium barbadense),the interaction between PGIP and PG is one of the crucial steps in the defense mechanism against major pathogens such as Xanthomonas citri pv.malvacearum and Alternaria mac-rospora,which are responsible for bacterial leaf blight and leaf spot diseases,respectively.Results To unravel the molecular mechanisms underlying these PR proteins,we conducted a comprehensive study involving molecular modeling,protein-protein docking,site-specific double mutation(E169G and F242K),and molec-ular dynamics simulations.Both wild-type and mutated cotton PGIPs were examined in the interaction with the PG enzyme of a bacterial and fungal pathogen.Our findings revealed that changes in conformations of double-mutated residues in the active site of PGIP lead to the inhibition of PG binding.The molecular dynamics simulation studies provide insights into the dynamic behaviour and stability of the PGIP-PG complexes,shedding light on the intricate details of the inhibitory and exhibitory mechanism against the major fungal and bacterial pathogens of G.barbadense,respectively.Conclusions The findings of this study not only enhance our understanding of the molecular interactions between PGs of Xanthomonas citri pv.malvacearum and Alternaria macrospora and PGIP of G.barbadense but also pre-sent a potential strategy for developing the disease-resistant cotton varieties.By variations in the binding affinities of PGs through specific mutations in PGIP,this research offers promising avenues for the development of enhanced resistance to cotton plants against bacterial leaf blight and leaf spot diseases.展开更多
The orientation effect of θʹ precipitates in stress-aged Al-Cu alloys has ambiguous interpretations. One view is that θʹ precipitates prefer to grow on the habit planes perpendicular to the applied compressive stres...The orientation effect of θʹ precipitates in stress-aged Al-Cu alloys has ambiguous interpretations. One view is that θʹ precipitates prefer to grow on the habit planes perpendicular to the applied compressive stress, while the other view suggests growth on habit planes parallel to the applied stress. In this study, stress-aged Al-4 wt.%Cu single crystal was sampled from three different <100>Al zone axes to observe the distribution of θʹ precipitates. A hybrid Monte-Carlo/ molecular dynamics simulations were used to investigate the orientation effect of θʹ precipitates. The simulation results are consistent with experimental observations and indicate that θʹ precipitates prefer to grow on the habit planes that are parallel to the direction of the applied compressive stress, not along the planes perpendicular to it. It is also found that 1/2<110> perfect dislocations are generated as θʹ precipitates plates grow thicker, and the reaction of 1/2<110> perfect dislocations decomposing into 1/6<112> Shockley dislocations lead to an increase in the length of θʹ precipitates. The former does not enhance the orientation effect, whereas the latter leads to a more significant orientation effect. Additionally, the degree of the orientation effect of θʹ precipitates is determined by the reduction of 1/2<110> dislocations with a positive correlation between them.展开更多
In the present study,molecular dynamic simulation(MD)was used to investigate the plastic deformation process of the Fe-Mn alloys with different Mn contents.The influences of Mn contents ranging from 10%to 30%(at%)on t...In the present study,molecular dynamic simulation(MD)was used to investigate the plastic deformation process of the Fe-Mn alloys with different Mn contents.The influences of Mn contents ranging from 10%to 30%(at%)on the deformation behavior and the controlling mechanism of the Fe-base alloys were analyzed.The results show that phase transformations and{112}<111>_(BCC)deformation twinning occur in all Fe-Mn alloys but follow different deformation paths.In the Fe-10%Mn alloy the deformation twinning mechanism obeys the FCC-related path,the Fe-20%Mn alloy involves both the FCC-and HCP-related paths,and the deformation of the Fe-30%Mn alloy is dominated by the HCP-related twinning path.The addition of Mn can increase the stacking fault energy and retard the activation of slip systems as well as the formation of stacking faults.Thus,a higher content of Mn can delay the FCC®ε-martensite and the subsequentε-martensite®BCC phase transition at the intersection of twoε-martensitic bands.Therefore,the addition of Mn alloying element increases the yield strength and reduces the elastic modulus of the Fe-Mn alloys.The formation of deformation twins will contribute to the work-hardening effect and delay the necking and fracture of alloys.It is expected that the results in the present study will provide theoretical reference for the design and optimization of high-performance steels.展开更多
Nanograins are characterized by a typical grain size from 1 to 100 nm. Molecular dynamics simulations have been carried out for the nanograin sphere with the diameters from 1.45 to 10.12 nm. We study the influence of ...Nanograins are characterized by a typical grain size from 1 to 100 nm. Molecular dynamics simulations have been carried out for the nanograin sphere with the diameters from 1.45 to 10.12 nm. We study the influence of grain size on structure and diffusion properties of the nanograins. The results reveal that as the grain size is reduced, the fraction of grain surface increases significantly, and the surface width is approximately constant; the mean atomic energy of the surface increases distinctly, but that of the grain interior varies insignificantly; the diffusion coefficient is increased sharply, and the relation of the diffusion coefficient and the grain size is close to exponential relation below 10 nm.展开更多
The distribution of Al （j） and the structural units distribution of Qi T in calcium aluminosilicate melts were studied by means of molecular dynamics simulation. The results show that provided there exists lower-fie...The distribution of Al （j） and the structural units distribution of Qi T in calcium aluminosilicate melts were studied by means of molecular dynamics simulation. The results show that provided there exists lower-field strength cation relative to Al3+, such as alkaline and alkaline earth metals, Al will be four-coordinated but not six-coordinated. Meanwhile, if there exist a large number of higher-field strength cations such as Si4+ and little lower-field strength cation, six-coordinated aluminum will be formed. The relation of structural units distribution of Qi T with chemical composition shift was also extracted, showing that as Ca2+ exists, the distributions of Qi Si, Qi Al or Qi T have the similar changing trend with the variation of component. Because of high-temperature effect, the Al-tetrahedral units in melts are greatly active and unstable and there exist dynamic transforming equilibria of Al（3）Al（4） and （Al（5））Al（4）. The three-coordinated oxygen and charge-compensated bridging oxygen are proposed to explain phenomena of the negative charge redundancy of AlO4 and location of network modifier with charge-compensated function in aluminosilicate melts.展开更多
The shock-induced reaction mechanism and characteristics of Ni/Al system,considering an Al nanoparticle-embedded Ni single crystal,are investigated through molecular dynamics simulation.For the shock melting of Al nan...The shock-induced reaction mechanism and characteristics of Ni/Al system,considering an Al nanoparticle-embedded Ni single crystal,are investigated through molecular dynamics simulation.For the shock melting of Al nanoparticle,interfacial crystallization and dissolution are the main characteristics.The reaction degree of Al particle first increases linearly and then logarithmically with time driven by rapid mechanical mixing and following dissolution.The reaction rate increases with the decrease of particle diameter,however,the reaction is seriously hindered by interfacial crystallization when the diameter is lower than 9 nm in our simulations.Meanwhile,we found a negative exponential growth in the fraction of crystallized Al atoms,and the crystallinity of B2-NiAl(up to 20%)is positively correlated with the specific surface area of Al particle.This can be attributed to the formation mechanism of B2-NiAl by structural evolution of finite mixing layer near the collapsed interface.For shock melting of both Al particle and Ni matrix,the liquid-liquid phase inter-diffusion is the main reaction mechanism that can be enhanced by the formation of internal jet.In addition,the enhanced diffusion is manifested in the logarithmic growth law of mean square displacement,which results in an almost constant reaction rate similar to the mechanical mixing process.展开更多
To improve the thermal properties of aluminum(Al)in the energetic system,a coated structure with ammonium perchlorate(AP)was prepared by a facile approach.And N,N-Dimethylformamide(DMF)was chosen as an ideal solvent b...To improve the thermal properties of aluminum(Al)in the energetic system,a coated structure with ammonium perchlorate(AP)was prepared by a facile approach.And N,N-Dimethylformamide(DMF)was chosen as an ideal solvent based on heterogeneous nucleation theory and molecular dynamics simulation.This coated structure could enlarge the contact area and improve the reaction environment to enhance the thermal properties.The addition of AP could accelerate oxidation temperature of Al with around 17.5°C.And the heat release of 85@15 composition rises to 26.13 k J/g and the reaction degree is97.6%with higher peak pressure(254.6 k Pa)and rise rate(1.397 MPa/s).An ideal ratio with 15 wt%AP was probed primarily.The high energy laser-induced shockwave experiment was utilized to simulate the reaction behavior in hot field.And the larger activated mixture of coated powder could release more energy to promote the growth of shockwave with higher speed up to 518.7±55.9 m/s.In conclusion,85@15 composition is expected to be applied in energetic system as a novel metal fuel.展开更多
基金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 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.
基金Projects(51901248,51828102)supported by the National Natural Science Foundation of ChinaProject(2018JJ3649)supported by the Natural Science Foundation of Hunan Province,ChinaProject(2019CX026)supported by the Innovation-driven Plan in Central South University,China。
文摘The effects of tensile temperatures ranging from 100 K to 900 K on the phase transition of hexagonal close-packed(HCP)zirconium were investigated by molecular dynamics simulations,which were combined with experimental observation under high resolution transmission electron microscopy.The results show that externally applied loading first induced the HCP to body-centered cubic(BCC)phase transition in the Pitsch-Schrader(PS)orientation relationship(OR).Then,the face-centered cubic(FCC)structure transformed from the BCC phase in the Bain path.However,the HCP-to-BCC transition was incomplete at 100 K and 300 K,resulting in a prismatic-type OR between the FCC and original HCP phase.Additionally,at the temperature ranging from 100 K to 600 K,the inverse BCC-to-HCP transition occurred locally following other variants of the PS OR,resulting in a basal-type relation between the newly generated HCP and FCC phases.A higher tensile temperature promoted the amount of FCC phase transforming into the BCC phase when the strain exceeded 45%.Besides,the crystal stretched at lower temperatures exhibits relatively higher strength but by the compromise of plasticity.This study reveals the deformation mechanisms in HCP-Zr at different temperatures,which may provide a better understanding of the deformation mechanism of zirconium alloys under different application environments.
基金supported by the National Natural Science Foundation of China(Grant No.40102005 and No.49725205).
文摘Molecular dynamics simulations are performed to observe the evolutions of 512 and 51262 cage-like water clusters filled with or without a methane molecule immersed in bulk liquid water at 250 K and 230 K. The lifetimes of these clusters are calculated according to their Lindemann index δ (t) using the criteria of δ≥0.07. For both the filled and empty clusters, we find the dynamics of bulk water determines the lifetimes of cage-like water clusters, and that the lifetime of 512 62 cage-like cluster is the same as that of 512 cage-like cluster. Although the methane molecule indeed makes the filled cage-like cluster more stable than the empty one, the empty cage-like cluster still has chance to be long-lived compared with the filled clusters. These observations support the labile cluster hypothesis on the formation mechanisms of gas hydrates.
文摘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.
基金supported by the fellowship of China Postdoctoral Science Foundation(Grant No.2021TQ0267)。
文摘The characteristic of easy sintering of aluminum nanoparticle(ANP)limits its application in solid propellants.Coating ANP with fluoropolymer could effectively improve its combustion performance.To find out how the coating layer inhibits sintering and promotes complete combustion of particles from an atomic view,a comparative study has been done for bare ANP and PTFE coated ANP by using reactive molecular dynamics simulations.The sintering process is quantified by shrinkage ratio and gyration radius.Our results show that,at the same heating rate and combustion temperatures,bare ANPs are sintered together after the temperature exceeds the melting point of aluminum but the decomposition of PTFE coating layer pushes particles away and increases reaction surface area by producing small Al-F clusters.The sintering of ANPs which are heated in PTFE is alleviated compared with particles heated in oxygen,but particles still sinter together due to the lack of intimate contact between PTFE and alumina surface.The effect of temperature on the combustion of PTFE coated ANPs is also studied from 1000 to3500 K.The number density analysis shows the particles will not be sintered at any temperature.Aluminum fluoride prefers diffusing to the external space and the remained particles are mainly composed of Al,C and O.Fast ignition simulations are performed by adopting micro canonical ensemble.With the expansion of aluminum core and the melting of alumina shell,bare ANPs are sintered into a liquid particle directly.For PTFE coated ANPs,the volatilization of gaseous aluminum fluoride products continually endows particles opposite momentum.
基金Projects(6512009004A,51908119,U1706222)supported by the National Natural Science Foundation of ChinaProject(BK20190367)supported by the Natural Science Foundation of Jiangsu Province,China。
文摘The interface properties between hydrated cement paste(hcp)and aggregates largely determine the various performances of concrete.In this work,molecular dynamics simulations were employed to explore the atomistic interaction mechanisms between the commonly used aggregate phase calcite/silica and calcium silicate hydrates(C-S-H),as well as the effect of moisture.The results suggest that the C-S-H/calcite interface is relatively strong and stable under both dry and moist conditions,which is caused by the high-strength interfacial connections formed between calcium ions from calcite and high-polarity non-bridging oxygen atoms from the C-S-H surface.Silica can be also adsorbed on the dry C-S-H surface by the H-bonds;however,the presence of water molecules on the interface may substantially decrease the affinities.Furthermore,the dynamics interface separation tests of C-S-H/aggregates were also implemented by molecular dynamics.The shape of the calculated stress-separation distance curves obeys the quasi-static cohesive law obtained experimentally.The moisture conditions and strain rates were found to affect the separation process of C-S-H/silica.A wetter interface and smaller loading rate may lead to a lower adhesion strength.The mechanisms interpreted here may shed new lights on the understandings of hcp/aggregate interactions at a nano-length scale and creation of high performance cementitious materials.
基金CABin grant(F.no.Agril.Edn.4-1/2013-A&P)Indian Council of Agricul-tural Research,Ministry of Agriculture and Farmers’Welfare,Govt.of India and Department of Biotechnology,Govt.of India for BIC project grant(BT/PR40161/BTIS/137/32/2021)。
文摘Background Polygalacturonase inhibiting proteins(PGIPs)play a pivotal role in plant defense against plant patho-gens by inhibiting polygalacturonase(PG),an enzyme produced by pathogens to degrade plant cell wall pectin.PGIPs,also known as leucine-rich repeat pathogenesis-related(PR)proteins,activate the host’s defense response upon interaction with PG,thereby reinforcing the host defense against plant pathogens attacks.In Egyptian or extra-long staple cotton(Gossypium barbadense),the interaction between PGIP and PG is one of the crucial steps in the defense mechanism against major pathogens such as Xanthomonas citri pv.malvacearum and Alternaria mac-rospora,which are responsible for bacterial leaf blight and leaf spot diseases,respectively.Results To unravel the molecular mechanisms underlying these PR proteins,we conducted a comprehensive study involving molecular modeling,protein-protein docking,site-specific double mutation(E169G and F242K),and molec-ular dynamics simulations.Both wild-type and mutated cotton PGIPs were examined in the interaction with the PG enzyme of a bacterial and fungal pathogen.Our findings revealed that changes in conformations of double-mutated residues in the active site of PGIP lead to the inhibition of PG binding.The molecular dynamics simulation studies provide insights into the dynamic behaviour and stability of the PGIP-PG complexes,shedding light on the intricate details of the inhibitory and exhibitory mechanism against the major fungal and bacterial pathogens of G.barbadense,respectively.Conclusions The findings of this study not only enhance our understanding of the molecular interactions between PGs of Xanthomonas citri pv.malvacearum and Alternaria macrospora and PGIP of G.barbadense but also pre-sent a potential strategy for developing the disease-resistant cotton varieties.By variations in the binding affinities of PGs through specific mutations in PGIP,this research offers promising avenues for the development of enhanced resistance to cotton plants against bacterial leaf blight and leaf spot diseases.
基金Project(2023YFB3710503) supported by the National Key R&D Program of ChinaProject(52305439) supported by the National Natural Science Foundation of China。
文摘The orientation effect of θʹ precipitates in stress-aged Al-Cu alloys has ambiguous interpretations. One view is that θʹ precipitates prefer to grow on the habit planes perpendicular to the applied compressive stress, while the other view suggests growth on habit planes parallel to the applied stress. In this study, stress-aged Al-4 wt.%Cu single crystal was sampled from three different <100>Al zone axes to observe the distribution of θʹ precipitates. A hybrid Monte-Carlo/ molecular dynamics simulations were used to investigate the orientation effect of θʹ precipitates. The simulation results are consistent with experimental observations and indicate that θʹ precipitates prefer to grow on the habit planes that are parallel to the direction of the applied compressive stress, not along the planes perpendicular to it. It is also found that 1/2<110> perfect dislocations are generated as θʹ precipitates plates grow thicker, and the reaction of 1/2<110> perfect dislocations decomposing into 1/6<112> Shockley dislocations lead to an increase in the length of θʹ precipitates. The former does not enhance the orientation effect, whereas the latter leads to a more significant orientation effect. Additionally, the degree of the orientation effect of θʹ precipitates is determined by the reduction of 1/2<110> dislocations with a positive correlation between them.
基金Project(51901248)supported by the National Natural Science Foundation of ChinaProject(2023JJ40742)supported by the Natural Science Foundation of Hunan Province,China。
文摘In the present study,molecular dynamic simulation(MD)was used to investigate the plastic deformation process of the Fe-Mn alloys with different Mn contents.The influences of Mn contents ranging from 10%to 30%(at%)on the deformation behavior and the controlling mechanism of the Fe-base alloys were analyzed.The results show that phase transformations and{112}<111>_(BCC)deformation twinning occur in all Fe-Mn alloys but follow different deformation paths.In the Fe-10%Mn alloy the deformation twinning mechanism obeys the FCC-related path,the Fe-20%Mn alloy involves both the FCC-and HCP-related paths,and the deformation of the Fe-30%Mn alloy is dominated by the HCP-related twinning path.The addition of Mn can increase the stacking fault energy and retard the activation of slip systems as well as the formation of stacking faults.Thus,a higher content of Mn can delay the FCC®ε-martensite and the subsequentε-martensite®BCC phase transition at the intersection of twoε-martensitic bands.Therefore,the addition of Mn alloying element increases the yield strength and reduces the elastic modulus of the Fe-Mn alloys.The formation of deformation twins will contribute to the work-hardening effect and delay the necking and fracture of alloys.It is expected that the results in the present study will provide theoretical reference for the design and optimization of high-performance steels.
基金supported by the National Naturl Science Foundation of China(No.10172088)the Potdoctoral Saience Foundation of China
文摘Nanograins are characterized by a typical grain size from 1 to 100 nm. Molecular dynamics simulations have been carried out for the nanograin sphere with the diameters from 1.45 to 10.12 nm. We study the influence of grain size on structure and diffusion properties of the nanograins. The results reveal that as the grain size is reduced, the fraction of grain surface increases significantly, and the surface width is approximately constant; the mean atomic energy of the surface increases distinctly, but that of the grain interior varies insignificantly; the diffusion coefficient is increased sharply, and the relation of the diffusion coefficient and the grain size is close to exponential relation below 10 nm.
文摘The distribution of Al （j） and the structural units distribution of Qi T in calcium aluminosilicate melts were studied by means of molecular dynamics simulation. The results show that provided there exists lower-field strength cation relative to Al3+, such as alkaline and alkaline earth metals, Al will be four-coordinated but not six-coordinated. Meanwhile, if there exist a large number of higher-field strength cations such as Si4+ and little lower-field strength cation, six-coordinated aluminum will be formed. The relation of structural units distribution of Qi T with chemical composition shift was also extracted, showing that as Ca2+ exists, the distributions of Qi Si, Qi Al or Qi T have the similar changing trend with the variation of component. Because of high-temperature effect, the Al-tetrahedral units in melts are greatly active and unstable and there exist dynamic transforming equilibria of Al（3）Al（4） and （Al（5））Al（4）. The three-coordinated oxygen and charge-compensated bridging oxygen are proposed to explain phenomena of the negative charge redundancy of AlO4 and location of network modifier with charge-compensated function in aluminosilicate melts.
基金supported by the State Key Program of National Natural Science Foundation of China(Grant No.12132003)State Key Laboratory of Explosion Science and Technology(Grant No.QNKT20-07)。
文摘The shock-induced reaction mechanism and characteristics of Ni/Al system,considering an Al nanoparticle-embedded Ni single crystal,are investigated through molecular dynamics simulation.For the shock melting of Al nanoparticle,interfacial crystallization and dissolution are the main characteristics.The reaction degree of Al particle first increases linearly and then logarithmically with time driven by rapid mechanical mixing and following dissolution.The reaction rate increases with the decrease of particle diameter,however,the reaction is seriously hindered by interfacial crystallization when the diameter is lower than 9 nm in our simulations.Meanwhile,we found a negative exponential growth in the fraction of crystallized Al atoms,and the crystallinity of B2-NiAl(up to 20%)is positively correlated with the specific surface area of Al particle.This can be attributed to the formation mechanism of B2-NiAl by structural evolution of finite mixing layer near the collapsed interface.For shock melting of both Al particle and Ni matrix,the liquid-liquid phase inter-diffusion is the main reaction mechanism that can be enhanced by the formation of internal jet.In addition,the enhanced diffusion is manifested in the logarithmic growth law of mean square displacement,which results in an almost constant reaction rate similar to the mechanical mixing process.
基金supported by National Natural Science Foundation of China[No.21975024]。
文摘To improve the thermal properties of aluminum(Al)in the energetic system,a coated structure with ammonium perchlorate(AP)was prepared by a facile approach.And N,N-Dimethylformamide(DMF)was chosen as an ideal solvent based on heterogeneous nucleation theory and molecular dynamics simulation.This coated structure could enlarge the contact area and improve the reaction environment to enhance the thermal properties.The addition of AP could accelerate oxidation temperature of Al with around 17.5°C.And the heat release of 85@15 composition rises to 26.13 k J/g and the reaction degree is97.6%with higher peak pressure(254.6 k Pa)and rise rate(1.397 MPa/s).An ideal ratio with 15 wt%AP was probed primarily.The high energy laser-induced shockwave experiment was utilized to simulate the reaction behavior in hot field.And the larger activated mixture of coated powder could release more energy to promote the growth of shockwave with higher speed up to 518.7±55.9 m/s.In conclusion,85@15 composition is expected to be applied in energetic system as a novel metal fuel.
