Zarif Ahmad Razin Bhuiyan

Materials with exceptionally low thermal conductivity are desirable for thermal insulation and waste heat recovery. While foams and aerogels boast ultralow thermal conductivities akin to air, lack of mechanical stiffness in these soft materials necessitates a paradigm shift in materials design that can offer thermal insulation and mechanical rigidity simultaneously. Here, we show that spun-cast layered hybrid organic-inorganic perovskite thin films, azobenzene ethyl ammonium lead iodides, exhibit a record-low thermal conductivity, down to ∼0.04 watts per meter per kelvin at room temperature, while maintaining mechanical rigidity with an elastic modulus of 7.7 gigapascals that surpasses that of most plastics, foams, and aerogels. This unusual combination of ultralow thermal conductivity and high mechanical rigidity is attributed to the specially engineered organic cations in the layered structure. Our finding highlights the potential of molecular engineering in hybrid layered structures to push the extreme of thermal insulation in dense, rigid solids.

Achieving low thermal conductivity in nonporous polymer materials without compromising mechanical integrity remains a longstanding challenge. Conventional strategies, such as introducing porosity, are inherently limited in dense systems. Here, we demonstrate that thermal transport in polymer–organic filler hybrids is closely associated with the engineered availability of thermally accessible vibrational modes, rather than individual interfacial resistance or porosity. Using PU/organic filler hybrids as a model system, we show that incorporating rigid organic motifs shifts the vibrational density of states toward higher frequencies, where vibrational modes are weakly thermally populated at ambient conditions. This limits the number of heat-carrying channels, leading to suppressed thermal conductivity despite negligible individual interfacial resistance (~10⁻⁹ m² K W⁻¹). These findings provide initial evidence that vibrational mode engineering may offer a promising molecular-level route for reducing thermal conductivity in dense polymer-organic hybrid systems.

Molecular dynamics (MD) simulations were conducted to investigate the impact of hybrid wettability on thermal properties, specifically focusing on the constant volume molar heat capacity and thermal conductivity of Nano-confined liquid (NCL). The simulation domain, maintained at a temperature of 100 K, consisted of a $\sim 3\text{ nm}$ thin film of liquid argon entrapped between two solid copper surfaces with hybrid wettability. Hybrid wettability surfaces were produced by varying the solid–liquid interaction parameter and applying two different wettability factors ($\mu$) to the same surface. In this study, key findings pertaining to the influence of hybrid wettability on heat capacity and thermal conductivity include: (i) The heat capacity of liquid confined within hybrid wettability surfaces surpasses the heat capacity of the liquid in its bulk form. …

Project Overview

When liquids are confined within nanoscale gaps such as in advanced micro/nanochannel heat sinks, thermal interface materials, and nanofluidic cooling devices, their thermophysical behavior deviates fundamentally from bulk fluid dynamics. Interfacial effects, solid-liquid interaction potentials, and atomic density layering dominate the energy storage and transport characteristics of the entrapped fluid.

While previous studies focused primarily on uniformly hydrophilic or uniformly hydrophobic surfaces, practical engineering surfaces inherently exhibit chemical heterogeneity. This research investigated the thermal properties of nano-confined liquid (NCL) entrapped between solid copper surfaces engineered with hybrid wettability, utilizing equilibrium molecular dynamics (EMD) simulations performed with LAMMPS.

Key Scientific Findings

1. Enhanced Specific Heat Capacity: Demonstrated that the constant volume molar heat capacity ($C_v$) of liquid argon entrapped between hybrid wettability surfaces is significantly higher than its bulk counterpart ($20\text{ J/(mol}\cdot\text{K)}$).
2. Tunability via Surface Fraction: Identified how the hydrophilic area fraction ($f$) governs specific heat capacity, thermal conductivity, interfacial fluid layering, localized mass density profiles, and total system energy fluctuations.

💻 Open-Source LAMMPS Code & GitHub Repository

To promote reproducible computational research, the complete LAMMPS simulation script, parametric loops, and data post-processing routines are openly shared on GitHub:

👉 GitHub Repository: zarifrazin/ncl_hybrid_wettability

The repository includes:

• in.ncl_hybrid: The master LAMMPS input deck containing the system initialization, lattice creation, region partitioning, and velocity distribution.
• Parametric Loops: Automated nested scripts looping over wettability combinations ($w_I, w_O$) and area fractions ($f$).
• Property Calculations: Fluctuation formulas for calculating $C_v$, mean squared displacement (MSD) for diffusion, radial distribution functions (RDF), and density profiles across the channel.
• Log & Output Parsers: Guidance on parsing thermodynamic outputs and extracting publication-grade datasets.

📊 Presentation Slides

Below are the presentation slides from my talk on this research, outlining the molecular dynamics framework, interfacial phenomena, and thermal property results: