Nonequilibrium Molecular Dynamics Simulations Research Articles

Thermal transport in disordered wurtzite ScAlN alloys using machine learning interatomic potentials

Disordered wurtzite ScAlN alloy is an emerging material with a strong piezoelectric coefficient that promotes its application in 5G communication devices. To extend the service lifetime of the device, it is an urgent task to understand the thermal transport properties of this alloy for thermal management. Yet a thorough understanding of the thermal transport in the disordered ScAlN alloy is hindered by the inaccuracy of empirical interatomic potentials adopted for atomistic simulations. Herein, we developed an accurate and efficient neuro-evolution potential and combined it with Allen and Feldman theory and molecular dynamics to investigate thermal transport properties of disordered ScAlN alloys with varying ratios. This potential is shown to have the same computational accuracy as first principles by comparing the phonon spectrum and radial distribution function. Based on the homogeneous non-equilibrium molecular dynamics simulations, the lattice thermal conductivity of the disordered alloys reaches its maximum value (less than 12W/m·K) at room temperature. Further investigation of the temperature effects on the vibrational modes shows that it affects the mode lifetimes for the propagons and the mode heat capacities for the diffusions. In addition, it was found that both mode lifetime and diffusivity decrease with increasing Sc concentration. This work aims to train a universal machine learning interatomic potential of disordered ScAlN alloys, seeking to accelerate the thermal management design of related devices.

Comparative Retention Analysis of Intercalated Cations Inside the Interlayer Gallery of Lamellar and Nonlamellar Graphene Oxide Membranes in Reverse Osmosis Process: A Molecular Dynamics Study.

Over the past decade, multilayered graphene oxide (GO) membranes have emerged as promising candidates for desalination applications. Despite their potential, a comprehensive understanding of separation mechanisms remains elusive due to the intricate morphology and structural arrangement of interlayer galleries. Moreover, a critical concern of multilayered GO membranes is their susceptibility to swelling within aqueous environments, which hinders their practical implementation. Therefore, this study introduces cation intercalation within GO laminates to elucidate the underlying factors governing swelling behavior and subsequently mitigate it. Moreover, this study performed nonequilibrium molecular dynamics simulations on the cation (Mg2+ or K+)-intercalated lamellar and nonlamellar GO membranes to understand the effect of the arrangement of GO sheets on the retention time of intercalated cations within GO layers, water permeance, and salt rejection mechanism in the reverse osmosis process using cation-intercalated GO membranes. Our results highlight that lamellar GO membranes exhibit higher water permeance, attributed to their well-defined interlayer gallery structure. On the other hand, nonlamellar GO membranes display superior salt rejection due to their complex interlayer gallery structure that impedes salt permeation. Moreover, the structural complexity of nonlamellar GO membranes contributes to greater stability by retention of the more intercalated cations for a longer time within the layers. Furthermore, it is observed that a higher percentage of Mg2+ cations remained inside the GO laminates as compared to K+ cations, hence resulting in the greater stability of the Mg2+-intercalated GO membrane in the aqueous environment.

Thermal transport of confined water molecules in quasi-one-dimensional nanotubes

Dimensions and molecular structures play pivotal roles in the principle of heat conduction. The dimensional characteristics of a solution within nanoscale systems depend on the degrees of confinement. However, the influence of such variations on heat transfer remains inadequately understood. Here, we perform quasi-one-dimensional non-equilibrium molecular dynamics simulations to calculate the thermal conductivity of water molecules confined in carbon nanotubes. The structure of water molecules is determined depending on the nanotube radius, forming a single-file, a single-layer, and a double-layer structure, corresponding to an increasing radius order. We reveal that the thermal conductivity of liquid water has a sublinear dependency on nanotube length exclusively when water molecules form a single file. A stronger confinement leads to behavioral and structural characteristics closely resembling a one-dimensional nature. Moreover, single-layer-structured water molecules exhibit enhanced thermal conductivity. We elucidate that this is due to the increase in the local water density and the absence of transitions to another layer, which typically occurs in systems with double-layer water structures within relatively large radius nanotubes.

