Upgrading Gel Electrolytes Through Electrostatic-Induced Dual-Salt Paradigm for Superior Zn-Ion Battery Performance.

Gel electrolytes are gaining attention for rechargeable Zn-ion batteries because of their high safety, high flexibility, and excellent comprehensive electrochemical performances. However, current gel electrolytes still perform at mediocre levels due to incomplete Zn salts dissociation and side reactions. Herein, an electrostatic-induced dual-salt strategy is proposed to upgrade gel electrolytes to tackle intrinsic issues of Zn metal anodes. The competitive coordination mechanism driven by electrostatic repulsion and steric hindrance of dual anions promotes zinc salt dissociation at low lithium salt addition levels, improving ion transport and mechanical properties of gel electrolytes. Li+ ions and gel components coordinate with H2O, reducing active H2O molecules and inhibiting associated side reactions. The dual-salt gel electrolyte enables excellent reversibility of Zn anodes at both room and low temperatures. Zn||Polyaniline cells using the dual-salt gel electrolyte exhibit a high discharge capacity of 180 mAh g-1 and long-term cycling stability over 180 cycles at -20 °C. The dual-salt strategy offers a cost-effective approach to improving gel electrolytes for high-performance flexible Zn-ion batteries.

Mechanical properties of amorphous CO2 hydrates: insights from molecular simulations.

Understanding physicochemical properties of amorphous gas hydrate systems is of great significance to reveal structural stabilities of polycrystalline gas hydrate systems. Furthermore, amorphous gas hydrates can occur ordinarily in the nucleation events of gas hydrate systems. Herein, the mechanical properties of amorphous carbon dioxide hydrates are examined by means of all-atom classical molecular dynamic simulations. Our molecular simulation results reveal that mechanical strengths of amorphous carbon dioxide hydrates are evidently governed by temperatures, confining pressures, and ratios of water to carbon dioxide molecules. Notably, under compressive loads, amorphous carbon dioxide hydrates firstly exhibit monotonic strain hardening, followed by an interesting distinct phenomenon characterized by a steady flow stress at further large deformation strains. Furthermore, structural evolutions of amorphous carbon dioxide hydrates are analyzed on the basis of the N-Hbond DOP order parameter. These important findings can not only contribute to our understanding of the structural stabilities of amorphous gas hydrate systems, but also help to develop fundamental understandings about grain boundaries of gas hydrate systems.

Thermal conductivity of fivefold twinned silicon-germanium heteronanowires.

The thermal transport properties of five-fold twinned (5FT) germanium-silicon (Ge-Si) heteronanowires (h-NWs) with varying cross-sectional areas, germanium (Ge) domain ratios and heterostructural patterns are investigated using homogeneous nonequilibrium molecular dynamics (HNEMD) simulations. The results demonstrate a distinctive behavior in the thermal conductivity (κ) of 5FT-NWs, characterized by a "flipped" trend at a critical cross-sectional area. This behavior is attributed to the hydrodynamic phonon flow, arising from the normal three-phonon scattering process in the low-frequency region. In addition, the composition ratio of 5FT-NWs has a significant impact on reducing the κ of 5FT-NWs and suppressing the hydrodynamic effect. Intriguingly, as the homogeneous element domains are separated, stronger phonon hydrodynamic flows are observed in comparison to the adjacent homogeneous element domains. By analyzing various phonon properties, including phonon dispersion, three-phonon scattering rate, and phonon mean free path, critical insights into the origin of the differential κ in different 5FT-NW structures are provided. The findings deepen the understanding of the thermal transport properties of nanomaterials and hold implications for the design and development of nanoelectronics and thermoelectric devices.

Accurate prediction of heat conductivity of water by a neuroevolution potential.

We propose an approach that can accurately predict the heat conductivity of liquid water. On the one hand, we develop an accurate machine-learned potential based on the neuroevolution-potential approach that can achieve quantum-mechanical accuracy at the cost of empirical force fields. On the other hand, we combine the Green-Kubo method and the spectral decomposition method within the homogeneous nonequilibrium molecular dynamics framework to account for the quantum-statistical effects of high-frequency vibrations. Excellent agreement with experiments under both isobaric and isochoric conditions within a wide range of temperatures is achieved using our approach.

Role of mechanical deformation in the thermal transport of sI-type methane hydrate.

