Inadequacy of the Lorentz-Berthelot combining rules for accurate predictions of equilibrium properties by molecular simulation

The wetting of H2O by CO2

Biphasic interfaces are complex but fascinating regimes that display a number of properties distinct from those of the bulk. The CO2-H2O interface, in particular, has been the subject of a number of studies on account of its importance for the carbon life cycle as well as carbon capture and sequestration schemes. Despite this attention, there remain a number of open questions on the nature of the CO2-H2O interface, particularly concerning the interfacial tension and phase behavior of CO2 at the interface. In this paper, we seek to address these ambiguities using ab initio-quality simulations. Harnessing the benefits of machine-learned potentials and enhanced statistical sampling methods, we present an ab initio-level description of the CO2-H2O interface. Interfacial tensions are predicted from 1 to 500 bars and found to be in close agreement with experiment at pressures for which experimental data are available. Structural analyses indicate the buildup of an adsorbed, saturated CO2 film forming at a low pressure (20 bars) with properties similar to those of the bulk liquid, but preferential perpendicular alignment with respect to the interface. The CO2 monolayer buildup coincides with a reduced structuring of water molecules close to the interface. This study highlights the predictive nature of machine-learned potentials for complex macroscopic properties of biphasic interfaces, and the mechanistic insight obtained into carbon dioxide aggregation at the water interface is of high relevance for geoscience, climate research, and materials science.

Machine-learned molecular mechanics force fields from large-scale quantum chemical data

The development of reliable and extensible molecular mechanics (MM) force fields-fast, empirical models characterizing the potential energy surface of molecular systems-is indispensable for biomolecular simulation and computer-aided drug design. Here, we introduce a generalized and extensible machine-learned MM force field, espaloma-0.3, and an end-to-end differentiable framework using graph neural networks to overcome the limitations of traditional rule-based methods. Trained in a single GPU-day to fit a large and diverse quantum chemical dataset of over 1.1 M energy and force calculations, espaloma-0.3 reproduces quantum chemical energetic properties of chemical domains highly relevant to drug discovery, including small molecules, peptides, and nucleic acids. Moreover, this force field maintains the quantum chemical energy-minimized geometries of small molecules and preserves the condensed phase properties of peptides and folded proteins, self-consistently parametrizing proteins and ligands to produce stable simulations leading to highly accurate predictions of binding free energies. This methodology demonstrates significant promise as a path forward for systematically building more accurate force fields that are easily extensible to new chemical domains of interest.

Influence of molecular parameters on the representativeness of interfacial properties of simple fluids

New parameterizations for the Lennard-Jones 12/6 potential capable of reproducing the vapor pressure and surface tension with sufficient precision, but not the liquid-vapor equilibrium densities for the case of simple fluids that include Ar, Kr, Xe, Ne, and CH4 are presented in this work. These results are compared with those derived from the family of Mie(n, 6) potentials, which adequately reproduce the coexistence curve and the vapor pressure, leaving aside the surface tension. In addition, a detailed analysis is presented on different parameterizations and methodologies, which have been developed in recent decades to estimate the interfacial properties of interest here for simple fluids, such as argon, which is a molecule that is, in principle, "simple" to study but that clearly reveals the enormous discrepancy between the results reported in the literature throughout these years. These facts undoubtedly reveal one of the fundamental problems in the context of molecular thermodynamics of fluids: reproducing different thermodynamic properties with sufficient precision from a single set of free parameters for some interaction potential. In order to show the scope of the parameterizations presented for the Lennard-Jones model, they were successfully applied to the case of binary mixtures, which included Ar-Kr, Ar-CH4, and Xe-Kr. Finally, and with the aim of showing a possible solution to the problem posed in this research, results of the same interfacial properties above mentioned for argon and methane were presented in this work by using a set of molecular interactions, called ANC2s, whose flexibility allowed to reproduce the experimental evidence with just one parameterization. The results reported in this work were generated using molecular dynamics simulations.

Stochastic computer experiments of the thermodynamic irreversibility of bulk nanobubbles in supersaturated and weak gas-liquid solutions

Spontaneous gas-bubble nucleation in weak gas-liquid solutions has been a challenging topic in theory, experimentation, and computer simulations. In analogy with recent advances in crystallization and droplet formation studies, the diffusive-shielding stabilization and thermodynamic irreversibility of bulk nanobubble (bNB) mechanisms are revisited and deployed to characterize nucleation processes in a stochastic framework of computer experiments using the large-scale atomic/molecular massively parallel simulator code. Theoretical bases, assumptions, and limitations underlying the irreversibility hypothesis of bNBs, and their computational counterparts, are extensively described and illustrated. In essence, it is established that the irreversibility hypothesis can be numerically investigated by converging the system volume (due to the finiteness of interatomic forces) and the initial dissolved-gas concentration in the solution (due to the single-bNB limitation). Helium nucleation in liquid Pb17Li alloy is selected as a representative case study, where it exhibits typical characteristics of noble-gas/liquid-metal systems. The proposed framework lays down the bases on which the stability of gas-bNBs in weak and supersaturated gas-liquid solutions can be inferred and explained from a novel perspective. In essence, it stochastically marches toward a unique irreversible state along out-of-equilibrium nucleation/growth trajectories. Moreover, it does not attempt to characterize the interface or any interface-related properties, neither theoretically nor computationally. It was concluded that bNBs of a few tens of He-atoms are irreversible when dissolved-He concentrations in the weak gas-liquid solution are at least ∼50 and ∼105 mol m-3 at 600 and 1000 K (and ∼80 MPa), respectively, whereas classical molecular dynamics -estimated solubilities are at least two orders of magnitude smaller.

