Estimating the sky background is critical for ground-based astronomical research. In the optical, scattered moonlight dominates the sky background, when the moon is above the horizon. The most uncertain component of a scattered moonlight model is the aerosol scattering. The current sky background model for Cerro Paranal uses an extrapolated aerosol extinction curve. With a set of X-Shooter sky observations, we have tested the current model as well as determined the aerosol extinction from the ultra-violet to near-infrared. To our knowledge, this is the first time that a scattered moonlight model has been used for this purpose. These observations were taken of blank sky, during three different lunar phases, and at six different angular distances from the moon for each night/lunar phase. Overall, the current model does reproduce the observations for average conditions decently well. Using a set of sky background models with varying aerosol distributions to compare with the observations, we found the most likely aerosol extinction curves, phase functions, and volume densities for the three nights of observations and compare them with the current model. While there were some degeneracies in the aerosol scattering properties, the extinction curves tend to flatten towards redder wavelengths and are overall less steep compared to the extrapolated curve used in the current model. Also, the current model had significantly less coarse particles compared to the favored volume densities from the X-Shooter data. For the three nights of sky observations, the aerosol size distributions differed, most likely reflecting the changes in atmospheric conditions and aerosol content. In short, the current sky background model is in fair agreement with the observations, and we have determined better aerosol extinction curves and phase functions for Cerro Paranal.

Many round or nearly roundish Planetary Nebulae (PNe) show multiple shells and halo structures during their evolutionary stage near the maximum temperature of their central star. Controversial debate is currently ongoing if these structures are recombination halos, as suggested by hydrodynamic modelling efforts, or ionized material. Recently we discovered a halo with even somewhat unusual structures around the sparsely studied PN IC~5148 and present for the first time spectroscopy going out to the halo of such a PN.} resolution spectroscopy is used to derive dust chemistry and mineralogy. We investigate the spatial distribution of material and its ionization state from the center of the nebula up to the very outskirts of the halo. We obtained long-slit low resolution spectroscopy (FORS2@VLT) of the nebula in two position angles, which we used to investigate the nebular structure and its halo in the optical range from 450 to 880\,nm. In addition we used medium resolution spectra taken with X-SHOOTER@VLT ranging from 320 nm to 2.4 mu to derive atmospheric parameters for the central star. We obtained the distance and position in the Galaxy from various methods combined with GAIA DR2 data. We also applied Cloudy models to the nebula in order to derive physical parameters of the various regions. We obtained spatially resolved structures and detailed descriptions of the outrunning shock front and a set of unusual halo structures denoted to further shock. The halo structures appears clearly as hot ionized material. Furthermore we derived a reliable photometric value for the central star at a GAIA distance of D=1.3kpc. Considering the large distance $z=1.0$\,kpc from the galactic plane together to its non-circular motion in the galaxy and, a metallicity only slightly below that of typical disk PNe, most likely IC 5148 originates from a thick disk population star.

MICADO will enable the ELT to perform diffraction limited near-infrared observations at first light. The instrument's capabilities focus on imaging (including astrometric and high contrast) as well as single object spectroscopy. This contribution looks at how requirements from the observing modes have driven the instrument design and functionality. Using examples from specific science cases, and making use of the data simulation tool, an outline is presented of what we can expect the instrument to achieve.

During the start of a survey program using FORS2 long slit spectroscopy on planetary nebulae (PN) and their haloes, we serendipitously discovered six background emission line galaxies (ELG) with redshifts of z = 0.2057, 0.3137, 0.37281, 0.4939, 0.7424 and 0.8668. Thus they clearly do not belong to a common cluster structure. We derived the major physical properties of the targets. Since the used long slit covers a sky area of only 570 arcsec^2, we discuss further potential of serendipitous discoveries in archival data, beside the deep systematic work of the ongoing and upcoming big surveys. We conclude that archival data provide a decent potential for extending the overall data on ELGs without any selection bias.

