We develop a data-driven model to map stellar parameters (effective temperature, surface gravity and metallicity) accurately and precisely to broad-band stellar photometry. This model must, and does, simultaneously constrain the passband-specific dust reddening vector in the Milky Way. The model uses a neural network to learn the (de-reddened) absolute magnitude in one band and colors across many bands, given stellar parameters from spectroscopic surveys and parallax constraints from Gaia. To demonstrate the effectiveness of this approach, we train our model on a dataset with spectroscopic labels from LAMOST, APOGEE and GALAH, Gaia parallaxes, and optical and near-infrared photometry from Gaia, Pan-STARRS 1, 2MASS and WISE. Testing the model on these datasets leads to an excellent fit and a precise -- and by construction accurate -- prediction of the color-magnitude diagrams in many bands. This flexible approach rigorously links spectroscopic and photometric surveys, and also results in an improved, temperature-dependent reddening vector. As such, it provides a simple and accurate method for predicting photometry in stellar evolutionary models. Our model will form a basis to infer stellar properties, distances and dust extinction from photometric data, which should be of great use in 3D mapping of the Milky Way. Our trained model may be obtained at https://doi.org/10.5281/zenodo.3902382.

In a cosmological setting, the disc of a galaxy is expected to continuously experience gravitational torques and perturbations from a variety of sources, which can cause the disc to wobble, flare and warp. Specifically, the study of galactic warps and their dynamical nature can potentially reveal key information on the formation history of galaxies and the mass distribution of their halos. Our Milky Way presents a unique case study for galactic warps, thanks to the detailed knowledge of its stellar distribution and kinematics. Using a simple model of how the warp's orientation is changing with time, we here measure the precession rate of the Milky Way's warp using 12 million giant stars from Gaia Data Release 2, finding that it is precessing at $10.86 \pm 0.03_{stat} \pm 3.20_{sys}$ km/s/kpc in the direction of Galactic rotation, about one third the angular velocity at the Sun's position in the Galaxy. The direction and magnitude of the warp's precession rate favour the scenario that the warp is the results of a recent or ongoing encounter with a satellite galaxy, rather than the relic of the ancient assembly history of the Galaxy.

We present a mock stellar catalog, matching in volume, depth and data model the content of the planned Gaia early data release 3 (Gaia EDR3). We have generated our catalog (GeDR3mock) using galaxia, a tool to sample stars from an underlying Milky Way (MW) model or from N-body data. We used an updated Besan\c{c}on Galactic model together with the latest PARSEC stellar evolutionary tracks, now also including white dwarfs. We added the Magellanic clouds and realistic open clusters with internal rotation. We empirically modelled uncertainties based on Gaia DR2 (GDR2) and scaled them according to the longer baseline in Gaia EDR3. The apparent magnitudes were reddened according to a new selection of 3D extinction maps. To help with the Gaia selection function we provide all-sky magnitude limit maps in G and BP for a few relevant GDR2 subsets together with the routines to produce these maps for user-defined subsets. We supplement the catalog with photometry and extinctions in non-Gaia bands. The catalog is available in the Virtual Observatory and can be queried just like the actual Gaia EDR3 will be. We highlight a few capabilities of the Astronomy Data Query Language (ADQL) with educative catalog queries. We use the data extracted from those queries to compare GeDR3mock to GDR2, which emphasises the importance of adding observational noise to the mock data. Since the underlying truth, e.g. stellar parameters, is know in GeDR3mock, it can be used to construct priors as well as mock data tests for parameter estimation. All code, models and data used to produce GeDR3mock are linked and contained in galaxia_wrap, a python package, representing a fast galactic forward model, able to project MW models and N-body data into realistic Gaia observables.

