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

# Create the solvent.
solvent = CosmoRealSolvent(
    298*Kelvin,
    water,
    parameters=CosmoRSParameters(),
)

# Calculate chemical potential of a methanol molecule in solvent water.
chemical_potential = solvent.chemicalPotential(methanol)
nlprint(f'The chemical potential of dissolving methanol in water is {chemical_potential}')

# Get the sigma potential of water.
charge_densities, sigma_potential = solvent.sigmaPotential()

# Get the sigma potential plot of water
plot = solvent.sigmaPotentialPlot()
Plot.show(plot)