The high-resolution echelle spectrograph UVES of the Very Large Telescope at Cerro Paranal in Chile has been regularly operated since April 2000. Thus, UVES archival data originally taken for astronomical projects but also including sky emission can be used to study airglow variations on a time scale longer than a solar cycle. Focusing on OH emission and observations until March 2015, we considered about 3,000 high-quality spectra from two instrumental set-ups centred on 760 and 860 nm, which cover about 380 nm each. These data allowed us to measure line intensities for several OH bands in order to derive band intensities and rotational temperatures for different upper vibrational levels as a function of solar activity and observing date. The results were compared with those derived from emission and temperature profile data of the radiometer SABER on the TIMED satellite taken in the Cerro Paranal area between 2002 and 2015. In agreement with the SABER data, the long-term variations in OH intensity and temperature derived from the UVES data are dominated by the solar cycle, whereas secular trends appear to be negligible. Combining the UVES and SABER results, the solar cycle effects for the OH intensity and temperature are about 12 to 17% and 4 to 5 K per 100 sfu and do not significantly depend on the selected OH band. The data also reveal that variations of the effective OH emission layer height and air density can cause significant changes in the OH rotational temperatures due to a varying ratio of OH thermalising collisions by air molecules and OH radiation, deactivation, and destruction processes which impede the rotational relaxation. However, this effect appears to be of minor importance for the explanation of the rotational temperature variations related to the solar activity cycle, which causes only small changes in the OH emission profile.

MICADO will be the first-light wide-field imager for the European Extremely Large Telescope (E-ELT) and will provide difiraction limited imaging (7mas at 1.2mm) over a ~53 arcsecond field of view. In order to support various consortium activities we have developed a first version of SimCADO: an instrument simulator for MICADO. SimCADO uses the results of the detailed simulation efforts conducted for each of the separate consortium-internal work packages in order to generate a model of the optical path from source to detector readout. SimCADO is thus a tool to provide scientific context to both the science and instrument development teams who are ultimately responsible for the final design and future capabilities of the MICADO instrument. Here we present an overview of the inner workings of SimCADO and outline our plan for its further development.

MICADO will equip the E-ELT with a first light capability for diffraction limited imaging at near-infrared wavelengths. The instrument's observing modes focus on various flavours of imaging, including astrometric, high contrast, and time resolved. There is also a single object spectroscopic mode optimised for wavelength coverage at moderately high resolution. This contribution provides an overview of the key functionality of the instrument, outlining the scientific rationale for its observing modes. The interface between MICADO and the adaptive optics system MAORY that feeds it is summarised. The design of the instrument is discussed, focussing on the optics and mechanisms inside the cryostat, together with a brief overview of the other key sub-systems.

Rotational temperatures Trot derived from lines of the same OH band are an important method to study the mesopause region near 87 km. To measure realistic temperatures, the rotational level populations have to be in local thermodynamic equilibrium (LTE). However, this might not be fulfilled, especially at high emission altitudes. In order to quantify possible non-LTE contributions to the OH Trot as a function of the upper vibrational level v', we studied a sample of 343 echelle spectra taken with the X-shooter spectrograph at the Very Large Telescope at Cerro Paranal in Chile. These data allowed us to analyse 25 OH bands in each spectrum. Moreover, we could measure lines of O2b(0-1), which peaks at 94 to 95 km, and O2a(0-0) with an emission peak at about 90 km. Since the radiative lifetimes are relatively long, the derived O2 Trot are not significantly affected by non-LTE contributions. For a comparison with OH, the differences in the emission profiles were corrected by using OH emission, O2a(0-0) emission, and CO2-based temperature profile data from the multi-channel radiometer SABER on the TIMED satellite. For a reference profile at 90 km, we found a good agreement of the O2 with the SABER-related temperatures, whereas the OH temperatures, especially for the high and even v', showed significant excesses with a maximum of more than 10 K for v' = 8. We could also find a nocturnal trend towards higher non-LTE effects, particularly for high v'.