We construct a supervised classifier based on Gaussian Mixture Models to probabilistically classify objects in Gaia data release 2 (GDR2) using only photometric and astrometric data in that release. The model is trained empirically to classify objects into three classes -- star, quasar, galaxy -- for G<=14.5 mag down to the Gaia magnitude limit of G=21.0 mag. Galaxies and quasars are identified for the training set by a cross-match to objects with spectroscopic classifications from the Sloan Digital Sky Survey. Stars are defined directly from GDR2. When allowing for the expectation that quasars are 500 times rarer than stars, and galaxies 7500 times rarer than stars (the class imbalance problem), samples classified with a threshold probability of 0.5 are predicted to have purities of 0.43 for quasars and 0.28 for galaxies, and completenesses of 0.58 and 0.72 respectively. The purities can be increased up to 0.60 by adopting a higher threshold. Not accounting for this expected low frequency of extragalactic objects (the class prior) would give both erroneously optimistic performance predictions and severely impure samples. Applying our model to all 1.20 billion objects in GDR2 with the required features, we classify 2.3 million objects as quasars and 0.37 million objects as galaxies (with individual probabilities above 0.5). The small number of galaxies is due to the strong bias of the satellite detection algorithm and on-ground data selection against extended objects. We infer the true number of quasars and galaxies -- as these classes are defined by our training set -- to be 690,000 and 110,000 respectively (+/- 50%). The aim of this work is to see how well extragalactic objects can be classified using only GDR2 data. Better classifications should be possible with the low resolution spectroscopy (BP/RP) planned for GDR3.

The mechanisms of the formation and evolution of the Milky Way are encoded in the orbits, chemistry and ages of its stars. With the 4MOST MIlky way Disk And BuLgE Low-Resolution Survey (4MIDABLE-LR) we aim to study kinematic and chemical substructures in the Milky Way disc and bulge region with samples of unprecedented size out to larger distances and greater precision than conceivable with Gaia alone or any other ongoing or planned survey. Gaia gives us the unique opportunity for target selection based almost entirely on parallax and magnitude range, hence increasing the efficiency in sampling larger Milky Way volumes with well-defined and effective selection functions. Our main goal is to provide a detailed chrono-chemo-kinematical extended map of our Galaxy and the largest Gaia follow-up down to $G = 19$ magnitudes (Vega). The complex nature of the disc components (for example, large target densities and highly structured extinction distribution in the Milky Way bulge and disc area), prompted us to develop a survey strategy with five main sub-surveys that are tailored to answer the still open questions about the assembly and evolution of our Galaxy, while taking full advantage of the Gaia data.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Using Gaia DR2 astrometry, we map the kinematic signature of the Galactic stellar warp out to a distance of 7 kpc from the Sun. Combining Gaia DR2 and 2MASS photometry, we identify, via a probabilistic approach, 599 494 upper main sequence stars and 12 616 068 giants without the need for individual extinction estimates. The spatial distribution of the upper main sequence stars clearly shows segments of the nearest spiral arms. The large-scale kinematics of both the upper main sequence and giant populations show a clear signature of the warp of the Milky Way, apparent as a gradient of 5-6 km/s in the vertical velocities from 8 to 14 kpc in Galactic radius. The presence of the signal in both samples, which have different typical ages, suggests that the warp is a gravitationally induced phenomenon.

For the vast majority of stars in the second Gaia data release, reliable distances cannot be obtained by inverting the parallax. A correct inference procedure must instead be used to account for the nonlinearity of the transformation and the asymmetry of the resulting probability distribution. Here we infer distances to essentially all 1.33 billion stars with parallaxes published in the second \gaia\ data release. This is done using a weak distance prior that varies smoothly as a function of Galactic longitude and latitude according to a Galaxy model. The irreducible uncertainty in the distance estimate is characterized by the lower and upper bounds of an asymmetric confidence interval. Although more precise distances can be estimated for a subset of the stars using additional data (such as photometry), our goal is to provide purely geometric distance estimates, independent of assumptions about the physical properties of, or interstellar extinction towards, individual stars. We analyse the characteristics of the catalogue and validate it using clusters. The catalogue can be queried on the Gaia archive using ADQL at http://gea.esac.esa.int/archive/ and downloaded from http://www.mpia.de/~calj/gdr2_distances.html .

