Source author record

Debashish Mukherji

Debashish Mukherji appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2025arXiv

Thermal conductivity of commodity polymers under high pressures

Understanding the thermal conductivity of polymers under high-pressure conditions is essential for a range of applications, from aerospace and deep-sea engineering to common lubricants. However, the complex relationship between pressure, $P$, the thermal transport coefficient, $κ$, and polymer architecture poses substantial challenges to both experimental and theoretical investigations. In this work, we study the pressur-dependent thermal transport properties of a widely used commodity polymer -- poly(methyl methacrylate) (PMMA) -- using a combination of all-atom molecular dynamics simulations and semi-analytical approaches. While we report both classical and quantum-corrected estimates of $κ$, the latter approach reveals that as the pressure increases from 1 atm to 10 GPa, $κ$ rises by up to a factor of four -- from 0.21 W m$^{-1}$ K$^{-1}$ to 0.80 W m$^{-1}$ K$^{-1}$. To better understand the mechanisms behind this increase, we disentangle the contributions from bonded and nonbonded monomer interactions. Our analysis shows that nonbonded energy-transfer rates increase by a factor of six over the pressure range, while bonded interactions show a more modest increase -- about a factor of three. This observation further consolidates the fact that the nonbonded interactions play the dominant role in dictating the microscopic heat flow in polymers. These individual energy-transfer rates are also incorporated into a simplified heat diffusion model to predict $κ$. The results obtained from different approaches show internal consistency and align well with available experimental data. Additionally, some data for polylactic acid (PLA) are presented.

preprint2022arXiv

Comparison of all atom and united atom models for thermal transport calculations of amorphous polyethylene

Polymer simulations routinely employ models with different molecular resolutions. United atom (UA) models are one such example, where groups of certain atoms in a molecule are clustered into superatoms. Although their computational simplicity makes them particularly attractive for studying a wide range of polymer properties, the missing degrees of freedom in UA models can impact certain properties that are intimately linked to localized vibrations, such as the heat capacity and the thermal transport coefficient $κ$. In contrast, the numerically exhausting all atom (AA) models produce results that better match experimental data. In this work, we systematically investigate and compare $κ$ obtained from an AA and a UA models for an amorphous polyethylene system. The results indicate that the UA description may not be a suitable model for evaluating thermal transport, since it underestimates $κ$ in comparison to an AA description and the experimental value. The coarse-graining leads to the softer interactions and its presence is highlighted in a weaker mechanical response from the UA model, thus also underestimates $κ$. We further consolidate our findings by extracting the bonded and the nonbonded contributions to $κ$ within the framework of the single chain energy transfer model.

preprint2022arXiv

Computational Indentation in Highly Cross-linked Polymer Networks

Indentation is a common experimental technique to study the mechanics of polymeric materials. The main advantage of using indentation is because this provides a direct correlation between the microstructure and the small-scale mechanical response, which is otherwise difficult within the standard tensile testing. Here, majority of studies have investigated hydrogels, microgels and/or elastomers. However, a lesser investigated system is the indentation in highly cross-linked polymer (HCP) networks, where the complex network structure plays a key role in dictating their physical properties. In this work, we investigate the structure-property relationship in HCP networks using the computational indentation of a generic model. We establish a correlation between the local bond breaking, the network rearrangement, and the small-scale mechanics. The results are compared with the elastic-plastic deformation model. HCPs harden upon indentation.

preprint2021arXiv

Correlating thermodynamics, morphology, mechanics and thermal transport in PMMA-PLA blends