Structural superlubricity at the interface of penta-BN2

Two-dimensional (2D) materials have been widely used as lubricants due to their weak interlayer interaction and low shear resistance for interlayer sliding. Composed entirely of five-membered rings, penta-BN2 monolayer has excellent thermal and mechanical stability, higher hardness and a negative Poisson's ratio. In this work, we investigate the frictional properties at both the commensurate and incommensurate contacting interfaces of penta-BN2 by adopting the molecular dynamics (MD) simulation method. Our calculations demonstrate robust superlubricity at the incommensurate contacting interface of penta-BN2. The ultra-low friction is explained by the potential energy surface (PES) fluctuations, interlayer binding energy and out-of-plane motion amplitude of the sliding layer. In addition, our calculations show that the anisotropy of friction at the commensurate contacting interface is more obvious compared with that at the incommensurate contacting interface. Finally, the influences of the size of the Moiré pattern, normal force, temperature and sliding velocity on the friction are examined. Our results show that 2D penta-BN2 is a promising solid lubricant, enriching the family of 2D lubrication materials.

Cluster-Assisted Mesoplasma Chemical Vapor Deposition for Fast Epitaxial Growth of SiGe/Si Heterostructures: A Molecular Dynamics Simulation Study

Co-condensation of mixed SiGe nanoclusters and impingement of SiGe nanoclusters on a Si substrate were applied using molecular dynamics (MD) simulation in this study to mimic the fast epitaxial growth of SiGe/Si heterostructures under mesoplasma chemical vapor deposition (CVD) conditions. The condensation dynamics and properties of the SiGe nanoclusters during the simulations were investigated first, and then the impingement of transient SiGe nanoclusters on both Si smooth and trench substrate surfaces under varying conditions was studied theoretically. The results show that the mixed nanoclusters as precursors demonstrate potential for enhancing epitaxial SiGe film growth at a high growth rate, owing to their loosely bound atomic structures and high mobility on the substrate surface. By varying cluster sizes and substrate temperatures, this study also reveals that smaller clusters and higher substrate temperatures contribute to faster structural ordering and smoother surface morphologies. Furthermore, the formed layers display a consistent SiGe composition, closely aligning with nominal values, and the cluster-assisted deposition method achieves the epitaxial bridging of heterostructures during cluster impingement, highlighting its additional distinctive characteristics. The implications of this work make it clear that the mechanism of fast alloyed epitaxial film growth by cluster-assisted mesoplasma CVD is critical for extending it as a versatile platform for synthesizing various epitaxial films.

Adsorption and Diffusion Properties of Functionalized MOFs for CO2 Capture: A Combination of Molecular Dynamics Simulation and Density Functional Theory Calculation

The capture of carbon dioxide (CO2) from fuel gases is a significant method to solve the global warming problem. Metal-organic frameworks (MOFs) are considered to be promising porous materials and have shown great potential for CO2 adsorption and separation applications. However, the adsorption and diffusion mechanisms of CO2 in functionalized MOFs from the perspective of binding energies are still not clear. Actually, the adsorption and diffusion mechanisms can be revealed more intuitively by the binding energies of CO2 with the functionalized MOFs. In this work, a combination of molecular dynamics simulation and density functional theory calculation was performed to study CO2 adsorption and diffusion mechanisms in five different functionalized isoreticular MOFs (IRMOF-1 through -5), considering the influence of functionalized linkers on the adsorption capacity of functionalized MOFs. The results show that the CO2 uptake is determined by two elements: the binding energy and porosity of MOFs. The porosity of the MOFs plays a dominant role in IRMOF-5, resulting in the lowest level of CO2 uptake. The potential of mean force (PMF) of CO2 is strongest at the CO2/functionalized MOFs interface, which is consistent with the maximum CO2 density distribution at the interface. IRMOF-3 with the functionalized linker -NH2 shows the highest CO2 uptake due to the higher porosity and binding energy. Although IRMOF-5 with the functionalized linker -OC5H11 exhibits the lowest diffusivity of CO2 and the highest binding energy, it shows the lowest CO2 uptake. Accordingly, among the five simulated functionalized MOFs, IRMOF-3 is an excellent CO2 adsorbent and IRMOF-5 can be used to separate CO2 from other gases, which will be helpful for the designing of CO2 capture devices. This work will contribute to the design and screening of materials for CO2 adsorption and separation in practical applications.

Universal Pairwise Interatomic van der Waals Potentials Based on Quantum Drude Oscillators

