# Create the custom CosmoRSParameters object
parameters = CosmoRSParameters(
    hydrogen_bond_donors=Hydrogen,
    hydrogen_bond_acceptors=(Nitrogen, Oxygen, Fluorine),
    r_average=0.334*Ang,
    alpha_prime=1428.3*kiloCaloriePerMol*Ang**2/elementary_charge**2,
    f_correlation=2.784,
    c_hydrogen_bonding=8800*kiloCaloriePerMol*Ang**2/elementary_charge**2,
    sigma_hydrogen_bonding=0.869*elementary_charge/nm**2,
    hydrogen_bond_temperature=False,
    a_effective=7.79*Ang**2,
    lambda_comb=0.13,
    omega_ring=-0.216*kiloCaloriePerMol,
    eta_gas=-10.05,
    lambda_reg=0.611,
    lambda_ortho=0.820,
    pka_scale=1.0,
    pka_offset=0.0,
    dispersion={
        Hydrogen: -3.59*kiloCaloriePerMol/nm**2,
        Carbon: -3.47*kiloCaloriePerMol/nm**2,
        Nitrogen: -2.26*kiloCaloriePerMol/nm**2,
        Oxygen: -3.17*kiloCaloriePerMol/nm**2,
        Chlorine: -5.00*kiloCaloriePerMol/nm**2,
    },
)

# Save the parameters to the database
database = CosmoRSParameterSets()
database.saveParameterSet(
    'custom_parameter_set',
    parameters,
    reference='DOI: 10.1139/V09-008',
    description=(
        'Taken from "C. C. Pye, T. Zeigler, E. van Lenthe and J. L. Louwen. '
        'Can. J. Chem., 87, 790-797 (2009).'
    )
)

# Print the names of the available parameter sets
parameter_names = database.parameterSets()
nlprint('Available parameter set names:')
for name in parameter_names:
    nlprint(name)
nlprint('')

# Remove the entry from the database again
database.deleteParameterSet('custom_parameter_set')