Estimaci on del desplazamiento al rojo para las galaxias en el Southern H-ATLAS Regions Ks-band Survey - Trabajo de Fin de Grado

Página creada Beltran Nou
 
SEGUIR LEYENDO
Estimaci on del desplazamiento al rojo para las galaxias en el Southern H-ATLAS Regions Ks-band Survey - Trabajo de Fin de Grado
Trabajo de Fin de Grado

 Estimación del desplazamiento al
rojo para las galaxias en el Southern
 H-ATLAS Regions Ks-band Survey

         Universidad de La Laguna
           Facultad de Ciencias

         David González Álvarez

                   Tutores:
            Aurelio Carnero Rosell
             Helmut Dannerbauer

                  Julio 2021
Estimaci on del desplazamiento al rojo para las galaxias en el Southern H-ATLAS Regions Ks-band Survey - Trabajo de Fin de Grado
Contents
1 Introducción                                                                                                                  3
  1.1 La astronomı́a    en el   infrarrojo cercano      .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .    3
  1.2 SHARKS . .        . . .   . . . . . . . . . . .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .    4
  1.3 Redshift . . .    . . .   . . . . . . . . . . .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .    6
  1.4 Objetivos . .     . . .   . . . . . . . . . . .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .    8

2 Datos                                                                                                                          9
  2.1 SHARKS-DR1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .                                                   9
  2.2 DES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .                                             10
  2.3 2MASS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .                                             11

3 Validación fotométrica de los datos de SHARKS-DR1                                                                           12
  3.1 Selección de la muestra . . . . . . . . . . . . . . . . . . . . . . . . . . .                                            12
  3.2 Pruebas de calibración . . . . . . . . . . . . . . . . . . . . . . . . . . .                                             13

4 Photo-z de los datos de SHARKS-DR1                                                                                            18
  4.1 Lephare . . . . . . . . . . . . . . . . .         .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   18
      4.1.1 Definición de librerı́a SED . . .          .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   19
      4.1.2 Definición de filtros . . . . . . .        .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   19
      4.1.3 Estimación del photo-z . . . . .           .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   21
  4.2 Preparación de los datos . . . . . . . .         .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   21
      4.2.1 Fotometrı́a de la muestra . . . .           .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   22
      4.2.2 Separación estrella-galaxia . . .          .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   22
  4.3 Validación de los datos . . . . . . . . .        .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   22
      4.3.1 Métricas usadas . . . . . . . . .          .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   23
      4.3.2 Muestra de calibración . . . . .           .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   24
  4.4 Resultados y discusión . . . . . . . . .         .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   24
      4.4.1 Validación frente a Spec-z . . .           .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   25
      4.4.2 Comparación de resultados . . .            .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   25

5 Conclusiones                                                                                                                  31

A Parámetros de Lephare                                                                                                        34

                                              1
Estimaci on del desplazamiento al rojo para las galaxias en el Southern H-ATLAS Regions Ks-band Survey - Trabajo de Fin de Grado
Abstract

    This project focuses on obtaining the photometric redshift (photo-z) for galax-
ies in the first data release (DR1) of the Southern H-ATLAS Regions Ks Band
Survey (SHARKS). In particular, we will study the effect of incorporating the
Ks -band in the photo-z estimates.

    A limiting factor in obtaining the redshift by photometric means is the lack
of spectral coverage, especially when the SHARKS data are only in the Ks-band.
Because of this, we will make use of an astronomical catalog that covers the same
region of the sky in the optical range, the Dark Energy Survey (DES), in partic-
ular DES-DR2.

    Before proceeding to the photo-z calculation, we study the calibration of the
SHARKS-DR1 data. For this purpose, we will proceed with a near-infrared as-
tronomical catalog that will act as a reference, 2MASS. In this section of the
project, using a linear regression between the magnitudes of both catalogs, we
will measure the calibration difference between them. In this case, we expect only
to measure the difference between the calibration in the Vega and AB system.

    As for the photometric redshift, it is noteworthy that since this is the first
time that these are calculated with SHARKS data, we have focused on a subsam-
ple where we have selected only galaxies with a distinct SHARKS-DES detection.
For the selection of galaxies, we will use color-color space separation that oc-
curs when using infrared bands. The calculation of the photometric redshift has
been performed using the Lephare software. Before the execution, we must define
the library of galaxy model templates that we will use (SEDs, Spectral Energy
Distribution), here we will use a library of 8 templates. For their subsequent con-
volution, we will also introduce the filters with which we took the measurements,
i.e. those of the DECam and VIRCAM cameras. Then, we proceed to calculate
the magnitudes as a function of the redshift for the templates, and finally, we
perform a minimization of the reshift for the templates. Finally, a minimization
by the 2 fitting method is carried out to obtain the redshift and the template that
best fits the observed galaxy.

    To validate the results, we calculate statistical metrics of the bias and dis-
persion of the photo-z, concerning a small calibration sample with spectroscopic
redshift, limited to low redshift. Next, we compare the overall redshift distribu-
tions obtained with and without using the Ks -band and look at some examples
of outliers.

    In general, we find that using SHARKS Ks -band improves the photo-z esti-
mates, compared to the similar estimates using only DES data, through a decrease
in the number of catastrophic results (outliers).

                                        2
1         Introducción
Este estudio se centra en la estimación del desplazamiento al rojo (redshift) de las
galaxias detectadas por el Southern H-ATLAS Regions Ks -band Survey (SHARKS), en
particular, en los datos del primer data release (DR1). El proyecto SHARKS1 (Oteo,
I., et al., 2016) es un proyecto público del Observatorio Europeo del Sur (ESO, del
inglés European Southern Observatory), donde las observaciones se realizaron durante
1200 horas y provienen del telescopio de 4 metros de diámetro VISTA (Visible and
Infrared Survey Telescope for Astronomy, Sutherland et al. 2015) ubicado en el Obser-
vatorio Paranal, Desierto de Atacama, Chile, con la cámara infrarroja VIRCAM. Los
300 grados cuadrados que observa están cubiertos por el proyecto Herschel-ATLAS2 (H-
ATLAS, Eales et al. 2010), siendo este uno de los mayores proyectos del Observatorio
Espacial Herschel.

    Uno de los objetivos principales de SHARKS es la medida en el infrarrojo cercano
(NIR) de, al menos, el 90% de las galaxias de H-ATLAS hasta redshifts ≈ 3, de ahı́ la
importancia de medir el redshifts de las galaxias de SHARKS que confirmen este hito.
    Cabe señalar aquı́, la importancia de los eventos ocurridos en el Universo a redshift
≈ 2, cuando el Universo tenı́a una edad de aproximadamente 4.75 Giga años. Según las
observaciones, en esta época se produce la formación de la mayor parte de la masa estelar
en las galaxias (Dickinson et al., 2003), alcanzando el máximo de la tasa formación
estelar (Behroozi et al., 2013). A esta época se la llama “el mediodı́a cósmico” (cosmic
noon en inglés). Es más, la mayor parte de la formación estelar de esta época tuvo lugar
en “galaxias polvorientas”, donde la luminosidad infrarroja predomina por encima de
la luminosidad ultravioleta (Reddy et al., 2010, 2011). Es por esto la importancia de
observar en el infrarrojo para estudiar esta época, para comprender qué condujo a este
máximo de formación estelar a redshift ≈ 2 y cómo se relaciona con la cosmologı́a
subyacente.

1.1        La astronomı́a en el infrarrojo cercano
El rango infrarrojo del espectro se define, generalmente, entre longitudes de onda que
van de ≈ 700 nm hasta 1 mm. A su vez, por conveniencia, el infrarrojo se divide en
tres franjas denominadas: infrarrojo cercano (NIR, del inglés near-infrared) entre 700
y 2500 nm, infrarrojo medio (MIR) hasta 25 µm e infrarrojo lejano (FIR) hasta 1 mm.
    Si partimos del extremo del rango óptico y continuamos aumentando la longitud