Repulsive short-range and attractive long-range van der Waals (vdW) forces play an appreciable role in the behavior of extended molecular systems. When using empirical force fields, the most popular computational methods applied to such systems, vdW forces are typically described by Lennard-Jones-like potentials, which unfortunately have a limited predictive power. Here, we present a universal parameterization of a quantum-mechanical vdW potential, which requires only two free-atom properties─the static dipole polarizability α1 and the dipole-dipole C6 dispersion coefficient. This is achieved by deriving the functional form of the potential from the quantum Drude oscillator (QDO) model, employing scaling laws for the equilibrium distance and the binding energy, and applying the microscopic law of corresponding states. The vdW-QDO potential is shown to be accurate for vdW binding energy curves, as demonstrated by comparing to the ab initio binding curves of 21 noble-gas dimers. The functional form of the vdW-QDO potential has the correct asymptotic behavior at both zero and infinite distances. In addition, it is shown that the damped vdW-QDO potential can accurately describe vdW interactions in dimers consisting of group II elements. Finally, we demonstrate the applicability of the atom-in-molecule vdW-QDO model for predicting accurate dispersion energies for molecular systems. The present work makes an important step toward constructing universal vdW potentials, which could benefit (bio)molecular computational studies.

Study on the Sliding Tribological Behavior of Oleic Acid-Modified MoS2 under Boundary Lubrication

The effects of MoS2 and MoS2 modified by adding oleic acid (OA) on the friction properties of lithium-based grease under boundary lubrication conditions are studied by molecular dynamics (MD) simulation and experiment. A rough wall boundary lubrication MD model with peaks and grooves is established to simulate the mechanical properties and lubrication effects of three lubrication systems on rough walls for the relative shear velocity between the two solid walls of 5 m/s at 500 MPa. The stress, wear amount, friction force, normal pressure, and friction heat of the friction surface are quantitively calculated. Simultaneously, a Retc friction and wear testing machine is used to measure the friction coefficient under different concentrations of additives and different pressures. The results show that the grease added with MoS2 can reduce friction, wear, and the temperature between friction pairs. However, under high pressure and shear, MoS2 can easily agglomerate and accumulate in the pits, reducing the lubricating effect. At the same time, since OA-modified MoS2 can reduce agglomeration, the modified MoS2 is adsorbed on the metal wall surface, forming a stable lubricant film. The main contributions of this article can be found in combining MD simulation and experimentation, establishing the connection between micronano structures and macroscopic properties, exploring the mechanism of the influence of wall roughness and particle size on the friction performance of lubricating oil, and providing a theory for predicting and developing high-performance new lubricating grease.

Desorption properties of the R600, R134a and their mixtures in several MOF structures: A molecular dynamics study

The efficiency of organic Rankine cycle is expected to be enhanced via using metal organic heat carriers. Thus, it is important to investigate the properties of organic working fluid in metal organic frameworks (MOFs). In this paper, molecular dynamics method was employed to study the desorption properties of R600, R134a, and their mixtures in MOF-5, IRMOF-16, and MOF-200 structures. The results show that the desorption capacity of pure working fluids are negatively correlated with their molecular size. In mixed refrigerant systems, the desorption of the working medium is competitive, and the higher the proportion of fluorine atoms in the system, the greater the desorption heat in IRMOF-16. It is noteworthy that the behavior of the desorption heat in the other two MOFs is diametrically opposed to this result. The desorption capacity is positively correlated with the average pore size, specific surface area, and porosity of the MOF. The self-diffusion coefficient of the working medium in MOF is inversely proportional to its molecular weight, but directly proportional to its molecular size and to the pore size of MOF.

Development, Validation, and Applications of Nonbonded Interaction Parameters between Coarse-Grained Amino Acid and Water Models

Interactions between amino acids and water play an important role in determining the stability and folding/unfolding, in aqueous solution, of many biological macromolecules, which affects their function. Thus, understanding the molecular-level interactions between water and amino acids is crucial to tune their function in aqueous solutions. Herein, we have developed nonbonded interaction parameters between the coarse-grained (CG) models of 20 amino acids and the one-site CG water model. The nonbonded parameters, represented using the 12-6 Lennard Jones (LJ) potential form, have been optimized using an artificial neural network (ANN)-assisted particle swarm optimization (PSO) (ANN-assisted PSO) method. All-atom (AA) molecular dynamics (MD) simulations of dipeptides in TIP3P water molecules were performed to calculate the Gibbs hydration free energies. The nonbonded force-field (FF) parameters between CG amino acids and the one-site CG water model were developed to accurately reproduce these energies. Furthermore, to test the transferability of these newly developed parameters, we calculated the hydration free energies of the analogues of the amino acid side chains, which showed good agreement with reported experimental data. Additionally, we show the applicability of these models by performing self-assembly simulations of peptide amphiphiles. Overall, these models are transferable and can be used to study the self-assembly of various biomaterials and biomolecules to develop a mechanistic understanding of these processes.

Dynamic phase transition induced by active molecules in a supercooled liquid

The purpose of this work is to use active particles to investigate the effect of facilitation on supercooled liquids. To this end we examine the behavior of a model supercooled liquid that is doped with a mixture of active particles and slowed particles. To simulate the facilitation mechanism, the activated particles are subjected to a force that follows the mobility of their most mobile neighboring molecule, while the slowed particles experience a friction force. Upon activation, we observe a fluidization of the entire medium along with a significant increase in dynamic heterogeneity. This effect is reminiscent of the fluidization observed experimentally when introducing molecular motors into soft materials. Interestingly, when the characteristic time τ_{μ}, used to define the mobility in the facilitation mechanism, matches the physical time t^{*} that characterizes the spontaneous cooperativity of the material, we observe a phase transition accompanied by structural aggregation of the active molecules. This transition is characterized by a sharp increase in fluidization and dynamic heterogeneity.

