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.

Ab initio molecular dynamics simulation reveals the influence of entropy effect on Co@BEA zeolite-catalyzed dehydrogenation of ethane

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.

A Self-Healing and Nonflammable Cross-Linked Network Polymer Electrolyte with the Combination of Hydrogen Bonds and Dynamic Disulfide Bonds for Lithium Metal Batteries

The self-healing solid polymer electrolytes(SHSPEs)can spontaneously eliminate mechanical damages or micro-cracks generated during the assembly or operation of lithium-ion batteries(LIBs),significantly improving cycling performance and extending service life of LIBs.Here,we report a novel cross-linked network SHSPE(PDDP)containing hydrogen bonds and dynamic disulfide bonds with excellent self-healing properties and nonflammability.The combination of hydrogen bonding between urea groups and the metathesis reaction of dynamic disulfide bonds endows PDDP with rapid self-healing capacity at 28°C without external stimulation.Furthermore,the addition of 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide(EMIMTFSI)improves the ionic conductivity(1.13×10^(−4)S cm^(−1)at 28°C)and non-flammability of PDDP.The assembled Li/PDDP/LiFePO_(4)cell exhibits excellent cycling performance with a discharge capacity of 137 mA h g^(−1)after 300 cycles at 0.2 C.More importantly,the self-healed PDDP can recover almost the same ionic conductivity and cycling performance as the original PDDP.

The bond graph model of planar flexible multibody mechanical systems and its dynamic principle

In order to increase the efficiency and reliability of the dynamic analysis for flexible planar linkage containing the coupling of multi-energy domains, a method based on bond graph is introduced. From the viewpoint of power conservation, the peculiar property of bond graph multiport element MTF is discussed. The procedure of modeling planar flexible muhibody mechanical systems by bond graphs and its dynamic principle are deseribed. To overcome the algebraic difficulty brought by differential causality anti nonlinear junction structure, the constraint forces at joints can be considered as unknown effort sources and added to the corresponding O-junctions of system bond graph model. As a result, the automatic modeling on a computer is realized. The validity of the procedure is illustrated by a practical example.

A New Three-dimensional Supramolecular Polymer Built from Non-covalent Bonding Interactions

A new complex, [Ni2（L）4（H2O）8]（1, L1 = 4-（1H-imidazol-4-yl）benzoic acid）, has been hydrothermally prepared and characterized by single-crystal X-ray diffraction, IR spectroscopy, elemental analysis and PXRD. Complex 1 crystallizes in monoclinic, space group P21/c with α = 22.281（2）, b = 7.3959（7）, c = 24.978（3） ？, β = 90.876（10）, V = 4115.6（7） ？3, Z = 8, C20H22N4O8Ni, Mr = 505.13, Dc = 1.630 g/cm3, μ = 1.001 mm-1, S = 1.080, F（000） = 2096, the final R = 0.452 and wR = 0.1152 for 9380 observed reflections （I 〉 2σ（I））. The result of X-ray diffraction analysis revealed three different kinds of Ni（II） centre mononuclear molecules in the asymmetric unit. The independent mononuclear units are bridged to form a three-dimensional supramolecular polymer by extensive hydrogen bonds and C–H… non-covalent bonding interactions.

Mechanism of microweld formation and breakage during Cu-Cu wire bonding investigated by molecular dynamics simulation

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.

Dissipative Particle Dynamics Simulation on Bonding Reaction Between Surface Modified Nanoparticles

A simulation study was carried out by using dissipative particle dynamics(DPD) method to explore the effects of properties of coating chains, such as length, density, rigidity of polymer chains, as well as the distance between nanoparticles on bonding reaction of coating chains grafted onto nanoparticles. The results show that bonding ratios of coated chains strongly depend on the length and density of coating chains. For nanoparticles with different coating densities, the optimum chain length for bonding reaction are varied. The rigidity of coating chains exhibits vigorous effects on bonding reaction that highly depends on chain lengths. DPD simulation can be used to study the bonding reaction between coated nanoparticles, which may help experimental synthesis of nanocomposites with excellent properties.

The dynamic stress intensity factor analysis of adhesively bonded material interface crack with damage under shear loading

