compute_hminus_bf_opacity

exo_skryer.opacity_special.compute_hminus_bf_opacity(state: Dict[str, Array], opac: Dict[str, Array], params: Dict[str, Array]) → Array

Compute H⁻ bound-free continuum mass opacity from the special registry.

This function uses the precomputed H⁻ bound-free cross-section table and applies the (n / ρ) normalization appropriate for a single-absorber continuum term.

Parameters:

statedict[str, ndarray]

Atmospheric state dictionary containing:

• wlndarray, shape (nwl,)

  Forward-model wavelength grid in microns.

• T_layndarray, shape (nlay,)

  Layer temperatures in Kelvin.

• nd_layndarray, shape (nlay,)

  Layer total number density in cm⁻³.

• rho_layndarray, shape (nlay,)

  Layer mass density in g cm⁻³.

• vmr_laydict[str, ndarray]

  Volume mixing ratios per species. Must include "H-" to enable this term. Values may be scalars or arrays with shape (nlay,).

• nlayint

  Number of atmospheric layers.

• nwlint

  Number of wavelength points.

paramsdict[str, ndarray]

Parameter dictionary (unused; kept for API compatibility).

Returns:

kappa_hminus_bfndarray, shape (nlay, nwl)

H⁻ bound-free continuum mass opacity in cm² g⁻¹. Returns zeros when the special registry is not loaded or when state["vmr_lay"] does not provide an "H-" mixing ratio.