Efficient desalination in fluorinated multilayered covalent organic framework: Tunning the energy landscape inside the nanochannel

The design and construction of advanced nanochannels can significantly enhance the performance of reverse osmosis (RO) membranes in desalination, offering potential solutions to the severe freshwater shortage crisis. In this study, we investigated the desalination performance of nanochannels formed by multilayered covalent organic frameworks (COF) that had undergone various types of chemical functionalization, through nonequilibrium molecular dynamics simulations. Our data reveals that fluorination of the nanochannel results in the most impressive desalination performance, as evidenced by significant improvements in both water permeability (increasing from 39.1 to 100.1 L cm−2 day−1 MPa−1) and salt rejection (rising from 91.9 % to 97.1 %). Especially, in multilayered systems, fluorinated nanochannel mitigates the typical decline in water permeability with increasing membrane thickness. For the working mechanism, the free energy analysis indicates that fluorination contributes to a smooth energy landscape for water molecules inside the nanochannel, reducing their interaction with the nanochannel, which results in a more dispersed arrangement of water molecules with fewer hydrogen bonds and a significantly enhanced water velocity, thus facilitating their fast and unobstructed permeation. Therefore, this study offers a promising semipermeable membrane for high-efficiency desalination, achieved through fluorination. It also demonstrates an indicative approach for advanced channel design through fluorination engineering.

Extreme near-field heat transfer between silica surfaces

Despite recent experiments exhibiting an impressive enhancement in radiative heat flux between parallel planar silica surfaces with gap sizes of about 10 nm, the exploration of sub-nanometric gap distances remains unexplored. In this work, by employing non-equilibrium molecular dynamics (NEMD) simulations, we study the heat transfer between two SiO2 plates in both their amorphous and crystalline forms. When the gap size is 2 nm, we find that the heat transfer coefficient experiences a substantial ∼30-fold increase compared to the experimental value at the gap size of 10 nm confirming the dependence on the distance inversely quadratic as predicted by the fluctuational electrodynamics (FE) theory. Comparative analysis between NEMD and FE reveals a generally good agreement, particularly for amorphous silica. Spectral heat transfer analysis demonstrates the profound influence of gap size on heat transfer, with peaks corresponding to the resonances of dielectric function. Deviations from the fluctuational electrodynamics theory at smaller gap sizes are interpreted in the context of acoustic phonon tunneling and the effects of a gradient of permittivity close to the surfaces.

Orientation-dependent deformation mechanisms of alpha-uranium single crystals under shock compression

Large-scale non-equilibrium molecular dynamics (NEMD) simulations were employed to investigate the dynamic deformations of alpha-uranium (α-U) single crystals subjected to varying shock strengths along low-index crystallographic orientations. The pronounced anisotropy of α-U gives rise to a complex microstructural evolution under shock loading. In-depth microstructural analysis of post-shock specimens reveals the identification of multiple dynamic deformation mechanisms. Notably, when the shock loading direction aligns with the a-axis, dynamic deformation of the α-U single crystals is primarily dominated by lattice instability, which attributes to a crystalline-to-amorphous transition serving as the dominant shear stress relaxation pathway. On the other hand, shock loading along the b-axis results in an abundance of deformation twins, with twinning planes identified as (130) and (13¯0). During the twinning event, the α-U matrix undergoes a transition to a metastable intermediate phase, subsequently decomposing into a composite structure comprising α-U twins and matrix. This unconventional twinning mechanism significantly deviates from classical theories. Furthermore, upon loading along the c-axis, twinning and a phase transition from α-U to body-centered tetragonal phase (bct-U) occur in α-U single crystal samples. Given that the pressure threshold of this phase transition predicted by ab initio calculations is as high as ∼270 GPa, the phase transition from α-U to bct-U might be implausible. An alternative interatomic potential of uranium with the higher pressure threshold was employed to reinvestigate the shock response of α-U single crystals along the c-axis. The phase transition of α-U to bct-U disappears, and twinning dominates the plastic deformation, with the twinning orientation conforming to the {112} twinning. The strong anisotropy of the α-U lattice triggers a wealth of orientation-dependent dynamic deformation mechanisms. The activation of the twinning system is evidently associated with the loading direction, constituting the potential cause for the discovery of multiple twinning variants during the deformation in polycrystalline uranium.