Molecular adsorption and self-diffusion of NO2, SO2, and their binary mixture in MIL-47(V) material

The loading dependence of self-diffusion coefficients (Ds) of NO2, SO2, and their equimolar binary mixture in MIL-47(V) have been investigated by using classical molecular dynamics (MD) simulations. The Ds of NO2 are found to be two orders of magnitude greater than SO2 at low loadings and temperatures, and its Ds decreases monotonically with loading. The Ds of SO2 exhibit two diffusion patterns, indicating the specific interaction between the gas molecules and the MIL-47(V) lattice. The maximum activation energy (Ea) in the pure component and in the mixture for SO2 are 16.43 and 18.35 kJ mol-1, and for NO2 are 2.69 and 1.89 kJ mol-1, respectively. It is shown that SO2 requires more amount of energy than NO2 to increase the diffusion rate. The radial distribution functions (RDFs) of gas-gas and gas-lattice indicate that the Oh of MIL-47(V) are preferential adsorption site for both NO2 and SO2 molecules. However, the presence of the hydrogen bonding (HB) interaction between the O of SO2 and the H of MIL-47(V) and also their binding angle (θ(OHC)) of 120° with the linkers of lattice indicate a stronger binding interaction between the SO2 and the MIL-47(V), but it does not occur with NO2. The jump-diffusion of SO2 between adsorption sites within the lattice has been confirmed by the 2D density distribution plots. Moreover, the extraordinarily high Sdiff for NO2/SO2 of 623.4 shows that NO2 can diffuse through the MIL-47(V) significantly faster than SO2, especially at low loading and temperature.

Insight into Oxygen Transport in Solid and High-Surface-Area Carbon Supports of Proton Exchange Membrane Fuel Cells

Understanding the oxygen transport mechanism through an ionomer film that covered the catalyst surface is essential for reducing local oxygen transport resistance and improving the low Pt-loading proton exchange membrane fuel cell performance. Besides the ionomer material, the carbon supports, upon which ionomers and catalyst particles are dispersed, also play a crucial role in local oxygen transport. Increasing attention has been paid to the effects of carbon supports on local transport, but the detailed mechanism is still unclear. Herein, the local oxygen transports based on conventional solid carbon (SC) and high-surface-area carbon (HSC) supports are investigated by molecular dynamics simulations. It is found that oxygen diffuses through the ionomer film that covered the SC supports via "effective diffusion" and "ineffective diffusion". The former denotes the process by which oxygen diffuses directly from the ionomer surface to the Pt upper surface through small and concentrated regions. In contrast, ineffective diffusion suffers more restrictions by both carbon- and Pt-dense layers, and thus, the oxygen pathways are long and tortuous. The HSC supports exhibit larger transport resistance relative to SC supports due to the existence of micropores. Also, the major transport resistance originates from the carbon-dense layer as it inhibits oxygen from diffusing downward and migrating toward the pore opening, while the oxygen transport inside the pore is facile along the pore's inner surface, which leads to a specific and short diffusion pathway. This work provides insight into oxygen transport behavior with SC and HSC supports, which is the basis for the development of high-performance electrodes with low local transport resistance.

Ion Separation Together with Water Purification via a New Type of Nanotube: A Molecular Dynamics Study

We propose a CNT-based concentric twin tube (CTT) as nanochannels for both water purification and ion separation at the nanoscale. In the model, a source reservoir dealing with the solution connects three containers via the CTT that has three subchannels for mass transfer. Before entering the three subchannels, the solution in the separating zone will form three layers (the aqua cations, water, and the aqua anions, respectively) by applying a charged capacitor with the two electrodes parallel to the flow direction of the solution. Under an electric field with moderate intensity, the three subchannels in the CTT have stable configurations for mass transfer. Since the water and the two types of aqua ions are collected by three different containers, the present model can realize both ion separation and water purification. The mass transfer in the subchannels will be sped up by an external pressure exerted on the solution in the source reservoir. The physical properties of the model, e.g., water purification speed, are analyzed with respect to the effects of the electric field, the size of CTT, and the concentration of solute, such as NaCl.

Nonreciprocal forces enable cold-to-hot heat transfer between nanoparticles

We study the heat transfer between two nanoparticles held at different temperatures that interact through nonreciprocal forces, by combining molecular dynamics simulations with stochastic thermodynamics. Our simulations reveal that it is possible to construct nano refrigerators that generate a net heat transfer from a cold to a hot reservoir at the expense of power exerted by the nonreciprocal forces. Applying concepts from stochastic thermodynamics to a minimal underdamped Langevin model, we derive exact analytical expressions predictions for the fluctuations of work, heat, and efficiency, which reproduce thermodynamic quantities extracted from the molecular dynamics simulations. The theory only involves a single unknown parameter, namely an effective friction coefficient, which we estimate fitting the results of the molecular dynamics simulation to our theoretical predictions. Using this framework, we also establish design principles which identify the minimal amount of entropy production that is needed to achieve a certain amount of uncertainty in the power fluctuations of our nano refrigerator. Taken together, our results shed light on how the direction and fluctuations of heat flows in natural and artificial nano machines can be accurately quantified and controlled by using nonreciprocal forces.

