A theoretical framework is presented for modeling the chemomechanical behavior of multiphase porous media, in general, and unsaturated soils, in particular, which can address skeletal deformation, fluid flow, heat conduction, solute diffusion, chemical reaction, and phase transition in a consistent and systematic way. A general expression is derived for the electrochemical potential of a fluid species with explicitly accounting for the effects of osmosis, capillarity, and adsorption. The equilibrium behavior of porous media is investigated, and the composition of pore water pressure is identified. Explicit formulations are developed for the effective stress and intergranular stress, with consideration of physicochemical effects. It is shown that the negative water pressure measured by a conventional transducer can be significantly different than the true pore water pressure. It is also theoretically revealed that, other than the soil water characteristic function, a new pressure (or potential) function accounting for the physicochemical effects is generally required in analyzing the coupled chemomechanical processes in unsaturated soils. The new theory is capable of effectively explaining many salient phenomena occurring in water-saturated porous media with a degree of saturation varying from an extremely low value to 100%, including Donnan osmosis, capillary fringe, air entry value, initial hydraulic head during seepage, and pressure solution. The new theory can be used to analyze the multiple coupled physical and chemical processes in the vadose zone.Abbreviations: AEV, air entry value; PWP, pore water pressure; REV, representative volume.Chemomechanical behavior of porous media with multiphases and multispecies is of great interest in many diverse fields of science and engineering. To mention a few, this includes the industries of nuclear, hazardous, and municipal waste isolation; petroleum and gas extraction; technologies of methane gas hydrates exploitation; CO 2 sequestration; and weak soil reinforcement, landform stability assessment, structural material durability, and weathering of rock masses. Understanding, controlling, and predicting the long-term effect of physicochemical processes on the mechanical performance of geomaterials is becoming an indispensable part of environmental impact assessment and performance assessment analysis. Thus, there is a clear need to develop comprehensive chemohydromechanical mathematical modeling capacities that are able to address realistically the reactive multispecies transport, multiphase flow, chemical reaction, phase transition, chemically induced deformation, and other related physicochemical processes in deformable soils.The coupling of multiple processes in porous media, including skeletal deformation, seepage, diffusion, and heat conduction, has been extensively investigated, and the poromechanic theory that is capable of describing these coupled physical processes has been very well developed (e.g., Coussy, 2004). Thus far, however, the behavior of porous medi...