Passing stars may play an important role in the evolution of our solar system. We search for close stellar encounters to the Sun among all 7.2 million stars in Gaia-DR2 that have six-dimensional phase space data. We characterize encounters by integrating their orbits through a Galactic potential and propagating the correlated uncertainties via a Monte Carlo resampling. After filtering to remove spurious data, we find 694 stars that have median (over uncertainties) closest encounter distances within 5 pc, all occurring within 15 Myr from now. 26 of these have at least a 50% chance of coming closer than 1 pc (and 7 within 0.5 pc), all but one of which are newly discovered here. We further confirm some and refute several other previously-identified encounters, confirming suspicions about their data. The closest encounter in the sample is Gl 710, which has a 95% probability of coming closer than 0.08 pc (17 000 AU). Taking mass estimates from Gaia astrometry and multiband photometry for essentially all encounters, we find that Gl 710 also has the largest impulse on the Oort cloud. Using a Galaxy model, we compute the completeness of the Gaia-DR2 encountering sample as a function of perihelion time and distance. Only 15% of encounters within 5 pc occurring within +/- 5 Myr of now have been identified, mostly due to the lack of radial velocities for faint and/or cool stars. Accounting for the incompleteness, we infer the present rate of encounters within 1 pc to be 19.7 +/- 2.2 per Myr, a quantity expected to scale quadratically with the encounter distance out to at least several pc. Spuriously large parallaxes in our sample from imperfect filtering would tend to inflate both the number of encounters found and this inferred rate. The magnitude of this effect is hard to quantify.

The second Gaia data release (Gaia-DR2) contains, beyond the astrometry, three-band photometry for 1.38 billion sources. We have used these three broad bands to infer stellar effective temperatures, Teff, for all sources brighter than G=17 mag with Teff in the range 3000-10 000 K (161 million sources). Using in addition the parallaxes, we infer the line-of-sight extinction, A_G, and the reddening, E[BP-RP], for 88 million sources. Together with a bolometric correction we derive luminosity and radius for 77 million sources. These quantities as well as their estimated uncertainties are part of Gaia-DR2. Here we describe the procedures by which these quantities were obtained, including the underlying assumptions, comparison with literature estimates, and the limitations of our results. Typical accuracies are of order 324 K (Teff), 0.46 mag (A_G), 0.23 mag (E[BP-RP]), 15% (luminosity), and 10% (radius). Being based on only a small number of observable quantities and limited training data, our results are necessarily subject to some extreme assumptions that can lead to strong systematics in some cases (not included in the aforementioned accuracy estimates). One aspect is the non-negativity contraint of our estimates, in particular extinction. Yet in several regions of parameter space our results show very good performance, for example for red clump stars and solar analogues. Large uncertainties render the extinctions less useful at the individual star level, but they show good performance for ensemble estimates. We identify regimes in which our parameters should and should not be used and we define a "clean" sample. Despite the limitations, this is the largest catalogue of uniformly-inferred stellar parameters to date. More precise and detailed astrophysical parameters based on the full BP/RP spectrophotometry are planned as part of the third Gaia data release.

We present a mock catalog of Milky Way stars, matching in volume and depth the content of the Gaia data release 2 (GDR2). We generated our catalog using Galaxia, a tool to sample stars from a Besancon Galactic model, together with a realistic 3D dust extinction map. The catalog mimicks the complete GDR2 data model and contains most of the entries in the Gaia source catalog: 5-parameter astrometry, 3-band photometry, radial velocities, stellar parameters, and associated scaled nominal uncertainty estimates. In addition, we supplemented the catalog with extinctions and photometry for non-Gaia bands. This catalog can be used to prepare GDR2 queries in a realistic runtime environment, and it can serve as a Galactic model against which to compare the actual GDR2 data in the space of observables. The catalog is hosted through the virtual observatory GAVO's Heidelberg data center service and thus can be queried using ADQL as for GDR2 data.