Aims: We studied the precipitable water vapour (PWV) content near Cerro Armazones and discuss the potential use of our technique of modelling the telluric absorbtion lines for the investigation of other molecular layers. The site is designated for the European Extremely Large Telescope (E-ELT) and the nearby planned site for the Cherenkov Telescope Array (CTA). Methods: Spectroscopic data from the Bochum Echelle Spectroscopic Observer (BESO) instrument were investigated by using line-by-line radiative transfer model (LBLRTM) radiative transfer models for the Earths atmosphere with the telluric absorption correction tool molecfit. All observations from the archive in the period from December 2008 to the end of 2014 were investigated. The dataset completely covers the El Nino event registered in the period 2009-2010. Models of the 3D Global Data Assimilation System (GDAS) were used for further comparison. Moreover, for those days with coincidence of data from a similar study with VLT/X-shooter and microwave radiometer LHATPRO data at Cerro Paranal, a direct comparison is presented. Results: This analysis shows that the site has systematically lower PWV values, even after accounting for the decrease in PWV expected from the higher altitude of the site with respect to Cerro Paranal, using the average atmosphere found with radiosondes. We found that GDAS data are not a suitable method for predicting of local atmospheric conditions - they usually systematically overestimate the PWV values. Due to the large sample, we were furthermore able to characterize the site with respect to symmetry across the sky and variation with the years and within the seasons. This kind of technique of studying the atmospheric conditions is shown to be a promising step into a possible monitoring equipment for CTA.

Context: Absorption by molecules in the Earth's atmosphere strongly affects ground-based astronomical observations. The resulting absorption line strength and shape depend on the highly variable physical state of the atmosphere, i.e. pressure, temperature, and mixing ratio of the different molecules involved. Usually, supplementary observations of so-called telluric standard stars (TSS) are needed to correct for this effect, which is expensive in terms of telescope time. We have developed the software package molecfit to provide synthetic transmission spectra based on parameters obtained by fitting narrow ranges of the observed spectra of scientific objects. These spectra are calculated by means of the radiative transfer code LBLRTM and an atmospheric model. In this way, the telluric absorption correction for suitable objects can be performed without any additional calibration observations of TSS. Aims: We evaluate the quality of the telluric absorption correction using molecfit with a set of archival ESO-VLT X-Shooter visible and near-infrared spectra. Methods: Thanks to the wavelength coverage from the U to the K band, X-Shooter is well suited to investigate the quality of the telluric absorption correction with respect to the observing conditions, the instrumental set-up, input parameters of the code, the signal-to-noise of the input spectrum, and the atmospheric profiles. These investigations are based on two figures of merit, I_off and I_res, that describe the systematic offsets and the remaining small-scale residuals of the corrections. We also compare the quality of the telluric absorption correction achieved with moelcfit to the classical method based on a telluric standard star. (Abridged)

Context: The interaction of the light from astronomical objects with the constituents of the Earth's atmosphere leads to the formation of telluric absorption lines in ground-based collected spectra. Correcting for these lines, mostly affecting the red and infrared region of the spectrum, usually relies on observations of specific stars obtained close in time and airmass to the science targets, therefore using precious observing time. Aims: We present molecfit, a tool for correcting for telluric absorption lines based on synthetic modelling of the Earth's atmospheric transmission. Molecfit is versatile and can be used with data obtained with various ground-based telescopes and instruments. Methods: Molecfit combines a publicly available radiative transfer code, a molecular line database, atmospheric profiles, and various kernels to model the instrument line spread function. The atmospheric profiles are created by merging a standard atmospheric profile representative of a given observatory's climate, of local meteorological data, and of dynamically retrieved altitude profiles for temperature, pressure, and humidity. We discuss the various ingredients of the method, its applicability, and its limitations. We also show examples of telluric line correction on spectra obtained with a suite of ESO Very Large Telescope (VLT) instruments. Results: Compared to previous similar tools, molecfit takes the best results for temperature, pressure, and humidity in the atmosphere above the observatory into account. As a result, the standard deviation of the residuals after correction of unsaturated telluric lines is frequently better than 2% of the continuum. Conclusion: Molecfit is able to accurately model and correct for telluric lines over a broad range of wavelengths and spectral resolutions. (Abridged)

While the near-infrared wavelength regime is becoming more and more important for astrophysics there is a marked lack of spectrophotometric standard star data that would allow the flux calibration of such data. Furthermore, flux calibrating medium- to high-resolution \'echelle spectroscopy data is challenging even in the optical wavelength range, because the available flux standard data are often too coarsely sampled. We will provide standard star reference data that allow users to derive response curves from 300nm to 2500nm for spectroscopic data of medium to high resolution, including those taken with \'echelle spectrographs. In addition we describe a method to correct for moderate telluric absorption without the need of observing telluric standard stars. As reference data for the flux standard stars we use theoretical spectra derived from stellar model atmospheres. We verify that they provide an appropriate description of the observed standard star spectra by checking for residuals in line cores and line overlap regions in the ratios of observed (X-shooter) spectra to model spectra. The finally selected model spectra are then corrected for remaining mismatches and photometrically calibrated using independent observations. The correction of telluric absorption is performed with the help of telluric model spectra.We provide new, finely sampled reference spectra without telluric absorption for six southern flux standard stars that allow the users to flux calibrate their data from 300 nm to 2500 nm, and a method to correct for telluric absorption using atmospheric models.

