Solvate structures of salicylic acid and its derivatives (acetylsalicylic acid, methyl salicylate, salicylamide, and sodium salicylate) in supercritical carbon dioxide modified by polar cosolvent (methanol, 0.03 mol fraction) at density of 0.7 g/cm(3) and temperature of 318 K were studied by means of molecular dynamics simulations. It was found that salicylic and acetylsalicylic acids form stable hydrogen-bonded complexes with methanol through hydrogen atom of the carboxyl group. In methyl salicylate molecule the carboxyl hydrogen is substituted by methyl radical, that is why strong hydrogen bonds with methanol are not typical for it. Salicylamide is able to form hydrogen bonds with methanol through several functional groups (preferably as proton donor), but their probability and lifetimes are lower than those for salicylic acid - methanol and acetylsalicylic acid - methanol hydrogen bonds. The continuous rearrangement of salicylamide solvation shell due to forming and breaking solute - cosolvent hydrogen bonds was observed. Sodium salicylate is in the form of undissociated ion pair in supercritical carbon dioxide, despite the presence of a polar co-solvent, but it is characterized by very high selective solvation. It is able to form hydrogen bonds with several (up to 6) methanol molecules simultaneously. Unlike other compounds considered, it forms highly stable hydrogen bonds as a proton acceptor. (C) 2016 Elsevier B.V. All rights reserved.