From exploratory studies and theoretical expectations it is known that simplifying approximations in spectroscopic analysis (LTE, 1D) lead to systematic biases of stellar parameters and abundances. These biases depend strongly on surface gravity, temperature, and, in particular, for LTE vs. non-LTE (NLTE) on metallicity of the stars. Here we analyse the [Mg/Fe] and [Fe/H] plane of a sample of 326 stars, comparing LTE and NLTE results obtained using 1D hydrostatic models and averaged <3D> models. We show that compared to the <3D>NLTE benchmark, all other three methods display increasing biases towards lower metallicities, resulting in false trends of [Mg/Fe] against [Fe/H], which have profound implications for interpretations by chemical evolution models. In our best <3D> NLTE model, the halo and disc stars show a clearer behaviour in the [Mg/Fe] - [Fe/H] plane, from the knee in abundance space down to the lowest metallicities. Our sample has a large fraction of thick disc stars and this population extends down to at least [Fe/H] ~ -1.6 dex, further than previously proven. The thick disc stars display a constant [Mg/Fe] ~ 0.3 dex, with a small intrinsic dispersion in [Mg/Fe] that suggests that a fast SN Ia channel is not relevant for the disc formation. The halo stars reach higher [Mg/Fe] ratios and display a net trend of [Mg/Fe] at low metallicities, paired with a large dispersion in [Mg/Fe]. These indicate the diverse origin of halo stars from accreted low-mass systems to stochastic/inhomogeneous chemical evolution in the Galactic halo.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

The most frequently proposed model for the origin of quasars holds that the high accretion rates seen in luminous active galactic nuclei are primarily triggered during major mergers between gas-rich galaxies. While plausible for decades, this model has only begun to be tested with statistical rigor in the past few years. Here we report on a Hubble Space Telescope study to test this hypothesis for $z=2$ quasars with high super-massive black hole masses ($M_\mathrm{BH}=10^9-10^{10}~M_\odot{}$), which dominate cosmic black hole growth at this redshift. We compare Wide Field Camera 3 $F160W$ (rest-frame $V$-band) imaging of 19 point source-subtracted quasar hosts to a matched sample of 84 inactive galaxies, testing whether the quasar hosts have greater evidence for strong gravitational interactions. Using an expert ranking procedure, we find that the quasar hosts are uniformly distributed within the merger sequence of inactive galaxies, with no preference for quasars in high-distortion hosts. Using a merger/non-merger cutoff approach, we recover distortion fractions of $f_\mathrm{m,qso}=0.39\pm{}0.11$ for quasar hosts and $f_\mathrm{m,gal}=0.30\pm{}0.05$ for inactive galaxies (distribution modes, 68% confidence intervals), with both measurements subjected to the same observational conditions and limitations. The slight enhancement in distorted fraction for quasar hosts over inactive galaxies is not significant, with a probability that the quasar fraction is higher of $P(f_\mathrm{m,qso}>f_\mathrm{m,gal}) = 0.78$ ($0.78\,\sigma{}$), in line with results for lower mass and lower $z$ AGN. We find no evidence that major mergers are the primary triggering mechanism for the massive quasars that dominate accretion at the peak of cosmic quasar activity.

The Gaia satellite will survey the entire celestial sphere down to 20th magnitude, obtaining astrometry, photometry, and low resolution spectrophotometry on one billion astronomical sources, plus radial velocities for over one hundred million stars. Its main objective is to take a census of the stellar content of our Galaxy, with the goal of revealing its formation and evolution. Gaia's unique feature is the measurement of parallaxes and proper motions with hitherto unparalleled accuracy for many objects. As a survey, the physical properties of most of these objects are unknown. Here we describe the data analysis system put together by the Gaia consortium to classify these objects and to infer their astrophysical properties using the satellite's data. This system covers single stars, (unresolved) binary stars, quasars, and galaxies, all covering a wide parameter space. Multiple methods are used for many types of stars, producing multiple results for the end user according to different models and assumptions. Prior to its application to real Gaia data the accuracy of these methods cannot be assessed definitively. But as an example of the current performance, we can attain internal accuracies (RMS residuals) on F,G,K,M dwarfs and giants at G=15 (V=15-17) for a wide range of metallicites and interstellar extinctions of around 100K in effective temperature (Teff), 0.1mag in extinction (A0), 0.2dex in metallicity ([Fe/H]), and 0.25dex in surface gravity (logg). The accuracy is a strong function of the parameters themselves, varying by a factor of more than two up or down over this parameter range. After its launch in November 2013, Gaia will nominally observe for five years, during which the system we describe will continue to evolve in light of experience with the real data.