Influence of the Lennard-Jones Combination Rules on the Simulated Properties of Organic Liquids at Optimal Force-Field Parametrization

We recently introduced the CombiFF scheme [Oliveira et al., J. Chem. Theory Comput. 2020, 16, 7525], an approach for the automated refinement of force-field parameters against experimental condensed-phase data for large compound families. Using this scheme, once the time-consuming task of target-data selection and curation has been performed, the force-field optimization itself is both straightforward and fast. As a result, CombiFF provides an ideal framework for evaluating the influence of functional-form decisions on the accuracy of a force field at an optimal level of parametrization. We already used this approach to assess the effect of using an all-atom representation compared to united-atom representations in the force field [Oliveira et al., J. Chem. Theory Comput. 2022, 18, 6757]. Here, CombiFF is applied to assess the effect of three Lennard-Jones combination rules, geometric-mean (GM), Lorentz-Berthelot (LB), or Waldman-Hagler (WH), on the simulated properties of organic liquids. The comparison is performed in terms of the experimental liquid density ρ_liq, vaporization enthalpy ΔH_vap, surface-tension coefficient γ, static relative dielectric permittivity ϵ, and self-diffusion coefficient D. The calibrations of the three force-field variants are carried out independently against 2044 experimental values for ρ_liq, and ΔH_vap concerning 1516 compounds. The resulting root-mean-square deviations from experiment are 30.0, 26.9, and 36.7 kg m^-3 for ρ_liq and 2.8, 2.8, and 2.9 kJ mol^-1 for ΔH_vap, when applying the GM, LB, and WH combination rules, respectively. In terms of these (and the other) properties, the three combination rules perform comparatively well, with the GM and LB results being more similar to each other and slightly more accurate compared to experiment. In contrast, the use of distinct combination rules for the parameter calibration and property calculation leads to much larger errors.

A study of how solid–liquid interactions affect flow resistance and heat transfer at different temperatures based on molecular dynamics simulations

The heat transfer performance is improved as surface temperature increases, but the rise of surface temperature increases the flow resistance when solid–liquid interaction is weak but decreases the flow resistance when solid–liquid interaction is strong.

Innovative Insight into O2/N2 Permeation Behavior through an Ionomer Film in Cathode Catalyst Layers of Polymer Electrolyte Membrane Fuel Cells

It is crucial to clarify the permeation behavior of O₂ through the ionomer film for enhancing local O₂ transport in cathodes of fuel cells. However, all existing studies mainly deal with pure O₂ rather than air. Herein, the permeation behavior of the O₂/N₂ mixture through the ionomer film has been well explored in view of molecular bond length variations by molecular dynamics simulations. The bond lengths for O₂ and N₂ are shortened under a low hydration level when permeating through a dense layer with small free voids while no obvious change occurs at higher hydration. In the bulk ionomer region, O₂ molecules residing in water domains are energetically unstable because the bond lengths deviate far from the equilibrium length; thus, O₂ diffuses through the interfacial or hydrophobic regions. N₂ molecules show similar properties with O₂. This study provides a novel perspective on the permeation behavior of O₂ and N₂ through the ionomer film, which definitely benefits enhancing local O₂ transport.

Mechanism of Hexane Displaced by Supercritical Carbon Dioxide: Insights from Molecular Simulations

Supercritical carbon dioxide (sCO₂) has great potential for displacing shale oil as a result of its high solubility and low surface tension and viscosity, but the underlying mechanisms have remained unclear up to now. By conducting equilibrium molecular dynamics (EMD) simulations, we found that the displacing process could be divided into three steps: the CO₂ molecules were firstly injected in the central region of shale slit, then tended to adsorb on the SiO₂-OH wall surface and mix with hexane, resulting in loose hexane layer on the shale surface, and finally displaced hexane from the wall due to strong interactions between CO₂ and wall. In that process, the displacing velocity and efficiency of hexane exhibit parabolic and increased trends with pressure, respectively. To gain deep insights into this phenomenon, we further performed non-equilibrium molecular dynamics (NEMD) simulations and found that both the Onsager coefficients of CO₂ and hexane were correlated to increase with pressure, until the diffusion rate of hexane being suppressed by the highly dense distribution of CO₂ molecules at 12 MPa. The rapid transportation of CO₂ molecules in the binary components (CO₂ and hexane) actually promoted the hexane diffusion, which facilitated hexane flowing out of the nanochannel and subsequently enhanced oil recovery efficiency. The displacing process could occur effectively at pressures higher than 7.5 MPa, after which the interaction energies of the CO₂-wall were stronger than that of the hexane-wall. Taking displacing velocity and efficiency and hexane diffusion rate into consideration, the optimal injection pressure was found at 10.5 MPa in this work. This study provides detailed insights into CO₂ displacing shale oil and is in favor of deepening the understanding of shale oil exploitation and utilization.