de onda, las estrellas frı́as empezarán a ser cada vez más brillantes. Pasado el NIR,
a menores temperaturas, empezaremos a ver el brillo térmico dado a la presencia de
granos de polvo.

    En el rango de frecuencias NIR, es interesante destacar 3 cosas: medimos temper-
aturas efectivas de 700 a 740 and 4000 Kelvin, el polvo es transparente a la luz en estas
frecuencias y además, galaxias a altos redshifts siguen siendo brillantes aunque sean
muy tenues en el óptico.

    1
        http://research.iac.es/proyecto/sharks/pages/en/home.php
    2
        https://www.h-atlas.org/

                                              3
Figure 1: Representación de la transmisión relativa de la atmósfera frente a la longi-
tud de onda en el infrarrojo cercano, y el nombre que reciben las 3 ventanas princi-
pales, J, H y K. Imagen de https://webs.ucm.es/info/Astrof/users/jaz/IA/IA_
14_infrarrojo.pdf

    Desde el punto de vista observacional, hay que tener en cuenta, también, que la
atmósfera contribuye significativamente a la radiación NIR, por lo que únicamente
podemos observar desde la Tierra, en ciertas “ventanas” en el NIR. Figure 1 muestra
la transmisión relativa de la atmósfera en función de la longitud de onda.

     Entre estas, cabe destacar la ventana K, que es en la que observa SHARKS. Esta
ventana corresponde al rango 2.008 − 2.427 µm (λc = 2.2 µm). Diferentes filtros se han
definido para observar en ella, por ejemplos los filtros K, K’ y Ks (de K short). Esta
última es la usada por el telescopio VISTA y por el proyecto SHARKS. Este filtro es
muy parecido al filtro K, pero recortado tal que esté centrado en λc = 2.15 µm (2.002 −
2.390 µm).
     La astronomı́a infrarroja es una ciencia que se renueva dı́a a dı́a y gracias a ella,
se pueden realizar importantes estudios astrofı́sicos, tales como explorar regiones de
formación estelar (Ojha et al., 2004), el estudio del universo joven (Franx et al., 2003) o
el estudio de objetos frı́os en la Vı́a Láctea, como las enanas marrones (Carnero Rosell
et al., 2019). Sin embargo, la astronomı́a infrarroja añade algunas complicaciones, por
ejemplo, los dispositivos usados en astronomı́a óptica, como los CCDs, dejan de ser
sensibles en estos rangos electromagnéticos, y necesitamos detectores especiales.

1.2    SHARKS
En este contexto nació el proyecto SHARKS. Este proyecto, cuya finalidad es el estudio
de la formación y evolución de galaxias en el Universo temprano, es una colaboración
de institutos de todo el mundo (Australia, Reino Unido, Estados Unidos, . . . ) liderada

                                             4
Figure 2: Resumen de uno de los objetivos más importantes de SHARKS, que es
buscar contrapartidas a alto redshifts para galaxias del H-ATLAS (izquierda, Oteo, I.,
et al. 2016) y para galaxias en observatorios en radio, como EMU (derecha, Carnero,
comunicación privada). En estas imágenes, la sensibilidad se compara con VIKING,
survey muy parecido a SHARKS, pero de menor profundidad.

por el Instituto de Astrofı́sica de Canarias (IAC, Invesigador Principal: Dr. Helmut
Dannerbauer).

    Los campos de SHARKS se han elegido para que coincidan con tres campos de
H-ATLAS. Este es el estudio más amplio llevado a cabo por el Observatorio Espacial
Herschel y ha cartografiado aproximadamente 600 grados cuadrados en cinco bandas
fotométricas del infrarrojo lejano, desde 160 µm hasta 500 µm (Eales et al., 2010).
SHARKS observará todos los campos H-ATLAS ecuatoriales y meridionales excepto
GAMA 9, que está muy contaminado por la emisión de cirros galácticos.

    Los principales objetivos de SHARKS son: proporcionar la mejor identificación posi-
ble (Figura 2) para aproximadamente el 90% de las fuentes detectadas a 0 < z < 3 por
H-ATLAS, ASKAP3 (Australian Square Kilometre Array Pathfinder), SKA4 (Square
Kilometer Array) y LOFAR 5 (Low Frequency Array); estudiar la evolución de las
estructuras más masivas del Universo y catalogar una muestra de más de mil lentes
gravitacionales para estudios de cosmografı́a (Oteo, I., et al., 2016).

    Además, SHARKS se complementará con varios proyectos, creando un legado as-
tronómico de gran valor. En la Figura 3 se puede observar la distribución espacial del
proyecto SHARKS (en rojo) en sus tres campos, GAMA12, GAMA15 y en el campo en
el Plano Sur Galáctico (South Galactic Plane, SGP), comparada con otros surveys de
galaxias pasados, presentes y futuros. Cabe destacar el solapamiento de la región SGP
con el Dark Energy Survey (DES), Euclid y LSST.
  3
    https://www.atnf.csiro.au/projects/askap/index.html
  4
    https://www.skatelescope.org/
  5
    https://www.astron.nl/telescopes/lofar/

                                           5
Figure 3: Distribución espacial del proyecto SHARKS junto con otros surveys
astronómicos. Imagen obtenida de http://research.iac.es/proyecto/sharks/
pages/en/home.php

1.3    Redshift
Entender cómo evoluciona la formación de las galaxias con el tiempo cósmico es uno de
los temas principales de la astrofı́sica, por ello, la medida de la distancia a las galaxias
es uno de los parámetros crı́ticos en astrofı́sica.

   El redshift es un fenómeno que se observa en todas las galaxias (a excepción de algu-
nas cercanas a la Vı́a Láctea) y es debido a la expansión del Universo. Este fenómeno se
percibe como un aumento de la longitud de onda de la luz emitida por una fuente, una
vez detectada en la Tierra, respecto a la misma frecuencia medida en un laboratorio.
Esta se define por medio de la siguiente expresión:
                                 λobservada − λemitida
                                z=                     ,                          (1)
                                        λemitida
 siendo λ la longitud de onda. Lo interesante del redshift, es que está directamente

relacionado con la distancia según la ley de Hubble. Cuando en 1929 Edwin Hubble
(tras trabajos previos de Georges Lemaitre) descubrió que las galaxias se alejan unas de
otras (Hubble, 1929), lo explicó suponiendo que el espacio que las separa va creciendo
de manera progresiva, es decir, que el universo está en expansión. Hubble y compañı́a
establecieron, de forma empı́rica, que el valor del redshift (z) es proporcional a la
distancia que nos separa de una galaxia, d, tal que:

                                       z = H0 d/c ,                                     (2)
donde la constante de proporcionalidad H0 se denomina la constante de Hubble y donde
c corresponde a la velocidad de la luz.

    Por tanto, para obtener la distancia a una galaxia, en primera aproximación, y según
la ley de Hubble, basta con medir el redshift. Dependiendo de la precisión que necesite-

                                             6
mos, del número de galaxias a las que queramos medir el redshift y de la profundidad
a la que queramos ir, existen 2 técnicas para medir el redshift.

    Por un lado, se encuentra el redshift espectroscópico (spec-z). Esta es la técnica
