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.展开更多
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.展开更多
Improving the energy conversion efficiency in metallic fuel(e.g.,Al)combustion is always desirable but challenging,which often involves redox reactions of aluminum(Al)with various mixed oxidizing environments.For inst...Improving the energy conversion efficiency in metallic fuel(e.g.,Al)combustion is always desirable but challenging,which often involves redox reactions of aluminum(Al)with various mixed oxidizing environments.For instance,Al-O reaction is the most common pathway to release limited energy while Al-F reaction has received much attentions to enhance Al combustion efficiency.However,microscopic understanding of the Al-O/Al-F reaction dynamics remains unsolved,which is fundamentally necessary to further improve Al combustion efficiency.In this work,for the first time,Al-O/Al-F reaction dynamic effects on the combustion of aluminum nanoparticles(n-Al)in oxygen/fluorine containing environments have been revealed via reactive molecular dynamics(RMD)simulations meshing together combustion experiments.Three RMD simulation systems of Al core/O_(2)/HF,n-Al/O_(2)/HF,and n-Al/O_(2)/CF4 with oxygen percentage ranging from 0%to 100%have been performed.The n-Al combustion in mixed O_(2)/CF_4 environments have been conducted by constant volume combustion experiments.RMD results show that Al-O reaction exhibits kinetic benefits while Al-F reaction owns thermodynamic benefits for n-Al combustion.In n-Al/O_(2)/HF,Al-O reaction gives faster energy release rate than Al-F reaction(1.1 times).The optimal energy release efficiency can be achieved with suitable oxygen percentage of 10%and 50%for n-Al/O_(2)/HF and n-Al/O_(2)/CF_4,respectively.In combustion experiments,90%of oxygen percentage can optimally enhance the peak pressure,pressurization rate and combustion heat.Importantly,Al-O reaction prefers to occur on the surface regions while Al-F reaction prefers to proceed in the interior regions of n-Al,confirming the kinetic/thermodynamic benefits of Al-O/Al-F reactions.The synergistic effect of Al-O/Al-F reaction for greatly enhancing n-Al combustion efficiency is demonstrated at atomicscale,which is beneficial for optimizing the combustion performance of metallic fuel.展开更多
As the simplest hydrogen-bonded alcohol,liquid methanol has attracted intensive experimental and theoretical interest.However,theoretical investigations on this system have primarily relied on empirical intermolecular...As the simplest hydrogen-bonded alcohol,liquid methanol has attracted intensive experimental and theoretical interest.However,theoretical investigations on this system have primarily relied on empirical intermolecular force fields or ab initio molecular dynamics with semilocal density functionals.Inspired by recent studies on bulk water using increasingly accurate machine learning force fields,we report a new machine learning force field for liquid methanol with a hybrid functional revPBE0 plus dispersion correction.Molecular dynamics simulations on this machine learning force field are orders of magnitude faster than ab initio molecular dynamics simulations,yielding the radial distribution functions,selfdiffusion coefficients,and hydrogen bond network properties with very small statistical errors.The resulting structural and dynamical properties are compared well with the experimental data,demonstrating the superior accuracy of this machine learning force field.This work represents a successful step toward a first-principles description of this benchmark system and showcases the general applicability of the machine learning force field in studying liquid systems.展开更多
In order to research the feasibility of using the selective adsorption principle to achieve automatic shaping of nano patterns,in this study,using the liquid gallium as the conductive ink and graphene as the printing ...In order to research the feasibility of using the selective adsorption principle to achieve automatic shaping of nano patterns,in this study,using the liquid gallium as the conductive ink and graphene as the printing plate surface,by changing the surface wettability of patterned areas on the nanoscale of graphene printed boards,the automatic formation of liquid gallium patterns on the graphene printed plate surface was simulated.The results indicated that liquid gallium can achieve automatic patterning on the surface of graphene patterned areas;the greater the interaction energy between gallium and carbon atoms,the clearer the pattern;gallium liquid is prone to remain in complex local positions of the pattern,making it difficult to shape the pattern;if the spacing between adjacent pattern lines is too large or too small,it will result in residual gallium liquid between the lines;increasing the thickness of the gallium film will cause the pattern to expand beyond the boundary,but increasing the thickness of the gallium film can also enhance the thickness and uniformity of the pattern lines.In summary,the principle of selective adsorption can be used to achieve the automatic formation of nano patterns,and the pattern formation effect is influenced by factors such as atomic interaction energy and pattern configuration.展开更多
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.展开更多
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 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.展开更多
An OpenMP approach was proposed to parallelize the sequential molecular dynamics(MD) code on shared memory machines. When a code is converted from the sequential form to the parallel form, data dependence is a main pr...An OpenMP approach was proposed to parallelize the sequential molecular dynamics(MD) code on shared memory machines. When a code is converted from the sequential form to the parallel form, data dependence is a main problem. A traditional sequential molecular dynamics code is anatomized to find the data dependence segments in it, and the two different methods, i.e., recover method and backward mapping method were used to eliminate those data dependencies in order to realize the parallelization of this sequential MD code. The performance of the parallelized MD code was analyzed by using some performance analysis tools. The results of the test show that the computing size of this code increases sharply form 1 million atoms before parallelization to 20 million atoms after parallelization, and the wall clock during computing is reduced largely. Some hot-spots in this code are found and optimized by improved algorithm. The efficiency of parallel computing is 30% higher than that of before, and the calculation time is saved and larger scale calculation problems are solved.展开更多
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.展开更多
Surfactant molecules, when dispersed in solution, have been shown to spontaneously form aggregates. Our previous studies on molecular dynamics(MD) calculations have shown that ionic sodium dodecyl sulfate molecules qu...Surfactant molecules, when dispersed in solution, have been shown to spontaneously form aggregates. Our previous studies on molecular dynamics(MD) calculations have shown that ionic sodium dodecyl sulfate molecules quickly aggregated even when the aggregation number is small. The aggregation rate, however, decreased for larger aggregation numbers. In addition, studies have shown that micelle formation was not completed even after a 100 ns-long MD run(Chem. Phys. Lett. 2016, 646, 36). Herein, we analyze the free energy change of micelle formation based on chemical species model combined with molecular dynamics calculations. First, the free energy landscape of the aggregation, ?G_(i+j)^+, where two aggregates with sizes i and j associate to form the(i + j)-mer, was investigated using the free energy of micelle formation of the i-mer, G_i^+, which was obtained through MD calculations. The calculated ?G_(i+j)^+ was negative for all the aggregations where the sum of DS ions in the two aggregates was 60 or less. From the viewpoint of chemical equilibrium, aggregation to the stable micelle is desired. Further, the free energy profile along possible aggregation pathways was investigated, starting from small aggregates and ending with the complete thermodynamically stable micelles in solution. The free energy profiles, G(l, k), of the aggregates at l-th aggregation path and k-th state were evaluated by the formation free energy ∑_in_i( l,k)G_i^+ and the free energy of mixing ∑_in_i( l,k)k_BTln( n_i( l,k)/n( l,k)), where ni(l, k) is the number of i-mer in the system at the l-th i aggregation path and k-th state, with n(l,k)= ∑_n_i( l,k). All the aggregation pathways were obtained from the initial i state of 12 pentamers to the stable micelle with i = 60. All the calculated G(l, k) values monotonically decreased with increasing k. This indicates that there are no free energy barriers along the pathways. Hence, the slowdown is not due to the thermodynamic stability of the aggregates, but rather the kinetics that inhibit the association of the fragments. The time required for a collision between aggregates, one of the kinetic factors, was evaluated using the fast passage time, t_(FPT). The calculated t_(FPT) was about 20 ns for the aggregates with N = 31. Therefore, if aggregation is a diffusion-controlled process, it should be completed within the 100 ns-simulation. However, aggregation does not occur due to the free energy barrier between the aggregates, that is, the repulsive force acting on them. This may be caused by electrostatic repulsions produced by the overlap of the electric double layers, which are formed by the negative charge of the hydrophilic groups and counter sodium ions on the surface of the aggregates.展开更多
The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this wo...The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this work.The following properties were determined:density,self-diffusion coefficient,excess molar volume,and radial distribution function.The results show that with an increase in the mole fraction of[HPy][BF_(4)],the self-diffusion coefficient decreases.Additionally,the excess molar volume initially decreases,reaches a minimum,and then increases.The rules of radial distribution functions(RDFs)of characteristic atoms are different.With increasing the mole fraction of[HPy][BF_(4)],the first peak of the RDFs of HA1-F decreases,while that of CT6-CT6 rises at first and then decreases.This indicates that the solvent molecules affect the polar and non-polar regions of[HPy][BF_(4)]differently.展开更多
This work uses the molecular dynamics approach to study the effects of functionalization of carbon nanotubes(CNTs)on the mechanical properties of Cu64Zr36 metallic glass(MG).Three types of functional groups,carboxylic...This work uses the molecular dynamics approach to study the effects of functionalization of carbon nanotubes(CNTs)on the mechanical properties of Cu64Zr36 metallic glass(MG).Three types of functional groups,carboxylic,vinyl and ester were used.The effect of CNT volume fraction(Vf)and the number of functional groups attached to CNT,on the mechanical properties and thermal conductivity of CNT-MG composites was analysed using Biovia Materials Studio.At lower values of Vf(from 0 to 5%),the percentage increase in Young’s modulus was approximately 66%.As the value of Vf was increased further(from 5 to 12%),the rate of increase in Young’s modulus was reduced to 16%.The thermal conductivity was found to increase from 1.52 W/mK at Vf?0%to 5.88 W/mK at Vf?12%,thus giving an increase of approximately 286%.Functionalization of SWCNT reduced the thermal conductivity of the SWCNT-MG composites.展开更多
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.展开更多
Based on the crystal engineering, six models of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine(HMX)/3-nitro-1,2,4-triazol-5-one(NTO) supramolecular explosive were designed. The probable formation of HMX/NTO supramol...Based on the crystal engineering, six models of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine(HMX)/3-nitro-1,2,4-triazol-5-one(NTO) supramolecular explosive were designed. The probable formation of HMX/NTO supramolecular explosive was investigated by the molecular dynamic (MD) method. Interaction between oxygen atoms in HMX and hydrogen atoms in NTO or between hydrogen atoms in HMX and oxygen atoms in NTO were studied by the radial distribution function (RDF). It shows that there are strong hydrogen bonds and Van Der Waals forces between HMX and NTO, in which the hydrogen bonds between oxygen atoms in the NTO and hydrogen atoms in HMX are the main host-guest interactions. The distributions of bond length, bond angle and dihedral angle were simulated by MD. It shows that the structure of HMX is seriously distorted. The binding energies and X-ray powder diffraction (XRD) patterns were calculated on the basis of the final HMX/NTO supramolecular structures. The results show that the binding energies of six supramolecular models are E binding (1 1 1-) >E binding (1 0 0)>E binding (0 2 0)>E binding (random)>E binding (1 0 2-)>Ebinding (0 1 1), and the XRD patterns of six supramolecular models are quite different from pure HMX or NTO. Based on the investigation for growth morphology, binding energies and RDF, the model of HMX supercell substituted by NTO along the (1 1 1-) surface of HMX is easier to form.展开更多
The feasibility study of the AlCl(g) generated by Al_2O-AlCl_2-C system under vacuum was carried out by thermodynamic analysis and CASTEP package of the Material Studio program which was based on density functional th...The feasibility study of the AlCl(g) generated by Al_2O-AlCl_2-C system under vacuum was carried out by thermodynamic analysis and CASTEP package of the Material Studio program which was based on density functional theory(DFT) formalism. Thermodynamic calculations indicate that Al Cl and CO molecules can be formed under conditions of temperature 1760 K and the pressure of 60 Pa. The interaction of Al_2O and AlCl_2 with C shows that the chemical adsorption of Al_2O and AlCl_2 does take place on C(001) crystal plane, and at the same time, new chemical bond is formed between Al atom in Al_2O and Cl atoms from one of the Al—Cl bonds in AlCl_2. The results, after 1.25 ps dynamics simulation, indicate that adsorbed Al Cl molecules are generated and CO molecule will be formed in this system, and they will escape from C(001) surface after a longer period of dynamic simulation time. It means that the reaction of Al_2O and AlCl_2 with C can be carried out under given constraint condition.展开更多
文摘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 (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.
基金support by the National Natural Science Foundation of China(NSFC,Grant Nos.12002324,12372341,12172342)。
文摘Improving the energy conversion efficiency in metallic fuel(e.g.,Al)combustion is always desirable but challenging,which often involves redox reactions of aluminum(Al)with various mixed oxidizing environments.For instance,Al-O reaction is the most common pathway to release limited energy while Al-F reaction has received much attentions to enhance Al combustion efficiency.However,microscopic understanding of the Al-O/Al-F reaction dynamics remains unsolved,which is fundamentally necessary to further improve Al combustion efficiency.In this work,for the first time,Al-O/Al-F reaction dynamic effects on the combustion of aluminum nanoparticles(n-Al)in oxygen/fluorine containing environments have been revealed via reactive molecular dynamics(RMD)simulations meshing together combustion experiments.Three RMD simulation systems of Al core/O_(2)/HF,n-Al/O_(2)/HF,and n-Al/O_(2)/CF4 with oxygen percentage ranging from 0%to 100%have been performed.The n-Al combustion in mixed O_(2)/CF_4 environments have been conducted by constant volume combustion experiments.RMD results show that Al-O reaction exhibits kinetic benefits while Al-F reaction owns thermodynamic benefits for n-Al combustion.In n-Al/O_(2)/HF,Al-O reaction gives faster energy release rate than Al-F reaction(1.1 times).The optimal energy release efficiency can be achieved with suitable oxygen percentage of 10%and 50%for n-Al/O_(2)/HF and n-Al/O_(2)/CF_4,respectively.In combustion experiments,90%of oxygen percentage can optimally enhance the peak pressure,pressurization rate and combustion heat.Importantly,Al-O reaction prefers to occur on the surface regions while Al-F reaction prefers to proceed in the interior regions of n-Al,confirming the kinetic/thermodynamic benefits of Al-O/Al-F reactions.The synergistic effect of Al-O/Al-F reaction for greatly enhancing n-Al combustion efficiency is demonstrated at atomicscale,which is beneficial for optimizing the combustion performance of metallic fuel.
基金supported by the CAS Project for Young Scientists in Basic Research(YSBR-005)the National Natural Science Foundation of China(22325304,22221003 and 22033007)We acknowledge the Supercomputing Center of USTC,Hefei Advanced Computing Center,Beijing PARATERA Tech Co.,Ltd.,for providing high-performance computing services。
文摘As the simplest hydrogen-bonded alcohol,liquid methanol has attracted intensive experimental and theoretical interest.However,theoretical investigations on this system have primarily relied on empirical intermolecular force fields or ab initio molecular dynamics with semilocal density functionals.Inspired by recent studies on bulk water using increasingly accurate machine learning force fields,we report a new machine learning force field for liquid methanol with a hybrid functional revPBE0 plus dispersion correction.Molecular dynamics simulations on this machine learning force field are orders of magnitude faster than ab initio molecular dynamics simulations,yielding the radial distribution functions,selfdiffusion coefficients,and hydrogen bond network properties with very small statistical errors.The resulting structural and dynamical properties are compared well with the experimental data,demonstrating the superior accuracy of this machine learning force field.This work represents a successful step toward a first-principles description of this benchmark system and showcases the general applicability of the machine learning force field in studying liquid systems.
文摘In order to research the feasibility of using the selective adsorption principle to achieve automatic shaping of nano patterns,in this study,using the liquid gallium as the conductive ink and graphene as the printing plate surface,by changing the surface wettability of patterned areas on the nanoscale of graphene printed boards,the automatic formation of liquid gallium patterns on the graphene printed plate surface was simulated.The results indicated that liquid gallium can achieve automatic patterning on the surface of graphene patterned areas;the greater the interaction energy between gallium and carbon atoms,the clearer the pattern;gallium liquid is prone to remain in complex local positions of the pattern,making it difficult to shape the pattern;if the spacing between adjacent pattern lines is too large or too small,it will result in residual gallium liquid between the lines;increasing the thickness of the gallium film will cause the pattern to expand beyond the boundary,but increasing the thickness of the gallium film can also enhance the thickness and uniformity of the pattern lines.In summary,the principle of selective adsorption can be used to achieve the automatic formation of nano patterns,and the pattern formation effect is influenced by factors such as atomic interaction energy and pattern configuration.
基金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.
基金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(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.
基金Project (50371026) supported by the National Natural Science Foundation of China
文摘An OpenMP approach was proposed to parallelize the sequential molecular dynamics(MD) code on shared memory machines. When a code is converted from the sequential form to the parallel form, data dependence is a main problem. A traditional sequential molecular dynamics code is anatomized to find the data dependence segments in it, and the two different methods, i.e., recover method and backward mapping method were used to eliminate those data dependencies in order to realize the parallelization of this sequential MD code. The performance of the parallelized MD code was analyzed by using some performance analysis tools. The results of the test show that the computing size of this code increases sharply form 1 million atoms before parallelization to 20 million atoms after parallelization, and the wall clock during computing is reduced largely. Some hot-spots in this code are found and optimized by improved algorithm. The efficiency of parallel computing is 30% higher than that of before, and the calculation time is saved and larger scale calculation problems are solved.
基金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.
基金This work was supported by FLAGSHIP2020,MEXT within Priority Study 5(Development of New Fundamental Technologies for High-Efficiency Energy Creation,Conversion/Storage and Use)Using Computational Resources of the K Computer Provided by the RIKEN Advanced
文摘Surfactant molecules, when dispersed in solution, have been shown to spontaneously form aggregates. Our previous studies on molecular dynamics(MD) calculations have shown that ionic sodium dodecyl sulfate molecules quickly aggregated even when the aggregation number is small. The aggregation rate, however, decreased for larger aggregation numbers. In addition, studies have shown that micelle formation was not completed even after a 100 ns-long MD run(Chem. Phys. Lett. 2016, 646, 36). Herein, we analyze the free energy change of micelle formation based on chemical species model combined with molecular dynamics calculations. First, the free energy landscape of the aggregation, ?G_(i+j)^+, where two aggregates with sizes i and j associate to form the(i + j)-mer, was investigated using the free energy of micelle formation of the i-mer, G_i^+, which was obtained through MD calculations. The calculated ?G_(i+j)^+ was negative for all the aggregations where the sum of DS ions in the two aggregates was 60 or less. From the viewpoint of chemical equilibrium, aggregation to the stable micelle is desired. Further, the free energy profile along possible aggregation pathways was investigated, starting from small aggregates and ending with the complete thermodynamically stable micelles in solution. The free energy profiles, G(l, k), of the aggregates at l-th aggregation path and k-th state were evaluated by the formation free energy ∑_in_i( l,k)G_i^+ and the free energy of mixing ∑_in_i( l,k)k_BTln( n_i( l,k)/n( l,k)), where ni(l, k) is the number of i-mer in the system at the l-th i aggregation path and k-th state, with n(l,k)= ∑_n_i( l,k). All the aggregation pathways were obtained from the initial i state of 12 pentamers to the stable micelle with i = 60. All the calculated G(l, k) values monotonically decreased with increasing k. This indicates that there are no free energy barriers along the pathways. Hence, the slowdown is not due to the thermodynamic stability of the aggregates, but rather the kinetics that inhibit the association of the fragments. The time required for a collision between aggregates, one of the kinetic factors, was evaluated using the fast passage time, t_(FPT). The calculated t_(FPT) was about 20 ns for the aggregates with N = 31. Therefore, if aggregation is a diffusion-controlled process, it should be completed within the 100 ns-simulation. However, aggregation does not occur due to the free energy barrier between the aggregates, that is, the repulsive force acting on them. This may be caused by electrostatic repulsions produced by the overlap of the electric double layers, which are formed by the negative charge of the hydrophilic groups and counter sodium ions on the surface of the aggregates.
文摘The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this work.The following properties were determined:density,self-diffusion coefficient,excess molar volume,and radial distribution function.The results show that with an increase in the mole fraction of[HPy][BF_(4)],the self-diffusion coefficient decreases.Additionally,the excess molar volume initially decreases,reaches a minimum,and then increases.The rules of radial distribution functions(RDFs)of characteristic atoms are different.With increasing the mole fraction of[HPy][BF_(4)],the first peak of the RDFs of HA1-F decreases,while that of CT6-CT6 rises at first and then decreases.This indicates that the solvent molecules affect the polar and non-polar regions of[HPy][BF_(4)]differently.
文摘This work uses the molecular dynamics approach to study the effects of functionalization of carbon nanotubes(CNTs)on the mechanical properties of Cu64Zr36 metallic glass(MG).Three types of functional groups,carboxylic,vinyl and ester were used.The effect of CNT volume fraction(Vf)and the number of functional groups attached to CNT,on the mechanical properties and thermal conductivity of CNT-MG composites was analysed using Biovia Materials Studio.At lower values of Vf(from 0 to 5%),the percentage increase in Young’s modulus was approximately 66%.As the value of Vf was increased further(from 5 to 12%),the rate of increase in Young’s modulus was reduced to 16%.The thermal conductivity was found to increase from 1.52 W/mK at Vf?0%to 5.88 W/mK at Vf?12%,thus giving an increase of approximately 286%.Functionalization of SWCNT reduced the thermal conductivity of the SWCNT-MG composites.
基金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.
基金Sponsored by the National Natural Science Foundation of China (61106078)NUST Research Funding (2011ZDJH28)
文摘Based on the crystal engineering, six models of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine(HMX)/3-nitro-1,2,4-triazol-5-one(NTO) supramolecular explosive were designed. The probable formation of HMX/NTO supramolecular explosive was investigated by the molecular dynamic (MD) method. Interaction between oxygen atoms in HMX and hydrogen atoms in NTO or between hydrogen atoms in HMX and oxygen atoms in NTO were studied by the radial distribution function (RDF). It shows that there are strong hydrogen bonds and Van Der Waals forces between HMX and NTO, in which the hydrogen bonds between oxygen atoms in the NTO and hydrogen atoms in HMX are the main host-guest interactions. The distributions of bond length, bond angle and dihedral angle were simulated by MD. It shows that the structure of HMX is seriously distorted. The binding energies and X-ray powder diffraction (XRD) patterns were calculated on the basis of the final HMX/NTO supramolecular structures. The results show that the binding energies of six supramolecular models are E binding (1 1 1-) >E binding (1 0 0)>E binding (0 2 0)>E binding (random)>E binding (1 0 2-)>Ebinding (0 1 1), and the XRD patterns of six supramolecular models are quite different from pure HMX or NTO. Based on the investigation for growth morphology, binding energies and RDF, the model of HMX supercell substituted by NTO along the (1 1 1-) surface of HMX is easier to form.
基金Projects(51104078,51264023)supported by the National Natural Science Foundation of ChinaProject(2010CD022)supported by Yunnan Province Applied Basic Research Fund,China+2 种基金Project(IRT1250)supported by the Program for Innovative Research Team in University of Ministry of Education of ChinaProject(U1202271)supported by the National Natural Science Foundation of China-Yunnan United FundProject(KKZ3201252020)supported by Kunming University of Science and Technology for Talent Training,China
文摘The feasibility study of the AlCl(g) generated by Al_2O-AlCl_2-C system under vacuum was carried out by thermodynamic analysis and CASTEP package of the Material Studio program which was based on density functional theory(DFT) formalism. Thermodynamic calculations indicate that Al Cl and CO molecules can be formed under conditions of temperature 1760 K and the pressure of 60 Pa. The interaction of Al_2O and AlCl_2 with C shows that the chemical adsorption of Al_2O and AlCl_2 does take place on C(001) crystal plane, and at the same time, new chemical bond is formed between Al atom in Al_2O and Cl atoms from one of the Al—Cl bonds in AlCl_2. The results, after 1.25 ps dynamics simulation, indicate that adsorbed Al Cl molecules are generated and CO molecule will be formed in this system, and they will escape from C(001) surface after a longer period of dynamic simulation time. It means that the reaction of Al_2O and AlCl_2 with C can be carried out under given constraint condition.
