Molecular dynamics simulation of the flow mechanism of shear-thinning fluids in a microchannel

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.

Study of the Relationship Between New Ionic Interaction Parameters and Salt Solubility in Electrolyte Solutions Based on Molecular Dynamics Simulation

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.

Molecular Dynamics Simulation of Shock Response of CL-20 Co-crystals Containing Void Defects

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.

Temperature-Induced Unfolding Pathway of Staphylococcal Enterotoxin B:Insights from Circular Dichroism and Molecular Dynamics Simulation

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.

Investigation of Projectile Impact Behaviors of Graphene Aerogel Using Molecular Dynamics Simulations

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.

Probing the interaction between asphaltene-wax and its effects on the crystallization behavior of waxes in heavy oil via molecular dynamics simulation

High content of asphaltenes and waxes leads to the high pour point and the poor flowability of heavy oil,which is adverse to its efficient development and its transportation in pipe.Understanding the interaction mechanism between asphaltene-wax is crucial to solve these problems,but it is still unclear.In this paper,molecular dynamics simulation was used to investigate the interaction between asphaltenewax and its effects on the crystallization behavior of waxes in heavy oil.Results show that molecules in pure wax are arranged in a paralleled geometry.But wax molecules in heavy oil,which are close to the surface of asphaltene aggregates,are bent and arranged irregularly.When the mass fraction of asphaltenes in asphaltene-wax system(ω_(asp))is 0-25 wt%,the attraction among wax molecules decreases and the bend degree of wax molecules increases with the increase ofω_(asp).Theω_(asp)increases from 0 to 25 wt%,and the attraction between asphaltene-wax is stronger than that among waxes.This causes that the wax precipitation point changes from 353 to 333 K.While theω_(asp)increases to 50 wt%,wax molecules are more dispersed owing to the steric hindrance of asphaltene aggregates,and the interaction among wax molecules transforms from attraction to repulsion.It causes that the ordered crystal structure of waxes can't be formed at normal temperature.Simultaneously,the asphaltene,with the higher molecular weight or the more hetero atoms,has more obvious inhibition to the formation of wax crystals.Besides,resins also have an obvious inhibition on the wax crystal due to the formation of asphalteneresin aggregates with a larger radius.Our results reveal the interaction mechanism between asphaltene-wax,and provide useful guidelines for the development of heavy oil.

Molecular dynamics simulation study of nitrogen vacancy color centers prepared by carbon ion implantation into diamond

Nitrogen vacancy(NV)color centers in diamond have useful applications in quantum sensing andfluorescent marking.They can be gen-erated experimentally by ion implantation,femtosecond lasers,and chemical vapor deposition.However,there is a lack of studies of the yield of NV color centers at the atomic scale.In the molecular dynamics simulations described in this paper,NV color centers are pre-pared by ion implantation in diamond with pre-doped nitrogen and subsequent annealing.The differences between the yields of NV color centers produced by implantation of carbon(C)and nitrogen(N)ions,respectively,are investigated.It is found that C-ion implantation gives a greater yield of NV color centers and superior location accuracy.The effects of different pre-doping concentrations(400–1500 ppm)and implantation energies(1.0–3.0 keV)on the NV color center yield are analyzed,and it is shown that a pre-doping concentra-tion of 1000 ppm with 2 keV C-ion implantation can produce a 13%yield of NV color centers after 1600 K annealing for 7.4 ns.Finally,a brief comparison of the NV color center identification methods is presented,and it is found that the error rate of an analysis utiliz-ing the identify diamond structure coordination analysis method is reduced by about 7%compared with conventional identification+methods.

Molecular dynamics simulations on the interactions between nucleic acids and a phospholipid bilayer

Recently,lipid nanoparticles(LNPs)have been extensively investigated as non-viral carriers of nucleic acid vaccines due to their high transport efficiency,safety,and straightforward production and scalability.However,the molecular mechanism underlying the interactions between nucleic acids and phospholipid bilayers within LNPs remains elusive.In this study,we employed the all-atom molecular dynamics simulation to investigate the interactions between single-stranded nucleic acids and a phospholipid bilayer.Our findings revealed that hydrophilic bases,specifically G in single-stranded RNA(ssRNA)and single-stranded DNA(ssDNA),displayed a higher propensity to form hydrogen bonds with phospholipid head groups.Notably,ssRNA exhibited stronger binding energy than ssDNA.Furthermore,divalent ions,particularly Ca2+,facilitated the binding of ssRNA to phospholipids due to their higher binding energy and lower dissociation rate from phospholipids.Overall,our study provides valuable insights into the molecular mechanisms underlying nucleic acidphospholipid interactions,with potential implications for the nucleic acids in biotherapies,particularly in the context of lipid carriers.