Natural gas hydrates (NGHs) are rising as an unconventional energy resource. The fundamental thermal characteristics of NGHs are of importance for natural gas exploitation from permafrost and oceanic sediments that are geomechanically deformed. Here, utilizing classic molecular dynamics simulations with all-atom (AA) and coarse-grained (CG) models of the methane guest molecule, the effects of mechanical strain on the thermal conductivity of sI-type methane hydrate are for the first time examined. Upon triaxial tension and compression, methane hydrate exhibits strong asymmetry in the stress responses. As the triaxial loads go from compression to tension, a reduction trend in the thermal conductivity is revealed for methane hydrate with both AA and CG models of methane, within a maximum reduction of over 44%. This reduction is because triaxial strain from compression to tension softens the phonon modes. Interestingly, there is a sudden rise in thermal conductivity at critical triaxial strain of 0.06, originating from that, at which, the phonon modes are hardened and the peaks of radial distribution functions are shifted back. This study provides important information on the thermal conductivity of methane hydrate, which is helpful for the practical production of natural gas from geo-deformed NGH-bearing sediments via a heating technique as well as evaluating their stability.

GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations.

We present our latest advancements of machine-learned potentials (MLPs) based on the neuroevolution potential (NEP) framework introduced in Fan et al. [Phys. Rev. B 104, 104309 (2021)] and their implementation in the open-source package gpumd. We increase the accuracy of NEP models both by improving the radial functions in the atomic-environment descriptor using a linear combination of Chebyshev basis functions and by extending the angular descriptor with some four-body and five-body contributions as in the atomic cluster expansion approach. We also detail our efficient implementation of the NEP approach in graphics processing units as well as our workflow for the construction of NEP models and demonstrate their application in large-scale atomistic simulations. By comparing to state-of-the-art MLPs, we show that the NEP approach not only achieves above-average accuracy but also is far more computationally efficient. These results demonstrate that the gpumd package is a promising tool for solving challenging problems requiring highly accurate, large-scale atomistic simulations. To enable the construction of MLPs using a minimal training set, we propose an active-learning scheme based on the latent space of a pre-trained NEP model. Finally, we introduce three separate Python packages, viz., gpyumd, calorine, and pynep, that enable the integration of gpumd into Python workflows.

Mechanical Anisotropy in Two-Dimensional Selenium Atomic Layers.

Two-dimensional (2D) trigonal selenium (t-Se) has become a new member in 2D semiconducting nanomaterial families. It is composed of well-aligned one-dimensional Se atomic chains bonded via van der Waals (vdW) interaction. The contribution of this unique anisotropic nanostructure to its mechanical properties has not been explored. Here, for the first time, we combine experimental and theoretical analyses to study the anisotropic mechanical properties of individual 2D t-Se nanosheets. It was found that its fracture strength and Young's modulus parallel to the atomic chain direction are much higher than along the transverse direction, which was attributed to the weak vdW interaction between Se atomic chains as compared to the covalent bonding within individual chains. Additionally, two distinctive fracture modes along two orthogonal loading directions were identified. This work provides important insights into the understanding of anisotropic mechanical behaviors of 2D semiconducting t-Se and opens new possibilities for future applications.

Self-Assembly of MoS2 Monolayer Sheets by Desulfurization.

Self-assembled structures of two-dimensional (2D) materials exhibit novel physical properties distinct from those of their parent materials. Herein, the critical role of desulfurization on the self-assembled structural morphologies of molybdenum disulfide (MoS2) monolayer sheets is explored using molecular dynamics (MD) simulations. MD results show that there are differences in the atomic energetics of MoS2 monolayer sheets with different desulfurization contents. Both free-standing and substrate-hosted MoS2 monolayer sheets show diversity in structural morphologies, for example, flat plane structures, wrinkles, nanotubes, and folds, depending on the desulfurization contents, planar dimensions, and ratios of length to width of MoS2 sheets. Particularly, at the critical desulfurization content, they can roll up into nanotubes, which is in good agreement with previous experimental observations. Importantly, these observed differences in the molecular structural morphologies between free-standing and substrate-hosted MoS2 monolayer sheets can be attributed to interatomic interactions and interlayer van der Waals interactions. Furthermore, MD results have demonstrated that the surface-driven stability of MoS2 structures can be indicated by the desulfurization contents on one surface of MoS2 monolayer sheets, and the self-assembly of MoS2 monolayer sheets by desulfurization can emerge to adjust their surface-driven stability. The study provides important atomic insights into tuning the self-assembling structural morphologies of 2D materials through defect engineering in the future science and engineering applications.

Mechanical creep instability of nanocrystalline methane hydrates.