que nos da el valor del redshift más preciso. En esta técnica, medimos el espectro a
la galaxia, lo que nos permite identificar las lı́neas de emisión y absorción y por tanto,
podemos comparar la longitud de onda a la que se encuentran con la que deberı́a tener
si proviniese de una muestra en reposo. Esta técnica nos permite obtener un desplaza-
miento al rojo exacto, sin embargo, se necesita un largo tiempo de exposición (alrededor
de 3 horas para un telescopio de 4 m para un objeto de magnitud aparente cercana a
22, Juncosa Serrano 2010). Hoy en dı́a, existen técnicas para obtener varios espectros
a la vez (hasta miles de ellos), no obstante, es difı́cil llegar al orden del millón con estas
técnicas.

    Por otro lado se encuentra la técnica que se usa en SHARKS, que es medir el redshift
a partir de la fotometrı́a, el redshift fotométrico (photo-z). En este método, se emplean
varios filtros a diferentes frecuencias para medir la intensidad de la luz en dichos rangos
espectrales. Estas medidas se comparan luego con observaciones o modelos de galax-
ias a distintos redshifts para encontrar el mejor ajuste y ası́ estimar el redshift. Esta
técnica tiene la ventaja de que podemos cubrir amplias regiones del cielo con tiempos
de exposición mucho más cortos que con espectroscopı́a, es decir, tiene la ventaja frente
al redshift espectroscópico de ser mucho más rápido, además de poder medir redshifts
simultáneamente a cientos de miles de galaxias a mayores distancias. Como desventaja,
podemos resaltar que nuestras medidas irán acompañadas de una mayor imprecisión
que con la espectroscopı́a. La resolución del photo-z será una función tanto del rango
espectroscópico cubierto (cobertura espectral), como del número de filtros usados.

    La cobertura espectral es un factor limitante para la aplicación de esta técnica. Si
no disponemos de un número suficiente de medidas que cubran regiones suficientemente
amplias del espectro electromagnético, no podremos medir con fiabilidad el redshift. Es
por esto que a veces es necesario combinar datos provenientes de diferentes instrumen-
tos que complementen los rangos espectrales que no han sido observados por nuestro
experimento. Como es el caso de este trabajo. Teniendo en cuenta que SHARKS solo
observa en la banda Ks , la cobertura espectral es poca, por lo tanto, necesariamente
necesitaremos complementar estas medidas con observaciones en otros rangos espec-
trales. Es por ello que haremos uso de datos del Dark Energy Survey, que cubren la
misma región del cielo en el rango óptico, ampliando ası́ la cobertura espectral.

    Para obtener el photo-z, durante décadas, se han desarrollado diferentes técnicas,
que pueden resumirse en dos: estimación del photo-z por métodos de ajuste a plantil-
las (template fitting, Bolzonella et al. 2000) y estimación del photo-z por métodos de
entrenamiento (Collister & Lahav, 2004). El primero consiste en comparar los flujos
(o magnitudes) observados, con una librerı́a de espectros de galaxias (SED, del inglés
spectral energy distribution) que se hacen evolucionar a diferentes redshifts. Dado el
sistema de filtros donde se observa, y una vez obtenidos los flujos para las SED en el
rango de redshift determinado, se busca el mejor ajuste a los datos, estimando ası́ el
photo-z. En el segundo caso, se construyen redes neuronales, en las que, una vez en-

                                              7
trenadas con una muestra conocida de galaxias con el redshift exacto medido (spec-z),
se usa para estimar el photo-z a la muestra requerida.

   En nuestro caso, dado que los datos de SHARKS analizados solo ocupan 7 grados
cuadrados, tendremos muy pocos redshifts espectroscópicos para construir redes neu-
ronales fiables, y por tanto, estimaremos el photo-z con técnicas de template fitting.

1.4    Objetivos
En nuestro estudio se pueden destacar dos objetivos claros. En primer lugar, hablamos
de la validación de la calibración de los datos de SHARKS-DR1 respecto al catálogo
de referencia 2MASS, donde se realizará una regrasión lineal entre las magnitudes de
ambos catálogos.

    En segundo lugar, y como objetivo principlal, medimos photo-z para galaxias del
SHARKS-DR1, en particular, para una submuestra de galaxias asociadas con datos
ópticos del Dark Energy Survey. Calculamos dos estimaciones, una usando los datos
ópticos de DES nada más (la muestra DES Only) y usando datos ópticos más los
datos de SHARKS (la muestra DES-SHARKS). Como primera conclusión del estudio,
comparamos los resultados con una muestra espectroscópica y a continuación, hacemos
un estudio de aquellas galaxias para las que la medida de photo-z es muy diferente
entre los datos de DES Only y DES-SHARKS. En ambas validaciones, encontramos
que los resultados de DES-SHARKS mejoran el desepeño de la muestra DES Only,
confirmando que el uso de la banda Ks mejora la definición del corrimiento al rojo de
las galaxias en SHARKS.

                                          8
2         Datos
In this chapter, we discuss the data used to obtain the photometric redshifts, where we
describe the main features of the SHARKS-DR1 catalog, of the 2MASS sample, that we
use as a reference for the validation of the SHARKS photometry, and also, of the Dark
Energy Survey catalog, that we use to increase the spectral coverage of our galaxies, in
order to obtain redshifts.

2.1        SHARKS-DR1
SHARKS es un survey fotométrico que observa en la banda Ks , que observará alrededor
de 300 grados cuadrados en varias regiones definidas por H-ATLAS, en particular, en
el South Galactic Plane (SGP), GAMA-12h (G12) and GAMA-15h (G15).
    En este trabajo hemos utilizado los datos de la primera publicación (release) de
datos del proyecto, denominada SHARKS-DR1. Estos datos consisten en imágenes
calibradas y catálogos de fuentes para alrededor de 20 grados cuadrados, divididos en
10 imágenes de ≈ 2 grados cuadrados cada uno. En este primer release, la reducción
de imágenes se ha afinado para detectar el mayor número posible de fuentes. En este
catálogo, el lı́mite de magnitud a 5σ, en el sistema AB, para las fuentes puntuales es
de 22.7 mag (con una desviación estándar de alrededor de 0.25 mag). SHARKS-DR1
ha sido producido en colaboración entre el Instituto de Astrofı́sica de Canarias (IAC)
y el Wide-Field Astronomy Unit (WFAU) del Royal Observatory of Edinburgh. Estos
datos serán publicados en las próximas semanas y se encontrarán a disposición de todo
el mundo a través de la página del IAC y también en los archivos de la ESO.
    SHARKS-DR1 está calibrado astrométrica y fotométricamente con respecto al catálogo
astronómico 2MASS (Skrutskie et al., 2006). Las imágenes del DR1 están distribuidas
entre los tres campos de SHARKS, donde seis de ellas están aisladas y cuatro de ellas
cubren una región contigua en el campo este del SGP (SGP-E, Figura 4). Es esta región
en la que centramos el análisis, debido al solapamiento con el proyecto DES6 (DES Col-
laboration, 2021). En esta región hay un solapamiento de ≈ 1 grado cuadrado (13%
del área) entre las imágenes, por tanto, el área final de estudio cubre un área de 7.25
grados cuadrados. En esta área, encontramos 623543 objetos.
    Aunque centramos nuestro análisis únicamente en esta área, cabe recordar que el
área total del DR1 es de 19.75 grados cuadrados. Cabe resaltar aquı́ que las regiones
GAMA12 y GAMA15, tiene solapamiento con datos ópticos del HSC-SSP 7 (Hyper
Suprime-Cam Subaru Strategic Program, Aihara et al. 2019), lo que nos amplı́a la
cobertura espectral para estos campos. Desafortunadamente, en el SGP, más allá de la
región con solapamiento con DES, no hay datos en el óptico suficientemente profundos.

    6
        https://des.ncsa.illinois.edu/releases/dr2
    7
        https://hsc-release.mtk.nao.ac.jp/doc/

                                              9
