INS （Inelastic Neutron Scattering） spectrum of methane hydrate was measured on MARI （a direct-geometry chopper spectrometer） at Ruther-ford Appleton Laboratory, UK. Compared with ice Ih, it is found that the whole...INS （Inelastic Neutron Scattering） spectrum of methane hydrate was measured on MARI （a direct-geometry chopper spectrometer） at Ruther-ford Appleton Laboratory, UK. Compared with ice Ih, it is found that the whole spectrum of methane hydrate moves toward high-energy by about 1.5 meV. Using lattice dynamical （LD） technique, computer simulations of methane hydrate were carried out. In the simulations, four potential models （BE TIP3P, TIP4P, MCY） were employed to calculate the phonon density of states （PDOS）. Comparing the calculated PDOS spectrum with the experimental spectrum, it is found that BF, TIP4P, and TIP3P potential lattices give out well-separated translational and librational bands while MCY potential lattice is unstable to do so and this model is not suitable to describe hydrate system.展开更多
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.展开更多
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.展开更多
In order to ensure the penetrability of double-cased perforation in offshore oil and gas fields and to maximize the capacity of perforation completion, This study establishes a dynamic model of double-cased perforatio...In order to ensure the penetrability of double-cased perforation in offshore oil and gas fields and to maximize the capacity of perforation completion, This study establishes a dynamic model of double-cased perforation using ANSYS/LS-DYNA simulation technology. The combination of critical perforation parameters for double casing is obtained by studying the influencing factors of the jet-forming process,perforation depth, diameter, and stress changes of the inner and outer casing. The single-target perforation experiments under high-temperature and high-pressure(HTHP) conditions and ground full-scale ring target perforation tests are designed to verify the accuracy of numerical simulation results. The reduced factor is adopted as the quantitative measure of perforation depth and diameter for different types of perforation charge under different conditions. The results show that the perforation depth reduction increases with temperature and pressure, and the reduced factor is between 0.67 and 0.87 under HTHP conditions of 130℃/44 MPa and 137℃/60 MPa. Comparing the results of the numerical simulation and the full-scale test correction, the maximum error is less than 8.91%, and this numerical simulation has strong reliability. This research provides a basis for a reasonable range of double-cased perforation parameters and their optimal selection.展开更多
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.展开更多
Recently,there has been a growing prevalence in the utilization of graphdiyne(GDY)in the field of biomedicine,attributed to its distinctive physical structure and chemical properties.Additionally,its biocompatibility ...Recently,there has been a growing prevalence in the utilization of graphdiyne(GDY)in the field of biomedicine,attributed to its distinctive physical structure and chemical properties.Additionally,its biocompatibility has garnered increasing attention.However,there is a lack of research on the biological effects and physical mechanisms of GDYprotein interactions at the molecular scale.In this study,the villin headpiece subdomain(HP35)served as a representative protein model.Molecular dynamics simulations were employed to investigate the interaction process between the HP35 protein and GDY,as well as the structural evolution of the protein.The data presented in our study demonstrate that GDY can rapidly adsorb HP35 protein and induce denaturation to one of the a-helix structures of HP35 protein.This implies a potential cytotoxicity concern of GDY for biological systems.Compared to graphene,GDY induced less disruption to HP35 protein.This can be attributed to the presence of natural triangular vacancies in GDY,which prevents p–p stacking action and the limited interaction of GDY with HP35 protein is not conducive to the expansion of protein structures.These findings unveil the biological effects of GDY at the molecular level and provide valuable insights for the application of GDY in biomedicine.展开更多
Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear...Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear-thinning fluid in a microchannel.We validated the feasibility of our simulation method by evaluating the mean square displacement and Reynolds number of the solution layers.The results show that the change rule of the fluid system's velocity profile and interaction energy can reflect the shear-thinning characteristics of the fluids.The velocity profile resembles a top-hat shape,intensifying as the fluid's power law index decreases.The interaction energy between the wall and the fluid decreases gradually with increasing velocity,and a high concentration of non-Newtonian fluid reaches a plateau sooner.Moreover,the velocity profile of the fluid is related to the molecule number density distribution and their values are inversely proportional.By analyzing the radial distribution function,we found that the hydrogen bonds between solute and water molecules weaken with the increase in velocity.This observation offers an explanation for the shear-thinning phenomenon of the non-Newtonian flow from a micro perspective.展开更多
Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced...Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced critical load and hardness compared to nc-Ni,where perfect,stair-rod and Shockley dislocations are activated at (111),(111) and (111) slip planes in nt-Ni compared to only SSockley dislocation nucleation at (111) and (111) slip planes of nc-Ni.The nt-Ni exhibits a less significant indentation size effect in comparison with nc-Ni due to the dislocation slips hindrance of the twin boundary.The atomic deformation associated with the indentation size effect is investigated during dislocation transmission.Different from the decreasing partial slips parallel to the indenter surface in nc-Ni with increasing temperature,the temperaturedependent atomic deformation of nt-Ni is closely related to the twin boundary:from the partial slips parallel to the twin boundary (~10 K),to increased confined layer slips and decreased twin migration(300 K–600 K),to decreased confined layer slips and increased dislocation interaction of dislocation pinning and dissociation (900 K–1200 K).Dislocation density and atomic structure types through quantitative analysis are implemented to further reveal the above-mentioned dislocation motion and atomic structure alteration.Our study is helpful for understanding the temperature-dependent plasticity of twin boundary in nanotwinned materials.展开更多
According to the complex differential accumulation history of deep marine oil and gas in superimposed basins,the Lower Paleozoic petroleum system in Tahe Oilfield of Tarim Basin is selected as a typical case,and the p...According to the complex differential accumulation history of deep marine oil and gas in superimposed basins,the Lower Paleozoic petroleum system in Tahe Oilfield of Tarim Basin is selected as a typical case,and the process of hydrocarbon generation and expulsion,migration and accumulation,adjustment and transformation of deep oil and gas is restored by means of reservoine-forming dynamics simulation.The thermal evolution history of the Lower Cambrian source rocks in Tahe Oilfield reflects the obvious differences in hydrocarbon generation and expulsion process and intensity in different tectonic zones,which is the main reason controlling the differences in deep oil and gas phases.The complex transport system composed of strike-slip fault and unconformity,etc.controlled early migration and accumulation and late adjustment of deep oil and gas,while the Middle Cambrian gypsum-salt rock in inner carbonate platform prevented vertical migration and accumulation of deep oil and gas,resulting in an obvious"fault-controlled"feature of deep oil and gas,in which the low potential area superimposed by the NE-strike-slip fault zone and deep oil and gas migration was conducive to accumulation,and it is mainly beaded along the strike-slip fault zone in the northeast direction.The dynamic simulation of reservoir formation reveals that the spatio-temporal configuration of"source-fault-fracture-gypsum-preservation"controls the differential accumulation of deep oil and gas in Tahe Oilfield.The Ordovician has experienced the accumulation history of multiple periods of charging,vertical migration and accumulation,and lateral adjustment and transformation,and deep oil and gas have always been in the dynamic equilibrium of migration,accumulation and escape.The statistics of residual oil and gas show that the deep stratum of Tahe Oilfield still has exploration and development potential in the Ordovician Yingshan Formation and Penglaiba Formation,and the Middle and Upper Cambrian ultra-deep stratum has a certain oil and gas resource prospect.This study provides a reference for the dynamic quantitative evaluation of deep oil and gas in the Tarim Basin,and also provides a reference for the study of reservoir formation and evolution in carbonate reservoir of paleo-craton basin.展开更多
The Brownian dynamics （BD） simulation of a dilute surfactant solution is conducted in a steady shear flow. The rodlike micelle is assumed as a rigid rod composed of lined-up beads. A novel intercluster potential mod...The Brownian dynamics （BD） simulation of a dilute surfactant solution is conducted in a steady shear flow. The rodlike micelle is assumed as a rigid rod composed of lined-up beads. A novel intercluster potential model is introduced for describing the interactions between, micelles. In the model, the Lennard-Jones and the soft-sphere potentials are used as inter-bead potentials for end-end and interior-interior beads, respectively. The micelles are combined at their ends to form a network structure at lower shear rates and are disconnected to become more and more parallel to the shear flow direction with increasing shear rate. The change of micellar microstructures with the variation of the shear rate results in shear thinning characteristics of the computed shear viscosities and first normal stress difference coefficients. The effects of surfactant solution concentration on the micellar structures and rheological properties are also investigated. Results show that the shear viscosities and the first normal stress difference coefficients increase with increasing the viscosity of the surfactant solution.展开更多
The track model used in the dynamic analysis and design system software is investigated. A home made tank is taken as an example to illustrate the method for modeling an integral tracked vehicle and perform the dynam...The track model used in the dynamic analysis and design system software is investigated. A home made tank is taken as an example to illustrate the method for modeling an integral tracked vehicle and perform the dynamic simulation. The obtained results have demonstrated that the simulation method has the advantage of high efficiency, more convenience and more insight into the dynamical behavior of the system.展开更多
A molecular dynamics simulation study has been performed for the microstructure evolution in a liquid metal Ni system during crystallization process at two cooling rates by adopting the embedded atom method （EAM） mo...A molecular dynamics simulation study has been performed for the microstructure evolution in a liquid metal Ni system during crystallization process at two cooling rates by adopting the embedded atom method （EAM） model potential. The bond-type index method of Honeycutt-Andersen （HA） and a new cluster-type index method （CTIM-2） have been used to detect and analyse the microstructures in this system. It is demonstrated that the cooling rate plays a critical role in the microstructure evolution： below the crystallization temperature Tc, the effects of cooling rate are very remarkable and can be fully displayed. At different cooling rates of 2.0 × 10^13 K·s^-1 and 1.0 × 10^12 K·s^-1, two different kinds of crystal structures are obtained in the system. The first one is the coexistence of the hcp （expressed by （12 0 0 0 6 6） in CTIM-2） and the fcc （12 0 0 0 12 0） basic clusters consisting of 1421 and 1422 bond-types, and the hcp basic cluster becomes the dominant one with decreasing temperature, the second one is mainly the fcc （12 0 0 0 12 0） basic clusters consisting of 1421 bond-type, and their crystallization temperatures Tc would be 1073 and 1173 K, respectively.展开更多
In this paper molecular dynamics simulations are performed to study the accumulation behaviour of N2 and H2 at water/graphite interface under ambient temperature and pressure. It finds that both N2 and H2 molecules ca...In this paper molecular dynamics simulations are performed to study the accumulation behaviour of N2 and H2 at water/graphite interface under ambient temperature and pressure. It finds that both N2 and H2 molecules can accumulate at the interface and form one of two states according to the ratio of gas molecules number to square of graphite surface from our simulation results： gas films （pancake-like） for a larger ratio and nanobubbles for a smaller ratio. In addition, we discuss the stabilities of nanobubbles at different environment temperatures. Surprisingly, it is found that the density of both kinds of gas states can be greatly increased, even comparable with that of the liquid N2 and liquid H2. The present results are expected to be helpful for the understanding of the stable existence of gas film （pancake-like） and nanobubbles.展开更多
Selectivity control is a difficult scientific and industrial challenge in methanol-to-olefins(MTO)conversion.It has been experimentally established that the topology of zeolite catalysts influenced the distribution of...Selectivity control is a difficult scientific and industrial challenge in methanol-to-olefins(MTO)conversion.It has been experimentally established that the topology of zeolite catalysts influenced the distribution of products.Besides the topology effect on reaction kinetics,the topology influences the diffusion of reactants and products in catalysts as well.In this work,by using COMPASS force-field molecular dynamics method,we investigated the intracrystalline diffusion of ethene and propene in four different zeolites,CHA,MFI,BEA and FAU,at different temperatures.The self-diffusion coefficients and diffusion activation barriers were calculated.A strong restriction on the diffusion of propene in CHA was observed because the self-diffusion coefficient ratio of ethene to propene is larger than 18 and the diffusion activation barrier of propene is more than 20 kJ/mol in CHA.This ratio decreases with the increase of temperature in the four investigated zeolites.The shape selectivity on products from diffusion perspective can provide some implications on the understanding of the selectivity difference between HSAPO-34 and HZSM-5 catalysts for the MTO conversion.展开更多
This paper performs molecular dynamics simulations to investigate the role of the monovalent cations K, Na and the divalent cation Ca on the stability and swelling of montmorillonite. The recently developed CLAYFF for...This paper performs molecular dynamics simulations to investigate the role of the monovalent cations K, Na and the divalent cation Ca on the stability and swelling of montmorillonite. The recently developed CLAYFF force field is used to predict the basal spacing as a function of the water content in the interlayer. The simulations reproduced the swelling pattern of these montmorillonites, suggesting a mechanism of their hydration different （K＋ 〈 Na＋ 〈 Ca2＋） from that of K＋-, Na＋-, and Ca2＋-montmorillonites. In particular, these results indicate that the valence of the cations has the larger impact on the behaviour of clay water systems. It also finds that the differences in size and hydration energy of K＋, Na＋ and Ca2＋ ions have strong implications for the structure of interlayer. This leads to the differences in the layer spacings of the simulated K＋-, Na＋-, and Ca2＋-montmorillonites. Furthermore, these simulations show that the K cations interact strongly with the clay sheets for the dehydrated clay sheets, but for the hydrated clays the Ca cations interact clearly strongly with the clay sheets.展开更多
Molecular dynamics simulations are performed to study the growth mechanism of CH4-CO2 mixed hydrate in xco2 = 75%, xco2 = 50%, and zco2 = 25% systems at T = 250 K, 255 K and 260 K, respectively. Our simulation results...Molecular dynamics simulations are performed to study the growth mechanism of CH4-CO2 mixed hydrate in xco2 = 75%, xco2 = 50%, and zco2 = 25% systems at T = 250 K, 255 K and 260 K, respectively. Our simulation results show that the growth rate of CH4-CO2 mixed hydrate increases as the CO2 concentration in the initial solution phase increases and the temperature decreases. Via hydrate formation, the composition of CO2 in hydrate phase is higher than that in initial solution phase and the encaging capacity of CO2 in hydrates increases with the decrease in temperature. By analysis of the cage occupancy ratio of CH4 molecules and CO2 molecules in large cages to small cages, we find that CO2 molecules are preferably encaged into the large cages of the hydrate crystal as compared with CH4 molecules. Interestingly, CH4 molecules and CO2 molecules frequently replace with each other in some particular cage sites adjacent to hydrate/solution interface during the crystal growth process. These two species of guest molecules eventually act to stabilize the newly formed hydrates, with CO2 molecules occupying large cages and CH4 molecules occupying small cages in hydrate.展开更多
We use non-equilibrium molecular dynamics simulations to calculate the self-diffusion coefficient, D, of a Lennard Jones fluid over a wide density and temperature range. The change in self-diffusion coefficient with t...We use non-equilibrium molecular dynamics simulations to calculate the self-diffusion coefficient, D, of a Lennard Jones fluid over a wide density and temperature range. The change in self-diffusion coefficient with temperature decreases by increasing density. For density ρ＊ = ρσ3 = 0.84 we observe a peak at the value of the self-diffusion coefficient and the critical temperature T＊ = kT/ε = 1.25. The value of the self-diffusion coefficient strongly depends on system size. The data of the self-diffusion coefficient are fitted to a simple analytic relation based on hydrodynamic arguments. This correction scales as N-α, where α is an adjustable parameter and N is the number of particles. It is observed that the values of a 〈 1 provide quite a good correction to the simulation data. The system size dependence is very strong for lower densities, but it is not as strong for higher densities. The self-diffusion coefficient calculated with non-equilibrium molecular dynamic simulations at different temperatures and densities is in good agreement with other calculations fronl the literature.展开更多
The understanding of the structural requirements and the intermolecular-interaction mechanism are important for discovering potent angiotensin-converting enzyme(ACE)inhibitory peptides.In this study,we modifi ed an eg...The understanding of the structural requirements and the intermolecular-interaction mechanism are important for discovering potent angiotensin-converting enzyme(ACE)inhibitory peptides.In this study,we modifi ed an egg-white derived peptide,LAPYK,using the amino acids with different properties to produce the LAPYK-modified peptides.The ACE inhibitory activities of the modified peptides were determined to explore the structural requirements of ACE inhibitory peptides(ACEIPs).Molecular simulation and isothermal titration calorimetry analysis were used to investigate interactions between the peptides and ACE.We found that hydrophobicity and the amino acids with ring structures were benefi cial for the ACE inhibitory activities of the peptides.The results of the molecular mechanics poisson boltzmann surface area(MMPBSA)binding free energy calculations indicated that the polar solvation free energy(ΔG_(polar))of the charged peptides(LAPYK,LAPYE)were unfavorable for binding to ACE.On the other hand,the results of isothermal titration calorimetry analyses suggested that the enthalpy-driven ACE-peptide interactions were more favorable than the entropy-driven ACE-peptide interaction counterparts.展开更多
Shale oil formations contain both inorganic and organic media.The organic matter holds both free oil in the pores and dissolved oil within the kerogen molecules.The free oil flow in organic pores and the dissolved oil...Shale oil formations contain both inorganic and organic media.The organic matter holds both free oil in the pores and dissolved oil within the kerogen molecules.The free oil flow in organic pores and the dissolved oil diffusion in kerogen molecules are coupled together.The molecular flow of free n-alkanes is an important process of shale oil accumulation and production.To study the dynamics of imbibition process of n-alkane molecules into kerogen slits,molecular dynamics(MD)simulations are conducted.Effects of slit width,temperature,and n-alkane types on the penetration speed,dynamic contact angle,and molecular conformations were analyzed.Results showed that molecular transportation of n-alkanes is dominated by molecular structure and molecular motion at this scale.The space-confinement conformational changes of molecules slow down the filling speeds in the narrow slits.The n-alkane molecules with long carbon chains require more time to undergo conformational changes.The high content of short-chain alkanes and high temperature facilitate the flow of alkane mixtures in kerogen slits.Results obtained from this study are useful for understanding the underlying nanoscale flow mechanism in shale formations.展开更多
Currently,wire bonding is the most popular first-level interconnection technology used between the die and package terminals,but even with its long-term and excessive usage,the mechanism of wire bonding has not been c...Currently,wire bonding is the most popular first-level interconnection technology used between the die and package terminals,but even with its long-term and excessive usage,the mechanism of wire bonding has not been completely evaluated.Therefore,fundamental research is still needed.In this study,the mechanism of microweld formation and breakage during Cu-Cu wire bonding was investigated by using molecular dynamics simulation.The contact model for the nanoindentation process between the wire and substrate was developed to simulate the contact process of the Cu wire and Cu substrate.Elastic contact and plastic instability were investigated through the loading and unloading processes.Moreover,the evolution of the indentation morphology and distributions of the atomic stress were also investigated.It was shown that the loading and unloading curves do not coincide,and the unloading curve exhibited hysteresis.For the substrate,in the loading process,the main force changed from attractive to repulsive.The maximum von Mises stress increased and shifted from the center toward the edge of the contact area.During the unloading process,the main force changed from repulsive to attractive.The Mises stress reduced first and then increased.Stress concentration occurs around dislocations in the middle area of the Cu wire.展开更多
基金supported by the National Natural Science Foundation of China (Grant No. 10474085)
文摘INS （Inelastic Neutron Scattering） spectrum of methane hydrate was measured on MARI （a direct-geometry chopper spectrometer） at Ruther-ford Appleton Laboratory, UK. Compared with ice Ih, it is found that the whole spectrum of methane hydrate moves toward high-energy by about 1.5 meV. Using lattice dynamical （LD） technique, computer simulations of methane hydrate were carried out. In the simulations, four potential models （BE TIP3P, TIP4P, MCY） were employed to calculate the phonon density of states （PDOS）. Comparing the calculated PDOS spectrum with the experimental spectrum, it is found that BF, TIP4P, and TIP3P potential lattices give out well-separated translational and librational bands while MCY potential lattice is unstable to do so and this model is not suitable to describe hydrate system.
文摘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.
文摘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.
基金the support of the Foundation of Natural Science Foundation of Shaanxi Province, Grant/ Award nos. 2023-JC-YB-361National Natural Science Foundation (Number 52104033)。
文摘In order to ensure the penetrability of double-cased perforation in offshore oil and gas fields and to maximize the capacity of perforation completion, This study establishes a dynamic model of double-cased perforation using ANSYS/LS-DYNA simulation technology. The combination of critical perforation parameters for double casing is obtained by studying the influencing factors of the jet-forming process,perforation depth, diameter, and stress changes of the inner and outer casing. The single-target perforation experiments under high-temperature and high-pressure(HTHP) conditions and ground full-scale ring target perforation tests are designed to verify the accuracy of numerical simulation results. The reduced factor is adopted as the quantitative measure of perforation depth and diameter for different types of perforation charge under different conditions. The results show that the perforation depth reduction increases with temperature and pressure, and the reduced factor is between 0.67 and 0.87 under HTHP conditions of 130℃/44 MPa and 137℃/60 MPa. Comparing the results of the numerical simulation and the full-scale test correction, the maximum error is less than 8.91%, and this numerical simulation has strong reliability. This research provides a basis for a reasonable range of double-cased perforation parameters and their optimal selection.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant No.52072132).
文摘Recently,there has been a growing prevalence in the utilization of graphdiyne(GDY)in the field of biomedicine,attributed to its distinctive physical structure and chemical properties.Additionally,its biocompatibility has garnered increasing attention.However,there is a lack of research on the biological effects and physical mechanisms of GDYprotein interactions at the molecular scale.In this study,the villin headpiece subdomain(HP35)served as a representative protein model.Molecular dynamics simulations were employed to investigate the interaction process between the HP35 protein and GDY,as well as the structural evolution of the protein.The data presented in our study demonstrate that GDY can rapidly adsorb HP35 protein and induce denaturation to one of the a-helix structures of HP35 protein.This implies a potential cytotoxicity concern of GDY for biological systems.Compared to graphene,GDY induced less disruption to HP35 protein.This can be attributed to the presence of natural triangular vacancies in GDY,which prevents p–p stacking action and the limited interaction of GDY with HP35 protein is not conducive to the expansion of protein structures.These findings unveil the biological effects of GDY at the molecular level and provide valuable insights for the application of GDY in biomedicine.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.51775077 and 51909023)。
文摘Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear-thinning fluid in a microchannel.We validated the feasibility of our simulation method by evaluating the mean square displacement and Reynolds number of the solution layers.The results show that the change rule of the fluid system's velocity profile and interaction energy can reflect the shear-thinning characteristics of the fluids.The velocity profile resembles a top-hat shape,intensifying as the fluid's power law index decreases.The interaction energy between the wall and the fluid decreases gradually with increasing velocity,and a high concentration of non-Newtonian fluid reaches a plateau sooner.Moreover,the velocity profile of the fluid is related to the molecule number density distribution and their values are inversely proportional.By analyzing the radial distribution function,we found that the hydrogen bonds between solute and water molecules weaken with the increase in velocity.This observation offers an explanation for the shear-thinning phenomenon of the non-Newtonian flow from a micro perspective.
基金supported by the National Natural Science Foundation of China (Grant No.11572090)。
文摘Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced critical load and hardness compared to nc-Ni,where perfect,stair-rod and Shockley dislocations are activated at (111),(111) and (111) slip planes in nt-Ni compared to only SSockley dislocation nucleation at (111) and (111) slip planes of nc-Ni.The nt-Ni exhibits a less significant indentation size effect in comparison with nc-Ni due to the dislocation slips hindrance of the twin boundary.The atomic deformation associated with the indentation size effect is investigated during dislocation transmission.Different from the decreasing partial slips parallel to the indenter surface in nc-Ni with increasing temperature,the temperaturedependent atomic deformation of nt-Ni is closely related to the twin boundary:from the partial slips parallel to the twin boundary (~10 K),to increased confined layer slips and decreased twin migration(300 K–600 K),to decreased confined layer slips and increased dislocation interaction of dislocation pinning and dissociation (900 K–1200 K).Dislocation density and atomic structure types through quantitative analysis are implemented to further reveal the above-mentioned dislocation motion and atomic structure alteration.Our study is helpful for understanding the temperature-dependent plasticity of twin boundary in nanotwinned materials.
基金Supported by the Sichuan Province Regional Innovation Cooperation Project(21QYCX0048)Sinopec Science and Technology Department Project(P21048-3)。
文摘According to the complex differential accumulation history of deep marine oil and gas in superimposed basins,the Lower Paleozoic petroleum system in Tahe Oilfield of Tarim Basin is selected as a typical case,and the process of hydrocarbon generation and expulsion,migration and accumulation,adjustment and transformation of deep oil and gas is restored by means of reservoine-forming dynamics simulation.The thermal evolution history of the Lower Cambrian source rocks in Tahe Oilfield reflects the obvious differences in hydrocarbon generation and expulsion process and intensity in different tectonic zones,which is the main reason controlling the differences in deep oil and gas phases.The complex transport system composed of strike-slip fault and unconformity,etc.controlled early migration and accumulation and late adjustment of deep oil and gas,while the Middle Cambrian gypsum-salt rock in inner carbonate platform prevented vertical migration and accumulation of deep oil and gas,resulting in an obvious"fault-controlled"feature of deep oil and gas,in which the low potential area superimposed by the NE-strike-slip fault zone and deep oil and gas migration was conducive to accumulation,and it is mainly beaded along the strike-slip fault zone in the northeast direction.The dynamic simulation of reservoir formation reveals that the spatio-temporal configuration of"source-fault-fracture-gypsum-preservation"controls the differential accumulation of deep oil and gas in Tahe Oilfield.The Ordovician has experienced the accumulation history of multiple periods of charging,vertical migration and accumulation,and lateral adjustment and transformation,and deep oil and gas have always been in the dynamic equilibrium of migration,accumulation and escape.The statistics of residual oil and gas show that the deep stratum of Tahe Oilfield still has exploration and development potential in the Ordovician Yingshan Formation and Penglaiba Formation,and the Middle and Upper Cambrian ultra-deep stratum has a certain oil and gas resource prospect.This study provides a reference for the dynamic quantitative evaluation of deep oil and gas in the Tarim Basin,and also provides a reference for the study of reservoir formation and evolution in carbonate reservoir of paleo-craton basin.
文摘The Brownian dynamics （BD） simulation of a dilute surfactant solution is conducted in a steady shear flow. The rodlike micelle is assumed as a rigid rod composed of lined-up beads. A novel intercluster potential model is introduced for describing the interactions between, micelles. In the model, the Lennard-Jones and the soft-sphere potentials are used as inter-bead potentials for end-end and interior-interior beads, respectively. The micelles are combined at their ends to form a network structure at lower shear rates and are disconnected to become more and more parallel to the shear flow direction with increasing shear rate. The change of micellar microstructures with the variation of the shear rate results in shear thinning characteristics of the computed shear viscosities and first normal stress difference coefficients. The effects of surfactant solution concentration on the micellar structures and rheological properties are also investigated. Results show that the shear viscosities and the first normal stress difference coefficients increase with increasing the viscosity of the surfactant solution.
文摘The track model used in the dynamic analysis and design system software is investigated. A home made tank is taken as an example to illustrate the method for modeling an integral tracked vehicle and perform the dynamic simulation. The obtained results have demonstrated that the simulation method has the advantage of high efficiency, more convenience and more insight into the dynamical behavior of the system.
基金Project supported by the National Natural Science Foundation of China （Grant Nos 50271026 and 50571037）.
文摘A molecular dynamics simulation study has been performed for the microstructure evolution in a liquid metal Ni system during crystallization process at two cooling rates by adopting the embedded atom method （EAM） model potential. The bond-type index method of Honeycutt-Andersen （HA） and a new cluster-type index method （CTIM-2） have been used to detect and analyse the microstructures in this system. It is demonstrated that the cooling rate plays a critical role in the microstructure evolution： below the crystallization temperature Tc, the effects of cooling rate are very remarkable and can be fully displayed. At different cooling rates of 2.0 × 10^13 K·s^-1 and 1.0 × 10^12 K·s^-1, two different kinds of crystal structures are obtained in the system. The first one is the coexistence of the hcp （expressed by （12 0 0 0 6 6） in CTIM-2） and the fcc （12 0 0 0 12 0） basic clusters consisting of 1421 and 1422 bond-types, and the hcp basic cluster becomes the dominant one with decreasing temperature, the second one is mainly the fcc （12 0 0 0 12 0） basic clusters consisting of 1421 bond-type, and their crystallization temperatures Tc would be 1073 and 1173 K, respectively.
基金supported in part by National Natural Science Foundation of China (Grant Nos 10474109 and 10674146)supported is part by the Shanghai Supercomputer Center of China
文摘In this paper molecular dynamics simulations are performed to study the accumulation behaviour of N2 and H2 at water/graphite interface under ambient temperature and pressure. It finds that both N2 and H2 molecules can accumulate at the interface and form one of two states according to the ratio of gas molecules number to square of graphite surface from our simulation results： gas films （pancake-like） for a larger ratio and nanobubbles for a smaller ratio. In addition, we discuss the stabilities of nanobubbles at different environment temperatures. Surprisingly, it is found that the density of both kinds of gas states can be greatly increased, even comparable with that of the liquid N2 and liquid H2. The present results are expected to be helpful for the understanding of the stable existence of gas film （pancake-like） and nanobubbles.
基金supported by the National Basic Research Program of China (2009CB623504)the National Science Foundation of China (21103231)Shanghai Science Foundation (11ZR1449700)
文摘Selectivity control is a difficult scientific and industrial challenge in methanol-to-olefins(MTO)conversion.It has been experimentally established that the topology of zeolite catalysts influenced the distribution of products.Besides the topology effect on reaction kinetics,the topology influences the diffusion of reactants and products in catalysts as well.In this work,by using COMPASS force-field molecular dynamics method,we investigated the intracrystalline diffusion of ethene and propene in four different zeolites,CHA,MFI,BEA and FAU,at different temperatures.The self-diffusion coefficients and diffusion activation barriers were calculated.A strong restriction on the diffusion of propene in CHA was observed because the self-diffusion coefficient ratio of ethene to propene is larger than 18 and the diffusion activation barrier of propene is more than 20 kJ/mol in CHA.This ratio decreases with the increase of temperature in the four investigated zeolites.The shape selectivity on products from diffusion perspective can provide some implications on the understanding of the selectivity difference between HSAPO-34 and HZSM-5 catalysts for the MTO conversion.
基金Project supported by the Key Laboratory of Mountain Hazards and Earth Surface Processes, the Chinese Academy of Sciences
文摘This paper performs molecular dynamics simulations to investigate the role of the monovalent cations K, Na and the divalent cation Ca on the stability and swelling of montmorillonite. The recently developed CLAYFF force field is used to predict the basal spacing as a function of the water content in the interlayer. The simulations reproduced the swelling pattern of these montmorillonites, suggesting a mechanism of their hydration different （K＋ 〈 Na＋ 〈 Ca2＋） from that of K＋-, Na＋-, and Ca2＋-montmorillonites. In particular, these results indicate that the valence of the cations has the larger impact on the behaviour of clay water systems. It also finds that the differences in size and hydration energy of K＋, Na＋ and Ca2＋ ions have strong implications for the structure of interlayer. This leads to the differences in the layer spacings of the simulated K＋-, Na＋-, and Ca2＋-montmorillonites. Furthermore, these simulations show that the K cations interact strongly with the clay sheets for the dehydrated clay sheets, but for the hydrated clays the Ca cations interact clearly strongly with the clay sheets.
基金supported by the National Natural Science Foundation of China(No.51176192)CAS Program(KGZD-EW-301)NOG Program(GHZ2012006003)
文摘Molecular dynamics simulations are performed to study the growth mechanism of CH4-CO2 mixed hydrate in xco2 = 75%, xco2 = 50%, and zco2 = 25% systems at T = 250 K, 255 K and 260 K, respectively. Our simulation results show that the growth rate of CH4-CO2 mixed hydrate increases as the CO2 concentration in the initial solution phase increases and the temperature decreases. Via hydrate formation, the composition of CO2 in hydrate phase is higher than that in initial solution phase and the encaging capacity of CO2 in hydrates increases with the decrease in temperature. By analysis of the cage occupancy ratio of CH4 molecules and CO2 molecules in large cages to small cages, we find that CO2 molecules are preferably encaged into the large cages of the hydrate crystal as compared with CH4 molecules. Interestingly, CH4 molecules and CO2 molecules frequently replace with each other in some particular cage sites adjacent to hydrate/solution interface during the crystal growth process. These two species of guest molecules eventually act to stabilize the newly formed hydrates, with CO2 molecules occupying large cages and CH4 molecules occupying small cages in hydrate.
基金supported by the National Natural Science Foundation of China (Grant No. 51076128)the National High Technology Research and Development Program of China (Grant No. 2009AA05Z107)
文摘We use non-equilibrium molecular dynamics simulations to calculate the self-diffusion coefficient, D, of a Lennard Jones fluid over a wide density and temperature range. The change in self-diffusion coefficient with temperature decreases by increasing density. For density ρ＊ = ρσ3 = 0.84 we observe a peak at the value of the self-diffusion coefficient and the critical temperature T＊ = kT/ε = 1.25. The value of the self-diffusion coefficient strongly depends on system size. The data of the self-diffusion coefficient are fitted to a simple analytic relation based on hydrodynamic arguments. This correction scales as N-α, where α is an adjustable parameter and N is the number of particles. It is observed that the values of a 〈 1 provide quite a good correction to the simulation data. The system size dependence is very strong for lower densities, but it is not as strong for higher densities. The self-diffusion coefficient calculated with non-equilibrium molecular dynamic simulations at different temperatures and densities is in good agreement with other calculations fronl the literature.
基金funded by the National Natural Science Foundation of China(31972096)Jilin Province Science and Technology Youth Talent Support Project(QT202021)Interdisciplinary Integration and Innovation Project of JLU(JLUXKJC2021QZ11)。
文摘The understanding of the structural requirements and the intermolecular-interaction mechanism are important for discovering potent angiotensin-converting enzyme(ACE)inhibitory peptides.In this study,we modifi ed an egg-white derived peptide,LAPYK,using the amino acids with different properties to produce the LAPYK-modified peptides.The ACE inhibitory activities of the modified peptides were determined to explore the structural requirements of ACE inhibitory peptides(ACEIPs).Molecular simulation and isothermal titration calorimetry analysis were used to investigate interactions between the peptides and ACE.We found that hydrophobicity and the amino acids with ring structures were benefi cial for the ACE inhibitory activities of the peptides.The results of the molecular mechanics poisson boltzmann surface area(MMPBSA)binding free energy calculations indicated that the polar solvation free energy(ΔG_(polar))of the charged peptides(LAPYK,LAPYE)were unfavorable for binding to ACE.On the other hand,the results of isothermal titration calorimetry analyses suggested that the enthalpy-driven ACE-peptide interactions were more favorable than the entropy-driven ACE-peptide interaction counterparts.
基金financially supported by the National Natural Science Foundation of China(Grant No.52004317,42090024)the Natural Science Foundation of Shandong Province of China(No.ZR2020ME091)+1 种基金the Fundamental Research Funds for the Central Universities(20CX06016A)the National Science and Technology Major Project(2017ZX05049-004)
文摘Shale oil formations contain both inorganic and organic media.The organic matter holds both free oil in the pores and dissolved oil within the kerogen molecules.The free oil flow in organic pores and the dissolved oil diffusion in kerogen molecules are coupled together.The molecular flow of free n-alkanes is an important process of shale oil accumulation and production.To study the dynamics of imbibition process of n-alkane molecules into kerogen slits,molecular dynamics(MD)simulations are conducted.Effects of slit width,temperature,and n-alkane types on the penetration speed,dynamic contact angle,and molecular conformations were analyzed.Results showed that molecular transportation of n-alkanes is dominated by molecular structure and molecular motion at this scale.The space-confinement conformational changes of molecules slow down the filling speeds in the narrow slits.The n-alkane molecules with long carbon chains require more time to undergo conformational changes.The high content of short-chain alkanes and high temperature facilitate the flow of alkane mixtures in kerogen slits.Results obtained from this study are useful for understanding the underlying nanoscale flow mechanism in shale formations.
基金the National Key R&D Program of China(Grant No.2019YFB1704600)the Hubei Provincial Natural Science Foundation of China(Grant No.2020CFA032).
文摘Currently,wire bonding is the most popular first-level interconnection technology used between the die and package terminals,but even with its long-term and excessive usage,the mechanism of wire bonding has not been completely evaluated.Therefore,fundamental research is still needed.In this study,the mechanism of microweld formation and breakage during Cu-Cu wire bonding was investigated by using molecular dynamics simulation.The contact model for the nanoindentation process between the wire and substrate was developed to simulate the contact process of the Cu wire and Cu substrate.Elastic contact and plastic instability were investigated through the loading and unloading processes.Moreover,the evolution of the indentation morphology and distributions of the atomic stress were also investigated.It was shown that the loading and unloading curves do not coincide,and the unloading curve exhibited hysteresis.For the substrate,in the loading process,the main force changed from attractive to repulsive.The maximum von Mises stress increased and shifted from the center toward the edge of the contact area.During the unloading process,the main force changed from repulsive to attractive.The Mises stress reduced first and then increased.Stress concentration occurs around dislocations in the middle area of the Cu wire.