Mechanical creep behaviors of natural gas hydrates are of importance for understanding the mechanical instability of gas hydrate-bearing sediments on Earth. Limited by the experimental challenges, intrinsic creep mechanisms of nanocrystalline methane hydrates remain largely unknown yet at the molecular scale. Herein, using large-scale molecular dynamics simulations, mechanical creep behaviors of nanocrystalline methane hydrates are investigated. It is revealed that mechanical creep responses are greatly dictated by internal microstructures of crystalline grain size and external conditions of temperature and static stress. Interestingly, a long steady-state creep is observed in nanocrystalline methane hydrates, which can be described by a modified constitutive Bird-Dorn-Mukherjee model. Microstructural analysis shows that deformations of crystalline grains, grain boundary diffusion and grain boundary sliding collectively govern the mechanical creep behaviors of nanocrystalline methane hydrates. Furthermore, structural transformation also appears to be important in their mechanical creep behaviors. This study provides new insights into understanding the mechanical creep scenarios of gas hydrates.

Mechanical ductile detwinning in CH3NH3PbI3 perovskite.

Twin boundaries (TBs) were identified to show conflicting positive/negative effects on the physical properties of CH3NH3PbI3 perovskites, but their effects on the mechanical properties are still unclear. Herein, the tensile characteristics of a variety of TB-dominated bicrystalline CH3NH3PbI3 perovskites are explored using molecular simulations. The results show that TB-containing CH3NH3PbI3 perovskites can be classified into four types based on their tensile ductile detwinning characteristics. Type I is characterized by smooth loading flow stress-strain responses, originating from relatively uniform stress distribution induced gradual amorphization in the TB region. Types II and III are represented by a sudden drop in loading stresses but then distinct ductile flow stress-strain curves, resulting from limited and large-area amorphizations of TB-involved structures, respectively. However, Type IV is highlighted by double apparent peaks in the loading curve, followed by a ductile flow response, originating from the stress-concentration of localization-to-globalization in the TB region, as well as amorphization. This study provides critical insights into the mechanical characteristics of CH3NH3PbI3 perovskites and indicates that TB engineering is a promising strategy to design mechanically robust hybrid organic-inorganic perovskite-based device systems.

Classical nucleation theory of ice nucleation: Second-order corrections to thermodynamic parameters.

Accurately estimating the nucleation rate is crucial in studying ice nucleation and ice-promoting and anti-freeze strategies. In classical nucleation theory, estimates of the ice nucleation rate are very sensitive to thermodynamic parameters, such as the chemical potential difference between water and ice Δµ and the ice-water interfacial free energy γ. However, even today, there are still many contradictions and approximations when estimating these thermodynamic parameters, introducing a large uncertainty in any estimate of the ice nucleation rate. Starting from basic concepts for a general solid-liquid crystallization system, we expand the Gibbs-Thomson equation to second order and derive second-order analytical formulas for Δµ, γ, and the nucleation barrier ΔG*, which are used in molecular dynamics simulations. These formulas describe well the temperature dependence of these thermodynamic parameters. This may be a new method of estimating Δµ, γ, and ΔG*.

Surface-topology-controlled mechanical characteristics of triply periodic carbon Schwarzite foams.

Bulky sp2-carbon Schwarzites with negative Gaussian curvature are promising structures for practical applications due to their unique properties such as high surface area, large porosity, and stability against graphitization. Herein, a comprehensive study on the tension, compression and shear mechanical characteristics of seven triply periodic carbon Schwarzite foams with distinct topologies is performed using reactive molecular dynamics (MD) simulations. All carbon Schwarzites exhibit unique thermal and mechanical properties that are markedly dictated by the topology. One of the structures presents a negative thermal expansion coefficient. Under uniaxial tension, the temperature is able to play a positive or negative role in the tensile stiffness, and there is no apparent positive relationship between tensile strength and mass density. Subjected to compression and shear loads, carbon Schwarzites can fail due to brittle fracture, and uniform and stepwise structural instabilities. Both compression- and tension-negative Poisson's ratios are revealed to originate from a curvature-flattening deformation mechanism. Analysis of the crush force efficiency, the stroke efficiency and the energy-absorption demonstrates that carbon Schwarzites are effective energy-absorbers. This study provides a fundamental understanding of the relationship between the topology and mechanical properties of carbon Schwarzites for designing 3D graphitic nanostructures with good mechanical performances.

CO2 wetting on pillar-nanostructured substrates.