Thermodynamics controls structure, function, stability and morphology of polymer blends. However, obtaining the precise information about their mixing thermodynamics is a challenging task, especially when dealing with complex macromolecules. This is partially because of a delicate balance between the local concentration/composition fluctuations and the monomer level (multi-body) interactions. In this context, the Kirkwood-Buff (KB) theory serves as a useful tool that connects the local pairwise fluid structure to the mixing thermodynamics. Using larger scale molecular dynamics simulations, within the framework of KB theory, we investigate a set of technologically relevant poly(methyl methacrylate)-poly(lactic acid) (PMMA-PLA) blends with the aim to elucidate the underlying microscopic picture of their phase behavior. Consistent with these experiments, we emphasize the importance of properly accounting for the entropic contribution, to the mixing Gibbs free energy change $Δ{\mathcal G}_{\rm mix}$, that controls the phase morphology. We further show how the relative microscopic interaction details and the molecular level structures between different mixing species can control the non-linear mechanics and ductility. As a direct consequence, we provide a correlation that links thermodynamics, phase behavior, mechanics, and thus also thermal transport in polymer blends. Therefore, this study provides a guiding principle for the design of light weight functional materials with extraordinary physical properties.

preprint2020arXiv

Reduced thermal conductivity in molecular forests

Heat propagation in quasi-one dimensional materials (Q1DMs) often appears paradoxical. While an isolated Q1DM, such as a nanowire, carbon nanotube, or polymer, can exhibit a high thermal conductivity \k{appa}, forests of the same materials show a reduction in \k{appa}. Here, the complex structures of these assemblies have hindered the emergence of a clear molecular picture of this intriguing phenomenon. We combine multiscale (coarse-grained) simulation with the concepts known from polymer physics and thermal transport to unveil a generic (microscopic) picture of \k{appa} reduction in molecular forests. We show that a delicate balance between the bond orientations, the persistence length of the Q1DM and the flexural vibrations govern the knock-down of \k{appa}.

preprint2020arXiv

Tuning thermal transport in highly cross-linked polymers by bond induced void engineering

Tuning the heat flow is fundamentally important for the design of advanced functional materials. Here, polymers are of particular importance because they provide different pathways for the energy transfer. More specifically, the heat flow between two covalently bonded monomers is over 100 times faster than between the two non-bonded monomers interacting via van der Waals (vdW) forces. Therefore, the delicate balance between these two contributions often provide a guiding tool for the tunability in thermal transport coefficient k of the polymeric materials. Traditionally most studies have investigated k in the linear polymeric materials, the recent interests have also been directed towards the highly cross-linked polymers (HCP). In this work, using the generic molecular dynamics simulations we investigate the factors effecting k of HCP. We emphasize on the importance of the cross-linking bond types and its influence on the network microstructure with a goal to provide a guiding principle for the tunability in k. While these simulation results are discussed in the context of the available experimental data, we also make predictions.

preprint2020arXiv

Why Do Elastin-Like Polypeptides Possibly Have Different Solvation Behaviors in Water-Ethanol and Water-Urea Mixtures?

The solvent quality determines the collapsed or the expanded state of a polymer. For example, a polymer dissolved in a poor solvent collapses, whereas in a good solvent it opens up. While this standard understanding is generally valid, there are examples when a polymer collapses even in a mixture of two good solvents. This phenomenon, commonly known as co-non-solvency, is usually associated with smart polymers. Moreover, recent experiments have shown that the elastin-like polypeptides (ELPs) show co-non-solvency behavior in aqueous-ethanol mixtures. In this study, we investigate the phase behavior of ELPs in aqueous binary mixtures using molecular dynamics simulations of all-atom and complementary explicit solvent generic models. The model is parameterized by mapping the solvation free energy obtained from the all-atom simulations onto the generic interaction parameters. For this purpose, we derive segment based generic parameters for four different peptides, namely proline (P), valine (V), glycine (G) and alanine (A). Here we compare the conformational behavior of two ELP sequences, namely VPGGG and VPGVG, in aqueous-ethanol and -urea mixtures. Consistent with recent experiments, we find that ELPs show co-non-solvency in aqueous-ethanol mixtures. Ethanol molecules have preferential binding with all ELP residues and thus driving the coil-to-globule transition. On the contrary, ELP conformations show weak variation in aqueous-urea mixtures. Our simulations suggest that the glycine residues dictate the overall behavior of ELPs in aqueous-urea, where urea molecules have a rather weak preferential binding with glycine, i.e., less than kT. While the validation of the latter findings will require more detailed experimental investigation, the results presented here may provide a new twist to the present understanding of cosolvent interactions with peptides and proteins.