This paper studies the dynamic stress intensity factor （DSIF） at the interface in an adhesive joint under shear loading. Material damage is considered. By introducing the dislocation density function and using the integral transform, the problem is reduced to algebraic equations and can be solved with the collocation dots method in the Laplace domain. Time response of DSIF is calculated with the inverse Laplace integral transform. The results show that the mode Ⅱ DSIF increases with the shear relaxation parameter, shear module and Poisson ratio, while decreases with the swell relaxation parameter. Damage shielding only occurs at the initial stage of crack propagation. The singular index of crack tip is -0.5 and independent on the material parameters, damage conditions of materials, and time. The oscillatory index is controlled by viscoelastic material parameters.

Effect of interlayer bonded bilayer graphene on friction

We study the friction properties of interlayer bonded bilayer graphene by simulating the movement of a slider on the surface of bilayer graphene using molecular dynamics.The results show that the presence of the interlayer covalent bonds due to the local sp^(3) hybridization of carbon atoms in the bilayer graphene seriously reduces the frictional coefficient of the bilayer graphene surface to 30%,depending on the coverage of interlayer sp^(3) bonds and normal loads.For a certain coverage of interlayer sp3bonds,when the normal load of the slider reaches a certain value,the surface of this interlayer bonded bilayer graphene will lose the friction reduction effect on the slider.Our findings provide guidance for the regulation and manipulation of the frictional properties of bilayer graphene surfaces through interlayer covalent bonds,which may be useful for applications of friction related graphene based nanodevices.

Machine-Learning-Assisted Design of Deep Eutectic Solvents Based on Uncovered Hydrogen Bond Patterns

Non-ionic deep eutectic solvents(DESs)are non-ionic designer solvents with various applications in catalysis,extraction,carbon capture,and pharmaceuticals.However,discovering new DES candidates is challenging due to a lack of efficient tools that accurately predict DES formation.The search for DES relies heavily on intuition or trial-and-error processes,leading to low success rates or missed opportunities.Recognizing that hydrogen bonds(HBs)play a central role in DES formation,we aim to identify HB features that distinguish DES from non-DES systems and use them to develop machine learning(ML)models to discover new DES systems.We first analyze the HB properties of 38 known DES and 111 known non-DES systems using their molecular dynamics(MD)simulation trajectories.The analysis reveals that DES systems have two unique features compared to non-DES systems:The DESs have①more imbalance between the numbers of the two intra-component HBs and②more and stronger inter-component HBs.Based on these results,we develop 30 ML models using ten algorithms and three types of HB-based descriptors.The model performance is first benchmarked using the average and minimal receiver operating characteristic(ROC)-area under the curve(AUC)values.We also analyze the importance of individual features in the models,and the results are consistent with the simulation-based statistical analysis.Finally,we validate the models using the experimental data of 34 systems.The extra trees forest model outperforms the other models in the validation,with an ROC-AUC of 0.88.Our work illustrates the importance of HBs in DES formation and shows the potential of ML in discovering new DESs.

Interfacial coordination bonds accelerate charge separation for unprecedented hydrogen evolution over S-scheme heterojunction

Inspired by natural photosynthesis,fabricating high-performance S-scheme heterojunction is regarded as a successful tactic to address energy and environmental issues.Herein,NH_(2)-MIL-125(Ti)/Zn_(0.5)Cd_(0.5)S/NiS(NMT/ZCS/NiS)S-scheme heterojunction with interfacial coordination bonds is successfully synthesized through in-situ solvothermal strategy.Notably,the optimal NMT/ZCS/NiS S-scheme heterojunction exhibits comparable photocatalytic H_(2)evolution(PHE)rate of about 14876.7μmol h^(−1)g^(−1)with apparent quantum yield of 24.2%at 420 nm,which is significantly higher than that of recently reported MOFs-based photocatalysts.The interfacial coordination bonds(Zn–N,Cd–N,and Ni–N bonds)accelerate the separation and transfer of photogenerated charges,and the NiS as cocatalyst can provide more catalytically active sites,which synergistically improve the photocatalytic performance.Moreover,theoretical calculation results display that the construction of NMT/ZCS/NiS S-scheme heterojunction also optimize the binding energy of active site-adsorbed hydrogen atoms to enable fast adsorption and desorption.Photoassisted Kelvin probe force microscopy,in-situ irradiation X-ray photoelectron spectroscopy,femtosecond transient absorption spectroscopy,and theoretical calculations provide sufficient evidence of the S-scheme charge migration mechanism.This work offers unique viewpoints for simultaneously accelerating the charge dynamics and optimizing the binding strength between the active sites and hydrogen adsorbates over S-scheme heterojunction.

