cubos de datos IFU basados ​​en longitud de onda VS cubos de datos IFU basados ​​en velocidad

Tengo 2 cubos de datos IFU simulados de la misma observación, digamos A y B. La tercera dimensión del cubo A está en longitudes de onda. Para crear el cubo BI, simplemente vuelva a agrupar el cubo A en longitudes de onda logarítmicas y luego convierta eso a velocidad (km/s). El rebinning a longitud de onda logarítmica fue necesario para convertir la tercera dimensión del cubo B en km/s. Para una explicación ver esto:

Conversión de la dimensión espectral de ångström a la velocidad de la galaxia (km/s)

Ahora quiero aplicar la ampliación instrumental en ambos cubos. Para hacer eso, para cada cubo, tengo que crear un perfil gaussiano con FWHM igual a la resolución del instrumento (o Función de dispersión de línea, o LSF) y luego convolucionar eso con cada espectro del cubo. Por supuesto, la LSF del cubo A estará en longitudes de onda, mientras que la LSF del cubo B estará en km/s. Aquí está mi pregunta ahora:

  • Dado que todos los desplazamientos de píxeles en el cubo A no corresponden al mismo valor de km/s, uno de los dos LSF no será un gaussiano perfecto, ¿verdad? ¿Debo asumir que sería la velocidad LSF?

  • ¿Cuál es la forma correcta de crear la velocidad LSF? ¿Simplemente creo un LSF de longitud de onda y luego lo vuelvo a agrupar en km/s, de la misma manera que volví a agrupar el cubo A en el cubo B?

  • ¿Vale la pena el esfuerzo? ¿Hará la diferencia? He visto artículos revisados ​​por pares que usan LSF de velocidad que son Gaussianos perfectos. ¿Es esto porque la diferencia es insignificante o todo mi pensamiento/pregunta es totalmente incorrecto?

Respuestas (1)

Tu lógica parece completamente correcta. Hagas lo que hagas es una aproximación, la diferencia entre un perfil gaussiano en el cubo B y cómo "debería" verse será pequeña para cualquier resolución sensible. Por sensato, me refiero a si su eje de longitud de onda/velocidad cubre muchos elementos de resolución.

Sin embargo, hay una forma "correcta". Generar un gaussiano en el espacio de longitud de onda F ( λ ) (si crees que eso es lo que es tu LSD, ¿dónde λ aquí hay una abreviatura para λ λ 0 , la separación de la longitud de onda central) y luego se convierte en espacio de velocidad conservando el flujo:

F ( λ )   d λ = gramo ( Iniciar sesión λ )   d Iniciar sesión λ
gramo ( Iniciar sesión λ ) = λ F ( λ )

¿Quizás una mejor manera de hacer esto es obtener algunos datos, como una calibración de arco y ver cómo se ven las líneas sin resolver en su espacio de observación?