The ability to predict multiphase fluid transport in nanoporous rocks is critical for many geoscience applications. For instance, in applications such as geologic carbon sequestration, and nuclear waste disposal (Altman et al., 2014;Zheng et al., 2008), the long-term performance of the selected sites depends on the multiphase physical and chemical interactions in nanoporous caprock or clay-based barriers. Similarly, finding effective strategies for sustainable hydrocarbon production from shale requires understanding of multiphase processes in the nanoporous shale matrix (Alexander et al., 2011;Falk et al., 2015;Tokunaga et al., 2017). When pore sizes approach nanoscales, the impact of molecular interaction forces between fluids and solids becomes increasingly important. The fluid-solid interaction forces, also known as surface forces (e.g., van der Waals, electrostatic, structural forces) (Churaev, 2000;Israelachvili, 2011), result from electrostatic and electromagnetic fields generated by charges and oscillating molecular dipoles. These forces can alter macroscopic fluid phase behavior and control related processes such as sorption, wetting, and transport (e.g.,