Polyamide Reverse Osmosis Membrane Compaction and Relaxation: Mechanisms and Implications for Desalination Performance

This work investigates the compaction and relaxation behavior of composite reverse osmosis (RO) polyamide (PA) selective layers, utilizing non-equilibrium molecular dynamic (NEMD) simulations and well-controlled permeation experiments. Composite PA-RO membranes are prepared by interfacial polymerization of different amine monomer ratios to vary PA film crosslinking degree (CD). Wet-testing results suggest that “tighter” (more highly crosslinked) composite RO membranes with a higher PA layer CD undergo 65% less compaction and recover 17% more of their initial permeability (termed “relaxation”) when the pressure is relieved compared to lower CD PA layers. NEMD simulations provide a visual and quantitative characterization of PA layer’s free volume changes in operando. NEMD simulations also corroborate experimental findings and elucidate the viscoelastic properties of the PA layer that govern compaction and relaxation behavior. The mechanisms of compaction under water permeation derived from NEMD simulations further support the notion of viscous flow of water through interconnected free volume elements (i.e., “pores”) within crosslinked PA films.

Enhancing interfacial thermal conductivity of copper-carbon nanotube array composite via metallic bonding: Molecular dynamics simulations

In this work, a novel strategy for improving the interface thermal conductivity (ITC) between the vertically aligned carbon nanotube (VACNT) arrays and copper was proposed and verified using the non-equilibrium molecular dynamics (NEMD) simulations. The introduced metal elements (Ag, Al, Ti, and Ni) enhanced the ITC of the Cu-VACNT interface region via a strong adhesive force. Notably, different from the physisorption between Ag/Al elements and VACNT, a strong interfacial interaction was introduced between the chemically bonded Ni/Ti and VACNT. In addition, the matching degree of the phonon state density (PDOS) at the Ni-VACNT interface increases with increasing contact area between the Ni atoms and Cu substrate. Nevertheless, it was found that over-strong interfacial interaction between the VACNT and Ti element deteriorates the matching degree of high-frequency phonons at the Ti-VACNT interface. These findings offer valuable criteria for designing and selecting materials to enhance interfacial thermal conduction in composite electronic packaging.

Analysis of the gas transport resistance of CO2 and CH4 through ultra-thin DD3R zeolite membrane

Ultra-thin DD3R zeolite membranes exhibit excellent separation properties for CO2 capture from natural gas, as well as the lower transmembrane resistance. However, the interfacial transport of permeation gas in such systems may dominate the whole separation process. In this work, we employed external force non-equilibrium molecular dynamic (EF-NEMD) simulation to predict the single-component permeabilities of CO2 and CH4 through a defect-free DD3R zeolite membrane at different pressured drops. By explicitly including the gas transport from bulk phase to zeolitic channels, the predicted results of EF-NEMD simulation showed good agreements with the reported experimental data. The contribution of interfacial resistance over the total transport resistance (Rinter/Rtotal) was further evaluated for the membranes with different thicknesses under temperatures between 273 K and 373 K. The results revealed a decrease in Rinter/Rtotal with increasing membrane thickness, leading to a reduction in CO2/CH4 selectivity, and the critical thickness (Rinter/Rtotal < 0.01) was determined to be approximately 300 nm at pressure of 1.0 MPa and temperature of 298 K. Similarly, the Rinter/Rtotal decreased with increasing temperature due to the augmentation in molecule kinetic energy. Furthermore, it was confirmed that the gas adsorption effect could expand the effective size of eight membered ring (8-MR) channels in DD3R zeolite membrane, thereby decreasing the CO2/CH4 selectivity, despite higher permeation fluxes. The above discoveries essentially benefit the design of the ultra-thin DD3R zeolite membrane through an understanding of CO2 and CH4 transport behaviors.