Airglow emission lines, which dominate the optical-to-near-IR sky radiation, show strong, line-dependent variability on various time scales. Therefore, the subtraction of the sky background in the affected wavelength regime becomes a problem if plain sky spectra have to be taken at a different time as the astronomical data. A solution of this issue is the physically motivated scaling of the airglow lines in the plain sky data to fit the sky lines in the object spectrum. We have developed a corresponding instrument-independent approach based on one-dimensional spectra. Our code skycorr separates sky lines and sky/object continuum by an iterative approach involving a line finder and airglow line data. The sky lines are grouped according to their expected variability. The line groups in the sky data are then scaled to fit the sky in the science data. Required pixel-specific weights for overlapping groups are taken from a comprehensive airglow model. Deviations in the wavelength calibration are corrected by fitting Chebyshev polynomials and rebinning via asymmetric damped sinc kernels. The scaled sky lines and the sky continuum are subtracted separately. VLT X-Shooter data covering time intervals from two minutes to about one year were selected to illustrate the performance. Except for short time intervals of a few minutes, the sky line residuals were several times weaker than for sky subtraction without fitting. Further tests show that skycorr performs consistently better than the method of Davies (2007) developed for VLT SINFONI data.

Correcting for the sky signature usually requires supplementary calibration data which are very expensive in terms of telescope time. In addition, the scheduling flexibility is restricted as these data have to be taken usually directly before/after the science observations due to the high variability of the telluric absorption which depends on the state and the chemical composition of the atmosphere at the time of observations. Therefore, a tool for sky correction, which does not require this supplementary calibration data, saves a significant amount of valuable telescope time and increases its efficiency. We developed a software package aimed at performing telluric feature corrections on the basis of synthetic absorption spectra.

The largest natural source of light at night is the Moon, and it is the major contributor to the astronomical sky background. Being able to accurately predict the sky background, including scattered moonlight is important for scheduling astronomical observations. We have developed an improved scattered moonlight model, in which the components are computed with a better physical understanding as opposed to the simple empirical fit in the frequently used photometric model of Krisciunas & Schaefer (1991). Our spectroscopic model can better trace the spectral trends of scattered moonlight for any position of the Moon and target observation. This is the first scattered moonlight model that we know of which is this physical and versatile. We have incorporated an observed solar spectrum, accurate lunar albedo fit, and elaborate scattering and absorption calculations that include scattering off of molecules and aerosols. It was designed for Cerro Paranal, but can be modified for any location with known atmospheric properties. Throughout the optical range, the uncertainty is less than 20%. This advanced scattered moonlight model can predict the amount of scattered moonlight for any given geometry of the Moon and target, and lunar phase for the entire optical spectrum.

We calculate harmonic frequencies of the three most abundant carbonic acid conformers. For this, different model chemistries are investigated with respect to their benefits and shortcomings. Based on these results we use perturbation theory to calculate anharmonic corrections at the {\omega}B97XD/aug-cc-pVXZ, X=D,T,Q, level of theory and compare them with recent experimental data and theoretical predictions. A discrete variable representation method is used to predict the large anharmonic contributions to the frequencies of the stretching vibrations in the hydrogen bonds in the carbonic acid dimer. Moreover, we re-investigate the energetics of the formation of the carbonic acid dimer from its constituents water and carbon dioxide using a high-level extrapolation method. We find that the {\omega}B97XD functional performs well in estimating the fundamental frequencies of the carbonic acid conformers. Concerning the reaction energetics, the accuracy of {\omega}B97XD is even comparable to the high-level extrapolation method. We discuss possibilities to detect carbonic acid in various natural environments such as Earth's and Martian atmospheres.