Figure 4: En azul se encuentran las 10 regiones que conforman el SHARKS-DR1. En la
región SGP-E se encuentran los 7.23 grados cuadrados donde se centra nuestro estudio,
en verde el campo cubierto por DES (Carnero, comunicación privada).

2.2    DES
Para ampliar la cobertura espectral de las galaxias observadas en SHARKS, usamos
datos del catálogo astronómico Dark Energy Survey (The Dark Energy Survey Collab-
oration, 2005), concretamente, datos del catálogo público Data Release 2.
    DES es una colaboración internacional dedicada a cartografiar cientos de millones
de galaxias y miles de supernovas que revelan la naturaleza de la misteriosa energı́a os-
cura, causante de la expansión acelerada del Universo. DES comenzó sus observaciones
en 2013. A lo largo de seis años (2013-2019), tomó imágenes de 5000 grados cuadrados
(aproximadamente 1/8 del área celeste) en cinco filtros ópticos (g,r,i,z,Y) para obtener
información detallada a ≈ 500 millones de galaxias. Para ello observó durante 525
noches utilizando la Cámara de Energı́a Oscura (DECam), montada en el telescopio
Blanco de 4 metros del Observatorio Interamericano de Cerro Tololo, situado en lo alto
de los Andes chilenos. Esto ha permitido registrar información de galaxias tan débiles,
que para muchas su luz es ≈ 1 millón de veces más tenue que la estrella más tenue que
puede verse a simple vista.

   En este proyecto hemos utilizado la segunda liberación de datos, el DES-DR2 que in-
cluyen los datos de los seis años observados. La mediana de la profundidad del catálogo,
para fuentes puntuales es g = 24.7, r = 24.4, i = 23.8, z = 23.1, e Y = 21.7 mag (10σ,
AB). Además, presenta una calidad de imagen excelente, con un FWHM mediano de g
= 1.11, r = 0.95, i = 0.88, z = 0.83, Y = 0.90 arcosegundos. DES DR2 incluye ≈ 691
millones de objetos astronómicos distintos. Tras una selección de calidad básica, las
muestras de galaxias y estrellas 543 millones y 145 millones de objetos, respectivamente.

    Esta segunda versión pública ha duplicado aproximadamente, frente a la primera
versión pública (DES DR1, Abbott et al. 2018), el tiempo de exposición en la mayorı́a
del área, aumentando la profundidad entre 0.45 y 0.7mag, y elevando el número total
de objetos catalogados de ≈ 400M a ≈ 700M. También se han realizado mejoras en la

                                           10
uniformidad de la calibración fotométrica (< 3mmag), la precisión astrométrica interna
(≈ 27mas) y la precisión de la separación estrella-galaxia. Por tanto, con la profundidad
y precisión fotométrica alcanzadas, DES DR2 constituye el mayor conjunto de datos
fotométricos hasta la fecha. Estos datos son accesibles a través de diferentes interfaces
abiertos al público en https://des.ncsa.illinois.edu/releases/dr2/dr2-access.

2.3    2MASS
El proyecto Two Micron All Sky Survey (2MASS, Skrutskie et al. 2006) es un survey
en el infrarrojo cercano terminado hace casi una década. Pese a ello, sigue siendo el
survey de referencia en el NIR al cubrir todo el cielo con gran homogeneidad. En su
origen, fue diseñado para proporcionar respuestas a preguntas sobre la estructura a
gran escala de la Vı́a Láctea y el Universo Local. Para lograr estos objetivos, 2MASS
escaneó uniformemente todo el cielo en tres bandas del infrarrojo cercano (J,H, Ks ), en
las tres ventanas mostradas en Figura 1, para detectar y caracterizar fuentes puntuales
más brillantes que aproximadamente 1 mJy en cada banda, con una relación señal-ruido
(SNR) superior a 10.

    2MASS utilizó dos telescopios de 1.3 m altamente automatizados, uno en el Monte
Hopkins, Arizona, Estados Unidos, y otro en el CTIO, Chile. Cada telescopio estaba
equipado con un conjunto de detectores capaces de observar el cielo simultáneamente
en J (1.25 micras), H (1.65 micras) y Ks (2.17 micras). Uno de los grandes hitos de este
proyecto fue ser el construir el primer censo completo de galaxias brillantes (Ks < 13.5
mag), incluyendo galaxias en la “Zona de Evitación” alrededor del disco de la Vı́a
Láctea (zona de 60 grados de alto, donde el polvo dentro de la Vı́a Láctea hace que los
estudios de galaxias ópticas sean sesgados).

    Igualmente importante es el censo de estrellas de 2MASS, con el cual se han calibrado
los datos de SHARKS. Este censo se encuentra en el llamado Point Source Catalog
(PSC, Cutri et al. 2003), que incluye alrededor de 500 millones de estrellas y galaxias
sobre el 99.998% del área del cielo. La calibración de este catálogo astronómico se ha
realizado respecto a la estrella Vega. En este sistema, alcanzamos de profundidad de
J < 15.8, H < 15.1 y Ks < 14.3 mag para 99% del área del cielo. Debido a que 2MASS
tiene un alto grado de consistencia interna, es posible calibrar el sistema de coordenadas
y la fotometrı́a de las imágenes de VISTA. Las magnitudes de 2MASS proporcionan
puntos cero fotométricos de los filtros de VIRCAM con una precisión del 1-2%.

                                            11
3         Validación fotométrica de los datos de SHARKS-
          DR1
In this section we present the validation that we have made of the calibration of the
SHARKS-DR1 data, for this, understanding that the DES catalog is already calibrated,
we will use 2MASS to check the calibration in SHARKS only, by performing a linear
regression between the magnitudes of SHARKS and 2MASS. As we will see in this
section, it is shown that the SHARKS-DR1 calibration is correct and that the ordinate
at the origin of the fit shows the difference between the Vega and AB calibrations.

    Entrando en uno de los objetivos principales de este trabajo, un aspecto fundamen-
tal a la hora de obtener redshift fotométricos es la correcta calibración fotométrica de los
datos, más aún cuando estos provienen de diferentes experimentos, como es el caso de
este trabajo donde usaremos datos de DES y SHARKS. En el caso de DES, confiamos
en lo explicado en DES Collaboration (2021), donde la fotometrı́a está correctamente
calibrada, con una precisión de ≈ 10 mmag, en el sistema AB, y por lo tanto, usaremos
los datos directamente.

    En el caso de SHARKS-DR1, al ser un catálogo que todavı́a no está publicado,
procederemos a comprobar si la calibración respecto a 2MASS es correcta. Según la
información disponible, SHARKS está calibrado respecto 2MASS en el sistema AB,
mientras que 2MASS es en el sistema Vega. Por tanto, para validar la calibración de
SHARKS-DR1, asociaremos estos datos con el PSC de 2MASS, y con diferentes con-
figuraciones de ajuste, probaremos que la calibración es la correcta.

    En general, para estimar la precisión fotométrica de un survey respecto a otro,
comparamos las magnitudes de las fuentes puntuales (estrellas) entre ambos. En un
diagrama magnitud-magnitud, si la calibración es correcta, se deberá obtener una pen-
diente de valor unidad (es decir, que la calibración es independiente de la magnitud) y
una ordenada en el origen que representa la diferencia entre la calibración entre ambos
(la diferencia en el punto cero). En nuestro caso, ya que SHARKS está calibrado re-
specto a 2MASS en el sistema AB, y 2MASS está calibrado respecto a Vega, esperamos
obtener una diferencia constante (asumiendo linealidad), cuyo valor viene recogido en
la Figura 5, donde mostramos también las diferencias entre AB y Vega para todos
los filtros de VIRCAM, además de otras propiedades de los filtros. En este caso, la
diferencias entre sistemas AB y Vega en Ks deberı́a ser ∆Ks = 1.827.