TTCF4LAMMPS: A toolkit for simulation of the non-equilibrium behaviour of molecular fluids at experimentally accessible shear rates

We present TTCF4LAMMPS, a toolkit for performing non-equilibrium molecular dynamics (NEMD) simulations to study the fluid behaviour at low shear rates using the LAMMPS software. By combining direct NEMD simulations and the transient-time correlation function (TTCF) technique, we study the fluid response to shear rates spanning 15 orders of magnitude. We present two examples for simple monatomic systems: one consisting of a bulk liquid and another with a liquid layer confined between two solid walls. The small bulk system is suitable for testing on personal computers, while the larger confined system requires high-performance computing (HPC) resources. We demonstrate that the TTCF formalism can successfully detect the system response for arbitrarily weak external fields. We provide a brief mathematical explanation for this feature. Although we showcase the method for simple monatomic systems, TTCF can be readily extended to study more complex molecular fluids. Moreover, in addition to shear flows, the method can be extended to investigate elongational or mixed flows as well as thermal or electric fields. The high computational cost needed for the method is offset by the two following benefits: i) the cost is independent of the magnitude of the external field, and ii) the simulations can be made highly efficient on HPC architectures by exploiting the parallel design of the algorithm. We expect the toolkit to be useful for computational researchers striving to study the nonequilibrium behaviour of fluids under experimentally-accessible conditions. Program summaryProgram title: TTCF4LAMMPSCPC Library link to program files:https://doi.org/10.17632/hh2rkcxbrf.1Developer's repository link:https://github.com/edwardsmith999/TTCF4LAMMPSLicensing provisions: GNU General Public License 3Programming language: Python 3Nature of problem: Measuring the nonequilibrium behaviour of bulk and confined fluids under experimentally accessible strain rates in non-equilibrium molecular dynamics (NEMD) simulations.Solution method: Creating a Python-based code that utilises the transient-time correlation function method and the LAMMPS software to enable the bulk fluid properties (e.g. viscosity) and confined fluid interfacial properties (e.g. shear stress and slip velocity) to be computed at low shear rates with NEMD.

Effects of shear strain on shock response in single crystal iron

With large-scale non-equilibrium molecular dynamics simulations and in situ x-ray diffraction analysis, we conducted a systematic investigation into the effects of pre-existing shear strain (γxy) on the shock response of single crystal iron. Our findings reveal significant effects of γxy on the deformation of the crystal structure during shock loading, leading to noticeable alterations in the propagation of shock waves. Specifically, during the elastic stage, the presence of γxy results in a reduction of shock strength, consequently diminishing the magnitude of elastic lattice strain (εe). In the plastic stage, γxy stimulates the α–ε phase transformation, and structure deformation undergoes a transition from the sequential activity of dislocation-to-transformation to the synchronous activity of dislocation and transformation. This transition inhibits the propagation of plastic waves and consequently broadens the elastic regime. Additionally, the introduction of γxy activates different slip systems, as it alters the corresponding resolved shear stress. Concurrently, the presence of γxy triggers the activation of different high-pressure phase variants. Our investigation sheds light on the fundamental physics of iron under shock compression and the influence of pre-existing shear strain on its behavior.

Boundary Lubrication Performance of Polyelectrolyte-Surfactant Complexes on Biomimetic Surfaces.