The optical light curves of many quasars show variations of tenths of a magnitude or more on time scales of months to years. This variation often cannot be described well by a simple deterministic model. We perform a Bayesian comparison of over 20 deterministic and stochastic models on 6304 QSO light curves in SDSS Stripe 82. We include the damped random walk (or Ornstein-Uhlenbeck [OU] process), a particular type of stochastic model which recent studies have focused on. Further models we consider are single and double sinusoids, multiple OU processes, higher order continuous autoregressive processes, and composite models. We find that only 29 out of 6304 QSO lightcurves are described significantly better by a deterministic model than a stochastic one. The OU process is an adequate description of the vast majority of cases (6023). Indeed, the OU process is the best single model for 3462 light curves, with the composite OU process/sinusoid model being the best in 1706 cases. The latter model is the dominant one for brighter/bluer QSOs. Furthermore, a non-negligible fraction of QSO lightcurves show evidence that not only the mean is stochastic but the variance is stochastic, too. Our results confirm earlier work that QSO light curves can be described with a stochastic model, but place this on a firmer footing, and further show that the OU process is preferred over several other stochastic and deterministic models. Of course, there may well exist yet better (deterministic or stochastic) models which have not been considered here.

Close encounters between galaxies are expected to be a viable mechanism, as predicted by numerical simulations, by which accretion onto supermassive black holes can be initiated. To test this scenario, we construct a sample of 562 galaxies (M*>2.5x10^10 M_sun) in kinematic pairs over the redshift range 0.25 < z < 1.05 that are more likely to be interacting than a well-matched control sample of 2726 galaxies not identified as being in a pair, both from the zCOSMOS 20k spectroscopic catalog. Galaxies that harbor an active galactic nucleus (AGN) are identified on the basis of their X-ray emission (L_x>2x10^42 erg s^-1) detected by Chandra. We find a higher fraction of AGN in galaxies in pairs relative to isolated galaxies of similar stellar mass. Our result is primarily due to an enhancement of AGN activity, by a factor of 1.9 (observed) and 2.6 (intrinsic), for galaxies in pairs of projected separation less than 75 kpc and line-of-sight velocity offset less than 500 km s^-1. This study demonstrates that close kinematic pairs are conducive environments for black hole growth either indicating a causal physical connection or an inherent relation, such as, to enhanced star formation. In the Appendix, we describe a method to estimate the intrinsic fractions of galaxies (either in pairs or the field) hosting an AGN with confidence intervals, and an excess fraction in pairs. We estimate that 17.8_{-7.4}^{+8.4}% of all moderate-luminosity AGN activity takes place within galaxies undergoing early stages of interaction that leaves open the question as to what physical processes are responsible for fueling the remaining ~80% that may include late-stage mergers.

Testing theories of angular-momentum acquisition of rotationally supported disc galaxies is the key to understand the formation of this type of galaxies. The tidal-torque theory tries to explain this acquisition process in a cosmological framework and predicts positive autocorrelations of angular-momentum orientation and spiral-arm handedness on distances of 1Mpc/h. This disc alignment can also cause systematic effects in weak-lensing measurements. Previous observations claimed discovering such correlations but did not account for errors in redshift, ellipticity and morphological classifications. We explain how to rigorously propagate all important errors. Analysing disc galaxies in the SDSS database, we find that positive autocorrelations of spiral-arm handedness and angular-momentum orientations on distances of 1Mpc/h are plausible but not statistically significant. This result agrees with a simple hypothesis test in the Local Group, where we find no evidence for disc alignment. Moreover, we demonstrate that ellipticity estimates based on second moments are strongly biased by galactic bulges, thereby corrupting correlation estimates and overestimating the impact of disc alignment on weak-lensing studies. Finally, we discuss the potential of future sky surveys. We argue that photometric redshifts have too large errors, i.e., PanSTARRS and LSST cannot be used. We also discuss potentials and problems of front-edge classifications of galaxy discs in order to improve estimates of angular-momentum orientation.

