Distribution of solvent molecules near transition-metal complex is key information to comprehend the functionality, reactivity, and so forth. However, polarizable continuum solvent models still are the standard and conventional partner of molecularorbital (MO) calculations in the solution system including transition-metal complex. In this study, we investigate the conformation, hydration, and ligand substitution reac-, and [Ru(NO)(OH)(NO 2 ) 3 (H 2 O)] À (C) using a combination method of MO theories and a state-of-the-art molecular solvation technique (NI-MC-MOZ-SCF). A dominant species is found in the complex B conformers and, as expected, different between the solvent models, which reveals that molecular solvation beyond continuum media treatment are required for a reliable description of solvation near transition-metal complex. In the stability constant evaluation of ligand substitution reaction, an assumption that considers the direct association between the dissociated NO 2 À and complex C is useful to obtain a reliable stability constant.