VECTOR BOND GRAPH REPRESENTATION OF FINITE ELEMENT METHOD IN STRUCTURAL DYNAMICS

In this paper we have shown that the invariance of energy(kinetic energy,potential energy)and virtual work is the common feature of vector bond graph and finite element method in struc-tural dynamics.Then we have discussed the vector bond graph representation of finite elementmethod in detail,there are:(1)the transformation of reference systems,(2)the transformation ofinertia matrices,stiffness matrices and vectors of joint force,(3)verctor bond graph representationof Lagrangian dynamic equation of structure.

Ab initio path-integral molecular dynamics and the quantum nature of hydrogen bonds

The hydrogen bond （HB） is an important type of intermolecular interaction, which is generally weak, ubiquitous, and essential to life on earth. The small mass of hydrogen means that many properties of HBs are quantum mechanical in nature. In recent years, because of the development of computer simulation methods and computational power, the influence of nuclear quantum effects （NQEs） on the structural and energetic properties of some hydrogen bonded systems has been intensively studied. Here, we present a review of these studies by focussing on the explanation of the principles underlying the simulation methods, i.e., the ab initio path-integral molecular dynamics. Its extension in combination with the thermodynamic integration method for the calculation of free energies will also be introduced. We use two examples to show how this influence of NQEs in realistic systems is simulated in practice.

Studies on the Structures and Hydrogen-bonding Interactions in Aqueous Glycine Solutions Using All-atom Molecular Dynamics Simulations and NMR Spectroscopy

All-atom molecular dynamics （MD） simulations and chemical shifts were used to study interactions and structures in the glycine-water system. Radial distribution functions and the hydrogen-bond network were applied in MD simulations. Aggregates in the aqueous glycine solution could be classified into different regions by analysis of the hydrogen-bonding network. Temperature-dependent NMR spectra and the viscosity of glycine in aqueous solutions were measured to compare with the results of MD simulations. The variation tendencies of the hydrogen atom chemical shifts and viscosity with concentration of glycine agree with the statistical results of hydrogen bonds in the MD simulations.

Molecular Dynamics Simulation Study of CIF in Water： Halogen Bonding Interaction in Liquid

Does the halogen bonding interaction co-exist in liquid when it competes with the hydrogen bonding interaction？ The classical molecular dynamics simulations for the solvation properties of CLF molecule in water are performed with the Lennard-Jones plus Coulomb electrostatic potential parameters that are optimized with ab initio interaction energy calculations for the pre-reactive H2O-CLF complex. We find that the halogen bonding interactions occur between O and CL atoms and have the comparable strength and population with respect to the hydrogen bonding interactions of C1...H.

Zero-point fluctuation of hydrogen bond in water dimer from ab initio molecular dynamics

Dynamic nature of hydrogen bond (H-bond) is central in molecular science of substance transportation, energy transfer, and phase transition in H-bonding networks diversely expressed as solution, crystal, and interfacial systems, thus attracting the state-of-the-art revealing of its phenomenological edges and sophisticated causes. However, the current understanding of the ground-state fluctuation from zero-point vibration (ZPV) lacks a firm quasi-classical base, concerning three basic dimensions as geometry, electronic structure, and interaction energy. Here, based on the ab initio molecular dynamics simulation of a ground-state water dimer, temporally separated fluctuation features in the elementary H-bond as the long-time weakening and the minor short-time strengthening are respectively assigned to two low-frequency intermolecular ZPV modes and two O–H stretching ones. Geometrically, the former modes instantaneously lengthen H-bond up to 0.2 Å whose time-averaged effect coverages to about 0.03 Å over 1-picosecond. Electronic-structure fluctuation crosses criteria' borders, dividing into partially covalent and noncovalent H-bonding established for equilibrium models, with a 370% amplitude and the district trend in interaction energy fluctuation compared with conventional dragging models using frozen monomers. Extended physical picture within the normal-mode disclosure further approaches to the dynamic nature of H-bond and better supports the upper-building explorations towards ultrafast and mode-specific manipulation.