Analysis of Gas Nanoclusters in Water Using All-Atom Molecular Dynamics

The Young-Laplace equation suggests that nanosized gas clusters would dissolve under the effects of perturbation. The fact that nanobubbles are observed raises questions as to the mechanism underlying their stability. In the current study, we used all-atom molecular dynamics simulations to investigate the gas-water interfacial properties of gas clusters. We employed the instantaneous coarse-graining method to define the fluctuating boundaries and analyze the deformation of gas clusters. Fourier transform analysis of the cluster morphology revealed that the radius and morphology deformation variations exhibit power law relationships with the vibrational frequency, indicating that the surface energy dissipated through morphology variations. Increasing pressure in the liquid region was found to alter the network of water molecules at the interface, whereas increasing pressure in the gas region did not exhibit this effect. The overall gas concentration was oversaturated and proportional to the gas density inside the clusters. However, the result of comparison with Henry's law reveals that the gas pressure at the interface reduced by the interfacial effects is much lower than that inside the gas region, thus reducing the demanding degree of oversaturation. Originating from the interfacial charge allocation, the magnitude of the electrostatic stress is greater than that of the gas pressure inside the cluster. However, the magnitude of the reversed tension induced by electrostatic stress is far below the value of interfacial tension. The potential of mean force (PMF) profiles revealed that a barrier potential at the interface hindered gas particles from escaping the cluster. Several effects contribute to stabilizing the gas clusters in water, including high-frequency morphological deformation, electrostatic stress, reduced interfacial tension, and gas oversaturation conditions. Our results suggest that gas clusters can exist in water under gas oversaturation conditions in the absence of hydrophobic contaminants or pinning charges at interfaces.

Self-Bending Behavior and Varying Bending Stiffness of Black Phosphorus/Molybdenum Disulfide (BP/MoS2) Heterostructure

Vertically-stacked black phosphorus/molybdenum disulfide (BP/MoS₂) heterostructures have broad prospects in flexible electronics. Bending is a common and highly concerned deformation for these flexible devices. However, the discrepancy in structures and properties among the components of 2D heterostructures often induces complex bending deformations. Here, the bending behaviors of BP, MoS₂ and BP/MoS₂ are investigated based on a molecular dynamics simulation. Compared with the constant bending stiffness of individual BP and MoS₂, that of BP/MoS₂ varies with the bending angle. Notably, a self-bending configuration induced by the lattice mismatch and size difference is found in BP/MoS₂. The corresponding self-bending amplitude depends on the degree of size difference of each component and the "soft/hard" competition between them. Moreover, the size difference leads to a weakened bending stiffness, which is ascribed to the reduction in interlayer interaction. A prediction formula is proposed to evaluate the bending stiffness of BP/MoS₂ with the size difference. This finding reveals novel ways for regulating the bending properties of 2D heterostructures, including the bending angle, characteristic size and stacking order. It offers an effective strategy for designing flexible devices with tunable bending performance.

Carbon Nanodots from an In Silico Perspective

Carbon nanodots (CNDs) are the latest and most shining rising stars among photoluminescent (PL) nanomaterials. These carbon-based surface-passivated nanostructures compete with other related PL materials, including traditional semiconductor quantum dots and organic dyes, with a long list of benefits and emerging applications. Advantages of CNDs include tunable inherent optical properties and high photostability, rich possibilities for surface functionalization and doping, dispersibility, low toxicity, and viable synthesis (top-down and bottom-up) from organic materials. CNDs can be applied to biomedicine including imaging and sensing, drug-delivery, photodynamic therapy, photocatalysis but also to energy harvesting in solar cells and as LEDs. More applications are reported continuously, making this already a research field of its own. Understanding of the properties of CNDs requires one to go to the levels of electrons, atoms, molecules, and nanostructures at different scales using modern molecular modeling and to correlate it tightly with experiments. This review highlights different in silico techniques and studies, from quantum chemistry to the mesoscale, with particular reference to carbon nanodots, carbonaceous nanoparticles whose structural and photophysical properties are not fully elucidated. The role of experimental investigation is also presented. Hereby, we hope to encourage the reader to investigate CNDs and to apply virtual chemistry to obtain further insights needed to customize these amazing systems for novel prospective applications.

Enrichment Effects Induced by Non-uniform Wettability Surfaces in the Presence of Non-condensable Gas: A Molecular Dynamics Simulation

For vapor condensation, the control of heterogeneous nucleation and spatial distribution of nuclei are crucial for regulating droplet dynamics and improving condensation efficiency. However, due to the complex characteristics of multicomponent, multiphase, and multiscale, the underlying mechanism of mixed vapor condensation remains unclear, especially at the nucleation stage. In this paper, we focus on the enrichment effects of non-uniform wettability surfaces by molecular dynamics simulation, which could intensify the droplet nucleation and growth processes in a water-air mixed system. The results clarify the inhibitory effect of non-condensable gas on droplet nucleation and prove that only 1% of non-condensable gas could reduce one half of the condensation performance from a molecular perspective. Furthermore, non-uniform surfaces are designed to promote the efficient enrichment of vapor molecules on nucleation sites, and the synergistic effect of hydrophilic and hydrophobic regions is proposed. In addition, the non-uniform wettability surfaces are characterized by varying the proportion and dispersion of hydrophilic regions. The results reveal that an optimal proportion of hydrophilic region (R = 5/6) could furnish the non-uniform surface with the best transfer performance. Moreover, the enhancement of condensation performance can also be achieved through the dispersed arrangement of hydrophilic regions. The results provide guidance for the optimized design of functionalized surfaces with enhanced mixed vapor condensation.

