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.展开更多
Studying the relationship between ionic interactions and salt solubility in seawater has implications for seawater desalination and mineral extraction.In this paper,a new method of expressing ion-to-ion interaction is...Studying the relationship between ionic interactions and salt solubility in seawater has implications for seawater desalination and mineral extraction.In this paper,a new method of expressing ion-to-ion interaction is proposed by using molecular dynamics simulation,and the relationship between ion-to-ion interaction and salt solubility in a simulated seawater water-salt system is investigated.By analyzing the variation of distance and contact time between ions in an electrolyte solution,from both spatial and temporal perspectives,new parameters were proposed to describe the interaction between ions:interaction distance(ID),and interaction time ratio(ITR).The best correlation between characteristic time ratio and solubility was found for a molar ratio of salt-to-water of 10:100 with a correlation coefficient of 0.96.For the same salt,a positive correlation was found between CTR and the molar ratio of salt and water.For type 1-1,type 2-1,type 1-2,and type 2-2 salts,the correlation coefficients between CTR and solubility were 0.93,0.96,0.92,and 0.98 for a salt-to-water molar ratio of 10:100,respectively.The solubility of multiple salts was predicted by simulations and compared with experimental values,yielding an average relative deviation of 12.4%.The new ion-interaction parameters offer significant advantages in describing strongly correlated and strongly hydrated electrolyte solutions.展开更多
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.展开更多
Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectil...Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectile impacts is thoroughly investigated using full-atomic molecular dynamics(MD)simulations.The study results show that the porous structure and density are key factors determining the mechanical response of GA under impact loading.Specifically,the impact-induced penetration of the projectile leads to the collapse of the pore structure,causing stretching and subsequent rupture of covalent bonds in graphene sheets.Moreover,the effects of temperature on the mechanical performance of GA have been proven to be minimal,thereby highlighting the mechanical stability of GA over a wide range of temperatures.Finally,the energy absorption density(EAD)and energy absorption efficiency(EAE)metrics are adopted to assess the energy absorption capacity of GA during projectile penetration.The research findings of this work demonstrate the significant potential of GA for energy absorption applications.展开更多
In this work, we performed in situ nanoindentation in TEM to capture the real-time dislocation and twinning activities in pure Mg during loading and unloading. We demonstrated that the screw component of dislocations ...In this work, we performed in situ nanoindentation in TEM to capture the real-time dislocation and twinning activities in pure Mg during loading and unloading. We demonstrated that the screw component of dislocations glides continuously, while the edge components rapidly become sessile during loading. The twin tip propagation is intermittent, whereas the twin boundary migration is more continuous. During unloading, we observed the elastic strain relaxation causes both dislocation retraction and detwinning. Moreover,we note that the plastic zone comprised of dislocations in Mg is well-defined, which contrasts with the diffused plastic zones observed in face-centered cubic metals under the nanoindentation impressions. Additionally, molecular dynamics simulations were performed to study the formation and evolution of deformation-induced crystallographic defects at the early stages of indentation. We observed that,in addition to dislocations, the I1stacking fault bounded with a <1/2c+p> Frank loop can be generated from the plastic zone ahead of the indenter, and potentially serve as a nucleation source for abundant dislocations observed experimentally. These new findings are anticipated to provide new knowledge on the deformation mechanisms of Mg, which are difficult to obtain through conventional ex situ approaches. These observations may serve as a baseline for simulation work that investigate the dynamics of dislocation slip and twinning in Mg and alloys.展开更多
To better understand the benzene alkylation with chloroaluminate ionic liquids(ILs) as catalyst, the interfacial properties between the benzene/butene binary reactants and chloroaluminate ILs with varying cation alkyl...To better understand the benzene alkylation with chloroaluminate ionic liquids(ILs) as catalyst, the interfacial properties between the benzene/butene binary reactants and chloroaluminate ILs with varying cation alkyl chain length and different anions were investigated using molecular dynamics(MD) simulations. The results indicate that ILs can obviously improve the interfacial width, solubility and diffusion of reactants compared to H_(2)SO_(4). The longer alkyl chains of cations present a density enrichment at the interface and protrude into the binary reactants phase. Furthermore, the ILs consisting of 1-octyl-3-methylimidazolium cations([Omim]^(+)) and the stronger acidity heptachlorodialuminate anions([Al_(2)Cl_(7)]^(-)) are more beneficial to promote the interfacial width and facilitate the dissolution and diffusion of benzene in both the IL bulk and the interfacial region in comparison to the ones with shorter alkyl chains cations and weaker acidity anions. The information gives us a better guideline for the design of ILs for benzene alkylation.展开更多
The interface properties of Fe(101)/zinc silicate modified by organo-siloxane(KH-570)was studied by using the method of molecular dynamics simulation.By calculating the temperature and energy fluctuation of equilibriu...The interface properties of Fe(101)/zinc silicate modified by organo-siloxane(KH-570)was studied by using the method of molecular dynamics simulation.By calculating the temperature and energy fluctuation of equilibrium state,equilibrium concentration distribution,MSD of layer and different groups,and interaction energy of two interface models,the influencing mechanism on the interface properties of adding organosiloxane into coating system was studied at the atomic scale.It shows that the temperature and energy of interface oscillated in a small range and it was exited in a state of dynamic equilibrium within the initial simulation stage(t<20 ps).It can be seen from the multiple peak states of concentration distribution that the iron substrate,organo-siloxane and zinc silicate are distributed in the form of a concentration gradient in the real environment.The rapid diffusion of free zinc powder in zinc silicate coating was the essential reason that affected the comprehensive properties of coating.The interface thickness decreased from 7.45 to 6.82Å,the MSD of free zinc powder was effectively reduced,and the interfacial energy was increased from 104.667 to 347.158 kcal/mol after being modified by organo-siloxane.展开更多
The molecular dynamics method is used to investigate decomposition of methane hydrate at different temperatures,pressures and concentrations of inhibitor.By analyzing the parameters of system conformation,mean square ...The molecular dynamics method is used to investigate decomposition of methane hydrate at different temperatures,pressures and concentrations of inhibitor.By analyzing the parameters of system conformation,mean square displacement and radial distribution function,the decomposition of hydrate in the presence of alcohol inhibitors ethylene glycol and glycerol is explored.The results show that the hydroxyl groups in alcohol molecules can destroy the cage structure of hydrate,and form hydrogen bonds with nearby water molecules to effectively prevent the reformation of hydrate.Therefore,ethylene glycol and glycerol serve as inhibitors of methane hydrate,furthermore,in terms of inhibition effect,glycerol is better than ethylene glycol by comparing rate of hydrate decomposition.展开更多
Background:To elucidate the molecular mechanisms of Curcuma longa(C.longa)in breast cancer treatment.Methods:Phytocompounds of C.longa were obtained from Dr.Duke’s Phytochemical and Ethnobotanical Database.Potential ...Background:To elucidate the molecular mechanisms of Curcuma longa(C.longa)in breast cancer treatment.Methods:Phytocompounds of C.longa were obtained from Dr.Duke’s Phytochemical and Ethnobotanical Database.Potential active targets were retrieved from Bindingdb,SEA and Swiss Target Prediction databases.Breast cancer targets were retrieved from the Therapeutic Target Database.Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were done using DAVID and KOBAS3.0 databases respectively.The Cytoscape software was used to construct the phytocompound-target-pathway network.The PyRx and Desmond software were utilized for molecular docking and molecular dynamics simulation respectively.Results:Out of one hundred and fifty-six phytocompounds,fifty-four modulated proteins involved in breast cancer.Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated C.longa exerts its therapeutic effect through regulating several key pathways.Molecular docking analysis revealed that most phytocompounds of C.longa had a good affinity with the key targets.Molecular dynamics simulation showed that ethinylestradiol formed stable ligand-protein complexes.Conclusion:The results of this study will enhance our understanding of the potential molecular mechanisms by which C.longa inhibits breast cancer and lay a foundation for future experimental studies.展开更多
The human serotonin transporter(SERT)terminates neurotransmission by removing serotonin from the synaptic cleft,which is an essential process that plays an important role in depression.In addition to natural substrate...The human serotonin transporter(SERT)terminates neurotransmission by removing serotonin from the synaptic cleft,which is an essential process that plays an important role in depression.In addition to natural substrate serotonin,SERT is also the target of the abused drug cocaine and,clinically used antidepressants,escitalopram,and paroxetine.To date,few studies have attempted to investigate the unbinding mechanism underlying the orthosteric and allosteric modulation of SERT.In this article,the conserved property of the orthosteric and allosteric sites(S1 and S2)of SERT was revealed by combining the high resolutions of x-ray crystal structures and molecular dynamics(MD)simulations.The residues Tyr95 and Ser438 located within the S1 site,and Arg104 located within the S2 site in SERT illustrate conserved interactions(hydrogen bonds and hydrophobic interactions),as responses to selective serotonin reuptake inhibitors.Van der Waals interactions were keys to designing effective drugs inhibiting SERT and further,electrostatic interactions highlighted escitalopram as a potent antidepressant.We found that cocaine,escitalopram,and paroxetine,whether the S1 site or the S2 site,were more competitive.According to this potential of mean force(PMF)simulations,the new insights reveal the principles of competitive inhibitors that lengths of trails from central SERT to an opening were~18A for serotonin and~22 A for the above-mentioned three drugs.Furthermore,the distance between the natural substrate serotonin and cocaine(or escitalopram)at the allosteric site was~3A.Thus,it can be inferred that the potent antidepressants tended to bind at deeper positions of the S1 or the S2 site of SERT in comparison to the substrate.Continuing exploring the processes of unbinding four ligands against the two target pockets of SERT,this study observed a broad pathway in which serotonin,cocaine,escitalopram(at the S1 site),and paroxetine all were pulled out to an opening between MT1b and MT6a,which may be helpful to understand the dissociation mechanism of antidepressants.展开更多
Protein adsorption preferentially occurs and significantly affects the physicochemical reactions once the biodegradable magnesium alloys as bone replacements have been implanted. To date, interactions mechanisms betwe...Protein adsorption preferentially occurs and significantly affects the physicochemical reactions once the biodegradable magnesium alloys as bone replacements have been implanted. To date, interactions mechanisms between Mg implants and proteins remain unclear at a molecular level. Thereby, a combination of molecular dynamic(MD) simulations and experimental exploration is used to investigate the adsorption behavior and conformational change of bovine serum albumin(BSA), a representative protein of blood plasma, upon the surface of microarc oxidation(MAO) coated Mg alloy AZ31. The influences of absorbed proteins on the cytocompatibility of MAO coating are evaluated by virtue of cytotoxicity assay. Results indicate that the negatively charged O atoms(BSA) exhibit strong interaction with Mg^(2+) ions of Mg(OH)_(2), revealing that BSA molecules are ionically adsorbed on the AZ31 surface. Interestingly, MD simulation reveals that MAO coating demonstrates superior ability to capture BSA molecules during the process of adsorption owing to strong electric attraction between the negatively charged O atoms in BSA molecules with Mg atoms of MgO in MAO coating. Moreover, the α-helix part of absorbed BSA molecules on AZ31 substrate and MAO coating markedly decreases with an increase in β-sheet, β-turn and unordered contents, which is attributed to the reduction in the number of hydrogen bonds in BSA molecules. Furthermore, the adsorbed BSA molecules improve the cytocompatibility of MAO coating since the positively charged-NH_(3)^(+) group and β-sheet content of absorbed BSA molecules mediate the cell adhesion by interacting with the negatively charged cell membrane.展开更多
The mechanical property and deformation mechanism of twinned gold nanowire with non-uniform distribution of twinned boundaries(TBs)are studied by the molecular dynamics(MD)method.It is found that the twin boundary spa...The mechanical property and deformation mechanism of twinned gold nanowire with non-uniform distribution of twinned boundaries(TBs)are studied by the molecular dynamics(MD)method.It is found that the twin boundary spacing(TBS)has a great effect on the strength and plasticity of the nanowires with uniform distribution of TBs.And the strength enhances with the decrease of TBS,while its plasticity declines.For the nanowires with non-uniform distribution of TBs,the differences in distribution among different TBSs have little effect on the Young's modulus or strength,and the compromise in strength appears.But the differences have a remarkable effect on the plasticity of twinned gold nanowire.The twinned gold nanowire with higher local symmetry ratio has better plasticity.The initial dislocations always form in the largest TBS and the fracture always appears at or near the twin boundaries adjacent to the smallest TBS.Some simulation results are consistent with the experimental results.展开更多
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.展开更多
Single crystal silicon carbide(SiC)is widely used for optoelectronics applications.Due to the anisotropic characteristics of single crystal materials,the C face and Si face of single crystal SiC have different physica...Single crystal silicon carbide(SiC)is widely used for optoelectronics applications.Due to the anisotropic characteristics of single crystal materials,the C face and Si face of single crystal SiC have different physical properties,which may fit for particular application purposes.This paper presents an investigation of the material removal and associated subsurface defects in a set of scratching tests on the C face and Si face of 4H-SiC and 6H-SiC materials using molecular dynamics simulations.The investigation reveals that the sample material deformation consists of plastic,amorphous transformations and dislocation slips that may be prone to brittle split.The results showed that the material removal at the C face is more effective with less amorphous deformation than that at the Si face.Such a phenomenon in scratching relates to the dislocations on the basal plane(0001)of the SiC crystal.Subsurface defects were reduced by applying scratching cut depths equal to integer multiples of a half molecular lattice thickness,which formed a foundation for selecting machining control parameters for the best surface quality.展开更多
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.展开更多
The mobility and aggregation behavior of macromolecular lubricant oxidation products and their influences on the performance of base stock were probed by molecular dynamics(MD)simulation.The mean square displacement(M...The mobility and aggregation behavior of macromolecular lubricant oxidation products and their influences on the performance of base stock were probed by molecular dynamics(MD)simulation.The mean square displacement(MSD)of molecules was calculated to explore the mobility of molecules.The distribution appearance of lubricant oxidation products in models was acquired to explore the aggregation of molecules.The results show that the mobility of macromolecular oxidation products is lower than that of base stock.The MSD of macromolecular oxidation products reduces with an increasing macromolecular weight.Macromolecular oxidation products can also decrease the mobility of base stock.The interaction energy between the macromolecules and the base stock soars with the increase of macromolecular weight.Macromolecules with a larger molecular weight can affect more base stock molecules with stronger restriction,which leads to lower mobility of base stock molecules.There are aggregates formed among macromolecular oxidation products,and the molecules in aggregates are connected by hydrogen bonds.The quantity of hydrogen bonds in aggregates is related to temperature.展开更多
A molecuIar dynamics simulation has been used to study zero-fluence and low- bombarding63Cu-65Cu sputtering. Calculations show that the isotopic ratios at both θ≤35° and energy 63Cu-65Cu sputtering. Calculation...A molecuIar dynamics simulation has been used to study zero-fluence and low- bombarding63Cu-65Cu sputtering. Calculations show that the isotopic ratios at both θ≤35° and energy 63Cu-65Cu sputtering. Calculations show that the isotopic ratios at both θ≤35° and θ >35°, and the total isotopic ratio increase when the bombarding energy decreases. This result might impty the existence of the bombarding-energy-dependent momentum asymmetry.展开更多
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.展开更多
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.展开更多
This work applied molecular dynamics(MD)simulation to calculate densities of natural gas mixtures at extremely high pressure(>138 MPa)and high temperature(>200℃)conditions(x HPHT)to bridge the knowledge and tec...This work applied molecular dynamics(MD)simulation to calculate densities of natural gas mixtures at extremely high pressure(>138 MPa)and high temperature(>200℃)conditions(x HPHT)to bridge the knowledge and technical gaps between experiments and classical theories.The experimental data are scarce at these conditions which are also out of assumptions for classical predictive correlations,such as the Dranchuk&Abou-Kassem(DAK)equation of state(EOS).Force fields of natural gas components were carefully chosen from literatures and the simulation results are validated with experimental data.The largest relative error is 2.67%for pure hydrocarbons,2.99%for C1/C3 mixture,7.85%for C1/C4 mixture,and 8.47%for pure H2S.These satisfactory predictions demonstrate that the MD simulation approach is reliable to predict natural-and acid-gases thermodynamic properties.The validated model is further used to generate data for the study of the EOS with pressure up to 276 MPa and temperature up to 573 K.Our results also reveal that the Dranchuk&Abou-Kassem(DAK)EOS is capable of predicting natural gas compressibility to a satisfactory accuracy at x HPHT conditions,which extends the confidence range of the DAK EOS.展开更多
基金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.
基金supported by the National Natural Science Foundation of China(No.21776264).
文摘Studying the relationship between ionic interactions and salt solubility in seawater has implications for seawater desalination and mineral extraction.In this paper,a new method of expressing ion-to-ion interaction is proposed by using molecular dynamics simulation,and the relationship between ion-to-ion interaction and salt solubility in a simulated seawater water-salt system is investigated.By analyzing the variation of distance and contact time between ions in an electrolyte solution,from both spatial and temporal perspectives,new parameters were proposed to describe the interaction between ions:interaction distance(ID),and interaction time ratio(ITR).The best correlation between characteristic time ratio and solubility was found for a molar ratio of salt-to-water of 10:100 with a correlation coefficient of 0.96.For the same salt,a positive correlation was found between CTR and the molar ratio of salt and water.For type 1-1,type 2-1,type 1-2,and type 2-2 salts,the correlation coefficients between CTR and solubility were 0.93,0.96,0.92,and 0.98 for a salt-to-water molar ratio of 10:100,respectively.The solubility of multiple salts was predicted by simulations and compared with experimental values,yielding an average relative deviation of 12.4%.The new ion-interaction parameters offer significant advantages in describing strongly correlated and strongly hydrated electrolyte solutions.
基金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(No.12102256).
文摘Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectile impacts is thoroughly investigated using full-atomic molecular dynamics(MD)simulations.The study results show that the porous structure and density are key factors determining the mechanical response of GA under impact loading.Specifically,the impact-induced penetration of the projectile leads to the collapse of the pore structure,causing stretching and subsequent rupture of covalent bonds in graphene sheets.Moreover,the effects of temperature on the mechanical performance of GA have been proven to be minimal,thereby highlighting the mechanical stability of GA over a wide range of temperatures.Finally,the energy absorption density(EAD)and energy absorption efficiency(EAE)metrics are adopted to assess the energy absorption capacity of GA during projectile penetration.The research findings of this work demonstrate the significant potential of GA for energy absorption applications.
基金supported by the National Science Foundation (Division of Materials Research,Program Manager: Dr.Jonathan Madison) under Grant No.2144973financially supported by the startup fund from Drexel University。
文摘In this work, we performed in situ nanoindentation in TEM to capture the real-time dislocation and twinning activities in pure Mg during loading and unloading. We demonstrated that the screw component of dislocations glides continuously, while the edge components rapidly become sessile during loading. The twin tip propagation is intermittent, whereas the twin boundary migration is more continuous. During unloading, we observed the elastic strain relaxation causes both dislocation retraction and detwinning. Moreover,we note that the plastic zone comprised of dislocations in Mg is well-defined, which contrasts with the diffused plastic zones observed in face-centered cubic metals under the nanoindentation impressions. Additionally, molecular dynamics simulations were performed to study the formation and evolution of deformation-induced crystallographic defects at the early stages of indentation. We observed that,in addition to dislocations, the I1stacking fault bounded with a <1/2c+p> Frank loop can be generated from the plastic zone ahead of the indenter, and potentially serve as a nucleation source for abundant dislocations observed experimentally. These new findings are anticipated to provide new knowledge on the deformation mechanisms of Mg, which are difficult to obtain through conventional ex situ approaches. These observations may serve as a baseline for simulation work that investigate the dynamics of dislocation slip and twinning in Mg and alloys.
基金The financial support by State Key Laboratory of Heavy Oil Processing, China University of Petroleum。
文摘To better understand the benzene alkylation with chloroaluminate ionic liquids(ILs) as catalyst, the interfacial properties between the benzene/butene binary reactants and chloroaluminate ILs with varying cation alkyl chain length and different anions were investigated using molecular dynamics(MD) simulations. The results indicate that ILs can obviously improve the interfacial width, solubility and diffusion of reactants compared to H_(2)SO_(4). The longer alkyl chains of cations present a density enrichment at the interface and protrude into the binary reactants phase. Furthermore, the ILs consisting of 1-octyl-3-methylimidazolium cations([Omim]^(+)) and the stronger acidity heptachlorodialuminate anions([Al_(2)Cl_(7)]^(-)) are more beneficial to promote the interfacial width and facilitate the dissolution and diffusion of benzene in both the IL bulk and the interfacial region in comparison to the ones with shorter alkyl chains cations and weaker acidity anions. The information gives us a better guideline for the design of ILs for benzene alkylation.
基金supported by the National Science Fundation of China(No.U1937601),and the National Natural Science Foundation of China(Grant No.NSFC51905471).
文摘The interface properties of Fe(101)/zinc silicate modified by organo-siloxane(KH-570)was studied by using the method of molecular dynamics simulation.By calculating the temperature and energy fluctuation of equilibrium state,equilibrium concentration distribution,MSD of layer and different groups,and interaction energy of two interface models,the influencing mechanism on the interface properties of adding organosiloxane into coating system was studied at the atomic scale.It shows that the temperature and energy of interface oscillated in a small range and it was exited in a state of dynamic equilibrium within the initial simulation stage(t<20 ps).It can be seen from the multiple peak states of concentration distribution that the iron substrate,organo-siloxane and zinc silicate are distributed in the form of a concentration gradient in the real environment.The rapid diffusion of free zinc powder in zinc silicate coating was the essential reason that affected the comprehensive properties of coating.The interface thickness decreased from 7.45 to 6.82Å,the MSD of free zinc powder was effectively reduced,and the interfacial energy was increased from 104.667 to 347.158 kcal/mol after being modified by organo-siloxane.
基金Sponsored by the National Natural Science Foundation of China(Grant No.51876032)the Natural Science Foundation of Heilongjiang Province of China(Grant No.ZD2019E002).
文摘The molecular dynamics method is used to investigate decomposition of methane hydrate at different temperatures,pressures and concentrations of inhibitor.By analyzing the parameters of system conformation,mean square displacement and radial distribution function,the decomposition of hydrate in the presence of alcohol inhibitors ethylene glycol and glycerol is explored.The results show that the hydroxyl groups in alcohol molecules can destroy the cage structure of hydrate,and form hydrogen bonds with nearby water molecules to effectively prevent the reformation of hydrate.Therefore,ethylene glycol and glycerol serve as inhibitors of methane hydrate,furthermore,in terms of inhibition effect,glycerol is better than ethylene glycol by comparing rate of hydrate decomposition.
文摘Background:To elucidate the molecular mechanisms of Curcuma longa(C.longa)in breast cancer treatment.Methods:Phytocompounds of C.longa were obtained from Dr.Duke’s Phytochemical and Ethnobotanical Database.Potential active targets were retrieved from Bindingdb,SEA and Swiss Target Prediction databases.Breast cancer targets were retrieved from the Therapeutic Target Database.Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were done using DAVID and KOBAS3.0 databases respectively.The Cytoscape software was used to construct the phytocompound-target-pathway network.The PyRx and Desmond software were utilized for molecular docking and molecular dynamics simulation respectively.Results:Out of one hundred and fifty-six phytocompounds,fifty-four modulated proteins involved in breast cancer.Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated C.longa exerts its therapeutic effect through regulating several key pathways.Molecular docking analysis revealed that most phytocompounds of C.longa had a good affinity with the key targets.Molecular dynamics simulation showed that ethinylestradiol formed stable ligand-protein complexes.Conclusion:The results of this study will enhance our understanding of the potential molecular mechanisms by which C.longa inhibits breast cancer and lay a foundation for future experimental studies.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11904036 and 12175081)Fundamental Research Funds for the Central Universities(Grant No.CCNU22QNOO4)。
文摘The human serotonin transporter(SERT)terminates neurotransmission by removing serotonin from the synaptic cleft,which is an essential process that plays an important role in depression.In addition to natural substrate serotonin,SERT is also the target of the abused drug cocaine and,clinically used antidepressants,escitalopram,and paroxetine.To date,few studies have attempted to investigate the unbinding mechanism underlying the orthosteric and allosteric modulation of SERT.In this article,the conserved property of the orthosteric and allosteric sites(S1 and S2)of SERT was revealed by combining the high resolutions of x-ray crystal structures and molecular dynamics(MD)simulations.The residues Tyr95 and Ser438 located within the S1 site,and Arg104 located within the S2 site in SERT illustrate conserved interactions(hydrogen bonds and hydrophobic interactions),as responses to selective serotonin reuptake inhibitors.Van der Waals interactions were keys to designing effective drugs inhibiting SERT and further,electrostatic interactions highlighted escitalopram as a potent antidepressant.We found that cocaine,escitalopram,and paroxetine,whether the S1 site or the S2 site,were more competitive.According to this potential of mean force(PMF)simulations,the new insights reveal the principles of competitive inhibitors that lengths of trails from central SERT to an opening were~18A for serotonin and~22 A for the above-mentioned three drugs.Furthermore,the distance between the natural substrate serotonin and cocaine(or escitalopram)at the allosteric site was~3A.Thus,it can be inferred that the potent antidepressants tended to bind at deeper positions of the S1 or the S2 site of SERT in comparison to the substrate.Continuing exploring the processes of unbinding four ligands against the two target pockets of SERT,this study observed a broad pathway in which serotonin,cocaine,escitalopram(at the S1 site),and paroxetine all were pulled out to an opening between MT1b and MT6a,which may be helpful to understand the dissociation mechanism of antidepressants.
基金supported by the National Natural Science Foundation of China (52071191)。
文摘Protein adsorption preferentially occurs and significantly affects the physicochemical reactions once the biodegradable magnesium alloys as bone replacements have been implanted. To date, interactions mechanisms between Mg implants and proteins remain unclear at a molecular level. Thereby, a combination of molecular dynamic(MD) simulations and experimental exploration is used to investigate the adsorption behavior and conformational change of bovine serum albumin(BSA), a representative protein of blood plasma, upon the surface of microarc oxidation(MAO) coated Mg alloy AZ31. The influences of absorbed proteins on the cytocompatibility of MAO coating are evaluated by virtue of cytotoxicity assay. Results indicate that the negatively charged O atoms(BSA) exhibit strong interaction with Mg^(2+) ions of Mg(OH)_(2), revealing that BSA molecules are ionically adsorbed on the AZ31 surface. Interestingly, MD simulation reveals that MAO coating demonstrates superior ability to capture BSA molecules during the process of adsorption owing to strong electric attraction between the negatively charged O atoms in BSA molecules with Mg atoms of MgO in MAO coating. Moreover, the α-helix part of absorbed BSA molecules on AZ31 substrate and MAO coating markedly decreases with an increase in β-sheet, β-turn and unordered contents, which is attributed to the reduction in the number of hydrogen bonds in BSA molecules. Furthermore, the adsorbed BSA molecules improve the cytocompatibility of MAO coating since the positively charged-NH_(3)^(+) group and β-sheet content of absorbed BSA molecules mediate the cell adhesion by interacting with the negatively charged cell membrane.
基金the National Natural Science Foundation of China(Grant No.51771033).
文摘The mechanical property and deformation mechanism of twinned gold nanowire with non-uniform distribution of twinned boundaries(TBs)are studied by the molecular dynamics(MD)method.It is found that the twin boundary spacing(TBS)has a great effect on the strength and plasticity of the nanowires with uniform distribution of TBs.And the strength enhances with the decrease of TBS,while its plasticity declines.For the nanowires with non-uniform distribution of TBs,the differences in distribution among different TBSs have little effect on the Young's modulus or strength,and the compromise in strength appears.But the differences have a remarkable effect on the plasticity of twinned gold nanowire.The twinned gold nanowire with higher local symmetry ratio has better plasticity.The initial dislocations always form in the largest TBS and the fracture always appears at or near the twin boundaries adjacent to the smallest TBS.Some simulation results are consistent with the experimental results.
基金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.
基金financial support from National Natural Science Foundation of China(Grant No.51835004 and 51575197)Huaqiao University International Cultivation Program for Outstanding Postgraduates and Subsidized Projec for Postgraduates’Innovative Fund in Scientific Research of Huaqiao University(No.18011080010)。
文摘Single crystal silicon carbide(SiC)is widely used for optoelectronics applications.Due to the anisotropic characteristics of single crystal materials,the C face and Si face of single crystal SiC have different physical properties,which may fit for particular application purposes.This paper presents an investigation of the material removal and associated subsurface defects in a set of scratching tests on the C face and Si face of 4H-SiC and 6H-SiC materials using molecular dynamics simulations.The investigation reveals that the sample material deformation consists of plastic,amorphous transformations and dislocation slips that may be prone to brittle split.The results showed that the material removal at the C face is more effective with less amorphous deformation than that at the Si face.Such a phenomenon in scratching relates to the dislocations on the basal plane(0001)of the SiC crystal.Subsurface defects were reduced by applying scratching cut depths equal to integer multiples of a half molecular lattice thickness,which formed a foundation for selecting machining control parameters for the best surface quality.
基金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.
基金The authors are grateful for the calculation support of the Key Laboratory of Molecular Oil Refining of the Research Institute of Petroleum Processing in SINOPEC and the financial supports from the Natural Science Foundation of China(NSFCNo.51671100)+2 种基金the State Key Laboratory of Metal Material for Marine Equipment and Application-School of Material and Metallurgy,University of Science and Technology Liaoning Co-project(No.SKLMEA-USTLN 201905)the University of Science and Technology Liaoning Talent Project(No.601010314)the University of Science and Technology Liaoning Young Teachers Fund(No.2019QN08).
文摘The mobility and aggregation behavior of macromolecular lubricant oxidation products and their influences on the performance of base stock were probed by molecular dynamics(MD)simulation.The mean square displacement(MSD)of molecules was calculated to explore the mobility of molecules.The distribution appearance of lubricant oxidation products in models was acquired to explore the aggregation of molecules.The results show that the mobility of macromolecular oxidation products is lower than that of base stock.The MSD of macromolecular oxidation products reduces with an increasing macromolecular weight.Macromolecular oxidation products can also decrease the mobility of base stock.The interaction energy between the macromolecules and the base stock soars with the increase of macromolecular weight.Macromolecules with a larger molecular weight can affect more base stock molecules with stronger restriction,which leads to lower mobility of base stock molecules.There are aggregates formed among macromolecular oxidation products,and the molecules in aggregates are connected by hydrogen bonds.The quantity of hydrogen bonds in aggregates is related to temperature.
文摘A molecuIar dynamics simulation has been used to study zero-fluence and low- bombarding63Cu-65Cu sputtering. Calculations show that the isotopic ratios at both θ≤35° and energy 63Cu-65Cu sputtering. Calculations show that the isotopic ratios at both θ≤35° and θ >35°, and the total isotopic ratio increase when the bombarding energy decreases. This result might impty the existence of the bombarding-energy-dependent momentum asymmetry.
基金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.
基金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.
基金partial financial support from Ballard Petroleum Holdings and Yangtze Universitythe Schooner Supercomputing from the University of Oklahomathe startup support from the University of Oklahoma。
文摘This work applied molecular dynamics(MD)simulation to calculate densities of natural gas mixtures at extremely high pressure(>138 MPa)and high temperature(>200℃)conditions(x HPHT)to bridge the knowledge and technical gaps between experiments and classical theories.The experimental data are scarce at these conditions which are also out of assumptions for classical predictive correlations,such as the Dranchuk&Abou-Kassem(DAK)equation of state(EOS).Force fields of natural gas components were carefully chosen from literatures and the simulation results are validated with experimental data.The largest relative error is 2.67%for pure hydrocarbons,2.99%for C1/C3 mixture,7.85%for C1/C4 mixture,and 8.47%for pure H2S.These satisfactory predictions demonstrate that the MD simulation approach is reliable to predict natural-and acid-gases thermodynamic properties.The validated model is further used to generate data for the study of the EOS with pressure up to 276 MPa and temperature up to 573 K.Our results also reveal that the Dranchuk&Abou-Kassem(DAK)EOS is capable of predicting natural gas compressibility to a satisfactory accuracy at x HPHT conditions,which extends the confidence range of the DAK EOS.
