Probing Rate-Dependent Liquid Shear Viscosity Using Combined Machine Learning and Non-Equilibrium Molecular Dynamics

Accurately measuring liquid dynamic viscosity across a wide range of shear rates, from the linear-response to shear-thinning regimes, presents significant experimental challenges due to limitations in resolving high shear rates and controlling thermal effects. In this study, we integrated machine learning (ML) with non-equilibrium molecular dynamics (NEMD) simulations to address these challenges. A supervised artificial neural network (ANN) model was developed to predict viscosity as a function of shear rate, normal pressure, and temperature, effectively capturing the complex interplay among these variables. The model reveals distinct trends in shear viscosity, characterized by the shear-thinning exponent, and highlights non-monotonic behavior in the radius of gyration components, reflecting molecular morphological changes driven by rate-dependent volume expansion. Notably, temperature effects diminish at higher shear rates, where molecular alignment and spacing dominate the response to shear. By implementing the 'fix npt/sllod' command in LAMMPS, we achieve precise constant-pressure control in NEMD simulations, ensuring accurate representation of system dynamics. This study demonstrates the potential of ML-enhanced NEMD for efficient and accurate viscosity prediction, providing a robust framework for future research in complex fluid dynamics and material design.