Parametrising galaxy morphologies is a challenging task, e.g., in shear measurements of weak lensing or investigations of galaxy evolution. The huge variety of morphologies requires an approach that is highly flexible, e.g., accounting for azimuthal structure. We revisit the method of sersiclets, where galaxy morphologies are decomposed into basis functions based on the Sersic profile. This approach is justified by the fact that the Sersic profile is the first-order Taylor expansion of any real light profile. We show that sersiclets overcome the modelling failures of shapelets. However, sersiclets implicate an unphysical relation between the steepness of the light profile and the spatial scale of azimuthal structures, which is not obeyed by real galaxy morphologies and can therefore give rise to modelling failures. Moreover, we demonstrate that sersiclets are prone to undersampling, which restricts sersiclet modelling to highly resolved galaxy images. Analysing data from the Great08 challenge, we demonstrate that sersiclets should not be used in weak-lensing studies. We conclude that although the sersiclet approach appears very promising at first glance, it suffers from conceptual and practical problems that severly limit its usefulness. The Sersic profile can be enhanced by higher-order terms in the Taylor expansion, which can drastically improve model reconstructions of galaxy images. If orthonormalised, these higher-order profiles can overcome the problems of sersiclets while preserving their mathematical justification.

Reduced chi-squared is a very popular method for model assessment, model comparison, convergence diagnostic, and error estimation in astronomy. In this manuscript, we discuss the pitfalls involved in using reduced chi-squared. There are two independent problems: (a) The number of degrees of freedom can only be estimated for linear models. Concerning nonlinear models, the number of degrees of freedom is unknown, i.e., it is not possible to compute the value of reduced chi-squared. (b) Due to random noise in the data, also the value of reduced chi-squared itself is subject to noise, i.e., the value is uncertain. This uncertainty impairs the usefulness of reduced chi-squared for differentiating between models or assessing convergence of a minimisation procedure. The impact of noise on the value of reduced chi-squared is surprisingly large, in particular for small data sets, which are very common in astrophysical problems. We conclude that reduced chi-squared can only be used with due caution for linear models, whereas it must not be used for nonlinear models at all. Finally, we recommend more sophisticated and reliable methods, which are also applicable to nonlinear models.

What is the relevance of major mergers and interactions as triggering mechanisms for active galactic nuclei (AGN) activity? To answer this longstanding question, we analyze 140 XMM-selected AGN host galaxies and a matched control sample of 1264 inactive galaxies over z~0.3-1.0 and log(M_*/M_sun)<11.7 with high-resolution HST/ACS imaging from the COSMOS field. The visual analysis of their morphologies by 10 independent human classifiers yields a measure of the fraction of distorted morphologies in the AGN and control samples, i.e. quantifying the signature of recent mergers which might potentially be responsible for fueling/triggering the AGN. We find that (1) the vast majority (>85%) of the AGN host galaxies do not show strong distortions, and (2) there is no significant difference in the distortion fractions between active and inactive galaxies. Our findings provide the best direct evidence that, since z~1, the bulk of black hole accretion has not been triggered by major galaxy mergers, therefore arguing that the alternative mechanisms, i.e., secular processes and minor interactions, are the leading triggers for the episodes of major black hole growth. We also exclude an alternative interpretation of our results: a significant time lag between merging and the observability of the AGN phase could wash out the most significant merging signatures, explaining the lack of enhancement of strong distortions on the AGN hosts. We show that this alternative scenario is unlikely due to: (1) recent major mergers being ruled out for the majority of sources due to the high fraction of disk-hosted AGN, (2) the lack of a significant X-ray signal in merging inactive galaxies as a signature of a potential buried AGN, and (3) the low levels of soft X-ray obscuration for AGN hosted by interacting galaxies, in contrast to model predictions.