3.1        Selección de la muestra
Para medir la calibración de SHARKS-DR1 con 2MASS empezamos bajando los datos
del 2MASS-PSC del servidor Vizier 8 , en la región de SHARKS-DR1 definida entre
8.27 < RA < 12.15 y −32.04 < Dec < −29.7. El número de objetos que obtenemos
de 2MASS son 10957. De esta muestra, aplicamos cortes de calidad, de tal forma
que seleccionamos aquellos objetos con fotometrı́a correcta en las tres bandas J,H,Ks
    8
        http://vizier.unistra.fr/

                                              12
Figure 5: Propiedades actuales de los filtros de la cámara VIRCAM de VISTA (v1.5).
Para nuestro análisis, destacamos los valores para la banda Ks , en particular, la difer-
encia de calibración entre el sistema Vega y AB. Imagen de http://casu.ast.cam.
ac.uk/surveys-projects/vista/technical/filter-set

con señal-ruido > 10 (Qflg = AAA), que no tienen pı́xeles saturados (Rflg = 222) y
que no está afectada por ningún artefacto cercano (Cflg = 000). Finalmente, en base
a la profundidad estimada (ver Sección 2.3), seleccionamos objetos únicamente con
Ks ≤ 14.3. El catálogo resultante de 2MASS contiene 3167 objetos puntuales.
    A continuación, realizamos la asociación con los datos de SHARKS-DR1 en un radio
de 1”, seleccionando estrellas según el parámetro CLASSSTAT (CLASSSTAT>= 0.95)
y limpiando también aquellos objetos afectados por problemas de reduccion (ERRBITS
≤ 3). Con esto obtenemos una muestra preliminar de 2481 objetos. Finalmente, en
la Figura 6 comparamos las magnitudes de SHARKS-DR1 (apertura de 2”) con las de
2MASS (magnitud “profile fitting”, la magnitud por defecto). Podemos apreciar como
en el rango más brillante, perdemos la linealidad entre las magnitudes, por tanto, para
que esto no afecte a nuestro estudio y siendo conservadores, seleccionamos aquellos
objetos con una magnitud 2MASS Ks ≥ 12.8, obteniendo 1633 objetos que usaremos
en las siguientes secciones.

3.2      Pruebas de calibración
A continuación, para comprobar la correcta calibración de SHARKS-DR1, realizamos
una serie de ajustes a los datos anteriores, teniendo en cuenta diferentes configura-
ciones. Estas son, y en orden de precisión: ajuste asumiendo que no tenemos errores
en las magnitudes, y sin quitar outliers (entendiendo los outliers como observaciones
que numéricamente están alejadas de la regresión lineal, en nuestro caso, definido como
aquellas que se alejan más de un 15% de la regresión lineal). Segundo, ajustaremos las
magnitudes asumiendo errores, pero todavı́a sin quitar outliers. Tercero, asumiendo er-
rores, y además, quitando los errores. Y por último, veremos la calibración de SHARKS,
ajustando por un polinomio de orden cero, es decir, asumiendo una linealidad perfecta
entre SHARKS y 2MASS. Los resultados pueden verse en las Figuras 7, 8, 9, 10, 11
respectivamente. En la Tabla 1 resumimos los resultados.
   Tras analizar los resultados, donde obtenemos siempre un valor superior a lo esper-
ado, revisamos la calibración realizada en los datos de SHARKS-DR1, y descubrimos,
a posteriori, que la calibración aplicada corresponde a la versión v1.3 (según la página
de VISTA9 ) y no la reflejada en la Figura 5 (versión v1.5), con un valor usado de
∆Ks = 1.839. Sabiendo esto, ahora encontramos un acuerdo muy bueno, especialmente
  9
      http://casu.ast.cam.ac.uk/surveys-projects/vista/technical/filter-set

                                            13
Figure 6: Comparación de las magnitudes entre SHARKS-DR1 y 2MASS, para la
muestra de calibración. En base a esta comparación realizamos un corte extra donde
eliminamos aquellos objetos con Ks < 12.8 en 2MASS (en gris), seleccionando la región
de mayor linealidad entre SHARKS-DR1 y 2MASS. La muestra final contiene 1633
objetos.

Figure 7: Regresión lineal entre la magnitud de 2MASS (Kmag) y SHARKS-DR1
(APERMAG3), sin considerar errores, y admitiendo outliers en el ajuste.

                                         14
Figure 8: Regresión lineal entre la magnitud de 2MASS (Kmag) y SHARKS-DR1
(APERMAG3), considerando errores en el ajuste, pero admitiendo outliers en el ajuste.

Figure 9: Regresión lineal entre la magnitud de 2MASS (Kmag) y SHARKS-DR1
(APERMAG3), ignorando los errores en el ajuste, pero esta vez eliminando los out-
liers en el ajuste.

                                         15
Figure 10: Regresión lineal entre la magnitud de 2MASS (Kmag) y SHARKS-DR1
(APERMAG3), considerando errores en el ajuste, pero admitiendo outliers en el ajuste.

Figure 11: Regresión lineal entre la magnitud de 2MASS (Kmag) y SHARKS-DR1
(APERMAG3), considerando errores en el ajuste, admitiendo outliers en el ajuste y
fijando la pendiente a 1 (asumiendo linealidad).

                                         16
Ordenada en el origen Pendiente
  Sin errores y con outliers                             1.839               1.001
  Con errores y con outliers                             1.898               0.996
  Sin errores y sin outliers                             1.861               0.999
  Con errores y sin outliers                             1.890               0.997
  Con errores y sin outliers, fijando Pendiente          1.846            1. (fijo)

Table 1: Resultados obtenidos de ajustar la magnitud de SHARKS (en sistema AB) y la
magnitud de 2MASS (en sistema Vega). La diferencia es siempre ligeramente superior
al valor dado por VISTA (v1.5), pero muy parecida a la calibración de la versión v1.3,
según la página de VISTA.

en el caso de no usar errores, y sin eliminar los outliers. Entendemos que este ha sido
el método usado en la calibración fotométrica realizada por el Wide-Field Astronomy
Unit (los encargados de la calibración del SHARKS-DR1).

    Por tanto, asumimos a continuación, que los resultados están correctamente cal-
ibrados, teniendo en cuenta una diferencia sistemática de al menos 0.012 entre las
magnitudes de 2MASS y SHARKS-DR1.

                                          17
4     Photo-z de los datos de SHARKS-DR1
After the validation of SHARKS-DR1, we proceed to obtain the photometric redshifts
of galaxies identified in the DES-SHARKS association. For this purpose, we introduce
Lephare, the software that we use to calculate photometric redshifts, using generalized
galaxy templates and the filters in which the data have been observed. Also in this
chapter we describe the spectroscopic sample we use to validate the results at low
redsfhit. Finally, we discuss the results obtained. We have validated against spectra
and also comparing the results obtained with DES-SHARKS and DES alone, where we
see the positive impact of the Ks band in our results.

    Una vez comprobada la correcta fotometrı́a de SHARKS-DR1, pasamos a estimar
los primeros photo-z para galaxias detectadas en SHARKS.
    Antes de pasar a describir el código usado y los resultados obtenidos, es oportuno
decir lo siguiente: al ser estos los primeros photo-z obtenidos en SHARKS, hemos
sido conservadores en la muestra seleccionada, además de simplificar lo más posible la
librerı́a de plantillas usadas, para facilitar la interpretación de los resultados. Cabe
resaltar, además, que no siempre un conjunto mayor de plantillas se traduce en una
mejor determinación de photo-z, como han demostrado en Abdalla et al. (2011).
    Una vez terminado este análisis, lo que serı́a el primer paso, la colaboración pro-