Mechanical Properties of Copper with Dendritic Silver Inclusions:Insights from Molecular Dynamics Simulations

This study explores the mechanical behavior of single-crystal copper with silver inclusions,focusing on the effects of dendritic and spherical geometries using molecular dynamics simulations.Uniaxial tensile tests reveal that dendritic inclusions lead to an earlier onset of plasticity due to the presence of high-strain regions at the complex inclusion/matrix interfaces,whereas spherical inclusions exhibit delayed plasticity associated with their symmetric geometry and homogeneous strain distribution.During the plastic regime,the dislocation density is primarily influenced by the volume fraction of silver inclusions rather than their shape,with spherical inclusions showing the highest densities due to their larger volume and higher silver content.Stacking faults,quantified via hexagonal closed-packed atom populations,are strongly correlated with dislocation activity but exhibit transient behavior,indicating that many faults are swept out or transformed during deformation.This transfient effect is observed in all cases,independently of the inclusion size.These findings highlight the complex interplay between inclusion geometry,dislocation activity,and stacking fault evolution in shaping the mechanical properties of copper.The study underscores the need to account for inclusion morphology and defect dynamics when designing advanced copper-based materials and suggests further investigations into the role of dendrite orientation and distribution to enhance material performance in engineering applications.

Probing the electric double layer structure at nitrogen-doped graphite electrodes by constant-potential molecular dynamics simulations

Electric double layer(EDL)is a critical topic in electrochemistry and largely determines the working performance of lithium batteries.However,atomic insights into the EDL structures on heteroatom-modified graphite anodes and EDL evolution with electrode potential are very lacking.Herein,a constant-potential molecular dynamics(CPMD)method is proposed to probe the EDL structure under working conditions,taking N-doped graphite electrodes and carbonate electrolytes as an example.An interface model was developed,incorporating the electrode potential and atom electronegativities.As a result,an insightful atomic scenario for the EDL structure under varied electrode potentials has been established,which unveils the important role of doping sites in regulating both the EDL structures and the following electrochemical reactions at the atomic level.Specifically,the negatively charged N atoms repel the anions and adsorb Li~+at high and low potentials,respectively.Such preferential adsorption suggests that Ndoped graphite can promote Li~+desolvation and regulate the location of Li~+deposition.This CPMD method not only unveils the mysterious function of N-doping from the viewpoint of EDL at the atomic level but also applies to probe the interfacial structure on other complicated electrodes.

Molecular Dynamics Numerical Simulation of Adsorption Characteristics and Exploitation Limits in Shale Oil Microscopic Pore Spaces

Microscopic pore structure in continental shale oil reservoirs is characterized by small pore throats and complex micro-structures.The adsorption behavior of hydrocarbons on the pore walls exhibits unique physical and chemical properties.Therefore,studying the adsorption morphology of hydrocarbon components in nanometer-sized pores and clarifying the exploitation limits of shale oil at the microscopic level are of great practical significance for the efficient development of continental shale oil.In this study,molecular dynamics simulations were employed to investigate the adsorption characteristics of various single-component shale oils in inorganic quartz fissures,and the influence of pore size and shale oil hydrocarbon composition on the adsorption properties in the pores was analyzed.The results show that different molecules have different adsorption capacities in shale oil pores,with lighter hydrocarbon components(C6H14)exhibiting stronger adsorption abilities.For the same adsorbed molecule,the adsorption amount linearly increases with the increase in pore diameter,but larger pores contribute more to shale oil adsorption.In shale pores,the thickness of the adsorption layer formed by shale oil molecules ranges from 0.4 to 0.5 nm,which is similar to the width of alkane molecules.Shale oil in the adsorbed state that is difficult to be exploited is mainly concentrated in the first adsorption layer.Among them,the volume fraction of adsorbed shale oil in 6 nm shale pores is 40.8%,while the volume fraction of shale oil that is difficult to be exploited is 16.2%.

Machine learning molecular dynamics simulations of liquid methanol

As the simplest hydrogen-bonded alcohol,liquid methanol has attracted intensive experimental and theoretical interest.However,theoretical investigations on this system have primarily relied on empirical intermolecular force fields or ab initio molecular dynamics with semilocal density functionals.Inspired by recent studies on bulk water using increasingly accurate machine learning force fields,we report a new machine learning force field for liquid methanol with a hybrid functional revPBE0 plus dispersion correction.Molecular dynamics simulations on this machine learning force field are orders of magnitude faster than ab initio molecular dynamics simulations,yielding the radial distribution functions,selfdiffusion coefficients,and hydrogen bond network properties with very small statistical errors.The resulting structural and dynamical properties are compared well with the experimental data,demonstrating the superior accuracy of this machine learning force field.This work represents a successful step toward a first-principles description of this benchmark system and showcases the general applicability of the machine learning force field in studying liquid systems.

