Structures and stability of benzoic acid dimer ((BA) 2 ) and benzoic acid–water (BA–H 2 O) m:n complexes, with m and n=1, 2, were studied in benzene solutions, using molecular dynamics (MD) simulations. It appeared that nearly all hydrogen bond (H-bond) complexes suggested from different partition experiments existed in MD simulations, with the probability depending on their size and temperature. The MD results revealed the probability of finding H-bonds between water molecules, as well as the non-self-association of water molecules in benzene solutions. Although the H-bonds in (BA) 2 are quite strong in the gas phase and pure benzene, they can be opened by water molecules, forming microhydrates in benzene solutions. It was shown that, in order to provide insights into the structures and stability of the BA–H 2 O complexes in benzene solution, solvent molecules as well as dynamic and temperature effects have to be included in theoretical investigation.