Aqueous mixtures of oppositely charged polyelectrolytes and surfactants are useful in many industrial applications, such as shampoos and hair conditioners. In this work, we investigate the friction between biomimetic hair surfaces in the presence of adsorbed complexes formed from cationic polyelectrolytes and anionic surfactants in an aqueous solution. We apply nonequilibrium molecular dynamics (NEMD) simulations using the coarse-grained MARTINI model. We first developed new MARTINI parameters for cationic guar gum (CGG), a functionalized, plant-derived polysaccharide. The complexation of CGG and the anionic surfactant sodium dodecyl sulfate (SDS) on virgin and chemically damaged biomimetic hair surfaces was studied using a sequential adsorption approach. We then carried out squeeze-out and sliding NEMD simulations to assess the boundary lubrication performance of the CGG-SDS complex compressed between two hair surfaces. At low pressure, we observe a synergistic friction behavior for the CGG-SDS complex, which gives lower shear stress than either pure CGG or SDS. Here, friction is dominated by viscous dissipation in an interfacial layer comprising SDS and water. At higher pressures, which are probably beyond those usually experienced during hair manipulation, SDS and water are squeezed out, and friction increases due to interdigitation. The outcomes of this work are expected to be beneficial to fine-tune and screen sustainable hair care formulations to provide low friction and therefore a smooth feel and reduced entanglement.

Viscosities of inhomogeneous systems from generalized entropy scaling

This study extends entropy scaling to inhomogeneous fluids by using the classical density functional theory together with a new viscosity reference that takes into account the influence of solid–fluid interactions on the fluid viscosity. The density functional theory uses a Helmholtz energy functional based on the perturbed-chain statistical associating fluid theory; the local residual entropy per particle is determined from the temperature derivative of the Helmholtz energy functional in combination with an appropriate weighted density profile. The weighted density calculation requires a single transferable parameter, which is adjusted to a reference molecular dynamics simulation. In particular, local viscosity values for fluids under nanoconfinement near solid–fluid interfaces are predicted using the same entropy scaling parameters as for homogeneous fluids. We validate the model by comparing viscosity and velocity profiles with results from non-equilibrium molecular dynamics simulations of a Couette flow in a slit pore. Good agreement is found between the entropy scaling model and the non-equilibrium molecular dynamics results for both the viscosity and velocity profiles of the Lennard–Jones truncated and shifted fluid. The proposed model extrapolates well to systems with different temperatures, fluid densities, and shear forces as well as to systems with different wetting behaviors. These results demonstrate that entropy scaling can be generalized to inhomogeneous fluids using an appropriate combination of residual entropy profile and viscosity reference.

Mechanical and thermal transport in nanoporous graphene-like BC6N

This work analyzes the thermal transport (TT) and mechanical properties (MP) of nanoporous graphene-like BC6N monolayer utilizing non-equilibrium molecular dynamics (NEMD) and molecular dynamic (MD) simulations. Herein, the heat transfer and mechanical of nanoporous graphene-like BC6N membrane with various circle-pore size defects are systematically reviewed and numerically examined. The results show that the BC6N sheet has a balance on mechanical properties in two directions under biaxial tension as changing pore size from D = 2.0 to D = 7.0 nm. Young's modulus and ultimate strength of BC6N monolayer with a circle pore decrease when the pore size increases. The trending laws are also shown to predict the mechanical parameters of the porous BC6N sheet. Besides, the thermal conductivity (TC) of the monolayer BC6N with a smaller pore size is much higher than that of the sample with a bigger pore size due to higher energy loss. Changing pore size can adjust the TC and 2D temperature distribution around the circle-pore. Furthermore, the pristine BC6N TC is more heightened than the defective ones. The study results confirm that 2D-BC6N has outstanding optical, electronic, thermal transport, and mechanical characteristics. It is particularly inspiring for further experimental and theoretical works that advance different technological applications, such as gas sensors, nanoelectronics, and optoelectronic devices. It helps for real-time monitoring, which is critical to improving safety and efficiency in marine and air transportation infrastructure. Keywords: 2D-BC6N, mechanical, thermal transport, energy.