Shear Deformation of DLC Based on Molecular Dynamics Simulation and Machine Learning

Shear deformation mechanisms of diamond-like carbon(DLC)are commonly unclear since its thickness of several micrometers limits the detailed analysis of its microstructural evolution and mechanical performance,which further influences the improvement of the friction and wear performance of DLC.This study aims to investigate this issue utilizing molecular dynamics simulation and machine learning(ML)techniques.It is indicated that the changes in the mechanical properties of DLC are mainly due to the expansion and reduction of sp3 networks,causing the stick-slip patterns in shear force.In addition,cluster analysis showed that the sp2-sp3 transitions arise in the stick stage,while the sp3-sp2 transitions occur in the slip stage.In order to analyze the mechanisms governing the bond breaking/re-formation in these transitions,the Random Forest(RF)model in ML identifies that the kinetic energies of sp3 atoms and their velocities along the loading direction have the highest influence.This is because high kinetic energies of atoms can exacerbate the instability of the bonding state and increase the probability of bond breaking/re-formation.Finally,the RF model finds that the shear force of DLC is highly correlated to its potential energy,with less correlation to its content of sp3 atoms.Since the changes in potential energy are caused by the variances in the content of sp3 atoms and localized strains,potential energy is an ideal parameter to evaluate the shear deformation of DLC.The results can enhance the understanding of the shear deformation of DLC and support the improvement of its frictional and wear performance.

Study on the Mechanism of Nanopatterning in Printed Electronics Based on Molecular Dynamics Simulation

In order to research the feasibility of using the selective adsorption principle to achieve automatic shaping of nano patterns,in this study,using the liquid gallium as the conductive ink and graphene as the printing plate surface,by changing the surface wettability of patterned areas on the nanoscale of graphene printed boards,the automatic formation of liquid gallium patterns on the graphene printed plate surface was simulated.The results indicated that liquid gallium can achieve automatic patterning on the surface of graphene patterned areas;the greater the interaction energy between gallium and carbon atoms,the clearer the pattern;gallium liquid is prone to remain in complex local positions of the pattern,making it difficult to shape the pattern;if the spacing between adjacent pattern lines is too large or too small,it will result in residual gallium liquid between the lines;increasing the thickness of the gallium film will cause the pattern to expand beyond the boundary,but increasing the thickness of the gallium film can also enhance the thickness and uniformity of the pattern lines.In summary,the principle of selective adsorption can be used to achieve the automatic formation of nano patterns,and the pattern formation effect is influenced by factors such as atomic interaction energy and pattern configuration.

Exploring the molecular mechanism of action of curcumin for the treatment of diabetic retinopathy,using network pharmacology,molecular docking,and molecular dynamics simulation

Background:Based on network pharmacology and molecular docking,the present study investigated the mechanism of curcumin(CUR)in diabetic retinopathy treatment.Methods:Based on the DisGeNET,Swiss TargetPrediction,GeneCards,Online Mendelian Inheritance in Man,Gene Expression Omnibus,and Comparative Toxicogenomics Database,the intersection core targets of CUR and diabetic retinopathy were identified.The intersection target was imported into the STRING database to obtain the protein-protein interaction map.According to the Database for Annotation,Visualization and Integrated Discovery database,the intersected targets were enriched in Gene Ontology(GO)and Kyoto Encyclopedia of Genes and Genomes pathways.Then Cytoscape 3.9.1 is used to make the drug-target-disease-pathway network.The mechanism of CUR and diabetic retinopathy was further verified by molecular docking and molecular dynamics simulation.Results:There were 203 intersecting targets of CUR and diabetic retinopathy identified.1320 GO entries were enriched for GO functions,which were primarily involved in the composition of cells such as identical protein binding,protein binding,enzyme binding,etc.It was found that 175 pathways were enriched using Kyoto Encyclopedia of Genes and Genomes pathway enrichment methods,which were mainly included in the lipid and atherosclerosis,AGE-RAGE signaling pathway in diabetic complications,pathways in cancer,etc.In the molecular docking analysis,CUR was found to have a good ability to bind to the core targets of albumin,IL-1B,and IL-6.The binding of albumin to CUR was further verified by molecular dynamics simulation.Conclusion:As a result of this study,CUR may exert a role in the treatment of diabetic retinopathy through multi-target and multi-pathway regulation,which indicates a possible direction of future research.

Molecular dynamics simulation of relationship between local structure and dynamics during glass transition of Mg_7Zn_3 alloy