The Earth's atmosphere affects ground-based astronomical observations. Scattering, absorption, and radiation processes deteriorate the signal-to-noise ratio of the data received. For scheduling astronomical observations it is, therefore, important to accurately estimate the wavelength-dependent effect of the Earth's atmosphere on the observed flux. In order to increase the accuracy of the exposure time calculator of the European Southern Observatory's (ESO) Very Large Telescope (VLT) at Cerro Paranal, an atmospheric model was developed as part of the Austrian ESO In-Kind contribution. It includes all relevant components, such as scattered moonlight, scattered starlight, zodiacal light, atmospheric thermal radiation and absorption, and non-thermal airglow emission. This paper focuses on atmospheric scattering processes that mostly affect the blue (< 0.55 mum) wavelength regime, and airglow emission lines and continuum that dominate the red (> 0.55 mum) wavelength regime. While the former is mainly investigated by means of radiative transfer models, the intensity and variability of the latter is studied with a sample of 1186 VLT FORS1 spectra. For a set of parameters such as the object altitude angle, Moon-object angular distance, ecliptic latitude, bimonthly period, and solar radio flux, our model predicts atmospheric radiation and transmission at a requested resolution. A comparison of our model with the FORS1 spectra and photometric data for the night-sky brightness from the literature, suggest a model accuracy of about 20%. This is a significant improvement with respect to existing predictive atmospheric models for astronomical exposure time calculators.

We present a search for gravitational arcs in a sample of X-ray luminous, medium redshift clusters of galaxies. The sample of clusters is called ARCRAIDER, is based on the ROSAT Bright Survey (RBS) and fulfills the following criteria: (a) X-ray luminosity Lx>=0.5x10^45erg/s (0.5-2keV band), (b) redshift range 0.1<=z<=0.52, (c) classified as clusters in the RBS, (d) not a member of the Abell catalogue and, finally, (e) visible from the ESO sites La Silla/Paranal (declination \delta<=20deg). In total we found more than 35 (giant) arc/arclet candidates, including a possible radial arc, one galaxy-galaxy lensing event and a possible quasar triple image in 14 of the 21 clusters of galaxies. Hence 66% of the sample members are possible lenses.

We present a detailed analysis of the medium redshift (z=0.2906) galaxy cluster Z3146 which is part of the ongoing ARCRAIDER project, a systematic search for gravitational arcs in massive clusters of galaxies. The analysis of Z3146 is based on deep optical wide field observations in the B, V and R bands obtained with the WFI@ESO2.2m, and shallow archival WFPC2@HST taken with the F606W filter, which are used for strong as well as weak lensing analyses. Additionally we have used publicly available XMM/Newton observations for a detailed X-ray analysis of Z3146. Both methods, lensing and X-ray, were used to determine the dynamical state and to estimate the total mass. We also identified four gravitational arc candidates. We find this cluster to be in a relaxed state, which is confirmed by a large cooling flow with nominal ~1600M_\odot per year, regular galaxy density and light distributions and a regular shape of the weak lensing mass reconstruction. The mass content derived with the different methods agrees well within 25% at r_{200}=1661 h_{70}^{-1}kpc indicating a velocity dispersion of \sigma_v=869^{+124}_{-153}km/s.

We present the first three galaxy clusters of a larger sample of the most X-ray luminous galaxy clusters selected from the ROSAT Bright Survey. This project, which is a systematic search for strong lensing, aims at arc statistics, mass determinations and studies of distant lensed galaxies. The three galaxy clusters presented here have been observed with the Wide Field Imager at the ESO2.2m in the R- and V-band. The images show lensing features like distinct distorted galaxies and arcs. Mass distributions of the lensing galaxy clusters and photometric properties of some arc candidates are presented. In addition we report the discovery of three giant arcs.

We present a detailed spectroscopic investigation, spanning four winters, of the asymptotic giant branch (AGB) star IRAS 02091+6333. Zijlstra & Weinberger (2002) found a giant wall of dust around this star and modelled this unique phenomenon. However their work suffered from the quality of the optical investigations of the central object. Our spectroscopic investigation allowed us to define the spectral type and the interstellar foreground extinction more precisely. Accurate multi band photometry was carried out. This provides us with the possibility to derive the physical parameters of the system. The measurements presented here suggest a weak irregular photometric variability of the target, while there is no evidence of a spectroscopic variability over the last four years.