A molecular dynamics simulation study of thermal transport in hydrazinium cyclo-pentazolate

Cyclo-pentazolate salts (CPSs) as a new type of high-energy-density materials (HEDMs) with high nitrogen content have attracted considerable research attention. In contrast to the extensive studies on their energy properties, the thermal transport process in CPSs has been less studied which relates closely to the thermal safety of this material. Concerning the hydrazinium cyclo-pentazolate (HCP), we conduct a computational study to estimate the thermal conductivity of HCP by means of the non-equilibrium molecular dynamics (NEMD) simulation. To achieve that, we have customized interaction parameters based on the default OPLS force field for the HCP, as benchmarked by its crystal structure. Our simulations have revealed surprisingly anisotropic thermal conductivity for the HCP, while the thermal conductivity becomes highest roughly in the direction perpendicular to cyclo-N 5− with its value notably higher than common high explosives. By modulating the interaction parameters within the HCP molecule, we have further captured a dominant role of the interaction between cyclo-N 5− in regulating the thermal conductivity of the HCP. The attractive Lennard–Jones (LJ) potential may restrict the relative motion between cyclo-N 5− , which forms a long-range order thus enhances thermal transport in the direction perpendicular to cyclo-N 5− . Our simulations result on the effect of (cyclo-N 5− )-(cyclo-N 5− ) interaction provide insights to engineer thermal transport in CPSs at the molecular level.

Binding Sites of Bicarbonate in Phosphoenolpyruvate Carboxylase.

Phosphoenolpyruvate carboxylase (PEPC) is used in plant metabolism for fruit maturation or seed development as well as in the C4 and crassulacean acid metabolism (CAM) mechanisms in photosynthesis, where it is used for the capture of hydrated CO2 (bicarbonate). To find the yet unknown binding site of bicarbonate in this enzyme, we have first identified putative binding sites with nonequilibrium molecular dynamics simulations and then ranked these sites with alchemical free energy calculations with corrections of computational artifacts. Fourteen pockets where bicarbonate could bind were identified, with three having realistic binding free energies with differences with the experimental value below 1 kcal/mol. One of these pockets is found far from the active site at 14 Å and predicted to be an allosteric binding site. In the two other binding sites, bicarbonate is in direct interaction with the magnesium ion; neither sequence alignment nor the study of mutant K606N allowed to discriminate between these two pockets, and both are good candidates as the binding site of bicarbonate in phosphoenolpyruvate carboxylase.

Identification of Potential Inhibitors of the SARS-CoV-2 NSP13 Helicase via Structure-Based Ligand Design, Molecular Docking and Nonequilibrium Alchemical Simulations.

We have assembled a computational pipeline based on virtual screening, docking techniques, and nonequilibrium molecular dynamics simulations, with the goal of identifying possible inhibitors of the SARS-CoV-2 NSP13 helicase, catalyzing by ATP hydrolysis the unwinding of double or single-stranded RNA in the viral replication process inside the host cell. The druggable sites for broad-spectrum inhibitors are represented by the RNA binding sites at the 5' entrance and 3' exit of the central channel, a structural motif that is highly conserved across coronaviruses. Potential binders were first generated using structure-based ligand techniques. Their potency was estimated by using four popular docking scoring functions. Common docking hits for NSP13 were finally tested using advanced nonequilibrium alchemical techniques for binding free energy calculations on a high-performing parallel cluster. Four potential NSP13 inhibitors with potency from submicrimolar to nanomolar were finally identified.

A Comparative Computational Study of the Solidification Kinetic Coefficients for the Soft-Sphere BCC-Melt and the FCC-Melt Interfaces

