A nonlinear dynamic simulation model based on coordinated control of speed and flow rate for the molten salt reactor and combined cycle systems is proposed here to ensure the coordination and stability between the mol...A nonlinear dynamic simulation model based on coordinated control of speed and flow rate for the molten salt reactor and combined cycle systems is proposed here to ensure the coordination and stability between the molten salt reactor and power system.This model considers the impact of thermal properties of fluid variation on accuracy and has been validated with Simulink.This study reveals the capability of the control system to compensate for anomalous situations and maintain shaft stability in the event of perturbations occurring in high-temperature molten salt tank outlet parameters.Meanwhile,the control system’s impact on the system’s dynamic characteristics under molten salt disturbance is also analyzed.The results reveal that after the disturbance occurs,the controlled system benefits from the action of the control,and the overshoot and disturbance amplitude are positively correlated,while the system power and frequency eventually return to the initial values.This simulation model provides a basis for utilizing molten salt reactors for power generation and maintaining grid stability.展开更多
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.展开更多
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 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.展开更多
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.展开更多
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.展开更多
The C–H bond activation in alkane dehydrogenation reactions is a key step in determining the reaction rate.To understand the impact of entropy,we performed ab initio static and molecular dynamics free energy simulati...The C–H bond activation in alkane dehydrogenation reactions is a key step in determining the reaction rate.To understand the impact of entropy,we performed ab initio static and molecular dynamics free energy simulations of ethane dehydrogenation over Co@BEA zeolite at different temperatures.AIMD simulations showed that a sharp decrease in free energy barrier as temperature increased.Our analysis of the temperature dependence of activation free energies uncovered an unusual entropic effect accompanying the reaction.The unique spatial structures around the Co active site at different temperatures influenced both the extent of charge transfer in the transition state and the arrangement of 3d orbital energy levels.We provided explanations consistent with the principles of thermodynamics and statistical physics.The insights gained at the atomic level have offered a fresh interpretation of the intricate long-range interplay between local chemical reactions and extensive chemical environments.展开更多
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...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.展开更多
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...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%.展开更多
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,GeneCa...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.展开更多
The launch dynamics theory for multibody systems emerges as an innovative and efficacious approach for the study of launch dynamics,capable of addressing the challenges of complex modeling,diminished computational eff...The launch dynamics theory for multibody systems emerges as an innovative and efficacious approach for the study of launch dynamics,capable of addressing the challenges of complex modeling,diminished computational efficiency,and imprecise analyses of system dynamic responses found in the dynamics research of intricate multi-rigid-flexible body systems,such as self-propelled artillery.This advancement aims to enhance the firing accuracy and launch safety of self-propelled artillery.Recognizing the shortfall of overlooking the band engraving process in existing theories,this study introduces a novel coupling calculation methodology for the launch dynamics of a self-propelled artillery multibody system.This method leverages the ABAQUS subroutine interface VUAMP to compute the dynamic response of the projectile and barrel during the launch process of large-caliber self-propelled artillery.Additionally,it examines the changes in projectile resistance and band deformation in relation to projectile motion throughout the band engraving process.Comparative analysis of the computational outcomes with experimental data evidences that the proposed method offers a more precise depiction of the launch process of self-propelled artillery,thereby enhancing the accuracy of launch dynamics calculations for self-propelled artillery.展开更多
The advantages of a flat-panel X-ray source(FPXS)make it a promising candidate for imaging applications.Accurate imaging-system modeling and projection simulation are critical for analyzing imaging performance and res...The advantages of a flat-panel X-ray source(FPXS)make it a promising candidate for imaging applications.Accurate imaging-system modeling and projection simulation are critical for analyzing imaging performance and resolving overlapping projection issues in FPXS.The conventional analytical ray-tracing approach is limited by the number of patterns and is not applicable to FPXS-projection calculations.However,the computation time of Monte Carlo(MC)simulation is independent of the size of the patterned arrays in FPXS.This study proposes two high-efficiency MC projection simulators for FPXS:a graphics processing unit(GPU)-based phase-space sampling MC(gPSMC)simulator and GPU-based fluence sampling MC(gFSMC)simulator.The two simulators comprise three components:imaging-system modeling,photon initialization,and physical-interaction simulations in the phantom.Imaging-system modeling was performed by modeling the FPXS,imaging geometry,and detector.The gPSMC simulator samples the initial photons from the phase space,whereas the gFSMC simulator performs photon initialization from the calculated energy spectrum and fluence map.The entire process of photon interaction with the geometry and arrival at the detector was simulated in parallel using multiple GPU kernels,and projections based on the two simulators were calculated.The accuracies of the two simulators were evaluated by comparing them with the conventional analytical ray-tracing approach and acquired projections,and the efficiencies were evaluated by comparing the computation time.The results of simulated and realistic experiments illustrate the accuracy and efficiency of the proposed gPSMC and gFSMC simulators in the projection calculation of various phantoms.展开更多
We develop a numerical method for the time evolution of Gaussian wave packets on flat-band lattices in the presence of correlated disorder.To achieve this,we introduce a method to generate random on-site energies with...We develop a numerical method for the time evolution of Gaussian wave packets on flat-band lattices in the presence of correlated disorder.To achieve this,we introduce a method to generate random on-site energies with prescribed correlations.We verify this method with a one-dimensional(1D)cross-stitch model,and find good agreement with analytical results obtained from the disorder-dressed evolution equations.This allows us to reproduce previous findings,that disorder can mobilize 1D flat-band states which would otherwise remain localized.As explained by the corresponding disorder-dressed evolution equations,such mobilization requires an asymmetric disorder-induced coupling to dispersive bands,a condition that is generically not fulfilled when the flat-band is resonant with the dispersive bands at a Dirac point-like crossing.We exemplify this with the 1D Lieb lattice.While analytical expressions are not available for the two-dimensional(2D)system due to its complexity,we extend the numerical method to the 2D a–T3 model,and find that the initial flat-band wave packet preserves its localization when a=0,regardless of disorder and intersections.However,when a̸=0,the wave packet shifts in real space.We interpret this as a Berry phase controlled,disorder-induced wave-packet mobilization.In addition,we present density functional theory calculations of candidate materials,specifically Hg1−xCdxTe.The flat-band emerges near the G point(α=0)in the Brillouin zone.展开更多
We investigate the Floquet spectrum and excitation properties of a two-ultracold-atom system with periodically driven interaction in a three-dimensional harmonic trap.The interaction between the atoms is changed by va...We investigate the Floquet spectrum and excitation properties of a two-ultracold-atom system with periodically driven interaction in a three-dimensional harmonic trap.The interaction between the atoms is changed by varying the s-wave scattering length in two ways,the cosine and the square-wave modulations.It is found that as the driving frequency increases,the Floquet spectrum exhibits two main features for both modulations,the accumulating and the spreading of the quasienergy levels,which further lead to different dynamical behaviors.The accumulation is associated with collective excitations and the persistent growth of the energy,while the spread indicates that the energy is bounded at all times.The initial scattering length,the driving frequency and amplitude can all significantly change the Floquet spectrum as well as the dynamics.However,the corresponding relation between them is valid universally.Finally,we propose a mechanism for selectively exciting the system to one specific state by using the avoided crossing of two quasienergy levels,which could guide preparation of a desired state in experiments.展开更多
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 mecha...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.展开更多
Evolution and interaction of plane waves of the multidimensional zero-pressure gas dynamics system leads to the study of the corresponding one dimensional system.In this paper,we study the initial value problem for on...Evolution and interaction of plane waves of the multidimensional zero-pressure gas dynamics system leads to the study of the corresponding one dimensional system.In this paper,we study the initial value problem for one dimensional zero-pressure gas dynamics system.Here the first equation is the Burgers equation and the second one is the continuity equation.We consider the solution with initial data in the space of bounded Borel measures.First we prove a general existence result in the algebra of generalized functions of Colombeau.Then we study in detail special solutions withδ-measures as initial data.We study interaction of waves originating from initial data concentrated on two point sources and interaction with classical shock/rarefaction waves.This gives an understanding of plane-wave interactions in the multidimensional case.We use the vanishing viscosity method in our analysis as this gives the physical solution.展开更多
In an effort to clarify the formation mechanism of LPSO structure in Mg-Y-Zn alloy,the chemical environment and structural ordering in liquid Mg-rich Mg-Y-Zn system are investigated with the aid of ab-initio molecular...In an effort to clarify the formation mechanism of LPSO structure in Mg-Y-Zn alloy,the chemical environment and structural ordering in liquid Mg-rich Mg-Y-Zn system are investigated with the aid of ab-initio molecular dynamics simulation.In liquid Mg-rich Mg-Y alloys,the strong Mg-Y interaction is determined,which promotes the formation of fivefold symmetric local structure.For Mg-Zn alloys,the weak Mg-Zn interaction results in the fivefold symmetry weakening in the liquid structure.Due to the coexistence of Y and Zn,the strong attractive interaction is introduced in liquid Mg-Y-Zn ternary alloy,and contributes to the clustering of Mg,Y,Zn launched from Zn.What is more,the distribution of local structures becomes closer to that in pure Mg compared with that in binary Mg-Y and Mg-Zn alloys.These results should relate to the origins of the Y/Zn segregation zone and close-packed stacking mode in LPSO structure,which provides a new insight into the formation mechanism of LPSO structure at atomic level.展开更多
Grain boundaries(GBs)play a crucial role on the structural stability and mechanical properties of Cu and its alloys.In this work,molecular dynamics(MD)simulations are employed to study the effects of Fe solutes on the...Grain boundaries(GBs)play a crucial role on the structural stability and mechanical properties of Cu and its alloys.In this work,molecular dynamics(MD)simulations are employed to study the effects of Fe solutes on the formation energy,excess volume,dislocations and melting behaviors of GBs in CuFe alloys.It is illustrated that Fe solute affects the structural stability of Cu GBs substantially,the formation energy of GBs is reduced,but the thickness and melting point of GBs are increased,that is,the structural stability of Cu GBs is significantly improved owing to the Fe solutes.A strong scaling law exists between the formation energy,excess volume,thickness and melting point of GBs.Therefore,Fe solid solute plays an important role in the characteristics of GBs in bi-crystal Cu.展开更多
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,...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.展开更多
基金This work was supported by the Chinese TMSR Strategic Pioneer Science and Technology Project(No.XDA02010300).
文摘A nonlinear dynamic simulation model based on coordinated control of speed and flow rate for the molten salt reactor and combined cycle systems is proposed here to ensure the coordination and stability between the molten salt reactor and power system.This model considers the impact of thermal properties of fluid variation on accuracy and has been validated with Simulink.This study reveals the capability of the control system to compensate for anomalous situations and maintain shaft stability in the event of perturbations occurring in high-temperature molten salt tank outlet parameters.Meanwhile,the control system’s impact on the system’s dynamic characteristics under molten salt disturbance is also analyzed.The results reveal that after the disturbance occurs,the controlled system benefits from the action of the control,and the overshoot and disturbance amplitude are positively correlated,while the system power and frequency eventually return to the initial values.This simulation model provides a basis for utilizing molten salt reactors for power generation and maintaining grid stability.
基金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(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.
文摘In this study,circular dichroism(CD)and molecular dynamics(MD)simulation were used to investigate the thermal unfolding pathway of staphylococcal enterotoxin B(SEB)at temperatures of 298–371 and 298–500 K,and the relationship between the experimental and simulation results were explored.Our computational findings on the secondary structure of SEB showed that at room temperature,the CD spectroscopic results were highly consistent with the MD results.Moreover,under heating conditions,the changing trends of helix,sheet and random coil obtained by CD spectral fitting were highly consistent with those obtained by MD.In order to gain a deeper understanding of the thermal stability mechanism of SEB,the MD trajectories were analyzed in terms of root mean square deviation(RMSD),secondary structure assignment(SSA),radius of gyration(R_(g)),free energy surfaces(FES),solvent-accessible surface area(SASA),hydrogen bonds and salt bridges.The results showed that at low heating temperature,domain Ⅰ without loops(omitting the mobile loop region)mainly relied on hydrophobic interaction to maintain its thermal stability,whereas the thermal stability of domain Ⅱ was mainly controlled by salt bridges and hydrogen bonds.Under high heating temperature conditions,the hydrophobic interactions in domain Ⅰ without loops were destroyed and the secondary structure was almost completely lost,while domain Ⅱ could still rely on salt bridges as molecular staples to barely maintain the stability of the secondary structure.These results help us to understand the thermodynamic and kinetic mechanisms that maintain the thermal stability of SEB at the molecular level,and provide a direction for establishing safer and more effective food sterilization processes.
基金supported by the National Natural Science Foundation of China(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.
基金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.
文摘The C–H bond activation in alkane dehydrogenation reactions is a key step in determining the reaction rate.To understand the impact of entropy,we performed ab initio static and molecular dynamics free energy simulations of ethane dehydrogenation over Co@BEA zeolite at different temperatures.AIMD simulations showed that a sharp decrease in free energy barrier as temperature increased.Our analysis of the temperature dependence of activation free energies uncovered an unusual entropic effect accompanying the reaction.The unique spatial structures around the Co active site at different temperatures influenced both the extent of charge transfer in the transition state and the arrangement of 3d orbital energy levels.We provided explanations consistent with the principles of thermodynamics and statistical physics.The insights gained at the atomic level have offered a fresh interpretation of the intricate long-range interplay between local chemical reactions and extensive chemical environments.
基金supported by the National Natural Science Foundation of China(Grant Nos.52035009 and 51761135106)the State Key Laboratory of Precision Measuring Technology and Instruments(Pilt1705)+1 种基金the Henan Key Laboratory of Intelligent Manufacturing Equipment Integration for Superhard Materials(JDKJ2022-01)the“111”project by the State Administration of Foreign Experts Affairs and the Ministry of Education of China(Grant No.B07014).
文摘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.
文摘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%.
基金supported by the Hubei Province Research Innovation Team Project(T2021022)Scientific Research Projects of Hubei Health Commission(WJ2023M119).
文摘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.
基金supported by the National Natural Science Foundation of China (Grant Number:12372093)。
文摘The launch dynamics theory for multibody systems emerges as an innovative and efficacious approach for the study of launch dynamics,capable of addressing the challenges of complex modeling,diminished computational efficiency,and imprecise analyses of system dynamic responses found in the dynamics research of intricate multi-rigid-flexible body systems,such as self-propelled artillery.This advancement aims to enhance the firing accuracy and launch safety of self-propelled artillery.Recognizing the shortfall of overlooking the band engraving process in existing theories,this study introduces a novel coupling calculation methodology for the launch dynamics of a self-propelled artillery multibody system.This method leverages the ABAQUS subroutine interface VUAMP to compute the dynamic response of the projectile and barrel during the launch process of large-caliber self-propelled artillery.Additionally,it examines the changes in projectile resistance and band deformation in relation to projectile motion throughout the band engraving process.Comparative analysis of the computational outcomes with experimental data evidences that the proposed method offers a more precise depiction of the launch process of self-propelled artillery,thereby enhancing the accuracy of launch dynamics calculations for self-propelled artillery.
文摘The advantages of a flat-panel X-ray source(FPXS)make it a promising candidate for imaging applications.Accurate imaging-system modeling and projection simulation are critical for analyzing imaging performance and resolving overlapping projection issues in FPXS.The conventional analytical ray-tracing approach is limited by the number of patterns and is not applicable to FPXS-projection calculations.However,the computation time of Monte Carlo(MC)simulation is independent of the size of the patterned arrays in FPXS.This study proposes two high-efficiency MC projection simulators for FPXS:a graphics processing unit(GPU)-based phase-space sampling MC(gPSMC)simulator and GPU-based fluence sampling MC(gFSMC)simulator.The two simulators comprise three components:imaging-system modeling,photon initialization,and physical-interaction simulations in the phantom.Imaging-system modeling was performed by modeling the FPXS,imaging geometry,and detector.The gPSMC simulator samples the initial photons from the phase space,whereas the gFSMC simulator performs photon initialization from the calculated energy spectrum and fluence map.The entire process of photon interaction with the geometry and arrival at the detector was simulated in parallel using multiple GPU kernels,and projections based on the two simulators were calculated.The accuracies of the two simulators were evaluated by comparing them with the conventional analytical ray-tracing approach and acquired projections,and the efficiencies were evaluated by comparing the computation time.The results of simulated and realistic experiments illustrate the accuracy and efficiency of the proposed gPSMC and gFSMC simulators in the projection calculation of various phantoms.
基金the National Natural Sci-ence Foundation of China(Grant No.61988102)the Key Research and Development Program of Guangdong Province(Grant No.2019B090917007)+5 种基金the Science and Technology Planning Project of Guangdong Province(Grant No.2019B090909011)Q.L.acknowledges Guangzhou Basic and Applied Basic Research Project(Grant No.2023A04J0018)Z.L.acknowledges the support of fund-ing from Chinese Academy of Sciences E1Z1D10200 and E2Z2D10200from ZJ project 2021QN02X159 and from JSPS(Grant Nos.PE14052 and P16027)We gratefully ac-knowledge HZWTECH for providing computation facilities.Z.-X.H.was supported by the National Natural Science Foun-dation of China(Grant Nos.11974064 and 12147102)the Fundamental Research Funds for the Central Universities(Grant No.2020CDJQY-Z003).
文摘We develop a numerical method for the time evolution of Gaussian wave packets on flat-band lattices in the presence of correlated disorder.To achieve this,we introduce a method to generate random on-site energies with prescribed correlations.We verify this method with a one-dimensional(1D)cross-stitch model,and find good agreement with analytical results obtained from the disorder-dressed evolution equations.This allows us to reproduce previous findings,that disorder can mobilize 1D flat-band states which would otherwise remain localized.As explained by the corresponding disorder-dressed evolution equations,such mobilization requires an asymmetric disorder-induced coupling to dispersive bands,a condition that is generically not fulfilled when the flat-band is resonant with the dispersive bands at a Dirac point-like crossing.We exemplify this with the 1D Lieb lattice.While analytical expressions are not available for the two-dimensional(2D)system due to its complexity,we extend the numerical method to the 2D a–T3 model,and find that the initial flat-band wave packet preserves its localization when a=0,regardless of disorder and intersections.However,when a̸=0,the wave packet shifts in real space.We interpret this as a Berry phase controlled,disorder-induced wave-packet mobilization.In addition,we present density functional theory calculations of candidate materials,specifically Hg1−xCdxTe.The flat-band emerges near the G point(α=0)in the Brillouin zone.
基金supported by the National Natural Science Foundation of China(Grant No.12004049)the Fund of State Key Laboratory of IPOC(BUPT)(Grant Nos.600119525 and 505019124).
文摘We investigate the Floquet spectrum and excitation properties of a two-ultracold-atom system with periodically driven interaction in a three-dimensional harmonic trap.The interaction between the atoms is changed by varying the s-wave scattering length in two ways,the cosine and the square-wave modulations.It is found that as the driving frequency increases,the Floquet spectrum exhibits two main features for both modulations,the accumulating and the spreading of the quasienergy levels,which further lead to different dynamical behaviors.The accumulation is associated with collective excitations and the persistent growth of the energy,while the spread indicates that the energy is bounded at all times.The initial scattering length,the driving frequency and amplitude can all significantly change the Floquet spectrum as well as the dynamics.However,the corresponding relation between them is valid universally.Finally,we propose a mechanism for selectively exciting the system to one specific state by using the avoided crossing of two quasienergy levels,which could guide preparation of a desired state in experiments.
基金financially supported by the Project funded by China Postdoctoral Science Foundation (NO.2022M723500)the National Natural Science Foundation of China (NO.52204069)the Sinopec Science and Technology Project of China (NO.P22015)。
文摘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.
基金supported by the TIFR-CAM Doctoral Fellowshipthe NISER Postdoctoral Fellowship (through the project “Basic research in physics and multidisciplinary sciences” with identification # RIN4001) during the preparation of this papersupported by the Raja Ramanna Fellowship
文摘Evolution and interaction of plane waves of the multidimensional zero-pressure gas dynamics system leads to the study of the corresponding one dimensional system.In this paper,we study the initial value problem for one dimensional zero-pressure gas dynamics system.Here the first equation is the Burgers equation and the second one is the continuity equation.We consider the solution with initial data in the space of bounded Borel measures.First we prove a general existence result in the algebra of generalized functions of Colombeau.Then we study in detail special solutions withδ-measures as initial data.We study interaction of waves originating from initial data concentrated on two point sources and interaction with classical shock/rarefaction waves.This gives an understanding of plane-wave interactions in the multidimensional case.We use the vanishing viscosity method in our analysis as this gives the physical solution.
基金supported by National Natural Science Foundation of China,China(No.51901117,51801116)Youth Innovation and Technology Support Program of Shandong Provincial Colleges and Universities,China(No.2020KJA002)+2 种基金Youth Fund of Shandong Academy of Sciences,China(2020QN0021)Innovation Pilot Project for Fusion of Science,Education and Industry(International Cooperation)from Qilu University of Technology(Shandong Academy of Sciences),China(No.2020KJC-GH03)Several Policies on Promoting Collaborative Innovation and Industrialization of Achievements in Universities and Research Institutes,China(No.2019GXRC030)。
文摘In an effort to clarify the formation mechanism of LPSO structure in Mg-Y-Zn alloy,the chemical environment and structural ordering in liquid Mg-rich Mg-Y-Zn system are investigated with the aid of ab-initio molecular dynamics simulation.In liquid Mg-rich Mg-Y alloys,the strong Mg-Y interaction is determined,which promotes the formation of fivefold symmetric local structure.For Mg-Zn alloys,the weak Mg-Zn interaction results in the fivefold symmetry weakening in the liquid structure.Due to the coexistence of Y and Zn,the strong attractive interaction is introduced in liquid Mg-Y-Zn ternary alloy,and contributes to the clustering of Mg,Y,Zn launched from Zn.What is more,the distribution of local structures becomes closer to that in pure Mg compared with that in binary Mg-Y and Mg-Zn alloys.These results should relate to the origins of the Y/Zn segregation zone and close-packed stacking mode in LPSO structure,which provides a new insight into the formation mechanism of LPSO structure at atomic level.
基金supported by National Key Research and Development Program of China(No.2021YFB3400800)National Natural Science Foundation of China(Grant No.52271136,51901177)Natural Science Foundation of Shaanxi Province(No.2021JC-06,2019TD-020).
文摘Grain boundaries(GBs)play a crucial role on the structural stability and mechanical properties of Cu and its alloys.In this work,molecular dynamics(MD)simulations are employed to study the effects of Fe solutes on the formation energy,excess volume,dislocations and melting behaviors of GBs in CuFe alloys.It is illustrated that Fe solute affects the structural stability of Cu GBs substantially,the formation energy of GBs is reduced,but the thickness and melting point of GBs are increased,that is,the structural stability of Cu GBs is significantly improved owing to the Fe solutes.A strong scaling law exists between the formation energy,excess volume,thickness and melting point of GBs.Therefore,Fe solid solute plays an important role in the characteristics of GBs in bi-crystal Cu.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.12222506,12347102,and 12174184).
文摘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.