preprint2019arXiv

Tuning morphology and thermal transport of asymmetric smart polymer blends by macromolecular engineering

A grand challenge in designing polymeric materials is to tune their properties by macromolecular engineering. In this context, one of the drawbacks that often limits broader applications under high temperature conditions is their poor thermal conductivity $κ$. Using molecular dynamics simulations, we establish a structure-property relationship in hydrogen bonded polymer blends for possible improvement of $κ$. For this purpose, we investigate two experimentally relevant hydrogen bonded systems -- one system consists of short poly({N}-acryloyl piperidine) (PAP) blended with longer chains of poly(acrylic acid) (PAA) and the second system is a mixture of PAA and short poly(acrylamide) (PAM) chains. Simulation results show that PAA-PAP blends are at the onset of phase separation over the full range of PAP monomer mole fraction $ϕ_{PAP}$, which intensifies even more for $ϕ_{PAP} > 0.3$. While PAA and PAP interact with preferential hydrogen bonding, phase separation is triggered by the dominant van der Waals attraction between the hydrophobic side groups of PAP. However, if PAP is replaced with PAM, which has a similar chemical structure as PAP without the hydrophobic side group, PAA-PAM blends show much improved solubility. Better solubility is due to the preferential hydrogen bonding between PAA and PAM. As a result, PAM oligomers act as cross-linking bridges between PAA chains resulting in a three dimensional highly cross-linked network. While $κ$ for PAA-PAP blends remain almost invariant with $ϕ_{PAP}$, PAA-PAM systems show improved $κ$ with increasing PAM concentration and also with respect to PAA-PAP blends. Consistent with the theoretical prediction for the thermal transport of amorphous polymers, we show that $κ$ is proportional to the materials stiffness, i.e., the bulk modulus K and sound velocity v of PAA-PAM blends.

preprint2016arXiv

Why does high pressure destroy co-non-solvency of PNIPAm in aqueous methanol?

It is well known that poly(N-isopropylacrylamide) (PNIPAm) exhibits an interesting, yet puzzling, phenomenon of co-non-solvency. Co-non-solvency occurs when two competing good solvents for PNIPAm, such as water and alcohol, are mixed together. As a result, the same PNIPAm collapses within intermediate mixing ratios. This complex conformational transition is driven by preferential binding of methanol with PNIPAm. Interestingly, co-non-solvency can be destroyed when applying high hydrostatic pressures. In this work, using a large scale molecular dynamics simulation employing high pressures, we propose a microscopic picture behind the suppression of the co-non-solvency phenomenon. Based on thermodynamic and structural analysis, our results suggest that the preferential binding of methanol with PNIPAm gets partially lost at high pressures, making the background fluid reasonably homogeneous for the polymer. This is consistent with the hypothesis that the co-non-solvency phenomenon is driven by preferential binding and is not based on depletion effects.

preprint2015arXiv

Co-non-solvency: Mean-field polymer theory does not describe polymer collapse transition in a mixture of two competing good solvents

Smart polymers are a modern class of polymeric materials that often exhibit unpredictable behavior in mixtures of solvents. One such phenomenon is co-non-solvency. Co-non-solvency occurs when two (perfectly) miscible and competing good solvents, for a given polymer, are mixed together. As a result, the same polymer collapses into a compact globule within intermediate mixing ratios. More interestingly, polymer collapses when the solvent quality remains good and even gets increasingly better by the addition of the better cosolvent. This is a puzzling phenomenon that is driven by strong local concentration fluctuations. Because of the discrete particle based nature of the interactions, Flory-Huggins type mean field arguments become unsuitable. In this work, we extend the analysis of the co-non-solvency effect presented earlier [Nature Communications 5, 4882 (2014)]. We explain why co-non-solvency is a generic phenomenon that can be understood by the thermodynamic treatment of the competitive displacement of (co)solvent components. This competition can result in a polymer collapse upon improvement of the solvent quality. Specific chemical details are not required to understand these complex conformational transitions. Therefore, a broad range of polymers are expected to exhibit similar reentrant coil-globule-coil transitions in competing good solvents.