Effect of Ion Pair on Contact Angle for Phosphonium Ionic Liquids

The wettability of ionic liquids (ILs) is relevant to their use in various applications. However, a mechanistic understanding of how the cation-anion pair affects wettability is still evolving. Here, focusing on phosphonium ILs, wettability was characterized in terms of contact angle using experiments and classical molecular dynamics simulations. Both experiments and simulations showed that the contact angle was affected by the anion and increased as benzoate < salicylate < saccharinate. Further, the simulations showed that the contact angle decreased with increasing cation alkyl chain length for these anions paired with five different tetra-alkyl-phosphonium cations. The trends were explained in terms of adhesive and cohesive energies in the simulations and then correlated to the atomic scale differences between the anions and the cations.

Quantitative molecular simulations

All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.

Study on liquid-like SiGe cluster growth during co-condensation from supersaturated vapor mixtures by molecular dynamics simulation

Based on the co-condensation processes in the Si-Ge system upon cooling, as determined by molecular dynamics (MD) simulation, we explored the mixed cluster growth dynamics and structural properties leading to the synthesis of liquid-like SiGe nanoclusters. The results indicated that the cluster size quickly increased to large clusters by the coalescence of transient small clusters in the growth stage during co-condensation. The transient clusters at different temperatures were verified to have slightly Si-rich compositions and liquid-like structures. The coalescence of such nanoclusters at high temperatures led to spherical clusters with homogeneous intermixing. However, irregularly shaped clusters with attached mixed parts were obtained owing to incomplete coalescence at low temperatures. Ge atoms tended to move to the cluster surface to exploit their energetically favorable state during the restructuring process, leading to slightly Ge-rich components on the cluster surface. The degree of intermixing for SiGe nanoclusters was related to cluster size. Generally, small clusters appeared to be more segregated during restructuring.

Isothermal-isobaric algorithm to study the effects of rotational degrees of freedom-Benz water model

We have developed isothermal-isobaric algorithm for non-equilibrium Monte Carlo simulations. As first we have shown that the new method correctly predict density by comparing it to the density determined in canonical Monte Carlo simulations through the virial pressure. The new method was then used to study the effect of translational and rotational degrees of freedom on the structural and thermodynamic properties of the simple Mercedes-Benz water model. By holding one of the temperatures constant and varying the other one, we investigated how the position of the density maxima changes. We have observed that upon increase of rotational temperature the fluid become more Lennard-Jones like and the density maxima disappears.

The generalized continuous multiple step (GCMS) potential: model systems and benchmarks

The generalized continuous multiple step (GCMS) potential is presented in this work. Its flexible form allows forrepulsiveand/orattractivecontributions to be encoded through adjustable energy and length scales. The GCMS interaction provides a continuous representation of square-well, square-shoulder potentials and their variants for implementation in computer simulations. A continuous and differentiable energy representation is required to derive forces in conventional simulation algorithms. Molecular dynamics simulations are of particular interest when considering the dynamic properties of a system. The GCMS potential can mimic other interactions with a judicious choice of parameters due to the versatile sigmoid form. In this study, our benchmarks for the GCMS representation include triangular, Yukawa, Franzese, and Lennard-Jones potentials. Comparisons made with published data on volumetric phase diagrams, liquid structure, and diffusivity from model systems are in excellent agreement.

Molecular Dynamics Study on the Combined Effects of the Nanostructure and Wettability of Solid Surfaces on Bubble Nucleation

In this paper, molecular dynamics (MD) simulations are conducted to investigate the bubble nucleation process of liquid argon on surfaces with a nanostructure of different wettabilities. To account for the combined effects of the nanostructure and surface wettability on bubble nucleation, the variation of the bubble volume, the nucleation starting time, as well as the heat flux between the solid surface and fluid are examined. It is found that the position of bubble nucleation depends on the pillar wettability. Bubble nucleation occurs in the bulk of fluid when the pillar is hydrophilic, while it occurs on the pillar surface when the pillar is hydrophobic. Under an integrated influence of the free-energy barrier of nucleation and heat transfer, the nucleation occurs later as the wettability of the pillar gets weaker over surfaces with the hydrophilic pillar, while it occurs earlier as the wettability of the pillar gets weaker over surfaces with the hydrophobic pillar. Moreover, the peak heat flux decreases with the decrease of the pillar wettability over surfaces with the hydrophilic pillar, while it increases with the decrease of the pillar wettability over surfaces with the hydrophobic pillar, which can be explained from the perspective of the heat transfer efficiency and the timing of phase change occurrence. Finally, a new surface with mixed-wettable pillars is proposed, which is verified to be conducive to both bubble nucleation and heat transfer.

Quasi-Casimir coupling can induce thermal resonance of adsorbed liquid layers in a nanogap

In a vacuum nanogap, phonon heat transfer can be induced by quasi-Casimir coupling in the absence of electromagnetic fields. However, it is unknown whether phonons can be transmitted across a...

Enhanced local viscosity around colloidal nanoparticles probed by equilibrium molecular dynamics simulations

Nanofluids-dispersions of nanometer-sized particles in a liquid medium-have been proposed for a wide variety of thermal management applications. It is known that a solid-like nanolayer of liquid of typical thicknesses of 0.5-1 nm surrounding the colloidal nanoparticles can act as a thermal bridge between the nanoparticle and the bulk liquid. Yet, its effect on the nanofluid viscosity has not been elucidated so far. In this article, we compute the local viscosity of the nanolayer using equilibrium molecular dynamics based on the Green-Kubo formula. We first assess the validity of the method to predict the viscosity locally. We apply this methodology to the calculation of the local viscosity in the immediate vicinity of a metallic nanoparticle for a wide range of solid-liquid interaction strength, where a nanolayer of thickness 1 nm is observed as a result of the interaction with the nanoparticle. The viscosity of the nanolayer, which is found to be higher than its corresponding bulk value, is directly dependent on the solid-liquid interaction strength. We discuss the origin of this viscosity enhancement and show that the liquid density increment alone cannot explain the values of the viscosity observed. Rather, we suggest that the solid-like structure of the distribution of the liquid atoms in the vicinity of the nanoparticle contributes to the nanolayer viscosity enhancement. Finally, we observe a failure of the Stokes-Einstein relation between viscosity and diffusion close to the wall, depending on the liquid-solid interaction strength, which we rationalize in terms of the hydrodynamic slip.

A hydrogen bond-modulated soft nanoscale water channel for ion transport through liquid-liquid interfaces

Ion transport through interfaces is of ubiquitous importance in many fields such as electrochemistry, emulsion stabilization, phase transfer catalysis, liquid-liquid extraction and enhanced oil recovery. However, the knowledge of interfacial structures that significantly affect ion transport through liquid-liquid interfaces is still lacking due to the difficulty of observing nanoscale interfaces. We studied here the evolution of interfacial structures during ion transport through the decane-water interface under different ionic concentrations and external forces using molecular dynamics simulations. The roles of hydrogen bonds in ion transport through interfaces are revealed. We identified a soft nanoscale channel during ion transport through liquid-liquid interfaces and the decane phase under specific external force. The stability of the water channel and the ion transport velocity both increase with ionic concentration due to the layered ordering structures of the water near the channel surface. We observed that the stability and connectivity of the water channel in the decane phase are remarkably improved both by the high increase of the number of hydrogen bonds in the water channel with increasing ionic concentration, and by the conformational change in water molecules near the water channel surface. Our discovery of a soft nanoscale water channel by molecular simulations implies that there is a potential stable passage for ion transport through liquid-liquid interfaces.

Atomistic Insight into the Effects of Depositional Nanoparticle on Nanoscale Liquid Film Evaporation

Nanoscale liquid film evaporation plays an essential role in many engineering applications. This study carries out molecular dynamics simulations on the effects of the depositional nanoparticle's wettability and volume in base fluid on the evaporation process to understand how the depositional nanoparticle affects the evaporation heat transfer. Increasing the nanoparticle's wettability can enhance the evaporation heat transfer process, and the enhancement effect of the hydrophobic surface is more remarkable than that of the hydrophilic surface. This because the increasing wettability causes more significant solid-liquid interaction. However, the potential energy of argon atoms at the liquid-vapor interface is almost unaffected by wettability. Moreover, when the depositional nanoparticle locates below the free liquid film, increasing the nanoparticle volume has a better heat transfer performance. As the volume increases, the heat transfer through the nanoparticle becomes more obvious, which effectively enhances the heat transfer at the solid-liquid interface and the liquid-vapor interface. The latent heat of phase change at the liquid-vapor interface is almost unchanged so that the evaporation can be enhanced. This research provides an understanding of the effects of depositional nanoparticles on nanoscale evaporation, which can impact several engineering applications, including devices' cooling and fluid transport.

Chemical computational approaches for optimization of effective surfactants in enhanced oil recovery

The surfactant flooding becomes an attractive method among several Enhanced Oil Recovery (EOR) processes to improve the recovery of residual oil left behind in the reservoir after secondary oil recovery process. The designing of a new effective surfactant is a comparatively complex and often time consuming process as well as cost-effective due to its dependency on the crude oil and reservoir properties. An alternative chemical computational approach is focused in this article to optimize the performance of effective surfactant system for EOR. The molecular dynamics (MD), dissipative particle dynamics (DPD) and density functional theory (DFT) simulations are mostly used chemical computational approaches to study the behaviour in multiple phase systems like surfactant/oil/brine. This article highlighted a review on the impact of surfactant head group structure on oil/water interfacial property like interfacial tensions, interface formation energy, interfacial thickness by MD simulation. The effect of entropy in micelle formation has also discussed through MD simulation. The polarity, dipole moment, charge distribution and molecular structure optimization have been illustrated by DFT. A relatively new coarse-grained method, DPD is also emphasized the phase behaviour of surfactant/oil/brine as well as polymer-surfactant complex system.