The rapid solidification process of Mg7Zn3 alloy was simulated by the molecular dynamics method. The relationship between the local structure and the dynamics during the liquid-glass transition was deeply investigated. It was found that the Mg-centered FK polyhedron and the Zn-centered icosahedron play a critical role in the formation of Mg7Zn3 metallic glass. The self-diffusion coefficients of Mg and Zn atoms deviate from the Arrhenius law near the melting temperature and then satisfy the power law. According to the time correlation functions of mean-square displacement, incoherent intermediate scattering function and non-Gaussian parameter, it was found that the β-relaxation in Mg7Zn3 supercooled liquid becomes more and more evident with decreasing temperature, and the α-relaxation time rapidly increases in the VFT law. Moreover, the smaller Zn atom has a faster relaxation behavior than the Mg atom. Some local atomic structures with short-range order have lower mobility, and they play a critical role in the appearance of cage effect in theβ-relaxation regime. The dynamics deviates from the Arrhenius law just at the temperature as the number of local atomic structures begins to rapidly increase. The dynamic glass transition temperature （Tc） is close to the glass transition point in structure （TgStr）.

Local structure of calcium silicate melts from classical molecular dynamics simulation and a newly constructed thermodynamic model

The distributions of local structural units of calcium silicate melts were quantified by means of classical molecular dynamics simulation and a newly constructed structural thermodynamic model. The distribution of five kinds of Si-O tetrahedra Qi from these two methods was compared with each other and also with the experimental Raman spectra, an excellent agreement was achieved. These not only displayed the panorama distribution of microstructural units in the whole composition range, but also proved that the thermodynamic model is suitable for the utilization as the subsequent application model of spectral experiments for the thermodynamic calculation. Meanwhile, the five refined regions mastered by different disproportionating reactions were obtained. Finally, the distributions of two kinds of connections between Qi were obtained, denoted as Qi-Ca-Qj and Qi-[Ob]-Qj, from the thermodynamic model, and a theoretical verification was given that the dominant connections for any composition are equivalent connections.

Pyrolysis of CL20-BTF Co-crystal via ReaxFF-lg Reactive Force Field Molecular Dynamics Simulations

To obtain detailed information on the potential energy, the evolution of species, the initial reaction paths, and thermal decomposition products, we conducted simulations on pyrolysis process of CL20/BTF co-crystal using the ReaxFF/lg reaction force field, with temperature set at 2000 K to 3000 K. With the analysis of evolution curves of potential energy based on exponential function, we obtain the overall characteristic time. Via a description of the total package reaction with classical Arrhenius law, we obtain the activation energy of CL20/BTF co-crystal： Ea=60.8 kcal/mol. Based on the initial path of CL20/BTF co-crystal thermal decomposition we studied, we conclude that N-NO2 bond of CL20 molecules breaks first, working as a dominant role in the initial stage of thermal decomposition under the condition of different temperatures, and that all CL20 molecules completely decompose before BTF molecular regardless of different temperatures. We also find that the main products of CL20/BTF co-crystal are NO2, NO, NO3, HNO, O2, N2, H2O, CO2, N2O, and HONO, etc., on which the temperature forms certain influence.

Simulating the Dynamics of Bimetallic Clusters Deposited onto a Surface Using Molecular Dynamics

This paper examines the interface development between a single crystalline Ag matrix and core-shell AgnCom nanoclusters that have been deposited with energies varying between 0.25 eV and 1.5 eV per atom using computer modeling techniques. Clusters undergo deformation as a result of the slowing down;they may also become epitaxial with the substrate and maintain their core-shell structure. A detailed analysis of the effects of the cluster-surface interaction is conducted over a realistic size and energy range, and a model is created to show how clusters accumulate. It is discovered that both the silver shells and the cobalt cluster cores exhibit limited epitaxy with the substrate, and that the contact produced is only a few atomic layers thick. The effect is higher for Ag shells than for Co cores, and it is not very energy dependent.

Anisotropy of Thermal-expansion for β-Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine： Quantum Chemistry Calculation and Molecular Dynamics Simulation

Molecular dynamics simulations on octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine （HMX） at 303-383 K and atmospheric pressure are carried out under NPT ensemble and COMPASS force field, the equilibrium structures at elevated temperatures were obtained and showed that the stacking style of molecules don＇t change. The coefficient of thermal expansion （CTE） values were calculated by linear fitting method. The results show that the CTE values are close to the experimental results and show anisotropy. The total energies of HMX cells with separately increasing expansion rates （100%-105%） along each crystallographic axis was calculated by periodic density functional theory method, the results of the energy change rates are anisotropic, and the correlation equations of energy change-CTE values are established. Thus the hypostasis of the anisotropy of HMX crystal＇s thermal expansion, the determinate molecular packing style, is elucidated.