CO2 capture by dropwise CO2 condensation on cold solid surfaces is a promising technology. Understanding the role of the nanoscale surface and topographical features of CO2 droplet wetting characteristics is of importance for CO2 capture by this technology, but this remains unexplored as of yet. Here, using large-scale molecular dynamics (MD) simulations, the contact angle and wetting behaviors of CO2 droplets on pillar-structured Cu-like surfaces are investigated for the first time. Dynamic wetting simulations show that, by changing the strength of the solid-liquid attraction [Formula: see text] a smooth Cu-like surface offers a transition from CO2-philic to CO2-phobic. By periodically pillared roughening of the Cu-like surfaces, however, a higher contact angle and a smaller spreading exponent of a liquid CO2 droplet are realized. Particularly, a critical crossover of CO2-philic to CO2-phobic can appear. The wetting of the pillared surfaces by a liquid CO2 droplet proceeds non-uniformly. A liquid CO2 droplet is capable of exhibiting a transition from the Cassie state to the Wenzel state with increasing [Formula: see text] increasing inter-pillar distance, and increasing pillar width. The wetting morphologies of the metastable Wenzel state of a CO2 droplet are very different from each other. The findings will inform the ongoing design of CO2-phobic solid surfaces for practical dropwise condensation-based CO2 capture applications.

Stability, deformation and rupture of Janus oligomer enabled self-emulsifying water-in-oil microemulsion droplets.

Microemulsions exist widely in nature, daily life and industrial manufacturing processes, including petroleum production, food processing, drug delivery, new material fabrication, sewage treatment, etc. The mechanical properties of microemulsion droplets and a correlation to their molecular structures are of vital importance to those applications. Despite studies on their physicochemical determinants, there are lots of challenges of exploring the mechanical properties of microemulsions by experimental studies. Herein, atomistic modelling was utilized to study the stability, deformation, and rupture of Janus oligomer enabled water-in-oil microemulsion droplets, aiming at revealing their intrinsic relationship with Janus oligomer based surfactants and oil structures. The self-emulsifying process from a water, oil and surfactant mixture to a single microemulsion droplet was modulated by the amphiphilicity and structure of the surfactants. Four microemulsion systems with an interfacial thickness in the range of 7.4-17.3 Å were self-assembled to explore the effect of the surfactant on the droplet morphology. By applying counter forces on the water core and the surfactant shell, the mechanical stability of the microemulsion droplets was probed at different ambient temperatures. A strengthening response and a softening regime before and after a temperature-dependent peak force were identified followed by the final rupture. This work demonstrates a practical strategy to precisely tune the mechanical properties of a single microemulsion droplet, which can be applied in the formation, de-emulsification, and design of microemulsions in oil recovery and production, drug delivery and many other applications.

How to characterize interfacial load transfer in spiral carbon-based nanostructure-reinforced nanocomposites: is this a geometry-dependent process?

There is a great deal of attention given to spiral carbon-based nanostructures (SCBNs) because of their unique mechanical, thermal and electrical properties along with fascinating morphology. Dispersing SCBNs inside a polymer matrix leads to extraordinary properties of nanocomposites in diverse fields. However, the role of the interfacial mechanical properties of these nanocomposites remains unknown. Here, using molecular dynamics simulations, the characteristics of interfacial load transfer of SCBN-polyethylene nanocomposites are explored. Considering the geometric characteristics of SCBNs, new insight into the separation behavior of nanoparticles in normal and sliding modes is addressed. Interestingly, the results show that the maximum force and the separation energy of the SCBNs are much larger than those of graphene because of interlocking of the coils and polymer. The heavy influence of changes in the geometric characteristics of SCBNs on the separation behavior is observed. Pullout tests reveal that the influence of parameters such as the length and number of polyethylene chains, temperature, and functionalization of the SCBNs on the interfacial mechanical properties is also significant. This study sheds new light in understanding the crucial effect of the interaction of SCBNs with polymer chains on the interfacial mechanical properties, which can lead to better performance of nanocomposites.

Grain-Size-Controlled Mechanical Properties of Polycrystalline Monolayer MoS2.

Pristine monocrystalline molybdenum disulfide (MoS2) possesses high mechanical strength comparable to that of stainless steel. Large-area chemical-vapor-deposited monolayer MoS2 tends to be polycrystalline with intrinsic grain boundaries (GBs). Topological defects and grain size skillfully alter its physical properties in a variety of materials; however, the polycrystallinity and its role played in the mechanical performance of the emerging single-layer MoS2 remain largely unknown. Here, using large-scale atomistic simulations, GB structures and mechanical characteristics of realistic single-layered polycrystalline MoS2 of varying grain size prepared by confinement-quenched method are investigated. Depending on misorientation angle, structural energetics of polar-GBs in polycrystals favor diverse dislocation cores, consistent with experimental observations. Polycrystals exhibit grain-size-dependent thermally induced global out-of-plane deformation, although defective GBs in MoS2 show planar structures that are in contrast to the graphene. Tensile tests show that presence of cohesive GBs pronouncedly deteriorates the in-plane mechanical properties of MoS2. Both stiffness and strength follow an inverse pseudo Hall-Petch relation to grain size, which is shown to be governed by the weakest link mechanism. Under uniaxial tension, transgranular crack propagates with small deflection, whereas upon biaxial stretching, the crack grows in a kinked manner with large deflection. These findings shed new light in GB-based engineering and control of mechanical properties of MoS2 crystals toward real-world applications in flexible electronics and nanoelectromechanical systems.