cederá a un análisis de mayor extensión. Sin embargo, ese objetivo queda fuera de
nuestro trabajo. Es por ello que en nuestro análisis seleccionaremos únicamente aquel-
los objetos con una detección clara, es decir, de hasta 5σ en las bandas i, z en DES y
5σ en la banda Ks en SHARKS. De esta forma, no calcularemos photo-z a todos los
datos de DR1, si no a una submuestra de la misma. También puntualizar que debido al
pequeño tamaño de la muestra, tenemos muy pocos objetos con spec-z, lo que dificulta
aún más la validación del photo-z obtenido. Igualmente, disponemos de una muestra
de 1000 objetos con los valores de sus redshift espectroscópicos (que describiremos más
adelante) hasta z < 0.7 que usaremos como muestra de validación.
    Finalmente, cabe recordar otro objetivo del trabajo, que es medir el efecto de usar
la banda Ks en la estimación de photo-z, en combinación con datos ópticos. ¿Mejorará
la medida de photo-z si incluimos la banda Ks respecto al Photo-z usando únicamente
datos ópticos?. En este sentido, estudios previos (Banerji et al., 2015) han estimado
que la contribución de bandas infrarrojas junto con datos ópticos reducen la dispersión
del photo-z para z < 0.5 y z > 1.

4.1    Lephare
Lephare (Arnouts & Ilbert, 2011) es un programa para calcular desplazamiento al rojo
fotométricos, utilizando la técnica de template fitting. La ejecución se compone de tres
partes:
    • Una fase preliminar, donde seleccionamos la librerı́a de plantillas (SED) y defini-
      mos el conjunto de filtros.
    • Con estos dos ingredientes, calculamos las magnitudes de las plantillas, en función
      del redshift. De esta forma, creamos un hiperespacio definido por los diferentes

                                            18
filtros observados, y el redshift. En el caso de Lephare, interpolamos en la di-
        mensión de redshift. Este paso define también información básica relativa a los
        filtros usados (longitudes de onda, diferencia de corrección AB a Vega...) y a
        las SEDs. Cabe también destacar que estos pasos pueden utilizarse para generar
        catálogos multicolores simulados.
      • Finalmente, realizamos una minimización por el método de ajuste χ2 entre las
        magnitudes observadas, y los valores definidos en el hiperespacio de magnitudes
        teórico para obtener el redshift y la plantilla que mejor se ajusta a cada galaxia.
        Adicionalmente (pero no explorado en este trabajo), podemos obtener la dis-
        tribución de probabilidad completa para cada objeto.

4.1.1      Definición de librerı́a SED
El procedimiento de ajuste para obtener photo-z se basa en el ajuste de la forma de
los espectros en función del redshift. Las magnitudes observadas se comparan con las
obtenidas a partir de un conjunto de espectros de referencia, utilizando el mismo sistema
fotométrico. Estas SED de plantilla pueden ser observadas o sintéticas. Por tanto, en
primer lugar, debemos seleccionar la lista de SEDs de referencia que vamos a emplear
en nuestro estudio, es decir, las plantillas de modelos de galaxias que nos ayudarán a
aproximar cada galaxia que observemos a una distribución espectral de energı́a. La
medida de redshift se beneficiará si las galaxias seleccionadas tienen caracterı́sticas
significativas en el espectro, como la ruptura de 4000 Å o la ruptura de Lyman a 912
Å.
    Como dijimos anteriormente, según los resultados de Abdalla et al. (2011), un
número grande de SED no siempre significa mejores resultados. Por tanto, para sim-
plificar la interpretación en este primer análisis, elegimos un conjunto con 8 SEDs,
sin incluir extinción interna. Usamos la librerı́a de Lephare llamada bpz, provenientes
de Coe et al. (2006). Las SEDs que hemos empleado son:
      • 1 plantilla de galaxia elı́ptica (Ell)
      • 2 plantillas de galaxias espirales (Sbc y Sbd)
      • 1 plantilla de galaxia irregular (Im)
      • 2 plantillas de galaxias espirales barradas (SB2 y SB3)
      • 2 plantillas de galaxias azules (ssp 5M y ssp 25M)
      En la Figura 12 mostramos las plantillas (a z = 0) usadas en nuestro trabajo.

4.1.2      Definición de filtros
En segunda instancia, también debemos aportar al programa los filtros en los que se han
observado los datos, para convolucionar con las plantillas a diferentes redshifts. Cada
uno de estos determina un rango especı́fico del espectro electromagnético. Los filtros
que hemos utilizado son los filtros g,r,i,z,Y de DECam10 y el filtro Ks de VIRCAM11 .
Estos se pueden ver en la Figura 13.
 10
      http://www.ctio.noao.edu/noao/content/DECam-filter-information
 11
      http://svo2.cab.inta-csic.es/theory/fps

                                                 19
Figure 12: Distribución espectral de energı́a para las 8 plantillas usadas en este trabajo,
provenientes de Coe et al. (2006).

Figure 13: Filtros empleados en las observaciones de nuestros datos, resultantes de la
asociación entre datos de DES y SHARKS.

                                            20
4.1.3   Estimación del photo-z
Una vez que tenemos nuestros SEDs y nuestros filtros definidos, procederemos a la
obtención de las magnitudes aparentes en un rango de redshifts definidos. Estas mag-
nitudes aparentes se obtienen al hacer la convolución del SED con el filtro correspondi-
ente. El resultado de esta será un valor de la magnitud asociado al filtro que estamos
empleando (g, r, i, z, Y, Ks ).
     Después de ese proceso, se comparan las magnitudes aparentes resultantes, con las
magnitudes de las galaxias en el catálogo astronómico de entrada, y se encuentra el
mejor ajuste, dándonos, como resultado el SED y el redshift que mejor ajusta a cada
lı́nea del catálogo astronómico.
     A continuación, para medir el photo-z utilizamos un procedimiento de ajuste de χ2 ,
comparando el flujo observado (Fobs ) y sus correspondientes incertidumbres (σ) con el
flujo de cada una de las plantillas (Ftemp ) definido como:
                                     X Fobs,i − sFtemp,i
                              χ2 =     [                 ]2                          (3)
                                     i
                                               σi

donde i hace referencia a la banda utilizada para el análisis y s es el factor de escala
que se elige para minimizar los ajustes de χ2 :
                                                     2
                              X Fobs,j Ftemp,j X Ftemp,j
                           s=   [             ]/   [     ]                           (4)
                              j
                                      σj2        j
                                                     σj2

donde j se refiere a la banda utilizada para el escalado, pudiendo ser esta distinta de
i. Lephare finalmente realiza este ajuste a través de un procedimiento estándar de
minimización de χ2 .

4.2     Preparación de los datos
Como se dijo anteriormente, hemos querido ser conservadores en la muestra seleccionada
al ser esta la primera vez que se calculan photo-z con datos de SHARKS.
    En ese caso, partimos del catálogo SHARKS en el campo SGP-E, que contiene
623543 objetos, que incluirán galaxias, estrellas, quásares y objetos espurios. A con-
tinuación, realizamos una asociación espacial con los datos del DES-DR2 (que contiene
en el mismo área, alrededor de 940000 objetos), usando un radio de 1”. Tras esta
asociación, obtenemos una muestra de 514946 objetos.
    A continuación, seleccionamos objetos con una detección clara en las bandas i,z de
DES y la banda Ks de SHARKS. En las bandas g,r,Y de DES, no ponemos ningún
requerimiento, permitiendo ası́ que la muestra contenga objetos extremadamente rojos
(tampoco en la banda Y, ya que esta es en media, más brillante que las otras). Por
tanto, el siguiente corte que hacemos es seleccionar objetos con SN R > 5 en i,z,Ks .
    Este corte se traduce, usando magnitudes, en seleccionar objetos que pasan el corte