Estimating errors is a crucial part of any scientific analysis. Whenever a parameter is estimated (model-based or not), an error estimate is necessary. Any parameter estimate that is given without an error estimate is meaningless. Nevertheless, many (undergraduate or graduate) students have to teach such methods for error estimation to themselves when working scientifically for the first time. This manuscript presents an easy-to-understand overview of different methods for error estimation that are applicable to both model-based and model-independent parameter estimates. These methods are not discussed in detail, but their basics are briefly outlined and their assumptions carefully noted. In particular, the methods for error estimation discussed are grid search, varying $\chi^2$, the Fisher matrix, Monte-Carlo methods, error propagation, data resampling, and bootstrapping. Finally, a method is outlined how to propagate measurement errors through complex data-reduction pipelines.

We demonstrate that morphological observables (e.g. steepness of the radial light profile, ellipticity, asymmetry) are intertwined and cannot be measured independently of each other. We present strong arguments in favour of model-based parametrisation schemes, namely reliability assessment, disentanglement of morphological observables, and PSF modelling. Furthermore, we demonstrate that estimates of the concentration and Sersic index obtained from the Zurich Structure & Morphology catalogue are in excellent agreement with theoretical predictions. We also demonstrate that the incautious use of the concentration index for classification purposes can cause a severe loss of the discriminative information contained in a given data sample. Moreover, we show that, for poorly resolved galaxies, concentration index and M_20 suffer from strong discontinuities, i.e. similar morphologies are not necessarily mapped to neighbouring points in the parameter space. This limits the reliability of these parameters for classification purposes. Two-dimensional Sersic profiles accounting for centroid and ellipticity are identified as the currently most reliable parametrisation scheme in the regime of intermediate signal-to-noise ratios and resolutions, where asymmetries and substructures do not play an important role. We argue that basis functions provide good parametrisation schemes in the regimes of high signal-to-noise ratios and resolutions. Concerning Sersic profiles, we show that scale radii cannot be compared directly for profiles of different Sersic indices. Furthermore, we show that parameter spaces are typically highly nonlinear. This implies that significant caution is required when distance-based classificaton methods are used.

Context: The huge and still rapidly growing amount of galaxies in modern sky surveys raises the need of an automated and objective classification method. Unsupervised learning algorithms are of particular interest, since they discover classes automatically. Aims: We briefly discuss the pitfalls of oversimplified classification methods and outline an alternative approach called "clustering analysis". Methods: We categorise different classification methods according to their capabilities. Based on this categorisation, we present a probabilistic classification algorithm that automatically detects the optimal classes preferred by the data. We explore the reliability of this algorithm in systematic tests. Using a small sample of bright galaxies from the SDSS, we demonstrate the performance of this algorithm in practice. We are able to disentangle the problems of classification and parametrisation of galaxy morphologies in this case. Results: We give physical arguments that a probabilistic classification scheme is necessary. The algorithm we present produces reasonable morphological classes and object-to-class assignments without any prior assumptions. Conclusions: There are sophisticated automated classification algorithms that meet all necessary requirements, but a lot of work is still needed on the interpretation of the results.

We seek to find a shapelet-based scheme for deconvolving galaxy images from the PSF which leads to unbiased shear measurements. Based on the analytic formulation of convolution in shapelet space, we construct a procedure to recover the unconvolved shapelet coefficients under the assumption that the PSF is perfectly known. Using specific simulations, we test this approach and compare it to other published approaches. We show that convolution in shapelet space leads to a shapelet model of order $n_{max}^h = n_{max}^g + n_{max}^f$ with $n_{max}^f$ and $n_{max}^g$ being the maximum orders of the intrinsic galaxy and the PSF models, respectively. Deconvolution is hence a transformation which maps a certain number of convolved coefficients onto a generally smaller number of deconvolved coefficients. By inferring the latter number from data, we construct the maximum-likelihood solution for this transformation and obtain unbiased shear estimates with a remarkable amount of noise reduction compared to established approaches. This finding is particularly valid for complicated PSF models and low $S/N$ images, which renders our approach suitable for typical weak-lensing conditions.