# Load the COSMO species
database = CosmoRSSpeciesDatabase()
water = database.exportSpecies('water')

# Create the gas.
gas = CosmoRealGas(
    298*Kelvin,
    parameters=CosmoRSParameters(),
)

# Calculate chemical potential of a methanol molecule in solvent water.
chemical_potential = gas.chemicalPotential(water)
nlprint(f'The chemical potential of gas-phase water with respect to ideal screened state is {chemical_potential}')