Mechanical properties of monocrystalline and polycrystalline monolayer black phosphorus.

The mechanical properties of monocrystalline and polycrystalline monolayer black phosphorus (MBP) are systematically investigated using classic molecular dynamic simulations. For monocrystalline MBP, it is found that the shear strain rate, sample dimensions, temperature, atomic vacancies and applied statistical ensemble affect the shear behaviour. The wrinkled morphology is closely connected with the direction of the in-plane shear, dimensions of the samples, and applied ensembles. Particularly, small samples subjected to loading/unloading of the shear deformation along the armchair direction demonstrate a clear mechanical hysteresis loop. For polycrystalline MBP, the maximum shear stress as a function of the average grain size follows an inverse pseudo Hall-Petch type relationship under an isothermal-isobaric (NPT) ensemble, whereas under a canonical (NVT) ensemble, the maximum shear stress of polycrystalline MBP exhibits a 'flipped' behaviour. Furthermore, polycrystalline MBP subjected to uniaxial tension also exhibits a strongly grain size-dependent mechanical response, and it can fail by brittle intergranular and transgranular fractures because of its weaker grain boundary structures and the direction-dependent edge energy, respectively. These findings provide useful insight into the mechanical design of BP for nanoelectronic devices.

Carbon nanotubes kirigami mechanical metamaterials.

Imparting elasticity and functionality to materials is one of the key objects of materials science research. Here, inspired by the art of kirigami, mechanical metamaterials comprising carbon nanotubes (CNTs) are hypothetically constructed. Using classical molecular dynamics (MD) simulations, a systematic study of the elastic limit, extensibility and yield stress of as-generated CNTs kirigami (CNT-k) is performed. Three designated kirigami patterns are employed to achieve high stretchability of CNTs. It is shown that CNT-k typically exhibits three distinct deformation stages, of which the first stage, which is referred to as geometric deformation, contributes quite a high proportion of the ductility. Various geometric parameters of CNT-k that influence the key mechanical properties of interest are respectively discussed. Three types of CNT-k with specifically identical geometric parameters exhibit distinct mechanical characteristics. This study provides an interesting example of interplay between the geometry, ductility, and mechanical characteristics of tubular materials.

Computer simulation of water desalination through boron nitride nanotubes.

Development of high-efficiency and low-cost seawater desalination technologies is critical to solving the global water crisis. Herein we report a fast water filtering method with high salt rejection by boron nitride nanotubes (BNNTs). The effect of the radius of BNNTs on water filtering and salt rejection was investigated by molecular dynamics (MD) simulation. Our simulation results demonstrate that fast water permeation and high salt rejection could be achieved by BNNT(7,7) under both high pressure and low pressure. The potential of mean force (PMF) of Na+ ion and water molecule through BNNT(7,7) further revealed the mechanism of seawater desalination by BNNT(7,7). Using BNNT(7,7) array, a 10 cm2 nanotube membrane with 1.5 × 1013 pores per cm2 will produce freshwater with a flow rate of 98 L per day per MPa under 100 MPa. Our study shows the potential application of BNNTs membrane for fast and efficient desalination.

Extraordinary deformation capacity of smallest carbohelicene springs.

The extraordinary deformation and loading capacity of nine different [∞]carbohelicene springs under uniaxial tension up to their fracture were computed using the density functional theory. The simulations comprised either the experimentally synthetized springs of hexagonal rings or the hypothetical ones that contained irregularities (defects) as, for example, pentagons replacing the hexagons. The results revealed that the presence of such defects can significantly improve mechanical properties. The maximum reversible strain varied from 78% to 222%, the maximum tensile force varied in the range of 5 nN to 7 nN and, moreover, the replacement of hexagonal rings by pentagons or heptagons significantly changed the location of double bonds in the helicenes. The fracture analysis revealed two different fracture mechanisms that could be related to the configurations of double and single bonds located at the internal atomic chain. Simulations performed with and without van der Waals interactions between intramolecular atoms showed that these interactions played an important role only in the first deformation stage.