1.086/MAGERR > 5 para i,z,Ks , donde MAGERR es el error que nos proporciona
el propio catálogo SHARKS de la magnitud. Tras aplicar este corte, obtenemos fi-
nalmente una muestra de 247522 objetos, aproximadamente, la mitad de la muestra
DES-SHARKS inicial.

                                            21
4.2.1      Fotometrı́a de la muestra
A la hora de calcular el photo-z con Lephare, usaremos las magnitudes denominadas
AUTO obtenidas por el código de SExtractor12 (Bertin & Arnouts, 1996). De este
código podemos comentar brevemente que aplica técnicas automatizadas para detec-
tar, medir y clasificar fuentes de imágenes astronómicas y es especialmente adecuado
para el análisis de grandes estudios. En cuanto a las magnitudes AUTO, se definen
como el parámetro que proporciona una estimación de la “magnitud total” mediante
la integración de los valores de los pı́xeles dentro de una apertura escalada adaptativa-
mente.
     Antes de proceder, debemos tener en cuenta el efecto del polvo interestelar y corregir
las magnitudes convenientemente. El polvo interestelar se encuentra presente en todas
las observaciones. El estudio y el conocimiento de la distribución del polvo cósmico por
el espacio es de suma importancia, dado que contamina u oscurece cualquier medida
astrofı́sica, especialmente, las medidas extragalácticas. Varios estudios han estudiado
la distribución de polvo, siendo el más referenciado la distribución de Schlegel et al.
(1998). Usamos la estimación del exceso de color, E(B-V), de esta referencia, para
corregir las magnitudes del catálogo DES-SHARKS.
     La influencia en nuestras medidas que produce polvo cósmico en el rango óptico
e infrarrojo cercano, disminuye a medida que aumenta la longitud de onda. Hoy en
dı́a, encontramos tabulados las correcciones que se tienen que realizar a las magnitudes
observadas. Las magnitudes se corrigen por una cantidad ∆mb = E(B − V )Rb , donde
los valores de Rb son, en el caso de las bandas g,r,i,z,Y de DES (DES Collaboration,
2021): Rg = 3.186, Rr = 2.140, Ri = 1.569, Rz = 1.196, RY = 1.048 y para la banda Ks ,
RKs = 0.005 (González-Fernández et al., 2017). Finalmente, las magnitudes corregidas
son las que usaremos en lephare, con: M AGb,CORREGIDA = M AGb − ∆mb .

4.2.2      Separación estrella-galaxia
Finalmente, antes de proceder a calcular el photo-z, hacemos la selección de galaxias.
Este paso es uno de los temas más difı́ciles de abordar en fotometrı́a, y siempre habrá
un nivel de confusión asociado a esta separación. Por suerte, se ha demostrado (ver,
por ejemplo Sevilla-Noarbe et al. 2021) que el uso de una combinación de datos ópticos
e infrarrojos juntos, es un separador estrella-galaxia simple y eficaz, capaz de producir
resultados comparables a la clasificación por métodos de entrenamiento. Con esto, rep-
resentamos los datos en el espacio de colores (i-z), (z-Ks ) en la Figura 14, donde vemos
claramente la separación de galaxias de la secuencia principal de estrellas. Finalmente,
decidimos seleccionar galaxias como aquellas fuentes donde:

                                 (z − Ks ) > −0.1 + (i − z)                               (5)
ası́ obtenemos la muestra final de 200645 galaxias para la cuál calcularemos los photo-z.

4.3       Validación de los datos
La validación de los resultados de photo-z es un paso crı́tico para evaluar si la estimación
de photo-z ha sido buena o no. Como dijimos anteriormente, debido al pequeño área de
  12
       https://sextractor.readthedocs.io/en/latest/

                                             22
Figure 14: Representación de la separación entre el cúmulo de galaxias (parte superior
respecto a la recta) y la secuencia principal de estrellas (parte inferior). Para separar
las galaxias de las estrellas representamos el polinomio y = −0.1 + x.

la muestra de análisis, tenemos pocos objetos con redshift espectroscópicos y por tanto,
nuestra validación puede no ser representativa. Igualmente, como prueba de concepto,
procederemos a validar nuestros resultados con la muestra de validación existente. A
continuación, describimos las métricas usadas, y a seguir, las muestras espectroscópicas
disponibles para calibrar.

4.3.1   Métricas usadas
Una vez obtenidos los valores del redshift fotométrico vamos a estimar cuánto de bueno
son estos datos. Para ello, compararemos nuestros valores obtenidos con los datos
de redshift espectroscópico (spec-z) que ya han sido determinados en otros surveys
astronómicos, como veremos en la sección 4.3.2. La cantidad de valores de redshift
espectroscópicos que tenemos no son muy grandes, por tanto los valores estadı́sticos
que obtengamos pueden no ser representativos, y además, limitados en redshift. En
nuestro caso, tenemos galaxias hasta spec-z ≈ 0.7. Las métricas que vamos a emplear
van a ser, en primer lugar el bias de redshift, zbias :
                                            N
                                         1 X i         i
                               zbias   =      (z    − zspec )                          (6)
                                         N i=1 phot
y la desviación media absoluta normalizada, σN M AD :

                                                   zspec − zphot
                            σN M AD = median                                           (7)
                                                     1 + zspec

                                              23
Figure 15: Distribución de redshifts espectroscópicos disponibles para validación, prove-
nientes de 2dF, 6dF y del 2dFLens.

   Por otro lado, también mediremos σ68 , definido como el valor tal que 68% de las
galaxias tienen |zphot −zspec | < σ68 , y la fracción de “outliers”, definido como la cantidad
de galaxias con ∆z > 0.15 entre el valor de zspec y zphot :

                                fout = N [ |∆z| > 0.15] /Ngals                             (8)

4.3.2   Muestra de calibración
Dado el área de SHARKS-DR1, identificamos las siguientes muestras espectroscópicas
en el área de análisis:

   • 2dFLens (Blake et al., 2016): el 2-degree Field Lensing Survey (2dFLenS) es un
     survey espectroscópico realizado desde el Anglo-Australian Telescope. Se diseñó
     para observar un área previamente observado por el Kilo-Degree Survey para la
     búsqueda de lentes gravitacionales. En total, se observaron 70079 redshifts en el
     rango z < 0.9, en un área de 731 grados cuadrados. Esta muestra es la que nos
     aporta la mayor cantidad de redshifts espectroscópicos a z > 0.2.

   • 2dF (Colless et al., 2001) y 6dF (Jones et al., 2009): Estos dos surveys, anteriores
     al 2dFLens, nos aportan unos redshifts extras a muy bajo redshift (z < 0.7).
     Fueron observados también desde el mismo telescopio del Anglo-Australian Tele-
     scope.

   En la Figura 15 mostramos la distribución de spec-z disponibles para validación de
nuestros datos.

4.4     Resultados y discusión
Corremos Lephare con los templates descritos anteriormente. En el Apéndice A de-
tallamos todos los parámetros usados a la hora de ejecutar Lephare. Internamente,
hicimos pruebas con otras combinaciones de plantillas, por ejemplo, usando las plan-
tillas de Poggianti como usan en Abdalla et al. (2011). Aunque las diferencias en los

                                              24
resultados son pequeñas, y continenen el mismo número de plantillas, las de bpz de-
scriben mejor la región del infrarrojo cercano.

   A la hora de correr lephare, ejecutamos en 2 modalidades, con los datos de DES nada
más (DES Only), y con los datos de DES + SHARKS (DES-SHARKS). Igualmente,
corremos únicamente para las galaxias, de acuerdo con la selección explicada en la
sección 4.2.