Using the non-equilibrium molecular dynamics (NEMD) simulations and the time-dependent Ginzburg–Landau (TDGL) theory for solidification kinetics, we study the crystal-melt interface (CMI) kinetic coefficients for both the soft-sphere (SS) BCC-melt and the FCC-melt interfaces, modeled with the inverse-power repulsive potential (n=8). The collective dynamics of the interfacial liquids at four equilibrium CMIs are calculated and employed to eliminate the discrepancy between the predictions of the kinetic coefficient using the NEMD simulations and the TDGL solidification theory. The speedup of the two modes of the interfacial liquid collective dynamics (at wavenumbers equal to the principal and the secondary reciprocal lattice vector of the grown crystal) at the equilibrium FCC CMI is observed. The calculated local collective dynamics of the SS BCC CMIs are compared with the previously reported data for the BCC Fe CMIs, validating a hypothesis proposed recently that the density relaxation times of the interfacial liquids at the CMIs are anisotropic and material dependent. With the insights provided by the improved application of the TDGL solidification theory, an attempt has been made to interpret the variation physics of the crystal-structure dependence of the solidification kinetic coefficient.

Molecular Dynamics Simulations of Thermal Transport of Carbon Nanotube Interfaces

In this paper, non-equilibrium molecular dynamics simulations are used to study the interfacial heat exchange capacity of one-dimensional carbon nanotube nested structures. When the radius of the CNT substrate is increased from 1.356 to 2.712 nm, the ITC has a great enhancement from 1.340 to 2.949 nw/k. After this, we investigate the effects of overlap length, CNT length, and van der Waals interaction strength on the thermal resistance of the interface between carbon nanotubes. Firstly, we found that the nesting depth can significantly increase the ITC, and the increase in ITC is more obvious at an overlap length of 40 Å than at 30 Å. After this, the effect of length on the interfacial thermal conductivity is investigated, and the interfacial thermal conductivity is enhanced by 33.8% when the length is increased to 30 nm. Finally, the effect of van der Waals interaction strength was investigated, and the ITC increased from 1.60 nW/K to 2.71 nW/K when the scale factor was increased from 1 to 2.

Mechanism of enhanced thermal conductivity of hybrid nanofluids by adjusting mixing ratio of nanoparticles

At certain mixing ratios, hybrid nanofluids exhibit superior thermal conductivity performance compared with mono nanofluids, without significant enhancements in viscosity. Studying the microstructures formed by nanoparticles in the base fluid is key to revealing the mechanism of enhanced heat transfer in nanofluidic systems. Here, effects of varying the nanoparticle volume fraction (0.5–2.0 vol%) and Cu:Al volume mixing ratio (20:80, 40:60, 50:50, 60:40, and 80:20) of Cu-Al/Ar hybrid nanofluids were investigated using non-equilibrium molecular dynamics (NEMD) simulations at a temperature of 85 K. The mechanism of thermal conductivity enhancement was elucidated through quantitative analysis of the nanolayer structure and atom trajectory. The results showed that the thermal conductivities of Cu-Al/Ar hybrid nanofluids with a mixing ratio of 20:80 were lower than those of the corresponding Cu/Ar mono nanofluids with volume fraction below 1.3 vol% due to the low quantity of Cu nanoparticles. However, when the quantity of Cu nanoparticles was increased by using mixing ratios from 40:60 to 80:20, the thermal conductivities were higher at all volume fractions. Compared to Cu/Ar nanofluids, the thermal conductivity of Cu-Al/Ar hybrid nanofluids improved from 0.149 to 0.166 W/(m·K) which was higher 2.76 % to 14.48 % as the mixing ratio increased from 20:80 to 80:20 at 2.0 vol%. Radial distribution function (g(r)), the number density and diffusion rate calculations indicated that the density and velocity of Ar atoms around Cu atoms were higher than those of Al atoms due to stronger interactions between Cu-Ar atoms. The main reason for the enhanced thermal conductivity of hybrid nanofluids with increasing mixing ratio was the higher unit nanolayer density and Brownian motion velocity.