redback.transient_models.afterglow_models.vegas_models.vegas_tophat

redback.transient_models.afterglow_models.vegas_models.vegas_tophat(time, redshift, thv, loge0, thc, lognism, loga, p, logepse, logepsb, g0, **kwargs)[source]

VegasAfterglow tophat jet with unified medium (ISM/Wind/Hybrid)

High-performance C++ implementation providing ~1ms per light curve evaluation.

Parameters:

  • time – time in days in observer frame

  • redshift – source redshift

  • thv – viewing angle in radians

  • loge0 – log10 isotropic equivalent energy [erg]

  • thc – jet core opening angle in radians

  • lognism – log10 ISM number density [cm^-3] (use small number e.g., -20) for pure wind)

  • loga – log10 wind parameter A_* [g/cm] (use small number e.g., -20 for pure ISM)

  • p – electron power-law index

  • logepse – log10 electron energy fraction

  • logepsb – log10 magnetic field energy fraction

  • g0 – initial Lorentz factor

  • kwargs – Additional parameters - frequency: frequency array in Hz (required for flux_density mode) - output_format: ‘flux_density’ or magnitude - bands: photometric bands for magnitude calculation - cosmology: astropy cosmology object (default: Planck18) - phiv: azimuthal viewing angle [rad] (default: 0.0) - xie: electron acceleration efficiency (default: 1.0) - spreading: enable jet spreading (default: False) - duration: shell duration [s] (default: 1.0) - wind_k: wind power-law index (default: 2.0) - wind_n0: wind transition radius density [cm^-3] (default: inf) - reverse_shock: enable reverse shock (default: False) - reverse_logepse: reverse shock electron fraction (default: logepse) - reverse_logepsb: reverse shock magnetic fraction (default: logepsb) - reverse_p: reverse shock electron index (default: p) - reverse_xie: reverse shock electron efficiency (default: 1.0) - ssc: enable synchrotron self-Compton (default: False) - cmb_cooling: enable CMB cooling (default: False) - kn: enable Klein-Nishina corrections (default: False) - resolutions: (phi_res, theta_res, time_res) tuple (default: (0.3, 1, 10)) - rtol: relative tolerance for ODE solver (default: 1e-6) - axisymmetric: assume axisymmetric jet (default: True) - magnetar_L0: magnetar luminosity [erg/s] (optional) - magnetar_t0: magnetar spin-down time [s] (optional) - magnetar_q: magnetar braking index (optional)

Returns:

flux density [mJy] or AB magnitude

Examples:

Pure ISM: lognism=0.0, loga=-20 Pure Wind: lognism=-20, loga=11.0 Hybrid: lognism=0.0, loga=11.0 (wind with ISM floor)