4.4.1   Validación frente a Spec-z
Una vez calculados los photo-z para las 2 configuraciones, comparamos los resultados
en bines de redshift, entre 0 y 0.7, respecto a la muestra de validación definida en la
sección 4.3.2. En la figura 17 mostramos la validación para las métricas explicadas
anteriormente, para la muestra DES Only y DES-SHARKS. Vemos como, en general,
los resultados de usar la banda Ks son mejores que sin usar la banda de SHARKS.
Especialmente, vemos una disminución del número de outliers a Spec − z ≈ 0.2, de
50% a 30%, y una mejora en la desviación media absoluta normalizada. En la figura 17
mostramos la relación entre spec-z y photo−z para las 2 muestras, que nos sirven para
entender por qué el z-bias parece mejor en la muestra DES Only; esto sucede por un
incremento en el número de “outliers” en la dirección contraria al valor medio del bias.

4.4.2   Comparación de resultados
En la figura 18 vemos la distribución global de photo-z para las dos muestras.
    La diferencia más significativa es el incremento de galaxias en DES-Only a photo−z ≈
1.5. La pregunta es, ¿Es real este incremento? En caso afirmativo, cabrı́a afirmar que
el uso de la banda Ks empeora la determinación de photo-z. Pero, en caso contrario,
podrı́amos estar viendo un efecto de corrección, en el sentido de una disminución de
“outliers”, como vimos también a bajo redshift.
    Para investigar este efecto, mostramos en la Figura 19 la comparación entre los dos
photo-z, en el eje x, mostramos el photo-z de DES Only, y en el eje y, el photo-z para
DES-SHARKS. Claramente, dado que los objetos que se encuentran cercanos a la recta
y = x, son similares, vamos a centrarnos en aquellos alejados de la linealidad, es decir los
“outliers”. En estos identificamos una población de objetos para los cuales el resultado
de uno y otro photo-z cambia mucho. En la Figura 20, resaltamos estos outliers a alto
redshift. Lo primero que vemos, es que los outliers de DES Only son bastantes más que
aquellos con DES-SHARKS, es decir, se confirma lo visto anteriormente en la Figura 18.
    Para intentar entender cuál está dando la medida más precisa, investigamos la dis-
tribución de magnitudes para las 2 muestras de outliers indicadas en la Figura 20,
mirando el histograma de las magnitudes en las bandas i,z,Ks . Esto se muestra en la
Figura 21. Claramente, sobretodo en las bandas ópticas de DES, las magnitudes son
más brillantes, siendo poco probable una distribución tan brillante a tan alto redshift.
Esto indicarı́a que la determinación del photo-z para DES Only, para estos outliers, está
mal determinada, y que seguramente sea la estimacion de DES-SHARKS la correcta.
    Para terminar, seleccionamos 3 objetos en cada una de las muestras de outliers a alto
redshift, para DES Only y para DES-SHARKS. Para estos, comparamos visualmente
las magnitudes observadas con la distribución espectral de energı́a para la plantilla

                                            25
Figure 16: Validación de photo-z para la muestra DES Only (en azul) y para la muestra
DES-SHARKS (naranja). Aunque ambas configuraciones dan resultados parecidos,
podemos apreciar una disminución en el número de “outliers” y una mejora en la
desviación media absoluta normalizada cuando usamos la banda Ks . En cuanto al bias,
aunque parezca que el caso DES Only sea mejor, esto se debe justamente al aumento
de “outliers”, como puede verse en la figura 17.

                                         26
Figure 17: Distribución de densidad, a la izquierda, para DES Only y a la derecha,
para DES-SHARKS. Las distribuciones son parecidas, aunque en el caso de DES Only
hay una mayor densidad de “outliers” a Spec − z ≈ 0.2, que resulta en una mejorı́a
artificial del bias (como puede verse en la figura 16).

Figure 18: Histograma que compara el photo-z obtenido con DES-SHARKS y solo de
DES

                                        27
Figure 19: Comparación del redshift fotométrico obtenido de SHARKS-DES (SHD, eje
y) frente al obtenido por DES Only (DOY, eje x)

                                       28
Figure 20: Seleccionamos outliers a alto redshift para DES-SHARKS (SHD, verde) y
outliers a alto redshift para DES Only (DOY, azul)

Figure 21: Distribución de magnitudes en las bandas i,z,Ks para los outliers selecciona-
dos en la Figura 20. Las magnitudes para el caso DES Only (azul) son más brillantes
de lo esperado para una muestra a tan alto redshift como indica su photo-z

                                           29
Figure 22: Comparación de las magnitudes observadas y las plantillas y redshifts del
mejor ajuste, para 3 galaxias con photo-z muy alto en la muestra DES Only. Incluir la
banda Ks rompe la degeneración, dando un valor de redshift más bajo, en lı́nea con lo
visto en la Figura 21.

Figure 23: Comparación de las magnitudes observadas y las plantillas y redshifts del
mejor ajuste, para 3 galaxias con photo-z muy alto en la muestra DES-SHARKS. Incluir
la banda Ks , igualmente, rompe la degeneración, dando un valor de redshift mayor.

y el photo-z del mejor ajuste. En el caso de DES Only, el ajuste se realiza sin la
banda Ks . Puede verse, claramente, como incluir la banda de SHARKS resuelve la
degeneración existente entre las dos plantillas a diferente redshift. En todos los casos
vistos, claramente, el uso de la banda Ks de SHARKS ayuda a una mejor determinación
del redshift.

                                           30
5    Conclusiones
To conclude, we comment about our main results, including the validation of the cal-
ibration with respect to 2MASS, and the results obtained with the comparison of the
spec-z with the phot-z through the photo-z metrics and the density distribution. Fi-
nally, we conclude that using the Ks improves the photo-z estimates.

    Habiendo expuesto ya los resultados finales del estudio, llegamos a diferentes conclu-
siones. Pero antes de abarcar las conclusiones de los redshifts fotométricos obtenidos,
cabe destacar también la validación de la calibración de los datos procedentes de
SHARKS-DR1. En este aspecto, el valor estimado de la diferencia entre el sistema
Vega y el sistema AB, con respecto al valor dado por VISTA (v1.5), es siempre ligera-
mente superior. Tras un estudio de los datos descubrimos, a posteriori, que los datos
del SHARKS-DR1 fueron calibrados, posiblemente por error, respecto a la version v1.3,
con una diferencia entre Vega y AB, ahora sı́, compatible con lo obtenido en este tra-
bajo. Además, todo indica a que la calibración se realizó sin incluir errores.

    En cuanto a los resultados del redshift fotométrico calculados en la Sección 4 de
este trabajo, usamos una pequeña muestra espectroscópica para validación. En esta, al
realizar la comparación entre el redshift fotométrico frente al espectroscópico mediante
una distribución de densidad (Figura 17) somos capaces ver la mejorı́a de nuestros re-
sultados al usar la banda Ks , lo cual era uno de los objetivos de este proyecto.

   Por otro lado, comparamos los resultados de usar o no la banda Ks en todo los ran-
gos de redshift. En esta comparación, identificamos galaxias para las cuales el photo-z
cambia mucho si usamos o no la banda de SHARKS. De esta comparación, como es-
perábamos en un principio, concluimos que la banda Ks sirve para romper degenera-
ciones entre mejores ajustes, y ası́, disminuir el número de resultados catastróficos.

    En definitiva, podemos concluir que usar datos obtenidos en el infrarrojo cercano,
particularmente en la banda Ks , ayuda a mejorar la determinación del desplazamiento
al rojo, respecto a la misma medida usando datos únicamente del óptico.

                                            31
También puede leer