We present a comparison of several Difference Image Analysis (DIA) techniques, in combination with Machine Learning (ML) algorithms, applied to the identification of optical transients associated with gravitational wave events. Each technique is assessed based on the scoring metrics of Precision, Recall, and their harmonic mean F1, measured on the DIA results as standalone techniques, and also in the results after the application of ML algorithms, on transient source injections over simulated and real data. This simulations cover a wide range of instrumental configurations, as well as a variety of scenarios of observation conditions, by exploring a multi dimensional set of relevant parameters, allowing us to extract general conclusions related to the identification of transient astrophysical events. The newest subtraction techniques, and particularly the methodology published in Zackay et al. (2016) are implemented in an Open Source Python package, named properimage, suitable for many other astronomical image analyses. This together with the ML libraries we describe, provides an effective transient detection software pipeline. Here we study the effects of the different ML techniques, and the relative feature importances for classification of transient candidates, and propose an optimal combined strategy. This constitutes the basic elements of pipelines that could be applied in searches of electromagnetic counterparts to GW sources.

We report the first confirmed detection of the galaxy cluster VVV-J144321-611754 at very low latitudes (l = 315.836$^{\circ}$, b = -1.650$^{\circ}$) located in the tile d015 of the VISTA Variables in the V\'ia L\'actea (VVV) survey. We defined the region of 30$\times$ 30 $arcmin^2$ centered in the brightest galaxy finding 25 galaxies. For these objects, extinction-corrected median colors of (H - K$_{s}$) = 0.34 $\pm$ 0.05 mag, (J - H) = 0.57 $\pm$ 0.08 mag and (J - K$_{s}$) = 0.87 $\pm$ 0.06 mag, and R$_{1/2}$ = 1.59 $\pm$ 0.16 $arcsec$; C = 3.01 $\pm$ 0.08; and Sersic index, n = 4.63 $\pm$ 0.39 were estimated. They were visually confirmed showing characteristics of early-type galaxies in the near-IR images. An automatic clustering analysis performed in the whole tile found that the concentration of galaxies VVV-J144321-611754 is a real, compact concentration of early-type galaxies. Assuming a typical galaxy cluster with low X-ray luminosity, the photometric redshift of the brightest galaxy is $z = $ 0.196 $\pm$ 0.025. Follow-up near-IR spectroscopy with FLAMINGOS-2 at the Gemini-South telescope revealed that the two brighter cluster galaxies have typical spectra of early-type galaxies and the estimated redshift for the brightest galaxy VVV-J144321.06-611753.9 is $z =$ 0.234$\pm$0.022 and for VVV-J144319.02-611746.1 is $z =$ 0.232$\pm$0.019. Finally, these galaxies clearly follow the cluster Red Sequence in the rest-frame near-IR color--magnitude diagram with the slope similar to galaxy cluster at redshift of 0.2. These results are consistent with the presence of a bona fide galaxy cluster beyond the Milky Way disk.

We present the results of the optical follow-up, conducted by the TOROS collaboration, of gravitational wave events detected during the Advanced LIGO-Virgo second observing run (Nov 2016 -- Aug 2017). Given the limited field of view ($\sim100\arcmin$) of our observational instrumentation we targeted galaxies within the area of high localization probability that were observable from our sites. We analyzed the observations using difference imaging, followed by a Random Forest algorithm to discriminate between real and bogus transients. For all three events that we respond to, except GW170817, we did not find any bona fide optical transient that was plausibly linked with the observed gravitational wave event. Our observations were conducted using telescopes at Estaci\'{o}n Astrof\'{\i}sica de Bosque Alegre, Cerro Tololo Inter-American Observatory, and the Dr. Cristina V. Torres Memorial Astronomical Observatory. Our results are consistent with the LIGO-Virgo detections of a binary black hole merger (GW170104) for which no electromagnetic counterparts were expected, as well as a binary neutron star merger (GW170817) for which an optical transient was found as expected.

Machine learning algorithms are highly useful for the classification of time series data in astronomy in this era of peta-scale public survey data releases. These methods can facilitate the discovery of new unknown events in most astrophysical areas, as well as improving the analysis of samples of known phenomena. Machine learning algorithms use features extracted from collected data as input predictive variables. A public tool called Feature Analysis for Time Series (FATS) has proved an excellent workhorse for feature extraction, particularly light curve classification for variable objects. In this study, we present a major improvement to FATS, which corrects inconvenient design choices, minor details, and documentation for the re-engineering process. This improvement comprises a new Python package called "feets", which is important for future code-refactoring for astronomical software tools.

We present the results of prompt optical follow-up of the electromagnetic counterpart of the gravitational-wave event GW170817 by the Transient Optical Robotic Observatory of the South Collaboration (TOROS). We detected highly significant dimming in the light curves of the counterpart (Delta g=0.17+-0.03 mag, Delta r=0.14+-0.02 mag, Delta i=0.10 +- 0.03 mag) over the course of only 80 minutes of observations obtained ~35 hr after the trigger with the T80-South telescope. A second epoch of observations, obtained ~59 hr after the event with the EABA 1.5m telescope, confirms the fast fading nature of the transient. The observed colors of the counterpart suggest that this event was a "blue kilonova" relatively free of lanthanides.

Data processing pipelines represent an important slice of the astronomical software library that include chains of processes that transform raw data into valuable information via data reduction and analysis. In this work we present Corral, a Python framework for astronomical pipeline generation. Corral features a Model-View-Controller design pattern on top of an SQL Relational Database capable of handling: custom data models; processing stages; and communication alerts, and also provides automatic quality and structural metrics based on unit testing. The Model-View-Controller provides concept separation between the user logic and the data models, delivering at the same time multi-processing and distributed computing capabilities. Corral represents an improvement over commonly found data processing pipelines in Astronomy since the design pattern eases the programmer from dealing with processing flow and parallelization issues, allowing them to focus on the specific algorithms needed for the successive data transformations and at the same time provides a broad measure of quality over the created pipeline. Corral and working examples of pipelines that use it are available to the community at https://github.com/toros-astro.

We present the results of the optical follow-up conducted by the TOROS collaboration of the first gravitational-wave event GW150914. We conducted unfiltered CCD observations (0.35-1 micron) with the 1.5-m telescope at Bosque Alegre starting ~2.5 days after the alarm. Given our limited field of view (~100 square arcmin), we targeted 14 nearby galaxies that were observable from the site and were located within the area of higher localization probability. We analyzed the observations using two independent implementations of difference-imaging algorithms, followed by a Random-Forest-based algorithm to discriminate between real and bogus transients. We did not find any bona fide transient event in the surveyed area down to a 5-sigma limiting magnitude of r=21.7 mag (AB). Our result is consistent with the LIGO detection of a binary black hole merger, for which no electromagnetic counterparts are expected, and with the expected rates of other astrophysical transients.

FRIDA is a diffraction-limited imager and integral-field spectrometer that is being built for the adaptive-optics focus of the Gran Telescopio Canarias. In imaging mode FRIDA will provide scales of 0.010, 0.020 and 0.040 arcsec/pixel and in IFS mode spectral resolutions of 1500, 4000 and 30,000. FRIDA is starting systems integration and is scheduled to complete fully integrated system tests at the laboratory by the end of 2017 and to be delivered to GTC shortly thereafter. In this contribution we present a summary of its design, fabrication, current status and potential scientific applications.