Ab Initio Study of the Dynamical Si-O Bond Breaking Event in α-Quartz

The Si-O bond breaking event in the a-quartz at the first triplet （T1） excitation state is studied by using ab initio molecular dynamics （AIMD） and nudged elastic band calculations. A meta-stable non-bridging oxygen hole center and E1 center （NBOHC-E） is observed in the AIMD which consists of a broken Si-O bond with a Si-O distance of 2.54A. By disallowing the re-bonding of the Si and 0 atoms, another defect configuration （lll- Si/V-Si） is obtained and validated to be stable at both ground and excitation states. The NBOHC-E is found to present on the minimal energy pathway of the initial to IlI-Si/V-Si transition, showing that the generating of the NBOHC-E is an important step of the excitation induced structure defect. The energy barriers to produce the NBQHC-E＇ and Ⅲ-Si/V-Si defects are calculated to be 1.19 and 1.28eV, respectively. The electronic structures of the two defects are calculated by the self-consistent GW calculations and the results show a clear electron transition from the bonding orbital to the non-bonding orbital.

MOLECULAR DYNAMICS SIMULATION ON VIERATIONAL SPECTROSCOPIC FEATURES OF HYDROGEN BONDS IN CRYSTALLINE POLYMERS

Introduction The molecular dynamics simulation technique has recently proved to be a suitable alternative approachfor simulation of vibrational spectroscopy. In this study, molecular dynamics was utilized to understandlow frequency vibrations in highly ordered poly(ρ-phenylene terephthalmide) (PPTA). A key structuralfeature of this polymer is the presence of hydrogen bonds. There is little question that this strong localized

Static and dynamic mechanical behaviour of ECO-RPC

Ecological reactive powder concrete (ECO-RPC) with small sized and differentvolume fraction steel fibers was prepared by substitution of ultra-fine industrial waste powder for50% to 60% cement by weight and replacement of ground fine quartz sand with natural fine aggregate.The effect of steel fiber volume fraction and curing ages on the static mechanical behaviour ofECO-RPC was studied. Using the split Hopkinson pressure bar technique, the dynamic mechanicalbehaviour of ECO-RPC was investigated under different strain rates. The results show that the staticmechanical behaviour of ECO-RPC increases with the increase of steel fiber volume fraction andcuring ages. The type of ECO-RPC with the substitution of 25% ultra-fine slag, 25% ultra-fine flyash and 10% silica fume is better than the others with compressive strength, flexural strength, andfracture energy more than 200 MPa, 60 MPa and 30 kJ/m^2, respectively. ECO-RPC has excellent strainrate stiffening effects under dynamic load. Its peak stress, peak strain and the area understrain-stress curve increase with the increase of strain rate. Its fracture pattern changes frombrittleness to toughness under high strain rates.

Synergistic Dual Dynamic Bonds in Covalent Adaptable Networks

Covalent adaptable networks(CANs),comprising polymer networks crosslinked by dynamic covalent bonds(DCBs),have garnered considerable attention as sustainable materials.Mastering the stress relaxation of CANs is essential for controlling their viscoelastic properties.An unexpected acceleration of stress relaxation has been observed in CANs containing dual dynamic bonds.The dynamic behavior of the second dynamic bonds can accelerate stress relaxation and lower the relaxation activation energy of dual dynamic CANs compared to analogous CANs that rely on only one type of DCB.These findings complement current approaches that utilize catalysts or adjust network parameters.In this minireview,we summarize the synergistic acceleration effects in various CANs containing dual dynamic bonds.We classify these effects based on the second dynamic bonds,including noncovalent bonds,mechanical bonds,and the second DCBs.We also discuss the mechanisms behind this synergy.Finally,we highlight the challenges and offer perspectives on harnessing the synergistic effects of these dual dynamic systems to expand the chemistry and applications of CANs.