Stream stars removed by tides from their progenitor satellite galaxy or globular cluster act as a group of test particles on neighboring orbits, probing the gravitational field of the Milky Way. While constraints from individual streams have been shown to be susceptible to biases, combining several streams from orbits with various distances reduces these biases. We fit a common gravitational potential to multiple stellar streams simultaneously by maximizing the clustering of the stream stars in action space. We apply this technique to members of the GD-1, Pal 5, Orphan and Helmi streams, exploiting both the individual and combined data sets. We describe the Galactic potential with a St\"ackel model, and vary up to five parameters simultaneously. We find that we can only constrain the enclosed mass, and that the strongest constraints come from the GD-1 and Pal 5 streams whose combined data set yields $M(< 20\ \mathrm{kpc}) = 3.47^{+0.95}_{-1.44} \times 10^{11} \ M_{\odot}$. When including the Orphan and Helmi stream in the data set, the mass uncertainty increases to $M(< 20\ \mathrm{kpc}) = 3.12^{+5.69}_{-1.07} \times 10^{11} \ M_{\odot}$, indicative of the hidden systematic errors in fits to individual streams. These systematics are likely due to insufficiencies in our St\"ackel model of the potential and to the limited phase space explored by the data set of 4 streams, so the larger uncertainty of the combined fit is the more robust measure of the actual uncertainty in the Milky Way potential.

We present the first spectroscopic measurements of the ATLAS and Aliqa Uma streams from the Southern Stellar Stream Spectroscopic Survey ($S^5$), in combination with the photometric data from the Dark Energy Survey and astrometric data from $Gaia$. From the coherence of spectroscopic members in radial velocity and proper motion, we find out that these two systems are physically one stream with discontinuity in morphology and density on the sky (the "kink" feature). We refer to this entire stream as ATLAS-Aliqa Uma stream, or AAU stream. We perform a comprehensive exploration of the effect of baryonic substructures and find that only an encounter with the Sagittarius dwarf $\sim 0.5$ Gyr ago can create a feature similar to the observed "kink". In addition, we also identify two gaps in the ATLAS component associated with the broadening in the stream width (the "broadening" feature). These gaps have likely been created by small mass perturbers, such as dark matter halos, as AAU stream is the most distant cold stream known with severe variations in both the stream surface density and the stream track on the sky. With the stream track, stream distance and kinematic information, we determine the orbit of the AAU stream and find that it has been affected by the Large Magellanic Cloud, resulting in a misalignment between the proper motion and stream track. Together with the Orphan-Chenab Stream, AAU is the second stream pair that has been found to be a single stream separated into two segments by external perturbation.

Large stellar surveys are revealing the chemodynamical structure of the Galaxy across a vast spatial extent. However, the many millions of low-resolution spectra observed to date are yet to be fully exploited. We employ The Cannon, a data-driven approach to estimating abundances, to obtain detailed abundances from low-resolution (R = 1800) LAMOST spectra, using the GALAH survey as our reference. We deliver five (for dwarfs) or six (for giants) estimated abundances representing five different nucleosynthetic channels, for 3.9 million stars, to a precision of 0.05 - 0.23 dex. Using wide binary pairs, we demonstrate that our abundance estimates provide chemical discriminating power beyond metallicity alone. We show the coverage of our catalogue with radial, azimuthal and dynamical abundance maps, and examine the neutron capture abundances across the disk and halo, which indicate different origins for the in-situ and accreted halo populations. LAMOST has near-complete Gaia coverage and provides an unprecedented perspective on chemistry across the Milky Way.

The Radial Velocity Experiment (RAVE) is a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A). The 6th and final data release (DR6 or FDR) is based on 518387 observations of 451783 unique stars. RAVE observations were taken between 12 April 2003 and 4 April 2013. Here we present the genesis, setup and data reduction of RAVE as well as wavelength-calibrated and flux-normalized spectra and error spectra for all observations in RAVE DR6. Furthermore, we present derived spectral classification and radial velocities for the RAVE targets, complemented by cross matches with Gaia DR2 and other relevant catalogs. A comparison between internal error estimates, variances derived from stars with more than one observing epoch and a comparison with radial velocities of Gaia DR2 reveals consistently that 68% of the objects have a velocity accuracy better than 1.4 km/s, while 95% of the objects have radial velocities better than 4.0 km/s. Stellar atmospheric parameters, abundances and distances are presented in subsequent publication. The data can be accessed via the RAVE Web (http://rave-survey.org) or the Vizier database.

Lithium depletion and enrichment in the cosmos is not yet well understood. To help tighten constraints on stellar and Galactic evolution models, we present the largest high-resolution analysis of Li abundances A(Li) to date, with results for over 100 000 GALAH field stars spanning effective temperatures $5900\,\mathrm{K} \lesssim \rm{T_{eff}} \lesssim7000\,\mathrm{K}$ and metallicities $-3 \lesssim \rm[Fe/H] \lesssim +0.5$. We separated these stars into two groups, on the warm and cool side of the so-called Li-dip, a localised region of the Kiel diagram wherein lithium is severely depleted. We discovered that stars in these two groups show similar trends in the A(Li)-[Fe/H] plane, but with a roughly constant offset in A(Li) of 0.4 dex, the warm group having higher Li abundances. At $\rm[Fe/H]\gtrsim-0.5$, a significant increasing in Li abundance with increasing metallicity is evident in both groups, signalling the onset of significant Galactic production. At lower metallicity, stars in the cool group sit on the Spite plateau, showing a reduced lithium of around 0.4 dex relative to the primordial value predicted from Big Bang nucleosynthesis (BBN). However, stars in the warm group between [Fe/H] = -1.0 and -0.5, form an elevated plateau that is largely consistent with the BBN prediction. This may indicate that these stars in fact preserve the primordial Li produced in the early Universe.

We present part 2 of the 6th and final Data Release (DR6 or FDR) of the Radial Velocity Experiment (RAVE), a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A) and span the complete time frame from the start of RAVE observations on 12 April 2003 to their completion on 4 April 2013. In the second of two publications, we present the data products derived from 518387 observations of 451783 unique stars using a suite of advanced reduction pipelines focussing on stellar atmospheric parameters, in particular purely spectroscopically derived stellar atmospheric parameters (Teff, log(g), and the overall metallicity), enhanced stellar atmospheric parameters inferred via a Bayesian pipeline using Gaia DR2 astrometric priors, and asteroseismically calibrated stellar atmospheric parameters for giant stars based on asteroseismic observations for 699 K2 stars. In addition, we provide abundances of the elements Fe, Al, and Ni, as well as an overall [alpha/Fe] ratio obtained using a new pipeline based on the GAUGUIN optimization method that is able to deal with variable signal-to-noise ratios. The RAVE DR6 catalogs are cross matched with relevant astrometric and photometric catalogs, and are complemented by orbital parameters and effective temperatures based on the infrared flux method. The data can be accessed via the RAVE Web site (http://rave-survey.org) or the Vizier database.

We present multi-scale and multi-wavelength observations of the Galactic HII region S305, which is excited by massive O8.5V and O9.5V stars. Infrared images reveal an extended sphere-like shell (extension ~7.5 pc; at T_d = 17.5-27 K) enclosing the S305 HII region (size ~5.5 pc; age ~1.7 Myr). The extended structure observed in the Herschel temperature map indicates that the molecular environment of S305 is heated by the massive O-type stars. Regularly spaced molecular condensations and dust clumps are investigated toward the edges of the infrared shell, where the PAH and H_2 emission is also observed. The molecular line data show a signature of an expanding shell of molecular gas in S305. GMRT 610 and 1280 MHz continuum maps reveal overdensities of the ionized emission distributed around two O-type stars, which are surrounded by the horseshoe envelope (extension ~2.3 pc). A molecular gas deficient region/cavity is identified toward the center of the horseshoe envelope, which is well traced with PAH, H_2, molecular, and dust emission. The edges of the infrared shell are found to be located in the front of the horseshoe envelope. All these outcomes provide the observational evidence of the feedback of O-type stars in S305. Moreover, non-thermal radio emission is detected in S305 with an average spectral index alpha ~-0.45. The variations in alpha, ranging from -1.1 to 1.3, are explained due to soft synchrotron emission and either optically-thicker thermal emission at high frequencies or a suppression of the low-frequency emission by the Razin-Tsytovich effect.

We detail the sensitivity of the liquid xenon (LXe) DARWIN observatory to solar neutrinos via elastic electron scattering. We find that DARWIN will have the potential to measure the fluxes of five solar neutrino components: $pp$, $^7$Be, $^{13}$N, $^{15}$O and $pep$. The precision of the $^{13}$N, $^{15}$O and $pep$ components is hindered by the double-beta decay of $^{136}$Xe and, thus, would benefit from a depleted target. A high-statistics observation of $pp$ neutrinos would allow us to infer the values of the weak mixing angle, $\sin^2\theta_w$, and the electron-type neutrino survival probability, $P_e$, in the electron recoil energy region from a few keV up to 200 keV for the first time, with relative precision of 5% and 4%, respectively, at an exposure of 300 ty. An observation of $pp$ and $^7$Be neutrinos would constrain the neutrino-inferred solar luminosity down to 0.2%. A combination of all flux measurements would distinguish between the high (GS98) and low metallicity (AGS09) solar models with 2.1-2.5$\sigma$ significance, independent of external measurements from other experiments or a measurement of $^8$B neutrinos through coherent elastic neutrino-nucleus scattering in DARWIN. Finally, we demonstrate that with a depleted target DARWIN may be sensitive to the neutrino capture process of $^{131}$Xe.

We investigate the properties of 1306 red giant stars with high photospheric abundances of lithium observed by the GALAH, K2-HERMES and TESS-HERMES surveys, and discuss them in the context of proposed mechanisms for lithium enrichment in giant stars. We confirm that Li-rich giants are rare, making up only 1.1 per cent of our giant star sample. We use GALAH+ DR3 stellar parameters and a Bayesian isochrone analysis to divide the sample into first-ascent red giant branch and red clump stars, and confirm these classifications using asteroseismic data from K2. We find that red clump stars are more than three times as likely to be lithium-rich as red giant branch stars, and the occurrence rate of lithium richness with metallicity is quite different between the two populations. The probability for a rapidly rotating giant to be lithium-rich is distinctly higher as compared to a non rapidly rotating giant, independent of its evolutionary state. There is also a clear correspondence on the red clump between rapid rotation and extreme levels of lithium enrichment. There are almost no stars on the secondary red clump with lithium abundances above the primordial level, indicating a mass dependence in the red clump lithium enrichment mechanism(s). The complex distribution of lithium-rich giants across evolutionary phase, metallicity, rotation rate and mass implies multiple independent mechanisms for producing lithium enrichment.

Accurate and precise radius estimates of transiting exoplanets are critical for understanding their compositions and formation mechanisms. To know the planet, we must know the host star in as much detail as possible. We present complete results for planet-candidate hosts from the K2-HERMES survey, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R$\sim$28,000 spectra for more than 30,000 K2 stars. We present complete host-star parameters and planet-candidate radii for 224 K2 candidate planets from C1-C13. Our results cast severe doubt on 30 K2 candidates, as we derive unphysically large radii, larger than 2 $R_{Jup}$. This work highlights the importance of obtaining accurate, precise, and self-consistent stellar parameters for ongoing large planet search programs - something that will only become more important in the coming years, as TESS begins to deliver its own harvest of exoplanets.

The $([\alpha/{\rm Fe}],[{\rm Fe/H}])$ distribution of Milky Way stars shows at least two distinct sequences, which have traditionally been associated with the thin and thick disc components. The abundance distribution varies systematically with location $R$ and $|z|$ across the Galaxy. Using an analytical chemodynamical model that includes the effects of radial migration and kinematic heating, we show that it is possible to reproduce the observed abundance distribution at different locations. Unlike some earlier models, our scheme does not require a distinct thick disc component emerging from a separate evolutionary path. The proposed model has a continuous star formation history and a continuous age velocity dispersion relation. Moreover, [$\alpha$/Fe] is constant for stellar ages less than 8 Gyr, but increases sharply for older stars over a time scale of 1.5 Gyr. The gap between the two sequences is due to this sharp transition. We show that the high-[$\alpha$/Fe] sequence at the low metallicity end is simply a pile-up of old stars. At the high metallicity end, we find a sequence of stars having different ages with a similar birth radius that originated in the inner disc. Our model successfully explains the uniformity of the locus of the high-[$\alpha$/Fe] sequence across different locations. The low-[$\alpha$/Fe] sequence contains stars with different birth radii that owes its existence to radial migration. For the low-[$\alpha$/Fe] sequence, angular momentum is anti-correlated with [Fe/H], and the model can reproduce this trend at different Galactic locations. If radial migration is not included, the model fails to generate the double sequence and instead shows only a single sequence. Our simple scheme has major advantages over earlier chemodynamical models, as we show.

Binary stellar systems form a large fraction of the Galaxy's stars. They are useful as laboratories for studying the physical processes taking place within stars, and must be correctly taken into account when observations of stars are used to study the structure and evolution of the Galaxy. We present a sample of 12760 well-characterised double-lined spectroscopic binaries that are appropriate for statistical studies of the binary populations. They were detected as SB2s using a t-distributed stochastic neighbour embedding (t-SNE) classification and a cross-correlation analysis of GALAH spectra. This sample consists mostly of dwarfs, with a significant fraction of evolved stars and several dozen members of the giant branch. To compute parameters of the primary and secondary star ($T_{\rm eff[1,2]}$, $\log g_{[1,2]}$, [Fe/H], $V_{r[1,2]}$, $v_{\rm mic[1,2]}$, $v_{\rm broad[1,2]}$, $R_{[1,2]}$, and $E(B-V)$), we used a Bayesian approach that includes a parallax prior from Gaia DR2, spectra from GALAH, and apparent magnitudes from APASS, Gaia DR2, 2MASS, and WISE. The derived stellar properties and their distributions show trends that are expected for a population of close binaries (a $<$ 10 AU) with mass ratios $0.5 \leq q \leq 1$. The derived metallicity of these binary stars is statistically lower than that of single dwarf stars from the same magnitude-limited sample.

Deep and wide-field optical photometric observations along with multiwavelength archival datasets have been employed to study the physical properties of the cluster NGC 6910. The study also examines the impact of massive stars to their environment. The age, distance and reddening of the cluster are estimated to be $\sim$4.5 Myr, $1.72\pm0.08$ kpc, and $ E(B-V)_{min}= 0.95$ mag, respectively. The mass function slope ($\Gamma = -0.74\pm0.15$ in the cluster region is found to be flatter than the Salpeter value (-1.35), indicating the presence of excess number of massive stars. The cluster also shows mass segregation towards the central region due to their formation processes. The distribution of warm dust emission is investigated towards the central region of the cluster, showing the signature of the impact of massive stars within the cluster region. Radio continuum clumps powered by massive B-type stars (age range $\sim$ 0.07-0.12 Myr) are traced, which are located away from the center of the stellar cluster NGC 6910 (age $\sim$ 4.5 Myr). Based on the values of different pressure components exerted by massive stars, the photoionized gas associated with the cluster is found to be the dominant feedback mechanism in the cluster. Overall, the massive stars in the cluster might have triggered the birth of young massive B-type stars in the cluster. This argument is supported with evidence of the observed age gradient between the cluster and the powering sources of the radio clumps.

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 present a detailed investigation of the Rastall gravity extension of the standard $\Lambda$CDM model. We review the model for two simultaneous modifications of different nature in the Friedmann equation due to the Rastall gravity: the new contributions of the material (actual) sources (considered as effective source) and the altered evolution of the material sources. We discuss the role/behavior of these modifications with regard to some low redshift tensions, including the $H_0$ tension, prevailing within the standard $\Lambda$CDM. We constrain the model at the level of linear perturbations, and obtain the first constraints through a robust and accurate analysis using the latest full Planck CMB data, with and without including BAO data. We find that the Rastall parameter $\epsilon$ (null for general relativity) is consistent with zero at 68% CL (with a tendency towards positive values, $-0.0001 < \epsilon < 0.0007$ (CMB+BAO) at 68% CL), which in turn implies no significant statistical evidence for deviation from general relativity, and also a precision of $\mathcal{O}(10^{-4})$ for the coefficient $-1/2$ of the term $g_{\mu\nu}R$ in the Einstein field equations of general relativity (guaranteeing the local energy-momentum conservation). We explore the consequences led by the Rastall gravity on the cosmological parameters in the light of the observational analyses. It turns out that the effective source dynamically screens the usual vacuum energy at high redshift, but this mechanism barely works due to the opposition by the altered evolution of CDM. Consequently, two simultaneous modifications of different nature in the Friedmann equation act against each other, and do not help to considerably relax the so called low redshift tensions. Our results may offer a guide for the research community that studies the Rastall gravity in various aspects of gravitation and cosmology.

We explore the fundamental relations governing the radial and vertical velocity dispersions of stars in the Milky Way, from combined studies of complementary surveys including GALAH, LAMOST, APOGEE, the NASA $Kepler$ and K2 missions, and $Gaia$ DR2. We find that different stellar samples, even though they target different tracer populations and employ a variety of age estimation techniques, follow the same set of fundamental relations. We provide the clearest evidence to date that, in addition to the well-known dependence on stellar age, the velocity dispersions of stars depend on orbital angular momentum $L_z$, metallicity and height above the plane $|z|$, and are well described by a multiplicatively separable functional form. The dispersions have a power-law dependence on age with exponents of 0.441$\pm 0.007$ and 0.251$\pm 0.006$ for $\sigma_z$ and $\sigma_R$ respectively, and the power law is valid even for the oldest stars. For the solar neighborhood stars, the apparent break in the power law for older stars, as seen in previous studies, is due to the anti-correlation of $L_z$ with age. The dispersions decrease with increasing $L_z$ until we reach the Sun's orbital angular momentum, after which $\sigma_z$ increases (implying flaring in the outer disc) while $\sigma_R$ flattens. The dispersions increase with decreasing metallicity, suggesting that the dispersions increase with birth radius. The dispersions also increase linearly with $|z|$. The same set of relations that work in the solar neighborhood also work for stars between $3<R/{\rm kpc}<20$. Finally, the high-[$\alpha$/Fe] stars follow the same relations as the low-[$\alpha$/Fe] stars.

With Gaia Data Release 2, the astronomical community is entering a new era of multidimensional surveys of the Milky Way. This new phase-space view of our Galaxy demands new tools for comparing observations to simulations of Milky-Way-mass galaxies in a cosmological context, to test the physics of both dark matter and galaxy formation. We present ananke, a framework for generating synthetic phase-space surveys from high-resolution baryonic simulations, and use it to generate a suite of synthetic surveys resembling Gaia DR2 in data structure, magnitude limits, and observational errors. We use three cosmological simulations of Milky-Way-mass galaxies from the Latte suite of the Feedback In Realistic Environments (FIRE) project, which feature self-consistent clustering of star formation in dense molecular clouds and thin stellar/gaseous disks in live cosmological halos with satellite dwarf galaxies and stellar halos. We select three solar viewpoints from each simulation to generate nine synthetic Gaia-like surveys. We sample synthetic stars by assuming each star particle (of mass 7070 $M_{\odot}$) represents a single stellar population. At each viewpoint, we compute dust extinction from the simulated gas metallicity distribution and apply a simple error model to produce a synthetic Gaia-like survey that includes both observational properties and a pointer to the generating star particle. We provide the complete simulation snapshot at $z = 0$ for each simulated galaxy. We describe data access points, the data model, and plans for future upgrades. These synthetic surveys provide a tool for the scientific community to test analysis methods and interpret Gaia data.

Chirality transfer between fermions and gauge fields plays a crucial role for understanding the dynamics of anomalous transport phenomena such as the Chiral Magnetic Effect. In this proceeding we present a first principles study of these processes based on classical-statistical real-time lattice simulations of strongly coupled QED $(e^2N_f=64)$. Our simulations demonstrate that a chirality imbalance in the fermion sector triggers chiral plasma instabilities in the gauge field sector, which ultimately lead to the generation of long range helical magnetic fields via a self-similar turbulent cascade of the magnetic helicity.

We present multi-epoch deep ($\sim$20 mag) $I_{c}$~band photometric monitoring of the Sh 2-170 star-forming region to understand the variability properties of pre-main-sequence (PMS) stars. We report identification of 47 periodic and 24 non-periodic variable stars with periods and amplitudes ranging from $\sim$4 hrs to 18 days and from $\sim$0.1 to 2.0 mag, respectively. We have further classified 49 variables as PMS stars (17 Class\,{\sc ii} and 32 Class\,{\sc iii}) and 17 as main-sequence (MS)/field star variables. A larger fraction of MS/field variables (88\%) show periodic variability as compared to the PMS variables (59\%). The ages and masses of the PMS variable stars are found to be comparable with those of T-Tauri stars. Their variability amplitudes show an increasing trend with the near-IR/mid-IR excess. The period distribution of the PMS variables shows two peaks, one near $\sim$1.5 days and the other near $\sim$4.5 days. It is found that the younger stars with thicker discs and envelopes seem to rotate slower than their older counterparts. These properties of the PMS variables support the disc-locking mechanism. Both the period and amplitude of PMS stars show decrease with increasing mass probably due to the effective dispersal of circumstellar discs in massive stars. Our results favour the notion that cool spots on weak line T-Tauri stars are responsible for most of their variations, while hot spots on classical T-Tauri stars resulting from variable mass accretion from an inner disc contribute to their larger amplitudes and irregular behaviours.

Using our deep optical and near-infrared photometry along with multiwavelength archival data, we here present a detailed study of the Galactic H II region Sh 2-305, to understand the star/star-cluster formation. On the basis of excess infra-red emission, we have identified 116 young stellar objects (YSOs) within a field of view of ~ 18.5 arcminute x 18.5 arcminute, around Sh 2-305. The average age, mass and extinction (A_V) for this sample of YSOs are 1.8 Myr, 2.9 solar mass and 7.1 mag, respectively. The density distribution of stellar sources along with minimal spanning tree calculations on the location of YSOs reveals at least three stellar sub-clusterings in Sh 2-305. One cluster is seen toward the center (i.e., Mayer 3), while the other two are distributed toward the north and south directions. Two massive O-type stars (VM2 and VM4; ages ~ 5 Myr) are located at the center of the Sh 2-305 H II region. The analysis of the infrared and radio maps traces the photon dominant regions (PDRs) in the Sh 2-305. Association of younger generation of stars with the PDRs is also investigated in the Sh 2-305. This result suggests that these two massive stars might have influenced the star formation history in the Sh 2-305. This argument is also supported with the calculation of various pressures driven by massive stars, slope of mass function/K-band luminosity function, star formation efficiency, fraction of Class I sources, and mass of the dense gas toward the sub-clusterings in the Sh 2-305.

New observations in the $VI$ bands along with archival data from the 2MASS and $WISE$ surveys have been used to generate a catalog of young stellar objects (YSOs) covering an area of about $6^\circ\times6^\circ$ in the Auriga region centered at $l\sim173^\circ$ and $b\sim1^\circ.5$. The nature of the identified YSOs and their spatial distribution are used to study the star formation in the region. The distribution of YSOs along with that of the ionized and molecular gas reveals two ring-like structures stretching over an area of a few degrees each in extent. We name these structures as Auriga Bubbles 1 and 2. The center of the Bubbles appears to be above the Galactic mid-plane. The majority of Class\,{\sc i} YSOs are associated with the Bubbles, whereas the relatively older population, i.e., Class\,{\sc ii} objects are rather randomly distributed. Using the minimum spanning tree analysis, we found 26 probable sub-clusters having 5 or more members. The sub-clusters are between $\sim$0.5 pc - $\sim$ 3 pc in size and are somewhat elongated. The star formation efficiency in most of the sub-cluster region varies between 5$\%$ - 20$\%$ indicating that the sub-clusters could be bound regions. The radii of these sub-clusters also support it.

We present the Bayesian Asteroseismology data Modeling (BAM) pipeline, an automated asteroseismology pipeline that returns global oscillation parameters and granulation parameters from the analysis of photometric time-series. BAM also determines if a star is likely to be a solar-like oscillator. We have designed BAM to specially process ${\it K2}$ light curves, which suffer from unique noise signatures that can confuse asteroseismic analysis, though it may be used on any photometric time series --- including those from ${\it Kepler}$ and ${\it TESS}$. We demonstrate the BAM oscillation parameters are consistent within $\sim 1.53\%\ (\mathrm{random}) \pm 0.2\%\ (\mathrm{systematic})$ and $1.51\%\ (\mathrm{random}) \pm 0.6\%\ (\mathrm{systematic})$ for $\nu_{\mathrm{max}}$ and $\Delta \nu$ with benchmark results for typical ${\it K2}$ red giant stars in the ${\it K2}$ Galactic Archaeology Program's (GAP) Campaign 1 sample. Application of BAM to $13016$ ${\it K2}$ Campaign 1 targets not in the GAP sample yields $104$ red giant solar-like oscillators. Based on the number of serendipitous giants we find, we estimate an upper limit on the average purity in dwarf selection among C1 proposals is $\approx 99\%$, which could be lower when considering incompleteness in BAM detection efficiency, and proper motion cuts specific to C1 Guest Observer proposals.

We present isochrone ages and initial bulk metallicities ($\rm [Fe/H]_{bulk}$, by accounting for diffusion) of 163,722 stars from the GALAH Data Release 2, mainly composed of main sequence turn-off stars and subgiants ($\rm 7000 K>T_{eff}>4000 K$ and $\rm log g>3$ dex). The local age-metallicity relationship (AMR) is nearly flat but with significant scatter at all ages; the scatter is even higher when considering the observed surface abundances. After correcting for selection effects, the AMR appear to have intrinsic structures indicative of two star formation events, which we speculate are connected to the thin and thick disks in the solar neighborhood. We also present abundance ratio trends for 16 elements as a function of age, across different $\rm [Fe/H]_{bulk}$ bins. In general, we find the trends in terms of [X/Fe] vs age from our far larger sample to be compatible with studies based on small ($\sim$ 100 stars) samples of solar twins but we now extend it to both sub- and super-solar metallicities. The $\alpha$-elements show differing behaviour: the hydrostatic $\alpha$-elements O and Mg show a steady decline with time for all metallicities while the explosive $\alpha$-elements Si, Ca and Ti are nearly constant during the thin disk epoch (ages $\lessapprox $ 12 Gyr). The s-process elements Y and Ba show increasing [X/Fe] with time while the r-process element Eu have the opposite trend, thus favouring a primary production from sources with a short time-delay such as core-collapse supernovae over long-delay events such as neutron star mergers.

We introduce the Southern Stellar Stream Spectroscopy Survey (${S}^5$), an on-going program to map the kinematics and chemistry of stellar streams in the Southern Hemisphere. The initial focus of ${S}^5$ has been spectroscopic observations of recently identified streams within the footprint of the Dark Energy Survey (DES), with the eventual goal of surveying streams across the entire southern sky. Stellar streams are composed of material that has been tidally striped from dwarf galaxies and globular clusters and hence are excellent dynamical probes of the gravitational potential of the Milky Way, as well as providing a detailed snapshot of its accretion history. Observing with the 3.9-m Anglo-Australian Telescope's 2-degree-Field fibre positioner and AAOmega spectrograph, and combining the precise photometry of DES DR1 with the superb proper motions from $Gaia$ DR2, allows us to conduct an efficient spectroscopic survey to map these stellar streams. So far ${S}^5$ has mapped 9 DES streams and 3 streams outside of DES; the former are the first spectroscopic observations of these recently discovered streams. In addition to the stream survey, we use spare fibres to undertake a Milky Way halo survey and a low-redshift galaxy survey. This paper presents an overview of the ${S}^5$ program, describing the scientific motivation for the survey, target selection, observation strategy, data reduction and survey validation. Finally, we describe early science results on stellar streams and Milky Way halo stars drawn from the survey. Updates on ${S}^5$, including future public data release, can be found at \url{http://s5collab.github.io}.

Using kinematics from Gaia and the large elemental abundance space of the second data release of the GALAH survey, we identify two new members of the Fimbulthul stellar stream, and chemically tag them to massive, multi-metallic globular cluster $\omega$ Centauri. Recent analysis of the second data release of Gaia had revealed the Fimbulthul stellar stream in the halo of the Milky Way. It had been proposed that the stream is associated with the $\omega$ Centauri, but this proposition relied exclusively upon the kinematics and metallicities of the stars to make the association. In this work, we find our two new members of the stream to be metal-poor stars that are enhanced in sodium and aluminium, typical of second population globular cluster stars, but not otherwise seen in field stars. Furthermore, the stars share the s-process abundance pattern seen in $\omega$ Centauri, which is rare in field stars. Apart from one star within 1.5 deg of $\omega$ Centauri, we find no other stars observed by GALAH spatially near $\omega$ Centauri or the Fimbulthul stream that could be kinematically and chemically linked to the cluster. Chemically tagging stars in the Fimbulthul stream to $\omega$ Centauri confirms the earlier work, and further links this tidal feature in the Milky Way halo to $\omega$ Centauri.

Chemical tagging seeks to identify unique star formation sites from present-day stellar abundances. Previous techniques have treated each abundance dimension as being statistically independent, despite theoretical expectations that many elements can be produced by more than one nucleosynthetic process. In this work we introduce a data-driven model of nucleosynthesis where a set of latent factors (e.g., nucleosynthetic yields) contribute to all stars with different scores, and clustering (e.g., chemical tagging) is modelled by a mixture of multivariate Gaussians in a lower-dimensional latent space. We use an exact method to simultaneously estimate the factor scores for each star, the partial assignment of each star to each cluster, and the latent factors common to all stars, even in the presence of missing data entries. We use an information-theoretic Bayesian principle to estimate the number of latent factors and clusters. Using the second Galah data release we find that six latent factors are preferred to explain N = 2,566 stars with 17 chemical abundances. We identify the rapid- and slow-neutron capture processes, as well as latent factors consistent with Fe-peak and \alpha-element production, and another where K and Zn dominate. When we consider N ~ 160,000 stars with missing abundances we find another 7 factors, as well as 16 components in latent space. Despite these components showing separation in chemistry that is explained through different yield contributions, none show significant structure in their positions or motions. We argue that more data, and joint priors on cluster membership that are constrained by dynamical models, are necessary to realise chemical tagging at a galactic-scale. We release software that allows for model parameters to be optimised in seconds given a fixed number of latent factors, components, and $10^7$ abundance measurements.

We present a first principles study of chiral plasma instabilities and the onset of chiral turbulence in QED plasmas far from equilibrium. By performing classical-statistical lattice simulations of the microscopic theory, we show that the generation of strong helical magnetic fields from a helicity imbalance in the fermion sector proceeds via three distinct phases. During the initial linear instability regime the helicity imbalance of the fermion sector causes an exponential growth(damping) of magnetic field modes with right(left) handed polarization, for which we extract the characteristic growth (damping) rates. Secondary growth of unstable modes accelerates the helicity transfer from fermions to gauge fields and ultimately leads to the emergence of a self-similar scaling regime characteristic of decaying turbulence, where magnetic helicity is efficiently transferred to macroscopic length scales. Within this turbulent regime the evolution of magnetic helicity spectrum can be described by an infrared power-spectrum with spectral exponent $\kappa$ and dynamical scaling exponents $\alpha,\beta$, which we determine from our simulations.

The Orion complex is arguably the most studied star-forming region in the Galaxy. While stars are still being born in the Orion nebula, the oldest part was believed to be no more than 13 Myr old. In order to study the full hierarchy of star formation across the Orion complex, we perform a clustering analysis of the Ori OB1a region using new stellar surveys and derive robust ages for each identified stellar aggregate. We use Gaia DR2 parameters supplemented with radial velocities from the GALAH and APOGEE surveys to perform clustering of the Ori OB1a association. Five overdensities are resolved in a six-dimensional parameter space (positions, distance, proper motions, and radial velocity). Most correspond to previously known structures (ASCC 16, 25 Orionis, ASCC 20, ASCC 21). We use Gaia DR2, Pan-STARRS1 and 2MASS photometry to fit isochrones to the colour-magnitude diagrams of the identified clusters. The ages of the clusters can thus be measured with ~10% precision. While four of the clusters have ages between 11 and 13 Myr, the ASCC 20 cluster stands out at an age of 21 $\pm$ 3 Myr. This is significantly greater than the age of any previously known component of the Orion complex. To some degree, all clusters overlap in at least one of the six phase-space dimensions. We argue that the formation history of the Orion complex, and its relation to the Gould belt, must be reconsidered. A significant challenge in reconstructing the history of the Ori OB1a association is to understand the impact of the newly discovered 21 Myr old population on the younger parts of the complex, including their formation.

Gaia DR2 has revealed new small-scale and large-scale patterns in the phase-space distribution of stars in the Milky Way. In cylindrical Galactic coordinates $(R,\phi,z)$, ridge-like structures can be seen in the \vphiR{} plane and asymmetric arch-like structures in the \vphivR{} plane. We show that the ridges are also clearly present when the third dimension of the \vphiR{} plane is represented by $\langle z \rangle$, $\langle V_z \rangle$, $\langle V_R \rangle$, $\langle$[Fe/H]$\rangle$ and $\langle[\alpha/{\rm Fe}]\rangle$. The maps suggest that stars along the ridges lie preferentially close to the Galactic midplane ($|z|<0.2$ kpc), and have metallicity and $\alpha$ elemental abundance similar to that of the Sun. We show that phase mixing of disrupting spiral arms can generate both the ridges and the arches. It also generates discrete groupings in orbital energy $-$ the ridges and arches are simply surfaces of constant energy. We identify 8 distinct ridges in the \gaia{} data: six of them have constant energy while two have constant angular momentum. Given that the signature is strongest for stars close to the plane, the presence of ridges in $\langle z \rangle$ and $\langle V_z \rangle$ suggests a coupling between planar and vertical directions. We demonstrate, using N-body simulations that such coupling can be generated both in isolated discs and in discs perturbed by an orbiting satellite like the Sagittarius dwarf galaxy.

This paper describes the design and implementation of Stingray, a library in Python built to perform time series analysis and related tasks on astronomical light curves. Its core functionality comprises a range of Fourier analysis techniques commonly used in spectral-timing analysis, as well as extensions for analyzing pulsar data, simulating data sets, and statistical modeling. Its modular build allows for easy extensions and incorporation of its methods into data analysis workflows and pipelines. We aim for the library to be a platform for the implementation of future spectral-timing techniques. Here, we describe the overall vision and framework, core functionality, extensions, and connections to high-level command-line and graphical interfaces. The code is well-tested, with a test coverage of currently 95%, and is accompanied by extensive API documentation and a set of step-by-step tutorials.

We consider the simplest anisotropic generalization, as a correction, to the standard $\Lambda$CDM model, by replacing the spatially flat Robertson-Walker metric by the Bianchi type-I metric, which brings in a new term $\Omega_{\sigma 0}a^{-6}$ (mimicking the stiff fluid) in the average expansion rate $H(a)$ of the Universe. From Hubble and Pantheon data, relevant to the late Universe ($z\lesssim 2.4$), we obtain the constraint $\Omega_{\sigma0}\lesssim10^{-3}$, in line with the model-independent constraints. When the baryonic acoustic oscillations and cosmic microwave background (CMB) data are included, the constraint improves by 12 orders of magnitude, i.e., $\Omega_{\sigma0}\lesssim10^{-15}$. We find that this constraint could alter neither the matter-radiation equality redshift nor the peak of the matter perturbations. Demanding that the expansion anisotropy has no significant effect on the standard big bang nucleosynthesis (BBN), we find the constraint $\Omega_{\sigma0}\lesssim10^{-23}$. We show explicitly that the constraint from BBN renders the expansion anisotropy irrelevant to make a significant change in the CMB quadrupole temperature, whereas the constraint from the cosmological data in our model provides the temperature change up to $\sim11\, \rm mK$, though it is much beyond the CMB quadrupole temperature.

We describe the catalogs assembled and the algorithms used to populate the revised TESS Input Catalog (TIC), based on the incorporation of the Gaia second data release. We also describe a revised ranking system for prioritizing stars for 2-minute cadence observations, and assemble a revised Candidate Target List (CTL) using that ranking. The TIC is available on the Mikulski Archive for Space Telescopes (MAST) server, and an enhanced CTL is available through the Filtergraph data visualization portal system at the URL http://filtergraph.vanderbilt.edu/tess_ctl.

We present results of a search for variable stars in a region of the globular cluster NGC 4147 based on photometric observations with 4k x 4k CCD imager mounted at the axial port of the recently installed 3.6 m Devasthal optical telescope at Aryabhatta Research Institute of Observational Sciences, Nainital, India. We performed time series photometry of NGC 4147 in V and R bands, and identified 42 periodic variables in the region of NGC 4147, 28 of which have been detected for the first time. Seventeen variable stars are located within the half light radius $\lesssim$ 0.48 arcmin, of which 10 stars are newly identified variables. Two of 10 variables are located within the core radius $\lesssim$ 0.09 arcmin. Based on the location in the $V/(V-R)$ colour magnitude diagram and variability characteristics, 7, 8, 5 and 1 newly identified probable member variables are classified as RRc, EA/E, EW and SX Phe, respectively. The metallicity of NGC 4147 estimated from light curves of RRab and RRc stars with the help of Fourier decomposition is found to be characteristics of Oosterhoff II. The distance derived using light curves of RRab stars is consistent with that obtained from the observed $V/(V-R)$ colour-magnitude diagram.

The latest Gaia data release enables us to accurately identify stars that are more luminous than would be expected on the basis of their spectral type and distance. During an investigation of the 329 best Solar twin candidates uncovered among the spectra acquired by the GALAH survey, we identified 64 such over-luminous stars. In order to investigate their exact composition, we developed a data-driven methodology that can generate a synthetic photometric signature and spectrum of a single star. By combining multiple such synthetic stars into an unresolved binary or triple system and comparing the results to the actual photometric and spectroscopic observations, we uncovered 6 definitive triple stellar system candidates and an additional 14 potential candidates whose combined spectrum mimics the Solar spectrum. Considering the volume correction factor for a magnitude limited survey, the fraction of probable unresolved triple stars with long orbital periods is ~2 %. Possible orbital configurations of the candidates were investigated using the selection and observational limits. To validate the discovered multiplicity fraction, the same procedure was used to evaluate the multiplicity fraction of other stellar types.

Asteroseismology is a promising tool to study Galactic structure and evolution because it can probe the ages of stars. Earlier attempts comparing seismic data from the {\it Kepler} satellite with predictions from Galaxy models found that the models predicted more low-mass stars compared to the observed distribution of masses. It was unclear if the mismatch was due to inaccuracies in the Galactic models, or the unknown aspects of the selection function of the stars. Using new data from the K2 mission, which has a well-defined selection function, we find that an old metal-poor thick disc, as used in previous Galactic models, is incompatible with the asteroseismic information. We show that spectroscopic measurements of [Fe/H] and [$\alpha$/Fe] elemental abundances from the GALAH survey indicate a mean metallicity of $\log (Z/Z_{\odot})=-0.16$ for the thick disc. Here $Z$ is the effective solar-scaled metallicity, which is a function of [Fe/H] and [$\alpha$/Fe]. With the revised disc metallicities, for the first time, the theoretically predicted distribution of seismic masses show excellent agreement with the observed distribution of masses. This provides an indirect verification of the asteroseismic mass scaling relation is good to within five percent. Using an importance-sampling framework that takes the selection function into account, we fit a population synthesis model of the Galaxy to the observed seismic and spectroscopic data. Assuming the asteroseismic scaling relations are correct, we estimate the mean age of the thick disc to be about 10 Gyr, in agreement with the traditional idea of an old $\alpha$-enhanced thick disc.

The recently published Kepler mission Data Release 25 (DR25) reported on ~197,000 targets observed during the mission. Despite this, no wide search for red giants showing solar-like oscillations have been made across all stars observed in Kepler's long-cadence mode. In this work, we perform this task using custom apertures on the Kepler pixel files and detect oscillations in 21,914 stars, representing the largest sample of solar-like oscillating stars to date. We measure their frequency at maximum power, $\nu_{\mathrm{max}}$, down to $\nu_{\mathrm{max}}\simeq4\mu$Hz and obtain $\log(g)$ estimates with a typical uncertainty below 0.05 dex, which is superior to typical measurements from spectroscopy. Additionally, the $\nu_{\mathrm{max}}$ distribution of our detections show good agreement with results from a simulated model of the Milky Way, with a ratio of observed to predicted stars of 0.992 for stars with 10$\mu$Hz $ <\nu_{\mathrm{max}}<270\mu$Hz. Among our red giant detections, we find 909 to be dwarf/subgiant stars whose flux signal is polluted by a neighbouring giant as a result of using larger photometric apertures than those used by the NASA Kepler Science Processing Pipeline. We further find that only 293 of the polluting giants are known Kepler targets. The remainder comprises over 600 newly identified oscillating red giants, with many expected to belong to the galactic halo, serendipitously falling within the Kepler pixel files of targeted stars.

We present the results based on the optical $R$-band polarization observations of 280 stars distributed towards the dark globule LDN\,1225. {\it Gaia} data release 2 parallaxes along with the polarization data of $\sim$200 stars have been used to (a) constrain the distance of LDN\,1225 as 830$\pm$83~pc, (b) determine the contribution of interstellar polarization (ISP), and (c) characterize the dust properties and delineate the magnetic field (B-field) morphology of LDN\,1225. We find that B-fields are more organized and exhibit a small dispersion of 12$\degr$. Using the $^{12}$CO molecular line data from the Purple Mountain Observatory (PMO), along with the column density, dispersion in B-fields, we estimate B-field strength to be $\sim$56\,$\pm$\,10\,$\mu$G, magnetic to turbulence pressure to be $\sim$3\,$\pm$\,2, and the mass-to-magnetic flux ratio (in units of critical value) to be~$<$\,1. These results indicate the dominant role of B-fields in comparison to turbulence and gravity in rendering the cloud support. B-fields are aligned parallel to the low-density parts (traced by $^{12}$CO map) of the cloud, in contrast they are neither parallel nor perpendicular to the high-density core structures (traced by $^{13}$CO and C$^{18}$O maps). LDN\,1225 hosts two 70\,$\mu$m sources which seem to be of low-mass Class 0 sources. The total-to-selective extinction derived using optical and near-infrared photometric data is found to be anomalous ($R_{V}$~$=$~3.4), suggesting dust grain growth in LDN\,1225. Polarization efficiency of dust grains follows a power-law index of $-$0.7 inferring that optical polarimetry traces B-fields in the outer parts of the cloud.

We present the chemodynamic structure of the solar neighbourhood using 62 814 stars within a 500 pc sphere of the Sun observed by GALAH and with astrometric parameters from Gaia DR2. We measure the velocity dispersion for all three components (vertical, radial, and tangential) and find that it varies smoothly with [Fe/H] and [$\alpha$/Fe] for each component. The vertical component is especially clean, with $\sigma_{v_z}$ increasing from a low of $8$ km s$^{-1}$ at solar-[$\alpha$/Fe] and [Fe/H] to a high of more than 50 km s$^{-1}$ for more metal-poor and [$\alpha$/Fe] enhanced populations. We find no evidence of a large decrease in the velocity dispersion of the higher-[$\alpha$/Fe] populations as claimed in analysis prior to Gaia DR2, although the trend of increasing velocity dispersion with [$\alpha$/Fe] for the same metallicity does significantly flatten at high-[$\alpha$/Fe]. The eccentricity distribution for local stars varies most strongly as a function of [$\alpha$/Fe], where stars with [$\alpha$/Fe]$<0.1$ dex having generally circular orbits ($e<0.15$), while the median eccentricity increases rapidly for more [$\alpha$/Fe]-enhanced stellar populations up to $e\sim0.35$. These [$\alpha$/Fe]-enhanced populations have guiding radii consistent with origins in the inner Galaxy. Of the stars with metallicities much higher than the local ISM ([Fe/H]>0.1 dex), we find that more than 70\% have $e<0.2$ and are likely observed in the solar neighbourhood through churning/migration rather than blurring effects, as the epicyclic motion for these stars is not large enough to reach the radii at which they were likely born based on their metallicity.

We use the second data releases of the ESA Gaia astrometric survey and the high-resolution GALAH spectroscopic survey to analyse the structure of our Galaxy's disc components. With GALAH, we separate the alpha-rich and alpha-poor discs (with respect to Fe), which are superposed in both position and velocity space, and examine their distributions in action space. We study the distribution of stars in the zV_z phase plane, for both V_phi and V_R, and recover the remarkable "phase spiral" discovered by Gaia. We identify the anticipated quadrupole signature in zV_z of a tilted velocity ellipsoid for stars above and below the Galactic plane. By connecting our work with earlier studies, we show that the phase spiral is likely to extend well beyond the narrow solar neighbourhood cylinder in which it was found. The phase spiral is a signature of corrugated waves that propagate through the disc, and the associated non-equilibrium phase mixing. The radially asymmetric distribution of stars involved in the phase spiral reveals that the corrugation, which is mostly confined to the alpha-poor disc, grows in z-amplitude with increasing radius. We present new simulations of tidal disturbance of the Galactic disc by the Sagittarius (Sgr) dwarf. The effect on the zV_z phase plane lasts >2 Gyr but a subsequent disc crossing wipes out the coherent structure. We find that the phase spiral was excited < 0.5 Gyr ago by an object like Sgr with total mass 3 x 10^10 Msun (stripped down from 5 x 10^10 Msun when it first entered the halo) passing through the plane.

It is well known that most of the stars form in rich clusters. However, recent $Spitzer$ observations have shown that a significant number of stars also form in distributed mode, origin of which is not well understood. In this work, we aim to investigate clustered and distributed mode of star formation in the W4 complex. To do so, we identified and characterized the young stellar population associated with the region using homogeneous infra-red data-sets obtained from 2MASS, GLIMPSE, MIPS and WISE surveys. We make stellar surface density and minimum spanning tree maps to identify young clusters, and use {\it Spitzer} images to identify irradiated structures, such as elephant trunk-like structures (ETLSs) and pillars in the region. The surface density distribution of the young stellar objects (YSOs) reveals three new clusterings and $\sim$ 50\% distributed protostars in the H{\sc ii} region. The clusters are of low-mass nature but significantly younger than the central cluster IC~1805. We identified $\sim$ 38 ETLSs in the region, a majority of which consist of one or a few stars at their tips. We find these stars are low-mass ($<$ 2~M$_\odot$) YSOs, located at the outskirts ($>$ 17 pc) of the cluster IC~1805 and are part of scattered distributed population. We argued that the star formation in the ETLSs of W4 is going on possibly due to triggering effect of expanding W4 bubble. Although high-resolution photometric and spectroscopic data would be required to confirm the scenario, nonetheless, we discuss the implications of this scenario for our understanding of distributed low-mass star formation in cloud complexes as opposed to other mechanisms such as turbulent fragmentation and dynamical ejection.

The overlap between the spectroscopic Galactic Archaeology with HERMES (GALAH) survey & $Gaia$ provides a high-dimensional chemodynamical space of unprecedented size. We present a first analysis of a subset of this overlap, of 7066 dwarf, turn-off, & sub-giant stars. [...] We investigate correlations between chemical compositions, ages, & kinematics for this sample. Stellar parameters & elemental abundances are derived from the GALAH spectra with the spectral synthesis code SME. [...] We report Li, C, O, Na, Mg, Al, Si, K, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, as well as Ba & we note that we employ non-LTE calculations for Li, O, Al, & Fe. We show that the use of astrometric & photometric data improves the accuracy of the derived spectroscopic parameters, especially $\log g$. [...] we recover the result that stars of the high-$\alpha$ sequence are typically older than stars in the low-$\alpha$ sequence, the latter spanning $-0.7<$[Fe/H]$<+0.5$. While these two sequences become indistinguishable in [$\alpha$/Fe] vs. [Fe/H] at the metal-rich regime, we find that age can be used to separate stars from the extended high-$\alpha$ & the low-$\alpha$ sequence even in this regime. [...] we find that the old stars ($>8$ Gyr have lower angular momenta $L_z$ than the Sun, which implies that they are on eccentric orbits & originate from the inner disk. Contrary to some previous smaller scale studies we find a continuous evolution in the high-$\alpha$-sequence up to super-solar [Fe/H] rather than a gap, which has been interpreted as a separate "high-$\alpha$ metal-rich" population. Stars in our sample that are younger than 10 Gyr, are mainly found on the low $\alpha$-sequence & show a gradient in $L_z$ from low [Fe/H] ($L_z>L_{z,\odot}$) towards higher [Fe/H] ($L_z<L_{z,\odot}$), which implies that the stars at the ends of this sequence are likely not originating from the close solar vicinity.

Swan bands - characteristic molecular absorption features of the C$_2$ molecule - are a spectroscopic signature of carbon-enhanced stars. They can also be used to identify carbon-enhanced metal-poor (CEMP) stars. The GALAH (GALactic Archaeology with Hermes) is a magnitude-limited survey of stars producing high-resolution, high signal-to-noise spectra. We used 627,708 GALAH spectra to search for carbon-enhanced stars with a supervised and unsupervised classification algorithm, relying on the imprint of the Swan bands. We identified 918 carbon-enhanced stars, including 12 already described in the literature. An unbiased selection function of the GALAH survey allows us to perform a population study of carbon-enhanced stars. Most of them are giants, out of which we find 28 CEMP candidates. A large fraction of our carbon-enhanced stars with repeated observations show variation in radial velocity, hinting that there is a large fraction of variables among them. 32 of the detected stars also show strong Lithium enhancement in their spectra.

If the Galaxy is axisymmetric and in dynamical equilibrium, we expect negligible fluctuations in the residual line-of-sight velocity field. Recent results using the \apg{} survey find significant fluctuations in velocity for stars in the midplane ($|z|<$0.25 kpc) out to 5 kpc, suggesting that the dynamical influence of non-axisymmetric features i.e., the Milky Way's bar, spiral arms and merger events extends out to the Solar neighborhood. Their measured power spectrum has a characteristic amplitude of 11 \kms{} on a scale of 2.5 kpc. The existence of such large-scale streaming motions has important implications for determining the Sun's motion about the Galactic Centre. Using Red Clump stars from \glh{} and \apg{}, we map the line-of-sight velocities around the Sun (d$<$5 kpc), and $|z|<$1.25 kpc from the midplane. By subtracting a smooth axisymmetric model for the velocity field, we study the residual fluctuations and compare our findings with mock survey generated by \glx{}. We find negligible large-scale fluctuations away from the plane. In the mid-plane, we reproduce the earlier \apg{} power spectrum but with 20\% smaller amplitude (9.3 \kms{}) after taking into account a few systematics (e.g., volume completeness). Using a flexible axisymmetric model the power-amplitude is further reduced to 6.3 \kms{}. Additionally, our simulations show that, in the plane, distances are underestimated for high-mass Red Clump stars which can lead to spurious power-amplitude of about 5.2 \kms{}. Taking this into account, we estimate the amplitude of real fluctuations to be $<$4.6 \kms{}, about a factor of three less than the \apg{} result.

GALAH is a large-scale magnitude-limited southern stellar spectroscopic survey. Its second data release (GALAH DR2) provides values of stellar parameters and abundances of 23 elements for 342,682 stars (Buder et al.). Here we add a description of the public release of radial velocities with a typical accuracy of 0.1 km/s for 336,215 of these stars, achievable due to the large wavelength coverage, high resolving power and good signal to noise ratio of the observed spectra, but also because convective motions in stellar atmosphere and gravitational redshift from the star to the observer are taken into account. In the process we derive medians of observed spectra which are nearly noiseless, as they are obtained from between 100 and 1116 observed spectra belonging to the same bin with a width of 50 K in temperature, 0.2 dex in gravity, and 0.1 dex in metallicity. Publicly released 1181 median spectra have a resolving power of 28,000 and trace the well-populated stellar types with metallicities between -0.6 and +0.3. Note that radial velocities from GALAH are an excellent match to the accuracy of velocity components along the sky plane derived by Gaia for the same stars. The level of accuracy achieved here is adequate for studies of dynamics within stellar clusters, associations and streams in the Galaxy. So it may be relevant for studies of the distribution of dark matter.

We present spatially resolved analysis of a lensed galaxy, SDSS1958+5950 at $z = 2.225$, from the Cambridge Sloan Survey of Wide Arcs in the Sky (CASSOWARY). We use our new high resolution imaging data to construct a robust lens model for the galaxy group at $z = 0.214$. We employ the updated lens model to combine the Integral Field Spectrographic observations on two highly distorted images of the lensed target. We adopt a forward-modeling approach to deconvolve the effects of point spread function from the combined source-plane reconstruction. The approach is adapted to the lens model magnification and enables a resolution of $\sim$170 pc in the galaxy-source plane. We propose an ongoing merger as the origin of the lensed system on the basis of its source-plane morphology, kinematics and rest-frame emission line ratios. Using our novel technique of adaptive coadded source plane reconstruction, we are able to detect different components in the velocity gradient that were not seen in previous studies of this object, plausibly belonging to different components in the merging system.

Open cluster members are coeval and share the same initial bulk chemical compositions. Consequently, differences in surface abundances between members of a cluster that are at different evolutionary stages can be used to study the effects of mixing and internal chemical processing. We carry out an abundance analysis of seven elements (Li, O, Na, Mg, Al, Si, Fe) in 66 stars belonging to the open cluster M67, based on high resolution GALAH spectra, 1D MARCS model atmospheres, and, for the first time for a large sample of M67 stars, non-local thermodynamic equilibrium (non-LTE) radiative transfer. From the non-LTE analysis, we find a typical star-to-star scatter in the abundance ratios of around 0.05 dex; this scatter is slightly but systematically larger when LTE is assumed instead. We find trends in the abundance ratios with effective temperature, indicating systematic differences in the surface abundances between turn-off and giant stars; these trends are more pronounced when LTE is assumed. However, in the non-LTE analysis, most of the element trends have been flattened. Two species are exceptions to this behaviour, namely Al and Si, which both clearly display remaining trends in the non-LTE analysis. We comment on the possible origin of these trends, by comparing them with recent stellar models that include atomic diffusion.

The stars within our Galactic halo presents a snapshot of its ongoing growth and evolution, probing galaxy formation directly. Here, we present our first analysis of the stellar halo from detailed maps of Blue Horizontal Branch (BHB) stars drawn from the SkyMapper Southern Sky Survey. To isolate candidate BHB stars from the overall population, we develop a machine-learning approach through the application of an Artificial Neural Network (ANN), resulting in a relatively pure sample of target stars. From this, we derive the absolute $u$ magnitude for the BHB sample to be $\sim2\ mag$, varying slightly with $(v-g)_0$ and $(u-v)_0$ colours. We examine the BHB number density distribution from 5272 candidate stars, deriving a double power-law with a break radius of $r_s = 11.8\pm0.3\ kpc$, and inner and outer slopes of $\alpha_{in} = -2.5\pm0.1$ and $\alpha_{out} = -4.5\pm0.3$ respectively. Through isochrone fitting of simulated BHB stars, we find a colour-age/metallicity correlation, with older/more metal-poor stars being bluer, and establish a parameter to indicate this age (or metallicity) variation. Using this, we construct the three-dimensional population distribution of BHB stars in the halo and identify significant substructure. Finally, in agreement with previous studies, we also identify a systemic age/metallicity shift spanning $\sim3\ kpc$ to $ \sim20\ kpc$ in galactocentric distance.

Sparse open clusters can be found at high galactic latitudes where loosely populated clusters are more easily detected against the lower stellar background. Because most star formation takes place in the thin disk, the observed population of clusters far from the Galactic plane is hard to explain. We combined spectral parameters from the GALAH survey with the Gaia DR2 catalogue to study the dynamics and chemistry of five old sparse high-latitude clusters in more detail. We find that four of them (NGC 1252, NGC 6994, NGC 7772, NGC 7826) - originally classified in 1888 - are not clusters but are instead chance projections on the sky. Member stars quoted in the literature for these four clusters are unrelated in our multi-dimensional physical parameter space; the published cluster properties are therefore irrelevant. We confirm the existence of morphologically similar NGC 1901 for which we provide a probabilistic membership analysis. An overdensity in three spatial dimensions proves to be enough to reliably detect sparse clusters, but the whole 6-dimensional space must be used to identify members with high confidence, as demonstrated in the case of NGC 1901.

We present a study using the second data release of the GALAH survey of stellar parameters and elemental abundances of 15 pairs of stars identified by Oh et al 2017. They identified these pairs as potentially co-moving pairs using proper motions and parallaxes from Gaia DR1. We find that 11 very wide (>1.7 pc) pairs of stars do in fact have similar Galactic orbits, while a further four claimed co-moving pairs are not truly co-orbiting. Eight of the 11 co-orbiting pairs have reliable stellar parameters and abundances, and we find that three of those are quite similar in their abundance patterns, while five have significant [Fe/H] differences. For the latter, this indicates that they could be co-orbiting because of the general dynamical coldness of the thin disc, or perhaps resonances induced by the Galaxy, rather than a shared formation site. Stars such as these, wide binaries, debris of past star formation episodes, and coincidental co-orbiters, are crucial for exploring the limits of chemical tagging in the Milky Way.

We report the serendipitous observations of 571 luminous supergiants in the Magellanic Clouds by the spectroscopic GALAH and TESS-HERMES surveys: 434 stars in the Large Magellanic Cloud and 137 in the Small Magellanic Cloud. We also find one star that appears associated with structured star formation in the Magellanic Bridge. Both of these surveys are aimed at the local volume of the Galaxy but have simple, magnitude-limited selection functions that mean they include some observations of luminous extra-Galactic stars. The surveys determine stellar parameter and abundances using The Cannon, a data-driven generative modelling approach. In this work, we explore the results from The Cannon when it is fed the spectra of these intrinsically luminous supergiants in the Magellanic Clouds, which are well outside the normal bounds of The Cannon's training set. We find that, although the parameters are astrophysically incorrect, the $v\sin i$ and the abundances of lithium, barium, and magnesium are excellent discriminants of these stars. It shows that in the future, with an expanded training set, it should be possible to determine accurate values for these types of stars.

Temporal and spectral characteristics of X-ray emission from 60 flares of intensity $\ge$C class observed by Solar X-ray Spectrometer (SOXS) during 2003-2011 are presented. We analyse the X-ray emission observed in four and three energy bands by the Si and CZT detectors, respectively. The number of peaks in the intensity profile of the flares varies between 1 and 3. We find moderate correlation (R$\simeq$0.2) between the rise time and the peak flux of the first peak of the flare irrespective to energy band, which is indicative of its energy-independent nature. Moreover, magnetic field complexity of the flaring region is found to be anti-correlated (R=0.61) with the rise time of the flares while positively correlated (R=0.28) with the peak flux of the flare. The time delay between the peak of the X-ray emission in a given energy band and that in the 25-30 keV decreases with increasing energy suggesting conduction cooling to be dominant in the lower energies. Analysis of 340 spectra from 14 flares reveals that the peak of Differential Emission Measure (DEM) evolution delays by 60-360 s relative to that of the temperature, and this time delay is inversely proportional to the peak flux of the flare. We conclude that temporal and intensity characteristics of flares are energy dependent as well as magnetic field configuration of the active region.

We present criteria for the photometric selection of M-dwarfs using all-sky photometry, with a view to identifying M-dwarf candidates for inclusion in the input catalogues of upcoming all-sky surveys, including TESS and FunnelWeb. The criteria are based on Gaia, WISE and 2MASS all-sky photometry, and deliberately do not rely on astrometric information. In the lead-up to the availability of truly distance-limited samples following the release of Gaia DR2, this approach has the significant benefit of delivering a sample unbiased with regard to space velocity. Our criteria were developed by using Galaxia synthetic galaxy model predictions to evaluate both M-dwarf completeness and false-positive detections (i.e. non-M-dwarf contamination rates). In addition to the previously known sensitivity of J-H colour for giant-dwarf discrimination at cool temperatures, we find the WISE W1-W2 colour is also effective at discriminating M-dwarfs from cool giants. We have derived two sets of Gaia G > 14.5 criteria - a "high-completeness" set that contains 78,340 stars, of which 30.7-44.4% are expected to be M-dwarfs and contains 99.3% of the total number of expected M-dwarfs; and a "low-contamination" set that prioritises the stars most likely to be M-dwarfs at a cost of a reduction in completeness. This subset contains 40,505 stars and is expected to be comprised of 58.7-64.1% M-dwarfs, with a completeness of 98%. Comparison of the high-completeness set with the TESS Input Catalogue has identified 234 stars not currently in that catalogue, which preliminary analysis suggests could be useful M-dwarf targets for TESS. We also compared the criteria to selection via absolute magnitude and a combination of both methods. We found that colour selection in combination with an absolute magnitude limit provides the most effective way of selecting M-dwarfs en masse.

The Galactic Archaeology with HERMES (GALAH) survey is a large-scale stellar spectroscopic survey of the Milky Way and designed to deliver chemical information complementary to a large number of stars covered by the $Gaia$ mission. We present the GALAH second public data release (GALAH DR2) containing 342,682 stars. For these stars, the GALAH collaboration provides stellar parameters and abundances for up to 23 elements to the community. Here we present the target selection, observation, data reduction and detailed explanation of how the spectra were analysed to estimate stellar parameters and element abundances. For the stellar analysis, we have used a multi-step approach. We use the physics-driven spectrum synthesis of Spectroscopy Made Easy (SME) to derive stellar labels ($T_\mathrm{eff}$, $\log g$, $\mathrm{[Fe/H]}$, $\mathrm{[X/Fe]}$, $v_\mathrm{mic}$, $v \sin i$, $A_{K_S}$) for a representative training set of stars. This information is then propagated to the whole survey with the data-driven method of $The~Cannon$. Special care has been exercised in the spectral synthesis to only consider spectral lines that have reliable atomic input data and are little affected by blending lines. Departures from local thermodynamic equilibrium (LTE) are considered for several key elements, including Li, O, Na, Mg, Al, Si, and Fe, using 1D MARCS stellar atmosphere models. Validation tests including repeat observations, Gaia benchmark stars, open and globular clusters, and K2 asteroseismic targets lend confidence in our methods and results. Combining the GALAH DR2 catalogue with the kinematic information from $Gaia$ will enable a wide range of Galactic Archaeology studies, with unprecedented detail, dimensionality, and scope.

The primary goal of Galactic archaeology is to learn about the origin of the Milky Way from the detailed chemistry and kinematics of millions of stars. Wide-field multi-fibre spectrographs are increasingly used to obtain spectral information for huge samples of stars. Some surveys (e.g. GALAH) are attempting to measure up to 30 separate elements per star. Stellar abundance spectroscopy is a subtle art that requires a very high degree of spectral uniformity across each of the fibres. However wide-field spectrographs are notoriously non-uniform due to the fast output optics necessary to image many fibre outputs onto the detector. We show that precise spectroscopy is possible with such instruments across all fibres by employing a photonic comb -- a device that produces uniformly spaced spots of light on the CCD to precisely map complex aberrations. Aberrations are parametrized by a set of orthogonal moments with $\sim100$ independent parameters. We then reproduce the observed image by convolving high resolution spectral templates with measured aberrations as opposed to extracting the spectra from the observed image. Such a forward modeling approach also trivializes some spectroscopic reduction problems like fibre cross-talk, and reliably extracts spectra with a resolution $\sim2.3$ times above the nominal resolution of the instrument. Our rigorous treatment of optical aberrations also encourages a less conservative spectrograph design in the future.

Using GALAH survey data of nearby stars, we look at how structure in the planar (u,v) velocity distribution depends on metallicity and on viewing direction within the Galaxy. In nearby stars, with distance d < 1 kpc, the Hercules stream is most strongly seen in higher metallicity stars [Fe/H] > 0.2. The Hercules stream peak v value depends on viewed galactic longitude, which we interpret as due to the gap between the stellar stream and more circular orbits being associated with a specific angular momentum value of about 1640 km/s kpc. The association of the gap with a particular angular momentum value supports a bar resonant model for the Hercules stream. Moving groups previously identified in Hipparcos observations are easiest to see in stars nearer than 250 pc, and their visibility and peak velocities in the velocity distributions depends on both viewing direction (galactic longitude and hemisphere) and metallicity. We infer that there is fine structure in local velocity distributions that varies over distances of a few hundred pc in the Galaxy.

Using data from the GALAH pilot survey, we determine properties of the Galactic thin and thick disks near the solar neighbourhood. The data cover a small range of Galactocentric radius ($7.9 \leq R_\mathrm{GC} \leq 9.5$ kpc), but extend up to 4 kpc in height from the Galactic plane, and several kpc in the direction of Galactic anti-rotation (at longitude $260 ^\circ \leq \ell \leq 280^\circ$). This allows us to reliably measure the vertical density and abundance profiles of the chemically and kinematically defined `thick' and `thin' disks of the Galaxy. The thin disk (low-$\alpha$ population) exhibits a steep negative vertical metallicity gradient, at d[M/H]/d$z=-0.18 \pm 0.01$ dex kpc$^{-1}$, which is broadly consistent with previous studies. In contrast, its vertical $\alpha$-abundance profile is almost flat, with a gradient of d[$\alpha$/M]/d$z$ = $0.008 \pm 0.002$ dex kpc$^{-1}$. The steep vertical metallicity gradient of the low-$\alpha$ population is in agreement with models where radial migration has a major role in the evolution of the thin disk. The thick disk (high-$\alpha$ population) has a weaker vertical metallicity gradient d[M/H]/d$z = -0.058 \pm 0.003$ dex kpc$^{-1}$. The $\alpha$-abundance of the thick disk is nearly constant with height, d[$\alpha$/M]/d$z$ = $0.007 \pm 0.002$ dex kpc$^{-1}$. The negative gradient in metallicity and the small gradient in [$\alpha$/M] indicate that the high-$\alpha$ population experienced a settling phase, but also formed prior to the onset of major SNIa enrichment. We explore the implications of the distinct $\alpha$-enrichments and narrow [$\alpha$/M] range of the sub-populations in the context of thick disk formation.

We use data from the Radial Velocity Experiment (RAVE) and the Tycho-Gaia astrometric solution catalogue (TGAS) to compute the velocity fields yielded by the radial (VR), azimuthal (Vphi) and vertical (Vz) components of associated Galactocentric velocity. We search in particular for variation in all three velocity components with distance above and below the disc midplane, as well as how each component of Vz (line-of-sight and tangential velocity projections) modifies the obtained vertical structure. To study the dependence of velocity on proper motion and distance we use two main samples: a RAVE sample including proper motions from the Tycho-2, PPMXL and UCAC4 catalogues, and a RAVE-TGAS sample with inferred distances and proper motions from the TGAS and UCAC5 catalogues. In both samples, we identify asymmetries in VR and Vz. Below the plane we find the largest radial gradient to be dVR / dR = -7.01+- 0.61 km\s kpc, in agreement with recent studies. Above the plane we find a similar gradient with dVR / dR= -9.42+- 1.77 km\s kpc. By comparing our results with previous studies, we find that the structure in Vz is strongly dependent on the adopted proper motions. Using the Galaxia Milky Way model, we demonstrate that distance uncertainties can create artificial wave-like patterns. In contrast to previous suggestions of a breathing mode seen in RAVE data, our results support a combination of bending and breathing modes, likely generated by a combination of external or internal and external mechanisms.

TIFR Near Infrared Imaging Camera-II is a closed-cycle Helium cryo-cooled imaging camera equipped with a Raytheon 512 x 512 pixels InSb Aladdin III Quadrant focal plane array having sensitivity to photons in the 1-5 microns wavelength band. In this paper, we present the performance of the camera on the newly installed 3.6-m Devasthal Optical Telescope (DOT) based on the calibration observations carried out during 2017 May 11-14 and 2017 October 7-31. After the preliminary characterization, the camera has been released to the Indian and Belgian astronomical community for science observations since 2017 May. The camera offers a field-of-view of ~86.5 arcsec x 86.5 arcsec on the DOT with a pixel scale of 0.169 arcsec. The seeing at the telescope site in the near-infrared bands is typically sub-arcsecond with the best seeing of ~0.45 arcsec realized in the near-infrared K-band on 2017 October 16. The camera is found to be capable of deep observations in the J, H and K bands comparable to other 4-m class telescopes available world-wide. Another highlight of this camera is the observational capability for sources up to Wide-field Infrared Survey Explorer (WISE) W1-band (3.4 microns) magnitudes of 9.2 in the narrow L-band (nbL; lambda_{cen} ~3.59 microns). Hence, the camera could be a good complementary instrument to observe the bright nbL-band sources that are saturated in the Spitzer-Infrared Array Camera ([3.6] <= 7.92 mag) and the WISE W1-band ([3.4] <= 8.1 mag). Sources with strong polycyclic aromatic hydrocarbon (PAH) emission at 3.3 microns are also detected. Details of the observations and estimated parameters are presented in this paper.

The spherical Jeans equation is a widely used tool for dynamical study of gravitating systems in astronomy. Here we test its efficacy in robustly weighing the mass of Milky Way analogues, given they need not be in equilibrium or even spherical. Utilizing Milky Way stellar halos simulated in accordance with $\Lambda{\rm CDM}$ cosmology by Bullock and Johnston (2005) and analysing them under the Jeans formalism, we recover the underlying mass distribution of the parent galaxy, within distance $r/{\rm kpc}\in[10,100]$, with a bias of $\sim12\%$ and a dispersion of $\sim14\%$. Additionally, the mass profiles of triaxial dark matter halos taken from the SURFS simulation, within scaled radius $0.2<r/r_{\rm max}<3$, are measured with a bias of $\sim-2.4\%$ and a dispersion of $\sim10\%$. The obtained dispersion is not because of Poisson noise due to small particle numbers as it is twice the later. We interpret the dispersion to be due to the inherent nature of the $\Lambda{\rm CDM}$ halos, for example being aspherical and out-of-equilibrium. Hence the dispersion obtained for stellar halos sets a limit of about $12\%$ (after adjusting for random uncertainty) on the accuracy with which the mass profiles of the Milky Way-like galaxies can be reconstructed using the spherical Jeans equation. This limit is independent of the quantity and quality of the observational data. The reason for a non zero bias is not clear, hence its interpretation is not obvious at this stage.

Our nearest large cosmological neighbour, the Andromeda galaxy (M31), is a dynamical system, and an accurate measurement of its total mass is central to our understanding of its assembly history, the life-cycles of its satellite galaxies, and its role in shaping the Local Group environment. Here, we apply a novel approach to determine the dynamical mass of M31 using high velocity Planetary Nebulae (PNe), establishing a hierarchical Bayesian model united with a scheme to capture potential outliers and marginalize over tracers unknown distances. With this, we derive the escape velocity run of M31 as a function of galacto-centric distance, with both parametric and non-parametric approaches. We determine the escape velocity of M31 to be $470\pm{40}$ km s$^{-1}$ at a galacto-centric distance of 15 kpc, and also, derive the total potential of M31, estimating the virial mass and radius of the galaxy to be $0.8\pm{0.1}\times10^{12}\,M_\odot$ and $240\pm{10}$ kpc, respectively. Our M31 mass is on the low-side of the measured range, this supports the lower expected mass of the M31-Milky Way system from the timing and momentum arguments, satisfying the HI constraint on circular velocity between $10\lesssim R/\textrm{kpc}<35$, and agreeing with the stellar mass Tully-Fisher relation. To place these results in a broader context, we compare them to the key predictions of the $\Lambda{\rm CDM}$ cosmological paradigm, including the stellar-mass-halo-mass and the dark matter halo concentration-virial mass correlation, and finding it to be an outlier to this relation.

Obtaining accurate and precise masses and ages for large numbers of giant stars is of great importance for unraveling the assemblage history of the Galaxy. In this paper, we estimate masses and ages of 6940 red giant branch (RGB) stars with asteroseismic parameters deduced from \emph{Kepler} photometry and stellar atmospheric parameters derived from LAMOST spectra. The typical uncertainties of mass is a few per cent, and that of age is $\sim$\,20 per cent. The sample stars reveal two separate sequences in the age -- [$\alpha$/Fe] relation -- a high--$\alpha$ sequence with stars older than $\sim$\,8\,Gyr and a low--$\alpha$ sequence composed of stars with ages ranging from younger than 1\,Gyr to older than 11\,Gyr. We further investigate the feasibility of deducing ages and masses directly from LAMOST spectra with a machine learning method based on kernel based principal component analysis, taking a sub-sample of these RGB stars as a training data set. We demonstrate that ages thus derived achieve an accuracy of $\sim$\,24 per cent. We also explored the feasibility of estimating ages and masses based on the spectroscopically measured carbon and nitrogen abundances. The results are quite satisfactory and significantly improved compared to the previous studies.

We present a multi-wavelength analysis of a $\sim$30$' \times$30$'$ area around the Sh 2-53 region (hereafter S53 complex), which is associated with at least three H II regions, two mid-infrared bubbles (N21 and N22), and infrared dark clouds. The $^{13}$CO line data trace the molecular content of the S53 complex in a velocity range of 36--60 km s$^{-1}$, and show the presence of at least three molecular components within the selected area along this direction. Using the observed radio continuum flux of the H II regions, the derived spectral types of the ionizing sources agree well with the previously reported results. The S53 complex harbors clusters of young stellar objects (YSOs) that are identified using the photometric 2--24 $\mu$m magnitudes. It also hosts several massive condensations (3000-30000 $M_\odot$) which are traced in the {\em Herschel} column density map. The complex is found at the junction of at least five molecular filaments, and the flow of gas toward the junction is evident in the velocity space of the $^{13}$CO data. Together, the S53 complex is embedded in a very similar "hub-filament" system to those reported in Myers, and the active star formation is evident towards the central "hub" inferred by the presence of the clustering of YSOs.

Accurate and precise radius estimates of transiting exoplanets are critical for understanding their compositions and formation mechanisms. To know the planet, we must know the host star in as much detail as possible. We present first results from the K2-HERMES project, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R$\sim$28,000 spectra of up to 360 stars in one exposure. This ongoing project aims to derive self-consistent spectroscopic parameters for about half of K2 target stars. We present complete stellar parameters and isochrone-derived masses and radii for 46 stars hosting 57 K2 candidate planets in Campaigns 1-3. Our revised host-star radii cast severe doubt on three candidate planets: EPIC\,201407812.01, EPIC\,203070421.01, and EPIC\,202843107.01, all of which now have inferred radii well in excess of the largest known inflated Jovian planets.

We present wide and deep photometry of the northwest part of the halo of the Andromeda galaxy (M31) using Hyper Suprime-Cam on the Subaru Telescope. The survey covers 9.2 deg$^{2}$ field in the $g$, $i$, and $NB515$ bands and shows a clear red giant branch (RGB) of M31's halo stars and a pronounced red clump (RC) feature. The spatial distribution of RC stars shows a prominent stream feature, the North Western (NW) Stream, and a diffuse substructure in the south part of our survey field. We estimate the distances based on the RC method and obtain $(m-M)$ = 24.63$\pm 0.191$(random)$\pm0.057$(systematic) and 24.29$\pm 0.211$(random)$\pm0.057$(systematic) mag for the NW stream and diffuse substructure, respectively, implying that the NW Stream is located behind M31, whereas the diffuse substructure is located in front. We also estimate line-of-sight distances along the NW Stream and find that the south part of the stream is $\sim$20 kpc closer to us relative to the north part. The distance to the NW Stream inferred from the isochrone fitting to the color-magnitude diagram favors the RC-based distance, but the TRGB-based distance estimated for $NB515$-selected RGB stars does not agree with it. The surface number density distribution of RC stars across the NW Stream is found to be approximately Gaussian with a FWHM of $\sim$25 arcmin (5.7 kpc), with a slight skew to the south-west side. That along the NW Stream shows a complicated structure including variations in number density and a significant gap in the stream.

The variability of young stellar objects is mostly driven by star-disk interactions. In long-term photometric monitoring of the accreting T Tauri star GI Tau, we detect extinction events with typical depths of $\Delta V \sim 2.5$ mag that last for days-to-months and often appear to occur stochastically. In 2014 - 2015, extinctions that repeated with a quasi-period of 21 days over several months is the first empirical evidence of slow warps predicted from MHD simulations to form at a few stellar radii away from the central star. The reddening is consistent with $R_V=3.85\pm0.5$ and, along with an absence of diffuse interstellar bands, indicates that some dust processing has occurred in the disk. The 2015 -- 2016 multi-band lightcurve includes variations in spot coverage, extinction, and accretion, each of which results in different traces in color-magnitude diagrams. This lightcurve is initially dominated by a month-long extinction event and return to the unocculted brightness. The subsequent light-curve then features spot modulation with a 7.03 day period, punctuated by brief, randomly-spaced extinction events. The accretion rate measured from $U$-band photometry ranges from $1.3\times10^{-8}$ to $1.1\times10^{-10}$ M$_\odot$ yr$^{-1}$ (excluding the highest and lowest 5% of high- and low- accretion rate outliers), with an average of $4.7 \times 10^{-9}$ M$_\odot$ yr$^{-1}$. A total of 50% of the mass is accreted during bursts of $>12.8\times10^{-9}$ M$_\odot$ yr${^{-1}}$, which indicates limitations on analyses of disk evolution using single-epoch accretion rates.

We report an integral-field spectroscopic (IFS) observation of a gravitationally lensed spiral galaxy A1689B11 at redshift $z=2.54$. It is the most ancient spiral galaxy discovered to date and the second kinematically confirmed spiral at $z\gtrsim2$. Thanks to gravitational lensing, this is also by far the deepest IFS observation with the highest spatial resolution ($\sim$ 400 pc) on a spiral galaxy at a cosmic time when the Hubble sequence is about to emerge. After correcting for a lensing magnification of 7.2 $\pm$ 0.8, this primitive spiral disk has an intrinsic star formation rate of 22 $\pm$ 2 $M_{\odot}$ yr$^{-1}$, a stellar mass of 10$^{9.8 \pm 0.3}$$M_{\odot}$ and a half-light radius of $r_{1/2}=2.6 \pm 0.7$ kpc, typical of a main-sequence star-forming (SF) galaxy at $z\sim2$. However, the H\alpha\ kinematics show a surprisingly tranquil velocity field with an ordered rotation ($V_{\rm c}$ = 200 $\pm$ 12 km/s) and uniformly small velocity dispersions ($V_{\rm \sigma, mean}$ = 23 $\pm$ 4 km/s and $V_{\rm \sigma, outer-disk}$ = 15 $\pm$ 2 km/s). The low gas velocity dispersion is similar to local spiral galaxies and is consistent with the classic density wave theory where spiral arms form in dynamically cold and thin disks. We speculate that A1689B11 belongs to a population of rare spiral galaxies at $z\gtrsim2$ that mark the formation epoch of thin disks. Future observations with JWST will greatly increase the sample of these rare galaxies and unveil the earliest onset of spiral arms.

The influence of magnetic fields (B-fields) in the formation and evolution of bipolar bubbles, due to the expanding ionization fronts (I-fronts) driven by the Hii regions that are formed and embedded in filamentary molecular clouds, has not been well-studied yet. In addition to the anisotropic expansion of I-fronts into a filament, B-fields are expected to introduce an additional anisotropic pressure which might favor expansion and propagation of I-fronts to form a bipolar bubble. We present results based on near-infrared polarimetric observations towards the central $\sim$8'$\times$8' area of the star-forming region RCW57A which hosts an Hii region, a filament, and a bipolar bubble. Polarization measurements of 178 reddened background stars, out of the 919 detected sources in the JHKs-bands, reveal B-fields that thread perpendicular to the filament long axis. The B-fields exhibit an hour-glass morphology that closely follows the structure of the bipolar bubble. The mean B-field strength, estimated using the Chandrasekhar-Fermi method, is 91$\pm$8 {\mu}G. B-field pressure dominates over turbulent and thermal pressures. Thermal pressure might act in the same orientation as those of B-fields to accelerate the expansion of those I-fronts. The observed morphological correspondence among the B-fields, filament, and bipolar bubble demonstrate that the B-fields are important to the cloud contraction that formed the filament, gravitational collapse and star formation in it, and in feedback processes. The latter include the formation and evolution of mid-infrared bubbles by means of B-field supported propagation and expansion of I-fronts. These may shed light on preexisting conditions favoring the formation of the massive stellar cluster in RCW57A.

The Transiting Exoplanet Survey Satellite (TESS) will provide high precision time-series photometry for millions of stars with at least a half-hour cadence. Of particular interest are the circular regions of 12-degree radius centered around the ecliptic poles that will be observed continuously for a full year. Spectroscopic stellar parameters are desirable to characterize and select suitable targets for TESS, whether they are focused on exploring exoplanets, stellar astrophysics, or Galactic archaeology. Here, we present spectroscopic stellar parameters ($T_{\rm eff}$, $\log g$, [Fe/H], $v \sin i$, $v_{\rm micro}$) for about 16,000 dwarf and subgiant stars in TESS' southern continuous viewing zone. For almost all the stars, we also present Bayesian estimates of stellar properties including distance, extinction, mass, radius, and age using theoretical isochrones. Stellar surface gravity and radius are made available for an additional set of roughly 8,500 red giants. All our target stars are in the range $10<V<13.1$. Among them, we identify and list 227 stars belonging to the Large Magellanic Cloud. The data were taken using the the High Efficiency and Resolution Multi-Element Spectrograph (HERMES, R $\sim 28,000$) at the Anglo-Australian Telescope as part of the TESS-HERMES survey. Comparing our results with the TESS Input Catalog (TIC) shows that the TIC is generally efficient in separating dwarfs and giants, but it has flagged more than hundred cool dwarfs ($T_{\rm eff}< 4800$ K) as giants, which ought to be high-priority targets for the exoplanet search. The catalog can be accessed via http://www.physics.usyd.edu.au/tess-hermes/ , or at MAST via https://archive.stsci.edu/prepds/tess-hermes/ .

The technique of chemical tagging uses the elemental abundances of stellar atmospheres to `reconstruct' chemically homogeneous star clusters that have long since dispersed. The GALAH spectroscopic survey --which aims to observe one million stars using the Anglo-Australian Telescope -- allows us to measure up to 30 elements or dimensions in the stellar chemical abundance space, many of which are not independent. How to find clustering reliably in a noisy high-dimensional space is a difficult problem that remains largely unsolved. Here we explore t-distributed stochastic neighbour embedding (t-SNE) -- which identifies an optimal mapping of a high-dimensional space into fewer dimensions -- whilst conserving the original clustering information. Typically, the projection is made to a 2D space to aid recognition of clusters by eye. We show that this method is a reliable tool for chemical tagging because it can: (i) resolve clustering in chemical space alone, (ii) recover known open and globular clusters with high efficiency and low contamination, and (iii) relate field stars to known clusters. t-SNE also provides a useful visualization of a high-dimensional space. We demonstrate the method on a dataset of 13 abundances measured in the spectra of 187,000 stars by the GALAH survey. We recover 7 of the 9 observed clusters (6 globular and 3 open clusters) in chemical space with minimal contamination from field stars and low numbers of outliers. With chemical tagging, we also identify two Pleiades supercluster members (which we confirm kinematically), one as far as 6$^\circ$ -- one tidal radius away from the cluster centre.

Deep photometric surveys of the Milky Way have revealed diffuse structures encircling our Galaxy far beyond the "classical" limits of the stellar disk. This paper reviews results from our own and other observational programs, which together suggest that, despite their extreme positions, the stars in these structures were formed in our Galactic disk. Mounting evidence from recent observations and simulations implies kinematic connections between several of these distinct structures. This suggests the existence of collective disk oscillations that can plausibly be traced all the way to asymmetries seen in the stellar velocity distribution around the Sun. There are multiple interesting implications of these findings: they promise new perspectives on the process of disk heating, they provide direct evidence for a stellar halo formation mechanism in addition to the accretion and disruption of satellite galaxies, and, they motivate searches of current and near-future surveys to trace these oscillations across the Galaxy. Such maps could be used as dynamical diagnostics in the emerging field of "Galactoseismology", which promises to model the history of interactions between the Milky Way and its entourage of satellites, as well examine the density of our dark matter halo. As sensitivity to very low surface brightness features around external galaxies increases, many more examples of such disk oscillations will likely be identified. Statistical samples of such features not only encode detailed information about interaction rates and mergers, but also about long sought-after dark matter halo densities and shapes. Models for the Milky Way's own Galactoseismic history will therefore serve as a critical foundation for studying the weak dynamical interactions of galaxies across the universe.

Thanks to modern sky surveys, over twenty stellar streams and overdensity structures have been discovered in the halo of the Milky Way. In this paper, we present an analysis of spectroscopic observations of individual stars from one such structure, "A13", first identified as an overdensity using the M giant catalog from the Two Micron All-Sky Survey. Our spectroscopic observations show that stars identified with A13 have a velocity dispersion of $\lesssim$ 40 $\mathrm{km~s^{-1}}$, implying that it is a genuine coherent structure rather than a chance super-position of random halo stars. From its position on the sky, distance ($\sim$15~kpc heliocentric), and kinematical properties, A13 is likely to be an extension of another low Galactic latitude substructure -- the Galactic Anticenter Stellar Structure (also known as the Monoceros Ring) -- towards smaller Galactic longitude and farther distance. Furthermore, the kinematics of A13 also connect it with another structure in the southern Galactic hemisphere -- the Triangulum-Andromeda overdensity. We discuss these three connected structures within the context of a previously proposed scenario that one or all of these features originate from the disk of the Milky Way.

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.

The first and second moments of stellar velocities encode important information about the formation history of the Galactic halo. However, due to the lack of tangential motion and inaccurate distances of the halo stars, the velocity moments in the Galactic halo have largely remained "known unknowns". Fortunately, our off-centric position within the Galaxy allows us to estimate these moments in the galacto-centric frame using the observed radial velocities of the stars alone. We use these velocities coupled with the Hierarchical Bayesian scheme, which allows easy marginalisation over the missing data (the proper-motion, and uncertainty-free distance and line-of-sight velocity), to measure the velocity dispersions, orbital anisotropy ($\beta$) and streaming motion ($v_{\rm rot}$) of the halo main-sequence turn-off (MSTO) and K-giant (KG) stars in the inner stellar halo (r $\lesssim 15$ kpc). We study the metallicity bias in kinematics of the halo stars and observe that the comparatively metal-rich ([Fe/H]$>-1.4$) and the metal-poor ([Fe/H]$\leq - 1.4$) MSTO samples show a clear systematic difference in $v_{\rm rot} \sim 20-40$ km s$^{-1}$, depending on how restrictive the spatial cuts to cull the disk contamination are. The bias is also detected in KG samples but with less certainty. Both MSTO and KG populations suggest that the inner stellar halo of the Galaxy is radially biased i.e. $\sigma_r>\sigma_\theta$ or $\sigma_\phi$ and $\beta \simeq 0.5$. The apparent metallicity contrariety in the rotation velocity among the halo sub-populations supports the co-existence of multiple populations in the galactic halo that may have formed through distinct formation scenarios, i.e. in-situ versus accretion.

Markov Chain Monte Carlo based Bayesian data analysis has now become the method of choice for analyzing and interpreting data in almost all disciplines of science. In astronomy, over the last decade, we have also seen a steady increase in the number of papers that employ Monte Carlo based Bayesian analysis. New, efficient Monte Carlo based methods are continuously being developed and explored. In this review, we first explain the basics of Bayesian theory and discuss how to set up data analysis problems within this framework. Next, we provide an overview of various Monte Carlo based methods for performing Bayesian data analysis. Finally, we discuss advanced ideas that enable us to tackle complex problems and thus hold great promise for the future. We also distribute downloadable computer software (available at https://github.com/sanjibs/bmcmc/ ) that implements some of the algorithms and examples discussed here.

We present a comparison of parallaxes and radii from asteroseismology and Gaia DR1 (TGAS) for 2200 Kepler stars spanning from the main sequence to the red giant branch. We show that previously identified offsets between TGAS parallaxes and distances derived from asteroseismology and eclipsing binaries have likely been overestimated for parallaxes <~ 5-10 mas (~ 90-98% of the TGAS sample). The observed differences in our sample can furthermore be partially compensated by adopting a hotter Teff scale (such as the infrared flux method) instead of spectroscopic temperatures for dwarfs and subgiants. Residual systematic differences are at the ~ 2% level in parallax across three orders of magnitude. We use TGAS parallaxes to empirically demonstrate that asteroseismic radii are accurate to ~ 5% or better for stars between ~ 0.8-8 Rsun. We find no significant offset for main-sequence (<~ 1.5 Rsun) and low-luminosity RGB stars (~ 3-8 Rsun), but seismic radii appear to be systematically underestimated by ~5% for subgiants (~ 1.5-3 Rsun). We find no systematic errors as a function of metallicity between [Fe/H] ~ -0.8 to +0.4 dex, and show tentative evidence that corrections to the scaling relation for the large frequency separation (Dnu) improve the agreement with TGAS for RGB stars. Finally, we demonstrate that beyond ~ 3 kpc asteroseismology will provide more precise distances than end-of-mission Gaia data, highlighting the synergy and complementary nature of Gaia and asteroseismology for studying galactic stellar populations.

The first transiting planetesimal orbiting a white dwarf was recently detected in K2 data of WD1145+017 and has been followed up intensively. The multiple, long, and variable transits suggest the transiting objects are dust clouds, probably produced by a disintegrating asteroid. In addition, the system contains circumstellar gas, evident by broad absorption lines, mostly in the u'-band, and a dust disc, indicated by an infrared excess. Here we present the first detection of a change in colour of WD1145+017 during transits, using simultaneous multi-band fast-photometry ULTRACAM measurements over the u'g'r'i'-bands. The observations reveal what appears to be 'bluing' during transits; transits are deeper in the redder bands, with a u'-r' colour difference of up to ~-0.05 mag. We explore various possible explanations for the bluing. 'Spectral' photometry obtained by integrating over bandpasses in the spectroscopic data in- and out-of-transit, compared to the photometric data, shows that the observed colour difference is most likely the result of reduced circumstellar absorption in the spectrum during transits. This indicates that the transiting objects and the gas share the same line-of-sight, and that the gas covers the white dwarf only partially, as would be expected if the gas, the transiting debris, and the dust emitting the infrared excess, are part of the same general disc structure (although possibly at different radii). In addition, we present the results of a week-long monitoring campaign of the system.

We characterize the selection function of RAVE using 2MASS as our underlying population, which we assume represents all stars which could have potentially been observed. We evaluate the completeness fraction as a function of position, magnitude, and color in two ways: first, on a field-by-field basis, and second, in equal-size areas on the sky. Then, we consider the effect of the RAVE stellar parameter pipeline on the final resulting catalogue, which in principle limits the parameter space over which our selection function is valid. Our final selection function is the product of the completeness fraction and the selection function of the pipeline. We then test if the application of the selection function introduces biases in the derived parameters. To do this, we compare a parent mock catalogue generated using Galaxia with a mock-RAVE catalogue where the selection function of RAVE has been applied. We conclude that for stars brighter than I = 12, between $4000 \rm K < T_{\rm eff} < 8000 \rm K$ and $0.5 < \rm{log}\,g < 5.0$, RAVE is kinematically and chemically unbiased with respect to expectations from Galaxia.

Here we report in depth reanalysis of a paper by Vats et al. (2001) [Astrophys. J. 548, L87] based on the measurements of differential rotation with altitude as a function of observing frequencies (as lower and higher frequencies indicate higher and lower heights, respectively) in the solar corona. The radial differential rotation of the solar corona is estimated from daily measurements of the disc-integrated solar radio flux at 11 frequencies: (275, 405, 670, 810, 925, 1080, 1215, 1350, 1620, 1755 MHz and 2800 MHz). We use the same data as were used in Vats et al. (2001), but instead of the 12th maxima of autocorrelograms used there, we use the 1st secondary maxima to derive the synodic rotation period. We estimate synodic rotation by Gaussian fit of the 1st secondary maxima. Vats et al. (2001) reported that the sidereal rotation period increases with increasing frequency. The variation found by them was from 23.6 to 24.15 days in this frequency range with a difference of only 0.55 days. The present study finds that sidereal rotation period increases with decreasing frequency. The variation range is from 24.4 to 22.5 days and difference is about three times larger (1.9 days). However, at 925 MHz both studies give similar rotation period. In Vats et al. (2001) the Pearson factor with trend line was 0.86 whereas present analysis obtained a ~ 0.97 Pearson factor with the trend line. Our study shows that the solar corona rotates slower at higher altitudes, which is in contradiction to the findings reported in Vats et al. (2001). Key words: Solar radio flux- flux modulation method- Gaussian fit- sidereal rotation period.

Early attempts to apply asteroseismology to study the Galaxy have already shown unexpected discrepancies for the mass distribution of stars between the Galactic models and the data; a result that is still unexplained. Here, we revisit the analysis of the asteroseismic sample of dwarf and subgiant stars observed by Kepler and investigate in detail the possible causes for the reported discrepancy. We investigate two models of the Milky Way based on stellar population synthesis, Galaxia and TRILEGAL. In agreement with previous results, we find that TRILEGAL predicts more massive stars compared to Galaxia, and that TRILEGAL predicts too many blue stars compared to 2MASS observations. Both models fail to match the distribution of the stellar sample in $(\log g,T_{\rm eff})$ space, pointing to inaccuracies in the models and/or the assumed selection function. When corrected for this mismatch in $(\log g,T_{\rm eff})$ space, the mass distribution calculated by Galaxia is broader and the mean is shifted toward lower masses compared to that of the observed stars. This behaviour is similar to what has been reported for the Kepler red giant sample. The shift between the mass distributions is equivalent to a change of 2\% in $\nu_{\rm max}$, which is within the current uncertainty in the $\nu_{\rm max}$ scaling relation. Applying corrections to the $\Delta \nu$ scaling relation predicted by the stellar models makes the observed mass distribution significantly narrower, but there is no change to the mean.

Deep optical photometric data on the NGC 7538 region were collected and combined with archival data sets from $Chandra$, 2MASS and {\it Spitzer} surveys in order to generate a new catalog of young stellar objects (YSOs) including those not showing IR excess emission. This new catalog is complete down to 0.8 M$_\odot$. The nature of the YSOs associated with the NGC 7538 region and their spatial distribution are used to study the star formation process and the resultant mass function (MF) in the region. Out of the 419 YSOs, $\sim$91\% have ages between 0.1 to 2.5 Myr and $\sim$86\% have masses between 0.5 to 3.5 M$_\odot$, as derived by spectral energy distribution fitting analysis. Around 24\%, 62\% and 2\% of these YSOs are classified to be the Class I, Class II and Class III sources, respectively. The X-ray activity in the Class I, Class II and Class III objects is not significantly different from each other. This result implies that the enhanced X-ray surface flux due to the increase in the rotation rate may be compensated by the decrease in the stellar surface area during the pre-main sequence evolution. Our analysis shows that the O3V type high mass star `IRS 6' might have triggered the formation of young low mass stars up to a radial distance of 3 pc. The MF shows a turn-off at around 1.5 M$_\odot$ and the value of its slope `$\Gamma$' in the mass range $1.5 <$M/M$_\odot < 6$ comes out to be $-1.76\pm0.24$, which is steeper than the Salpeter value.

We present a set of 87 RAVE stars with detected solar like oscillations, observed during Campaign 1 of the K2 mission (RAVE K2-C1 sample). This dataset provides a useful benchmark for testing the gravities provided in RAVE Data Release 4 (DR4), and is key for the calibration of the RAVE Data Release 5 (DR5). In the present work, we use two different pipelines, GAUFRE (Valentini et al. 2013) and Sp_Ace (Boeche & Grebel 2016), to determine atmospheric parameters and abundances by fixing log(g) to the seismic one. Our strategy ensures highly consistent values among all stellar parameters, leading to more accurate chemical abundances. A comparison of the chemical abundances obtained here with and without the use of seismic log(g) information has shown that an underestimated (overestimated) gravity leads to an underestimated (overestimated) elemental abundance (e.g. [Mg/H] is underestimated by ~0.25 dex when the gravity is underestimated by 0.5 dex). We then perform a comparison between the seismic gravities and the spectroscopic gravities presented in the RAVE DR4 catalogue, extracting a calibration for log(g) of RAVE giants in the colour interval 0.50<(J - Ks)<0.85. Finally, we show a comparison of the distances, temperatures, extinctions (and ages) derived here for our RAVE K2-C1 sample with those derived in RAVE DR4 and DR5.DR5 performs better than DR4 thanks to the seismic calibration, although discrepancies can still be important for objects for which the difference between DR4/DR5 and seismic gravities differ by more than ~0.5 dex. The method illustrated in this work will be used for analysing RAVE targets present in the other K2 campaigns, in the framework of Galactic Archaeology investigations.

Galah is an ongoing high-resolution spectroscopic survey with the goal of disentangling the formation history of the Milky Way, using the fossil remnants of disrupted star formation sites which are now dispersed around the Galaxy. It is targeting a randomly selected, magnitude limited ($V \leq 14$) sample of stars, with the goal of observing one million objects. To date, 300,000 spectra have been obtained. Not all of them are correctly processed by parameter estimation pipelines and we need to know about them. We present a semi-automated classification scheme which identifies different types of peculiar spectral morphologies, in an effort to discover and flag potentially problematic spectra and thus help to preserve the integrity of the survey's results. To this end we employ a recently developed dimensionality reduction technique t-SNE (t-distributed Stochastic Neighbour Embedding), which enables us to represent the complex spectral morphology in a two-dimensional projection map while still preserving the properties of the local neighbourhoods of spectra. We find that the majority (178,483) of the 209,533 Galah spectra considered in this study represents normal single stars, whereas 31,050 peculiar and problematic spectra with very diverse spectral features pertaining to 28,579 stars are distributed into 10 classification categories: Hot stars, Cool metal-poor giants, Molecular absorption bands, Binary stars, H$\alpha$/H$\beta$ emission, H$\alpha$/H$\beta$ emission superimposed on absorption, H$\alpha$/H$\beta$ P-Cygni, H$\alpha$/H$\beta$ inverted P-Cygni, Lithium absorption, and Problematic. Classified spectra with supplementary information are presented in the catalogue, indicating candidates for follow-up observations and population studies of the short-lived phases of stellar evolution.

NASA's K2 mission is observing tens of thousands of stars along the ecliptic, providing data suitable for large scale asteroseismic analyses to inform galactic archaeology studies. Its first campaign covered a field near the north galactic cap, a region never covered before by large asteroseismic-ensemble investigations, and was therefore of particular interest for exploring this part of our Galaxy. Here we report the asteroseismic analysis of all stars selected by the K2 Galactic Archaeology Program during the mission's "North Galactic Cap" campaign 1. Our consolidated analysis uses six independent methods to measure the global seismic properties, in particular the large frequency separation, and the frequency of maximum power. From the full target sample of 8630 stars we find about 1200 oscillating red giants, a number comparable with estimates from galactic synthesis modeling. Thus, as a valuable by-product we find roughly 7500 stars to be dwarfs, which provide a sample well suited for galactic exoplanet occurrence studies because they originate from our simple and easily reproducible selection function. In addition, to facilitate the full potential of the data set for galactic archaeology we assess the detection completeness of our sample of oscillating red giants. We find the sample is at least near complete for stars with 40 < numax/microHz < 270, and numax_detec < 2.6*1e6 * 2e-Kp microHz. There is a detection bias against helium core burning stars with numax ~ 30 microHz, affecting the number of measurements of DeltaNu and possibly also numax. Although we can detect oscillations down to Kp = 15, our campaign 1 sample lacks enough faint giants to assess the detection completeness for stars fainter than Kp ~ 14.5.

We study the topological susceptibility in 2+1 flavor QCD above the chiral crossover transition temperature using Highly Improved Staggered Quark action and several lattice spacings, corresponding to temporal extent of the lattice, $N_\tau=6,8,10$ and $12$. We observe very distinct temperature dependences of the topological susceptibility in the ranges above and below $250$ MeV. While for temperatures above $250$ MeV, the dependence is found to be consistent with dilute instanton gas approximation, at lower temperatures the fall-off of topological susceptibility is milder. We discuss the consequence of our results for cosmology wherein we estimate the bounds on the axion decay constant and the oscillation temperature if indeed the QCD axion is a possible dark matter candidate.

Data Release 5 (DR5) of the Radial Velocity Experiment (RAVE) is the fifth data release from a magnitude-limited (9< I < 12) survey of stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra ($R\sim7500$) covering the Ca-triplet region (8410-8795\AA) span the complete time frame from the start of RAVE observations in 2003 to their completion in 2013. Radial velocities from 520,781 spectra of 457,588 unique stars are presented, of which 255,922 stellar observations have parallaxes and proper motions from the Tycho-Gaia astrometric solution (TGAS) in Gaia DR1. For our main DR5 catalog, stellar parameters (effective temperature, surface gravity, and overall metallicity) are computed using the RAVE DR4 stellar pipeline, but calibrated using recent K2 Campaign 1 seismic gravities and Gaia benchmark stars, as well as results obtained from high-resolution studies. Also included are temperatures from the Infrared Flux Method, and we provide a catalogue of red giant stars in the dereddened color $(J-Ks)_0$ interval (0.50,0.85) for which the gravities were calibrated based only on seismology. Further data products for sub-samples of the RAVE stars include individual abundances for Mg, Al, Si, Ca, Ti, Fe, and Ni, and distances found using isochrones. Each RAVE spectrum is complemented by an error spectrum, which has been used to determine uncertainties on the parameters. The data can be accessed via the RAVE Web site or the Vizier database.

We present a quantitative measurement of the amount of clustering present in the inner $\sim30$ kpc of the stellar halo of the Andromeda galaxy (M31). For this we analyse the angular positions and radial velocities of the carefully selected Planetary Nebulae (PNe) in the M31 stellar halo. We study the cumulative distribution of pair-wise distances in angular position and line of sight velocity space, and find that the M31 stellar halo contains substantially more stars in the form of close pairs as compared to that of a featureless smooth halo. In comparison to a smoothed/scrambled distribution we estimate that the clustering excess in the M31 inner halo is roughly $40\%$ at maximum and on average $\sim 20\%$. Importantly, comparing against the 11 stellar halo models of \cite{2005ApJ...635..931B}, which were simulated within the context of the $\Lambda{\rm CDM}$ cosmological paradigm, we find that the amount of substructures in the M31 stellar halo closely resembles that of a typical $\Lambda{\rm CDM}$ halo.

We have recorded three lunar occultations of Aldebaran (AlfTau) at different telescopes and using various band-passes, from the ultraviolet to the far red. The data have been analyzed using both model-dependent and model-independent methods. The derived uniform-disc angular diameter values have been converted to limb-darkened values using model atmosphere relations, and are found in broad agreement among themselves and with previous literature values. The limb-darkened diameter is about 20.3 milliarcseconds on average . However, we have found indications that the photospheric brightness profile of Aldebaran may have not been symmetric, a finding already reported by other authors for this and for similar late-type stars. At the sampling scale of our brightness profile, between one and two milliarcsecond, the uniform and limb-darkened disc models may not be a good description for Aldebaran. The asymmetries appear to differ with wavelength and over the 137 days time span of our measurements. Surface spots appear as a likely explanation for the differences between observations and the models.

The Galactic Archaeology with HERMES (GALAH) Survey is a massive observational project to trace the Milky Way's history of star formation, chemical enrichment, stellar migration and minor mergers. Using high-resolution (R$\simeq$28,000) spectra taken with the High Efficiency and Resolution Multi-Element Spectrograph (HERMES) instrument at the Anglo-Australian Telescope (AAT), GALAH will determine stellar parameters and abundances of up to 29 elements for up to one million stars. Selecting targets from a colour-unbiased catalogue built from 2MASS, APASS and UCAC4 data, we expect to observe dwarfs at 0.3 to 3 kpc and giants at 1 to 10 kpc. This enables a thorough local chemical inventory of the Galactic thin and thick disks, and also captures smaller samples of the bulge and halo. In this paper we present the plan, process and progress as of early 2016 for GALAH survey observations. In our first two years of survey observing we have accumulated the largest high-quality spectroscopic data set at this resolution, over 200,000 stars. We also present the first public GALAH data catalogue: stellar parameters (Teff, log(g), [Fe/H], [alpha/Fe]), radial velocity, distance modulus and reddening for 10680 observations of 9860 Tycho-2 stars that may be included in the first Gaia data release.

Assuming the Central Limit Theorem, experimental uncertainties in any data set are expected to follow the Gaussian distribution with zero mean. We propose an elegant method based on Kolmogorov-Smirnov statistic to test the above; and apply it on the measurement of Hubble constant which determines the expansion rate of the Universe. The measurements were made using Hubble Space Telescope. Our analysis shows that the uncertainties in the above measurement are non-Gaussian.

We present the data reduction procedures being used by the GALAH survey, carried out with the HERMES fibre-fed, multi-object spectrograph on the 3.9~m Anglo-Australian Telescope. GALAH is a unique survey, targeting 1 million stars brighter than magnitude V=14 at a resolution of 28,000 with a goal to measure the abundances of 29 elements. Such a large number of high resolution spectra necessitates the development of a reduction pipeline optimized for speed, accuracy, and consistency. We outline the design and structure of the Iraf-based reduction pipeline that we developed, specifically for GALAH, to produce fully calibrated spectra aimed for subsequent stellar atmospheric parameter estimation. The pipeline takes advantage of existing Iraf routines and other readily available software so as to be simple to maintain, testable and reliable. A radial velocity and stellar atmospheric parameter estimator code is also presented, which is used for further data analysis and yields a useful verification of the reduction quality. We have used this estimator to quantify the data quality of GALAH for fibre cross-talk level ($\lesssim0.5$%) and scattered light ($\sim5$ counts in a typical 20 minutes exposure), resolution across the field, sky spectrum properties, wavelength solution reliability (better than $1$ $\mathrm{km\ s^{-1}}$ accuracy) and radial velocity precision.

The nature of random errors in any data set that is Gaussian is a well established fact according to the Central Limit Theorem. Supernovae type Ia data have played a crucial role in major discoveries in cosmology. Unlike in laboratory experiments, astronomical measurements can not be performed in controlled situations. Thus, errors in astronomical data can be more severe in terms of systematics and non-Gaussianity compared to those of laboratory experiments. In this paper, we use the Kolmogorov-Smirnov statistic to test non-Gaussianity in high-z supernovae data. We apply this statistic to four data sets, i.e., Gold data(2004), Gold data(2007), Union2 catalogue and the Union2.1 data set for our analysis. Our results shows that in all four data sets the errors are consistent with the Gaussian distribution.

We present a multiwavelength investigation of the young stellar population and star formation activities around the \hii region Sharpless 311. Using our deep near-infrared observations and archival {\it Spitzer}-IRAC observations, we have detected a total of 125 young stellar objects (YSOs) in an area of $\sim$86 arcmin$^2$. The YSOs sample include 8 Class I and 117 Class II candidate YSOs. The mass completeness of the identified YSOs sample is estimated to be 1.0 \msun. The ages and masses of the majority of the candidate YSOs are estimated to be in the range of $\sim$0.1$-$5 Myr and $\sim$0.3$-$6 \msun, respectively. The 8 \mum image of S311 displays an approximately spherical cavity around the ionizing source which is possibly created due to the expansion of the \hii region. The spatial distribution of the candidate YSOs reveals that a significant number of them are distributed systematically along the 8 $\mu$m emission with a majority clustered around the eastern border of the \hii region. Four clumps/compact \hii regions are detected in the radio continuum observations at 1280 MHz, which might have been formed during the expansion of the \hii region. The estimated dynamical age of the region, main-sequence lifetime of the ionizing source, the spatial distribution and ages of the candidate YSOs indicate triggered star formation in the complex.

The K2 Mission uses the Kepler spacecraft to obtain high-precision photometry over ~80 day campaigns in the ecliptic plane. The Ecliptic Plane Input Catalog (EPIC) provides coordinates, photometry and kinematics based on a federation of all-sky catalogs to support target selection and target management for the K2 mission. We describe the construction of the EPIC, as well as modifications and shortcomings of the catalog. Kepler magnitudes (Kp) are shown to be accurate to ~0.1 mag for the Kepler field, and the EPIC is typically complete to Kp~17 (Kp~19 for campaigns covered by SDSS). We furthermore classify 138,600 targets in Campaigns 1-8 (~88% of the full target sample) using colors, proper motions, spectroscopy, parallaxes, and galactic population synthesis models, with typical uncertainties for G-type stars of ~3% in Teff, ~0.3 dex in log(g), ~40% in radius, ~10% in mass, and ~40% in distance. Our results show that stars targeted by K2 are dominated by K-M dwarfs (~41% of all selected targets), F-G dwarfs (~36%) and K giants (~21%), consistent with key K2 science programs to search for transiting exoplanets and galactic archeology studies using oscillating red giants. However, we find a significant variation of the fraction of cool dwarfs with galactic latitude, indicating a target selection bias due to interstellar reddening and the increased contamination by giant stars near the galactic plane. We discuss possible systematic errors in the derived stellar properties, and differences to published classifications for K2 exoplanet host stars. The EPIC is hosted at the Mikulski Archive for Space Telescopes (MAST): http://archive.stsci.edu/k2/epic/search.php.

We present an analysis of the relative throughputs of the 3.9-metre Anglo-Australian Telescope's 2dF/HERMES system, based upon spectra acquired during the first two years of the GALAH survey. Averaged spectral fluxes of stars were compared to their photometry to determine the relative throughputs of fibres for a range of fibre position and atmospheric conditions. We find that overall the throughputs of the 771 usable fibres have been stable over the first two years of its operation. About 2.5 per cent of fibres have throughputs much lower than the average. There are also a number of yet unexplained variations between the HERMES bandpasses, and mechanically & optically linked fibre groups known as retractors or slitlets related to regions of the focal plane. These findings do not impact the science that HERMES will produce.

With current space-based missions it is now possible to obtain age-sensitive asteroseismic information for tens of thousands of red giants. This provides a promising opportunity to study the Galactic structure and evolution. We use asteroseismic data of red giants, observed by Kepler, to test the current theoretical framework of modelling the Galaxy based on population synthesis modeling and the use of asteroseismic scaling relations for giants. We use the open source code Galaxia to model the Milky Way and find the distribution of the masses predicted by Galaxia to be systematically offset with respect to the seismically-inferred observed masses. The Galactic model overestimates the number of low mass stars, and these stars are predominantly old and of low metallicity. Using corrections to the $\Delta \nu$ scaling relation suggested by stellar models (available for download) significantly reduces the disagreement between predicted and observed masses. For a few cases where non-seismic mass estimates are available, the corrections to $\Delta \nu$ also improve the agreement between seismic and non-seismic mass estimates. The disagreement between predictions of the Galactic model and the observations is most pronounced for stars with ${\rm [Fe/H]}<-0.5$ and ${\rm [Fe/H]}>0$ or for $T_{\rm eff}>4700$ K. Altering the star formation rate in order to suppress stars older than 10 Gyr improves the agreement for mass but leads to inconsistent color distributions. We also tested the predictions of the TRILEGAL Galactic model. However, unlike {\sl Galaxia}, it had difficulties in reproducing the photometric properties of the Kepler Input Catalog because it overestimates the number of blue stars. We conclude that either the scaling relations and/or the Galactic models need to be revised to reconcile predictions of theory with asteroseismic observations.

We carried out deep and wide-field near- and mid-infrared observations for a sample of 8 bright-rimmed clouds (BRCs). Supplemented with the $Spitzer$ archival data, we have identified and classified 44 to 433 young stellar objects (YSOs) associated with these BRCs. The Class I sources are generally located towards the places with higher extinction and are relatively closer to each other than the Class II sources, confirming that the young protostars are usually found in regions having denser molecular material. On the other hand the comparatively older population, Class II objects, are more randomly found throughout the regions, which can be due to their dynamical evolution. Using the minimal sampling tree analyses, we have extracted 13 stellar cores of 8 or more members, which contains 60\% of the total YSOs. The typical core is $\sim$0.6 pc in radii and somewhat elongated (aspect ratio of 1.45), of relatively low stellar density (surface density 60 pc$^{-2}$), consisting of a small (35) number of YSOs of relatively young sources (66\% Class I ), and partially embedded (median $A_K$ =1.1 mag). But the cores show a wide range in their mass distribution ($\sim$20 to 2400 M$_\odot$) with a median value of around 130 M$_\odot$. We have found the star formation efficiencies in the cores to be between 3\% and 30\% with an average of $\sim$14\%, which agree with the efficiencies needed to link the core mass function to the initial mass function. We also found a linear relation between the density of the clouds and the number of YSOs. The peaked nearest neighbor spacing distributions of the YSOs and the ratio of Jeans lengths to the YSOs separations indicates a significant degree of non-thermally driven fragmentation in these BRCs.

We report the identification of extended tidal debris potentially associated with the globular cluster NGC 3201, using the RAVE catalogue. We find the debris stars are located at a distance range of 1-7 kpc based on the forthcoming RAVE distance estimates. The derived space velocities and integrals of motion show interesting connections to NGC 3201, modulo uncertainties in the proper motions. Three stars, which are among the 4 most likely candidates for NGC 3201 tidal debris, are separated by 80 degrees on the sky yet are well matched by the 12 Gyr, [Fe/H] = -1.5 isochrone appropriate for the cluster. This is the first time tidal debris around this cluster has been reported over such a large spatial extent, with implications for the cluster$'$s origin and dynamical evolution.

We have selected a sample of 80 candidates for obscured radio-loud active galactic nuclei and presented their basic optical/near-infrared (NIR) properties in Paper 1. In this paper, we present both high-resolution radio continuum images for all of these sources and HI 21cm absorption spectroscopy for a few selected sources in this sample. A-configuration 4.9 and 8.5 GHz VLA continuum observations find that 52 sources are compact or have substantial compact components with size <0.5" and flux density >0.1 Jy at 4.9 GHz. The most compact 36 sources were then observed with the VLBA at 1.4 GHz. One definite and 10 candidate Compact Symmetric Objects (CSOs) are newly identified, a detection rate of CSOs ~3 times higher than the detection rate previously found in purely flux-limited samples. Based on possessing compact components with high flux densities, 60 of these sources are good candidates for absorption-line searches. Twenty seven sources were observed for HI 21cm absorption at their photometric or spectroscopic redshifts with only 6 detections made (one detection is tentative). However, five of these were from a small subset of six CSOs with pure galaxy optical/NIR spectra and for which accurate spectroscopic redshifts place the redshifted 21cm line in a RFI-free spectral window. It is likely that the presence of ubiquitous RFI and the absence of accurate spectroscopic redshifts preclude HI detections in similar sources (only one detection out of the remaining 22 sources observed, 14 of which have only photometric redshifts). Future searches for highly-redshifted HI and molecular absorption can easily find more distant CSOs among bright, blank field' radio sources but will be severely hampered by an inability to determine accurate spectroscopic redshifts for them due to their lack of rest-frame UV continuum.

The High Efficiency and Resolution Multi Element Spectrograph, HERMES, is a facility-class optical spectrograph for the Anglo-Australian Telescope (AAT). It is designed primarily for Galactic Archaeology, the first major attempt to create a detailed understanding of galaxy formation and evolution by studying the history of our own galaxy, the Milky Way. The goal of the GALAH survey is to reconstruct the mass assembly history of the Milky Way through a detailed chemical abundance study of one million stars. The spectrograph is based at the AAT and is fed by the existing 2dF robotic fiber positioning system. The spectrograph uses volume phase holographic gratings to achieve a spectral resolving power of 28,000 in standard mode and also provides a high-resolution mode ranging between 40,000 and 50,000 using a slit mask. The GALAH survey requires an SNR greater than 100 for a star brightness of V ?= 14 in an exposure time of one hour. The total spectral coverage of the four channels is about 100 nm between 370 and 1000 nm for up to 392 simultaneous targets within the 2-degree field of view. HERMES has been commissioned over three runs, during bright time in October, November, and December 2013, in parallel with the beginning of the GALAH pilot survey, which started in November 2013. We present the first-light results from the commissioning run and the beginning of the GALAH survey, including performance results such as throughput and resolution, as well as instrument reliability.

NASA's re-purposed Kepler mission -- dubbed K2 -- has brought new scientific opportunities that were not anticipated for the original Kepler mission. One science goal that makes optimal use of K2's capabilities, in particular its 360-degree ecliptic field of view, is galactic archaeology -- the study of the evolution of the Galaxy from the fossil stellar record. The thrust of this research is to exploit high-precision, time-resolved photometry from K2 in order to detect oscillations in red giant stars. This asteroseismic information can provide estimates of stellar radius (hence distance), mass and age of vast numbers of stars across the Galaxy. Here we present the initial analysis of a subset of red giants, observed towards the North Galactic Gap, during the mission's first full science campaign. We investigate the feasibility of using K2 data for detecting oscillations in red giants that span a range in apparent magnitude and evolutionary state (hence intrinsic luminosity). We demonstrate that oscillations are detectable for essentially all cool giants within the $\log g$ range $\sim 1.9-3.2$. Our detection is complete down to $\mathit{Kp}\sim 14.5$, which results in a seismic sample with little or no detection bias. This sample is ideally suited to stellar population studies that seek to investigate potential shortcomings of contemporary Galaxy models.

The GALAH survey targets one million stars in the southern hemisphere down to a limiting magnitude of V = 14 at the Anglo- Australian Telescope. The project aims to measure up to 30 elemental abundances and radial velocities (~1 km/s accuracy) for each star at a resolution of R = 28000. These elements fall into 8 independent groups (e.g. alpha, Fe peak, r-process). For all stars, Gaia will provide distances to 1% and transverse velocities to 1 km/s or better, giving us a 14D set of parameters for each star, i.e. 6D phase space and 8D abundance space. There are many scientic applications but here we focus on the prospect of chemically tagging the thick disk and making a direct measurement of how stellar migration evolves with cosmic time.

The GALAH survey is a large high-resolution spectroscopic survey using the newly commissioned HERMES spectrograph on the Anglo-Australian Telescope. The HERMES spectrograph provides high-resolution (R ~28,000) spectra in four passbands for 392 stars simultaneously over a 2 degree field of view. The goal of the survey is to unravel the formation and evolutionary history of the Milky Way, using fossil remnants of ancient star formation events which have been disrupted and are now dispersed throughout the Galaxy. Chemical tagging seeks to identify such dispersed remnants solely from their common and unique chemical signatures; these groups are unidentifiable from their spatial, photometric or kinematic properties. To carry out chemical tagging, the GALAH survey will acquire spectra for a million stars down to V~14. The HERMES spectra of FGK stars contain absorption lines from 29 elements including light proton-capture elements, alpha-elements, odd-Z elements, iron-peak elements and n-capture elements from the light and heavy s-process and the r-process. This paper describes the motivation and planned execution of the GALAH survey, and presents some results on the first-light performance of HERMES.

Bulk flows of galaxies moving with respect to the cosmic microwave background are well established observationally and seen in the most recent LCDM simulations. With the aid of an idealised Gadget-2 simulation, we show that void asymmetries in the cosmic web can exacerbate local bulk flows of galaxies. The Cosmicflows-2 survey, which has mapped in detail the 3D structure of the Local Universe, reveals that the Local Group resides in a "local sheet" of galaxies that borders a "local void" with a diameter of about 40 Mpc. The void is emptying out at a rate of 16 km/s/Mpc. In a co-moving frame, the Local Sheet is found to be moving away from the Local Void at ~ 260 km/s. Our model shows how asymmetric collapse due to unbalanced voids on either side of a developing sheet or wall can lead to a systematic movement of the sheet. We conjectured that asymmetries could lead to a large-scale separation of dark matter and baryons, thereby driving a dependence of galaxy properties with environment, but we do not find any evidence for this effect.

The diffuse interstellar bands (DIBs) are absorption lines observed in visual and near infrared spectra of stars. Understanding their origin in the interstellar medium (ISM) is one of the oldest problems in astronomical spectroscopy, as DIBs ahave been known since 1922. In a completely new approach to understanding DIBs, we combined information from nearly 500,000 stellar spectra obtained by the massive spectroscopic survey RAVE (Radial Velocity Experiment) to produce the first pseudo three-dimensional map of the strength of the DIB at 8620~{\AA}ngstroms covering the nearest 3 kiloparsecs from the Sun, and show that it follows our independently constructed spatial distribution of extinction by interstellar dust along the Galactic plane. Despite having a similar distribution in the Galactic plane, the DIB 8620 carrier has a significantly larger vertical scale height than the dust. Even if one DIB may not represent the general DIB population, our observations outline the future direction of DIB research.

Halo stars orbit within the potential of the Milky Way and hence their kinematics can be used to understand the underlying mass distribution. However, the inferred mass distribution depends sensitively upon assumptions made on the density and the velocity anisotropy profiles of the tracers. Also, there is a degeneracy between the parameters of the halo and that of the disk or bulge. Here, we decompose the Galaxy into bulge, disk and dark matter halo and then model the kinematic data of the halo BHB and K-giants from the SEGUE. Additionally, we use the gas terminal velocity curve and the Sgr A$^*$ proper motion. With $R_\odot = 8.5$kpc, our study reveals that the density of the stellar halo has a break at $17.2^{+1.1}_{-1.0}$ kpc, and an exponential cut-off in the outer parts starting at $97.7^{+15.6}_{-15.8}$kpc. Also, we find the velocity anisotropy is radially biased with $\beta_s= 0.4\pm{0.2}$ in the outer halo. We measure halo virial mass $M_{\text{vir}} = 0.80^{+0.31}_{-0.16} \times 10^{12} M_{\odot}$, concentration $c=21.1^{+14.8}_{-8.3}$, disk mass of $0.95^{+0.24}_{-0.30}\times10^{11} M_{\odot}$, disk scale length of $4.9^{+0.4}_{-0.4}$ kpc and bulge mass of $0.91^{+0.31}_{-0.38} \times 10^{10} M_{\odot}$. The mass of halo is found to be small and this has important consequences. The giant stars reveal that the outermost halo stars have low velocity dispersion interestingly suggesting a truncation of the stellar halo density rather than a small overall mass of the Galaxy. Our estimates of local escape velocity $v_{\rm esc} = 550.9^{+32.4}_{-22.1}$ kms$^{-1}$ and dark matter density $\rho^{\rm DM}_{\odot} = 0.0088^{+0.0024}_{-0.0018} M_{\odot} {\rm pc^{-3}} $ ($0.35^{+0.08}_{-0.07}$ GeV cm$^{-3}$) are in good agreement with recent estimates. Some of the above estimates are depended on the adopted value of $R_\odot$ and of outer power-law index of the tracer number density.

We provide new constraints on the chemo-dynamical models of the Milky Way by measuring the radial and vertical chemical gradients for the elements Mg, Al, Si, Ti, and Fe in the Galactic disc and the gradient variations as a function of the distance from the Galactic plane ($Z$). We selected a sample of giant stars from the RAVE database using the gravity criterium 1.7$<$log g$<$2.8. We created a RAVE mock sample with the Galaxia code based on the Besan\c con model and selected a corresponding mock sample to compare the model with the observed data. We measured the radial gradients and the vertical gradients as a function of the distance from the Galactic plane $Z$ to study their variation across the Galactic disc. The RAVE sample exhibits a negative radial gradient of $d[Fe/H]/dR=-0.054$ dex kpc$^{-1}$ close to the Galactic plane ($|Z|<0.4$ kpc) that becomes flatter for larger $|Z|$. Other elements follow the same trend although with some variations from element to element. The mock sample has radial gradients in fair agreement with the observed data. The variation of the gradients with $Z$ shows that the Fe radial gradient of the RAVE sample has little change in the range $|Z|\lesssim0.6$ kpc and then flattens. The iron vertical gradient of the RAVE sample is slightly negative close to the Galactic plane and steepens with $|Z|$. The mock sample exhibits an iron vertical gradient that is always steeper than the RAVE sample. The mock sample also shows an excess of metal-poor stars in the [Fe/H] distributions with respect to the observed data. These discrepancies can be reduced by decreasing the number of thick disc stars and increasing their average metallicity in the Besan\c con model.

VISTA Variables in the V\'ia L\'actea (VVV) is one of six ESO Public Surveys using the 4 meter Visible and Infrared Survey Telescope for Astronomy (VISTA). The VVV survey covers the Milky Way bulge and an adjacent section of the disk, and one of the principal objectives is to search for new star clusters within previously unreachable obscured parts of the Galaxy. The primary motivation behind this work is to discover and analyze obscured star clusters in the direction of the inner Galactic disk and bulge. Regions of the inner disk and bulge covered by the VVV survey were visually inspected using composite JHKs color images to select new cluster candidates on the basis of apparent overdensities. DR1, DR2, CASU, and PSF photometry of 10x10 arcmin fields centered on each candidate cluster were used to construct color-magnitude and color-color diagrams. Follow-up spectroscopy of the brightest members of several cluster candidates was obtained in order to clarify their nature. We report the discovery of 58 new infrared cluster candidates. Fundamental parameters such as age, distance, and metallicity were determined for 20 of the most populous clusters.

The VISTA Variables in the V\'ia L\'actea (VVV) survey is one of six public ESO surveys, and is now in its 4th year of observing. Although far from being complete, the VVV survey has already delivered many results, some directly connected to the intended science goals (detection of variables stars, microlensing events, new star clusters), others concerning more exotic objects, e.g. novae. Now, at the end of the fourth observing period, and comprising roughly 50% of the proposed observations, the actual status of the survey, as well some of the results based on the VVV data, are presented.

The stellar content in and around Sh2-311 region have been studied using the deep optical observations as well as near-infrared (NIR) data from 2MASS. The region contains three clusters, viz. NGC 2467, Haffner 18 and Haffner 19. We have made an attempt to distinguish the stellar content of these individual regions as well as to re-determine their fundamental parameters such as distance, reddening, age, onto the basis of a new and more extended optical and infrared photometric data set. NGC 2467 and Haffner 19 are found to be located in the Perseus arm at the distances of 5.0 $\pm$ 0.4 kpc and 5.7 $\pm$ 0.4 kpc, respectively, whereas Haffner 18 is located at the distance of 11.2 $\pm$ 1.0 kpc. The clusters NGC 2467 and Haffner 19 might have formed from the same molecular cloud, whereas the cluster Haffner 18 is located in the outer galactic arm, i.e. the Norma-Cygnus arm. We identify 8 class II young stellar objects (YSOs) using the NIR $(J - H)/(H - K)$ two colour diagram. We have estimated the age and mass of the YSOs identified in the present work and those by Snider et al. (2009) using the $V/(V - I)$ colour-magnitude diagram. The estimated ages and mass range of the majority of the YSOs are $\lesssim$1 Myr and $\sim$0.4 - 3.5 \msun, respectively, indicating that these sources could be T-Tauri stars or their siblings. Spatial distribution of the YSOs shows that some of the YSOs are distributed around the H II region Sh2-311, suggesting a triggered star formation at its periphery.

We obtained deep $UBVRI$ H$\alpha$ photometric data of the field situated to the west of the main Carina nebula and centered on WR~22. Medium-resolution optical spectroscopy of a subsample of X-ray selected objects along with archival data sets from $Chandra$, $XMM-Newton$ and 2MASS surveys were used for the present study. Different sets of color-color and color-magnitude diagrams are used to determine reddening for the region and to identify young stellar objects (YSOs) and estimate their age and mass. Our spectroscopic results indicate that the majority of the X-ray sources are late spectral type stars. The region shows a large amount of differential reddening with minimum and maximum values of $E(B-V)$ as 0.25 and 1.1 mag, respectively. Our analysis reveals that the total-to-selective absorption ratio $R \rm_{V}$ is $\sim$3.7 $\pm$ 0.1, suggesting an abnormal grain size in the observed region. We identified 467 YSOs and studied their characteristics. The ages and masses of the 241 optically identified YSOs range from $\sim$0.1 to 10 Myr and $\sim$0.3 to 4.8M$_\odot$, respectively. However, the majority of them are younger than 1 Myr and have masses below 2 M$_\odot$. The high mass star WR 22 does not seem to have contributed to the formation of YSOs in the CrW region. The initial mass function slope, $\Gamma$, in this region is found to be -1.13 $\pm$ 0.20 in the mass range of 0.5 < M/M$_\odot$ < 4.8. The $K$-band luminosity function slope ($\alpha$) is also estimated as 0.31 $\pm$ 0.01. We also performed minimum spanning tree analysis of the YSOs in this region, which reveals that there are at least ten YSO cores associated with the molecular cloud, and that leads to an average core radius of 0.43 pc and a median branch length of 0.28 pc.

We investigate the kinematic parameters of the Milky Way disc using the RAVE and GCS stellar surveys. We do this by fitting a kinematic model to the data taking the selection function of the data into account. For stars in the GCS we use all phase-space coordinates, but for RAVE stars we use only $(l,b,v_{\rm los})$. Using MCMC technique, we investigate the full posterior distributions of the parameters given the data. We investigate the `age-velocity dispersion' relation for the three kinematic components ($\sigma_R,\sigma_{\phi},\sigma_z$), the radial dependence of the velocity dispersions, the Solar peculiar motion ($U_{\odot},V_{\odot}, W_{\odot} $), the circular speed $\Theta_0$ at the Sun and the fall of mean azimuthal motion with height above the mid-plane. We confirm that the Besan\c{c}on-style Gaussian model accurately fits the GCS data, but fails to match the details of the more spatially extended RAVE survey. In particular, the Shu distribution function (DF) handles non-circular orbits more accurately and provides a better fit to the kinematic data. The Gaussian distribution function not only fits the data poorly but systematically underestimates the fall of velocity dispersion with radius. We find that correlations exist between a number of parameters, which highlights the importance of doing joint fits. The large size of the RAVE survey, allows us to get precise values for most parameters. However, large systematic uncertainties remain, especially in $V_{\odot}$ and $\Theta_0$. We find that, for an extended sample of stars, $\Theta_0$ is underestimated by as much as $10\%$ if the vertical dependence of the mean azimuthal motion is neglected. Using a simple model for vertical dependence of kinematics, we find that it is possible to match the Sgr A* proper motion without any need for $V_{\odot}$ being larger than that estimated locally by surveys like GCS.

In this paper we present time series photometry of 104 variable stars in the cluster region NGC 1893. The association of the present variable candidates to the cluster NGC 1893 has been determined by using $(U-B)/(B-V)$ and $(J-H)/(H-K)$ two colour diagrams, and $V/(V-I)$ colour magnitude diagram. Forty five stars are found to be main-sequence variables and these could be B-type variable stars associated with the cluster. We classified these objects as $\beta$ Cep, slowly pulsating B stars and new class variables as discussed by Mowlavi et al. (2013). These variable candidates show $\sim$0.005 to $\sim$0.02 mag brightness variations with periods of $<$ 1.0 d. Seventeen new class variables are located in the $H-R$ diagram between the slowly pulsating B stars and $\delta$ Scuti variables. Pulsation could be one of the causes for periodic brightness variations in these stars. The X-ray emission of present main-sequence variables associated with the cluster lies in the saturated region of X-ray luminosity versus period diagram and follows the general trend by Pizzolato et al. (2003).

Shortly after the seminal paper {\sl "Self-Organized Criticality: An explanation of 1/f noise"} by Bak, Tang, and Wiesenfeld (1987), the idea has been applied to solar physics, in {\sl "Avalanches and the Distribution of Solar Flares"} by Lu and Hamilton (1991). In the following years, an inspiring cross-fertilization from complexity theory to solar and astrophysics took place, where the SOC concept was initially applied to solar flares, stellar flares, and magnetospheric substorms, and later extended to the radiation belt, the heliosphere, lunar craters, the asteroid belt, the Saturn ring, pulsar glitches, soft X-ray repeaters, blazars, black-hole objects, cosmic rays, and boson clouds. The application of SOC concepts has been performed by numerical cellular automaton simulations, by analytical calculations of statistical (powerlaw-like) distributions based on physical scaling laws, and by observational tests of theoretically predicted size distributions and waiting time distributions. Attempts have been undertaken to import physical models into the numerical SOC toy models, such as the discretization of magneto-hydrodynamics (MHD) processes. The novel applications stimulated also vigorous debates about the discrimination between SOC models, SOC-like, and non-SOC processes, such as phase transitions, turbulence, random-walk diffusion, percolation, branching processes, network theory, chaos theory, fractality, multi-scale, and other complexity phenomena. We review SOC studies from the last 25 years and highlight new trends, open questions, and future challenges, as discussed during two recent ISSI workshops on this theme.

Probability density functions are determined from new stellar parameters for the distance moduli of stars for which the RAdial Velocity Experiment (RAVE) has obtained spectra with S/N>=10. Single-Gaussian fits to the pdf in distance modulus suffice for roughly half the stars, with most of the other half having satisfactory two-Gaussian representations. As expected, early-type stars rarely require more than one Gaussian. The expectation value of distance is larger than the distance implied by the expectation of distance modulus; the latter is itself larger than the distance implied by the expectation value of the parallax. Our parallaxes of Hipparcos stars agree well with the values measured by Hipparcos, so the expectation of parallax is the most reliable distance indicator. The latter are improved by taking extinction into account. The effective temperature absolute-magnitude diagram of our stars is significantly improved when these pdfs are used to make the diagram. We use the method of kinematic corrections devised by Schoenrich, Binney & Asplund to check for systematic errors for general stars and confirm that the most reliable distance indicator is the expectation of parallax. For cool dwarfs and low-gravity giants <pi> tends to be larger than the true distance by up to 30 percent. The most satisfactory distances are for dwarfs hotter than 5500 K. We compare our distances to stars in 13 open clusters with cluster distances from the literature and find excellent agreement for the dwarfs and indications that we are over-estimating distances to giants, especially in young clusters.

Perturbation of FRW spacetime is carried out in NP formalism. The equation governing the scalar, vector and tensor modes take on a very simple and transparent form. All of them can be combined in one master equation for all helicities. The solutions for the closed, flat and open FRW are analytic continuation of the same functions, so only the solutions in the closed model are described. The scalar equation is the same as that of the conformally coupled massless Klein-Gordon field, the vectorial ones are the same as Maxwell equations, and the tensorial ones are for spin-2 fields. The corresponding eigen-functions are all determined, and in particular, the Green's function and the Lienard-Wiechert type potential also solved. These solutions reduce to the familiar form in flat space.

We explore a second order Hamiltonian vertical resonance model for X-shaped or peanut-shaped galactic bulges. The X-shape is caused by the 2:1 vertical Lindblad resonance with the bar, with two vertical oscillation periods per orbital period in the bar frame. We examine N-body simulations and find that due to the bar slowing down and disk thickening during bar buckling, the resonance and associated peanut-shape moves outward. The peanut-shape is consistent with the location of the vertical resonance, independent of whether the bar buckled or not. We estimate the resonance width from the potential m=4 Fourier component and find that the resonance is narrow, affecting orbits over a narrow range in the angular momentum distribution, dL/L ~ 0.05. As the resonance moves outward, stars originally in the mid plane are forced out of the mid plane into orbits just within the resonance separatrix. The height of the separatrix orbits, estimated from the Hamiltonian model, is approximately consistent with the peanut-shape height. The X-shape is comprised of stars in the vicinity of the resonance separatrix. The velocity distributions from the simulations illustrate that low inclination orbits are depleted within resonance. Within resonance, the vertical velocity distribution is broad, consistent with resonant heating caused by the passage of the resonance through the disk. In the Milky Way bulge we relate the azimuthally averaged mid-plane mass density near the vertical resonance to the rotation curve and bar pattern speed. At an estimated vertical resonance galactocentric radius of ~1.3 kpc, we confirm a mid-plane density of ~5x10^8 Msol/kpc^3, consistent with recently estimated mass distributions. We find that the rotation curve, bar pattern speed, 2:1 vertical resonance location, X-shape tips, and mid-plane mass density, are all self-consistent in the Milky Way galaxy bulge.

The RAVE survey, combined with proper motions and distance estimates, can be used to study in detail stellar kinematics in the extended solar neighbourhood (solar suburb). Using the red clump, we examine the mean velocity components in 3D between an R of 6 and 10 kpc and a Z of -2 to 2 kpc, concentrating on North-South differences. Simple parametric fits to the R, Z trends for VPHI and the velocity dispersions are presented. We confirm the recently discovered gradient in mean Galactocentric radial velocity, VR, finding that the gradient is more marked below the plane, with a Z gradient also present. The vertical velocity, VZ, also shows clear structure, with indications of a rarefaction-compression pattern, suggestive of wave-like behaviour. We perform a rigorous error analysis, tracing sources of both systematic and random errors. We confirm the North-South differences in VR and VZ along the line-of-sight, with the VR estimated independent of the proper motions. The complex three-dimensional structure of velocity space presents challenges for future modelling of the Galactic disk, with the Galactic bar, spiral arms and excitation of wave-like structures all probably playing a role.

We present the stellar atmospheric parameters (effective temperature, surface gravity, overall metallicity), radial velocities, individual abundances and distances determined for 425 561 stars, which constitute the fourth public data release of the RAdial Velocity Experiment (RAVE). The stellar atmospheric parameters are computed using a new pipeline, based on the algorithms of MATISSE and DEGAS. The spectral degeneracies and the 2MASS photometric information are now better taken into consideration, improving the parameter determination compared to the previous RAVE data releases. The individual abundances for six elements (magnesium, aluminum, silicon, titanium, iron and nickel) are also given, based on a special-purpose pipeline which is also improved compared to that available for the RAVE DR3 and Chemical DR1 data releases. Together with photometric information and proper motions, these data can be retrieved from the RAVE collaboration website and the Vizier database.

Aim: We aim at measuring the chemical gradients of the elements Mg, Al, Si, and Fe along the Galactic radius to provide new constraints on the chemical evolution models of the Galaxy and Galaxy models such as the Besancon model. Methods: We analysed three different samples selected from three independent datasets: a sample of 19,962 dwarf stars selected from the RAVE database, a sample of 10,616 dwarf stars selected from the Geneva-Copenhagen Survey (GCS) dataset, and a mock sample (equivalent to the RAVE sample) created by using the GALAXIA code, which is based on the Besancon model. We measured the chemical gradients as functions of the guiding radius (Rg) at different distances from the Galactic plane reached by the stars along their orbit (Zmax). Results: The chemical gradients of the RAVE and GCS samples are negative and show consistent trends, although they are not equal: at Zmax<0.4 kpc and 4.5<Rg(kpc)<9.5, the iron gradient for the RAVE sample is d[Fe/H]/dRg=-0.065 dex kpc^{-1}, whereas for the GCS sample it is d[Fe/H]/dRg=-0.043 dex kpc^{-1} with internal errors +-0.002 and +-0.004 dex kpc^{-1}, respectively. The gradients of the RAVE and GCS samples become flatter at larger Zmax. Conversely, the mock sample has a positive iron gradient of d[Fe/H]/dRg=+0.053+-0.003 dex kpc^{-1} at Zmax<0.4 kpc and remains positive at any Zmax. These positive and unrealistic values originate from the lack of correlation between metallicity and tangential velocity in the Besancon model. The discrepancies between the observational samples and the mock sample can be reduced by i) decreasing the density, ii) decreasing the vertical velocity, and iii) increasing the metallicity of the thick disc in the Besancon model.

Some solutions of the Maxwell equations with Dirac particles for the source in FRW spacetime are discussed. The Green's function of the equation for the radial component of the Maxwell fields, F_{r\eta} and F_{\theta\phi} is solved. Green's function is found to reduce to that of Minkowskian spacetime in the appropriate limit. Also, the Lienard-Wiechert type solution is derived. Also, the solutions with the Dirac particle current is also presented. It is found that the F_{r\eta} is composed of even angular momentum states while the odd states constitue F_{\theta\phi} .

An empirical formula for a Shu distribution function that reproduces a thin disc with exponential surface density to good accuracy is presented. The formula has two free parameters that specify the functional form of the velocity dispersion. Conventionally, this requires the use of an iterative algorithm to produce the correct solution, which is computationally taxing for applications like Markov Chain Monte Carlo (MCMC) model fitting. The formula has been shown to work for flat, rising and falling rotation curves. Application of this methodology to one of the Dehnen distribution functions is also shown. Finally, an extension of this formula to reproduce velocity dispersion profiles that are an exponential function of radius is also presented. Our empirical formula should greatly aid the efficient comparison of disc models with large stellar surveys or N-body simulations.

$UBVRI$ photometry of the five open clusters Czernik 4, Berkeley 7, NGC 2236, NGC 7226 and King 12 has been carried out using ARIES 104 cm telescope, Nainital. Fundamental cluster parameters such as foreground reddening $E(B-V)$, distance, and age have been derived by means of the observed two colour and colour-magnitude diagrams, coupled to comparisons with theoretical models. $E(B-V)$ values range from 0.55 to 0.74 mag, while ages derived for these clusters range from $\sim$10 to $\sim$500 Myr. We have also studied the spatial structure, mass function and mass segregation effects. The present study shows that evaporation of low mass stars from the halo of the clusters increases as they evolve.

We present a sample of ~5,000 RR Lyrae stars selected from the recalibrated LINEAR dataset and detected at heliocentric distances between 5 kpc and 30 kpc over ~8,000 deg^2 of sky. The coordinates and light curve properties, such as period and Oosterhoff type, are made publicly available. We find evidence for the Oosterhoff dichotomy among field RR Lyrae stars, with the ratio of the type II and I subsamples of about 1:4. The number density distribution of halo RRab stars as a function of galactocentric distance can be described as an oblate ellipsoid with the axis ratio q=0.63 and with either a single or a double power law with a power-law index in the range -2 to -3. Using a group-finding algorithm EnLink, we detected seven candidate halo groups, only one of which is statistically spurious. Three of these groups are near globular clusters (M53/NGC 5053, M3, M13), and one is near a known halo substructure (Virgo Stellar Stream); the remaining three groups do not seem to be near any known halo substructures or globular clusters, and seem to have a higher ratio of Oosterhoff type II to Oosterhoff type I RRab stars than what is found in the halo. The extended morphology and the position (outside the tidal radius) of some of the groups near globular clusters is suggestive of tidal streams possibly originating from globular clusters. Spectroscopic followup of detected halo groups is encouraged.

We simulate the evolution of a 10^9 Msun dark matter halo in a cosmological setting with an adaptive-mesh refinement code as an analogue to local low luminosity dwarf irregular and dwarf spheroidal galaxies. The primary goal of our study is to investigate the roles of reionization and supernova feedback in determining the star formation histories of low mass dwarf galaxies. We include a wide range of physical effects, including metal cooling, molecular hydrogen formation and cooling, photoionization and photodissociation from a metagalactic background, a simple prescription for self-shielding, star formation, and a simple model for supernova driven energetic feedback. We carry out simulations excluding each major effect in turn. We find that reionization is primarily responsible for expelling most of the gas in our simulations, but that supernova feedback is required to disperse the dense, cold gas in the core of the halo. Moreover, we show that the timing of reionization can produce an order of magnitude difference in the final stellar mass of the system. For our full physics run with reionization at z=9, we find a stellar mass of about 10^5 Msun at z=0, and a mass-to-light ratio within the half-light radius of approximately 130 Msun/Lsun, consistent with observed low-luminosity dwarfs. However, the resulting median stellar metallicity is 0.06 Zsun, considerably larger than observed systems. In addition, we find star formation is truncated between redshifts 4 and 7, at odds with the observed late time star formation in isolated dwarf systems but in agreement with Milky Way ultrafaint dwarf spheroidals. We investigate the efficacy of energetic feedback in our simple thermal-energy driven feedback scheme, and suggest that it may still suffer from excessive radiative losses, despite reaching stellar particle masses of about 100 Msun, and a comoving spatial resolution of 11 pc.

Aims: We study the relations between stellar kinematics and chemical abundances of a large sample of RAVE giants in search for selection criteria needed for disentangling different Galactic stellar populations. Methods: We select a sample of 2167 giant stars with signal-to-noise per spectral measurements above 75 from the RAVE chemical catalogue and follow the analysis performed by Gratton and colleagues on 150 subdwarf stars spectroscopically observed at high-resolution. We then use a larger sample of 9131 giants (with signal-to-noise above 60) to investigate the chemo-kinematical characteristics of our stars by grouping them into nine subsamples with common eccentricity ($e$) and maximum distance achieved above the Galactic plane ($Z_max$). Results: The RAVE kinematical and chemical data proved to be reliable by reproducing the results by Gratton et al. obtained with high-resolution spectroscopic data. Our analysis, based on the $e$-$Z_max$ plane combined with additional orbital parameters and chemical information, provides an alternative way of identifying different populations of stars. In addition to extracting canonical thick- and thin-disc samples, we find a group of stars in the Galactic plane ($Z_max<1$ kpc and 0.4 $< e < $0.6), which show homogeneous kinematics but differ in their chemical properties. We interpret this as a clear sign that some of these stars have experienced the effects of heating and/or radial migration, which have modified their original orbits. The accretion origin of such stars cannot be excluded.

Here we study the metallicity bias in the velocity dispersions, the derived quantity called anisotropy and the mean azimuthal velocity profiles of the Milky Way stellar halo using Blue Horizontal Branch (BHB) stars taken from SDSS/SEGUE survey. The comparatively metal-rich sample ([Fe/H]>-2) has prograde motion and is found to have an offset of 40 km/s in the mean azimuthal velocity with respect to a metal-poor sample ([Fe/H]<=-2) which has retrograde motion. The difference in rotation between the most metal-poor and most metal-rich population was found to be around 65 km/s. For galactocentric distances r<16 kpc, an offset in velocity dispersion profiles and anisotropy can also be seen. In the inner regions, the metal-poor population is in average tangential orbit; however, anisotropy is found to decrease monotonically with radius independent of metallicity. Beyond r = 16 kpc, both the metal-rich and the metal-poor samples are found to have tangential motion. The metallicity bias in the kinematics of the halo stars qualitatively supports the co-existence of at least two-components in the halo having different formation history e.g. in-situ formation and formation by accretion.

We present optical photometric and polarimetric observations of stars towards NGC 1931 with the aim to derive the cluster parameters such as distance, reddening, age and luminosity/mass function as well as to understand the dust properties and star formation in the region. The distance to the cluster is found to be 2.3$\pm$0.3 kpc and the reddening E(B-V) in the region is found to be variable. The stellar density contours reveal two clustering in the region. The observations suggest differing reddening law within the cluster region. Polarization efficiency of the dust grains towards the direction of the cluster is found to be less than that for the general diffuse interstellar medium (ISM). The slope of the mass function (-0.98$\pm$0.22) in the southern region in the mass range 0.8 \textless $M/M_{\sun}$ $\textless$ 9.8 is found to be shallower in comparison to that in the northern region (-1.26$\pm$0.23), which is comparable to the Salpeter value (-1.35). The K-band luminosity function (KLF) of the region is found to be comparable to the average value of slope (~0.4) for young clusters obtained by Lada & Lada (2003), however, the slope of the KLF is steeper in the northern region as compared to the southern region. The region is probably ionized by two B2 main-sequence type stars. The mean age of the young stellar objects (YSOs) is found to be 2$\pm$1 Myr which suggests that the identified YSOs could be younger than the ionizing sources of the region. The morphology of the region, the distribution of the YSOs as well as ages of the YSOs and ionizing sources indicate a triggered star formation in the region.

Here we present a kinematic study of the Galactic halo out to a radius of $\sim$ 60 kpc, using 4664 blue horizontal branch (BHB) stars selected from the SDSS/SEGUE survey, to determine key dynamical properties. Using a maximum likelihood analysis, we determine the velocity dispersion profiles in spherical coordinates ($\sigma_{r}$, $\sigma_{\theta}$, $\sigma_{\phi}$) and the anisotropy profile ($\beta$). The radial velocity dispersion profile ($\sigma_{r}$) is measured out to a galactocentric radius of $r \sim 60$ kpc, but due to the lack of proper-motion information, $\sigma_{\theta}$, $\sigma_{\phi}$ and $\beta$ could only be derived directly out to $r \sim25$ kpc. From a starting value of $\beta\approx 0.5$ in the inner parts ($9<r/\kpc<12$), the profile falls sharply in the range $r \approx 13-18$ kpc, with a minimum value of $\beta=-1.2$ at $r=17$ kpc, rising sharply at larger radius. In the outer parts, in the range $25<r/\kpc<56$, we predict the profile to be roughly constant with a value of $\beta\approx 0.5$. The newly discovered kinematic anomalies are shown not to arise from halo substructures. We also studied the anisotropy profile of simulated stellar halos formed purely by accretion and found that they cannot reproduce the sharp dip seen in the data. From the Jeans equation, we compute the stellar rotation curve ($v_{\rm circ}$) of the Galaxy out to $r \sim 25$ kpc. The mass of the Galaxy within $r \lesssim 25$ kpc is determined to be $2.1 \times 10^{11}$ $M_{\sun}$, and with a 3-component fit to $v_{\rm circ}(r)$, we determine the virial mass of the Milky Way dark matter halo to be $M_{\rm vir} = 0.9 ^{+0.4}_{-0.3} \times 10^{12}$ $M_{\sun}$ ($R_{\rm vir} = 249^{+34}_{-31}$ kpc).

The first generation of large-scale chemical tagging surveys, in particular the HERMES/GALAH million star survey, promises to vastly expand our understanding of the chemical and dynamical evolution of the Galaxy. This, however, is contingent on our ability to confidently perform chemical tagging on such a large data-set. Chemical homogeneity has been observed across a range of elements within several Galactic open clusters, yet the level to which this is the case globally, and particularly in comparison to the scatter across clusters themselves, is not well understood. The patterns of elements in coeval cluster members, occupying a complex chemical abundance space, are rooted in the evolution, ultimately the nature of the very late stages, of early generations of stars. The current astrophysical models of such stages are not yet sufficient to explain all observations, combining with our significant gaps in the understanding of star formation, makes this a difficult arena to tackle theoretically. Here, we describe a robust pair-wise metric used to gauge the chemical difference between two stellar components. This metric is then applied to a database of high-resolution literature abundance sources to derive a function describing the probability that two stars are of common evolutionary origin. With this cluster probability function, it will be possible to report a confidence, grounded in empirical observational evidence, with which clusters are detected, independent of the group finding methods. This formulation is also used to probe the role of chemical dimensionality, and that of individual chemical species, on the ability of chemical tagging to differentiate coeval groups of stars.

We present a comparative study of recent works on merger-timescales with dynamical friction and find a strong contrast between idealized/isolated mergers (Boylan-Kolchin et al. 2008) and mergers from a cosmological volume (Jiang et al. 2008). Our study measures the duration of mergers in a cosmological N-body simulation of dark matter, with emphasis on higher redshifts (z < 10) and a lower mass range. In our analysis we consider and compare two merger definitions; tidal disruption and coalescence. We find that the merger-time formula proposed by Jiang et al. (2008) describes our results well and conclude that cosmologically motivated merger-time formulae provide a more versatile and statistically robust approximation for practical applications such as semi-analytic/hybrid models.

We present a multiwavelength study of the NGC 281 complex which contains the young cluster IC 1590 at the center, using deep wide-field optical UBVI_c photometry, slitless spectroscopy along with archival data sets in the near-infrared (NIR) and X-ray. The extent of IC 1590 is estimated to be ~6.5 pc. The cluster region shows a relatively small amount of differential reddening. The majority of the identified young stellar objects (YSOs) are low mass PMS stars having age <1-2 Myr and mass 0.5-3.5 M_\odot. The slope (\Gamma) of the mass function for IC 1590, in the mass range 2 < M/M_\odot \le 54, is found to be -1.11+-0.15. The slope of the K-band luminosity function (0.37+-0.07) is similar to the average value (~0.4) reported for young clusters. The distribution of gas and dust obtained from the IRAS, CO and radio maps indicates clumpy structures around the central cluster. The radial distribution of the young stellar objects, their ages, \Delta(H-K) NIR-excess, and the fraction of classical T Tauri stars suggest triggered star formation at the periphery of the cluster region. However, deeper optical, NIR and MIR observations are needed to have a conclusive view of star formation scenario in the region. The properties of the Class 0/I and Class II sources detected by using the Spitzer mid-infrared observations indicate that a majority of the Class II sources are X-ray emitting stars, whereas X-ray emission is absent from the Class 0/I sources. The spatial distribution of Class 0/I and Class II sources reveals the presence of three sub-clusters in the NGC 281 West region.

We perform a set of non-radiative hydrodynamical simulations of merging spherical halos in order to understand the angular momentum (AM) properties of the galactic halos seen in cosmological simulations. The universal shape of AM distributions seen in simulations is found to be generically produced as a result of mergers. The universal shape is such that it has an excess of low AM material and hence cannot explain the exponential structure of disk galaxies. A resolution to this is suggested by the spatial distribution of low AM material which is found to be in the centre and a conical region close to the axis of rotation. A mechanism that preferentially discards the material in the centre and prevents the material along the poles from falling onto the disc is proposed as a solution. We implement a simple geometric criteria for selective removal of low AM material and show that in order for 90% of halos to host exponential discs one has to reject at least 40% of material. Next, we explore the physical mechanisms responsible for distributing the AM within the halo during a merger. For dark matter there is an inside-out transfer of AM, whereas for gas there is an outside-in transfer, which is due to differences between collisionless and gas dynamics. We also explain the apparent high spin of dark matter halos undergoing mergers and show that a criteria stricter than what is currently used, would be required to detect such unrelaxed halos. Finally, we demonstrate that the misalignment of AM between gas and dark matter only occurs when the intrinsic spins of the merging halos are not aligned with the orbital AM of the system. The self-misalignment (orientation of AM when measured in radial shells not being constant), which could be the cause of warps and anomalous rotation in disks galaxies, also occurs under similar conditions.

(Abridged) This paper presents the first connections made between two local features in velocity-space found in a survey of M giant stars and stellar spatial inhomogeneities on global scales. Comparison to cosmological, chemodynamical stellar halo models confirm that the M giant population is particularly sensitive to rare, recent and massive accretion events. These events can give rise to local observed velocity sequences - a signature of a small fraction of debris from a common progenitor, passing at high velocity through the survey volume, near the pericenters of their eccentric orbits. The majority of the debris is found in much larger structures, whose morphologies are more cloud-like than stream-like and which lie at the orbital apocenters. Adopting this interpretation, the full-space motions represented by the observed velocity features are derived under the assumption that the members within each sequence share a common velocity. Orbit integrations are then used to trace the past and future trajectories of these stars across the sky revealing plausible associations with large, previously-discovered, cloud-like structures. The connections made between nearby velocity structures and these distant clouds represent preliminary steps towards developing coherent maps of such giant debris systems. These maps promise to provide new insights into the origin of debris clouds, new probes of Galactic history and structure, and new constraints on the high-velocity tails of the local dark matter distribution that are essential for interpreting direct detection experiments.

In a series of two recent papers, the frequency and size distribution dependence of extinction spectra for astronomical silicate and graphite grains was analyzed by us in the context of MRN type interstellar dust models. These grains were taken to be homogeneous spheres following the power law $(a^{-3.5})$ size distribution which is very much in use. The analytic formulas were obtained for the graphite and silicate grains in wavelength range 1000\AA - 22,500\AA and their utility was demonstrated. In this paper of the series, we present analytic formulas for the scattering and absorption spectrum of another important constituent of interstellar dust models, namely, the polycyclic aromatic hydrocarbons (PAHs). Relative contribution of the PAHs to extinction {\it vis a vis} carbonaceous classical grains has been examined.

Jeans mass calculated with different combination of parameters involved has shown interesting variation with remarkable shifting of position of peak value from x = m/T = 0.5 to 5.5. The standard deviation is 2.217. In particular, using the harmonic mean square velocity, shifts the peak Jeans mass to x = m/T ~ 2, which is remarkably less than previously reported value of 4.2. Different scales of neutrino structures including virialized moments have also been compared.

The ESO Public Survey VISTA Variables in the Via Lactea (VVV) started in 2010. VVV targets 562 sq. deg in the Galactic bulge and an adjacent plane region and is expected to run for ~5 years. In this paper we describe the progress of the survey observations in the first observing season, the observing strategy and quality of the data obtained. The observations are carried out on the 4-m VISTA telescope in the ZYJHKs filters. In addition to the multi-band imaging the variability monitoring campaign in the Ks filter has started. Data reduction is carried out using the pipeline at the Cambridge Astronomical Survey Unit. The photometric and astrometric calibration is performed via the numerous 2MASS sources observed in each pointing. The first data release contains the aperture photometry and astrometric catalogues for 348 individual pointings in the ZYJHKs filters taken in the 2010 observing season. The typical image quality is ~0.9-1.0". The stringent photometric and image quality requirements of the survey are satisfied in 100% of the JHKs images in the disk area and 90% of the JHKs images in the bulge area. The completeness in the Z and Y images is 84% in the disk, and 40% in the bulge. The first season catalogues contain 1.28x10^8 stellar sources in the bulge and 1.68x10^8 in the disk area detected in at least one of the photometric bands. The combined, multi-band catalogues contain more than 1.63x10^8 stellar sources. About 10% of these are double detections due to overlapping adjacent pointings. These overlapping multiple detections are used to characterise the quality of the data. The images in the JHKs bands extend typically ~4 mag deeper than 2MASS. The magnitude limit and photometric quality depend strongly on crowding in the inner Galactic regions. The astrometry for Ks=15-18 mag has rms ~35-175 mas.

The properties of young stellar clusters (YSCs) in M33, identified from the center out to about twice the size of the bright star-forming disk,are investigated. We find 915 discrete MIR sources as far as the extent of the warped HI disk, i.e. 16 kpc from the galaxy center. Their surface density has a steep radial decline beyond 4.5 kpc, and flattens out beyond the optical radius at 8.5 kpc. We are able to identify YSCs out to 12 kpc. At large galactocentric radii, the paucity of very luminous clusters and the relevance of hot dust emission become evident from the analysis of the bolometric and MIR luminosity functions. The YSC mass and size are correlated with a log-log slope of 2.09, similar to that measured for giant molecular clouds in M33 and the Milky Way, which represent the protocluster environment. Most of the YSCs in our sample have low extinction and ages between 3 and 10 Myr. In the inner regions of M33 the clusters span a wide range of mass (10^2<M<3 10^5 msun) and luminosity 10^38<L{bol}<3 10^{41}erg/s, while at galactocentric radii larger than 4 kpc we find a deficiency of massive clusters. Beyond 7 kpc, where the Halpha surface brightness drops significantly, the dominant YSC population has M<10^3 msun and a slightly older age (10 Myrs). This implies the occurrence of star formation events about 10 Myr ago as far as 10-12 kpc from the center of M33. The cluster L{FUV}--L{Halpha} relation is non-linear for L{FUV}<10^{39}erg/s, in agreement with randomly sampled models of the IMF which, furthermore, shows no appreciable variation throughout the M33 disk.

We uncover the H2 flows in the Corona Australis molecular cloud and in particular identify the flows from the Coronet cluster. Near-infrared H2 v=1--0 S(1), 2.12micron-line, narrow-band imaging survey of the R CrA cloud core was carried out. We identify the best candidate-driving source for each outflow by comparing the flow properties, available proper motions, and the known/estimated properties of the driving sources. We also adopted the thumbrule of outflow power as proportional to source luminosity and inversely proportional to the source age to reach a consensus. Results: Continuum-subtracted, narrow-band images reveal several new Molecular Hydrogen emission-line Objects (MHOs). Together with previously known MHOs and Herbig-Haro objects we catalog at least 14 individual flow components of which 11 appear to be driven by the RCrA aggregate members. The flows originating in the Coronet cluster have lengths of ~0.1-0.2 pc. Eight out of nine submillimeter cores mapped in the Coronet cluster region display embedded stars driving an outflow component. Roughly 80% of the youngest objects in the Coronet are associated with outflows. The MHO flows to the west of the Coronet display lobes moving to the west and vice-versa, resulting in nondetections of the counter lobe in our deep imaging. We speculate that these counterflows may be experiencing a stunting effect in penetrating the dense central core. Conclusions:Although this work has reduced the ambiguities for many flows in the Coronet region, one of the brightest H2 feature (MHO2014) and a few fainter features in the region remain unassociated with a clear driving source. The flows from Coronet, therefore, continue to be interesting targets for future studies.

We apply the method of Burnett & Binney (2010) for the determination of stellar distances and parameters to the internal catalogue of the Radial Velocity Experiment (Steinmetz et al. 2006). Subsamples of stars that either have Hipparcos parallaxes or belong to well-studied clusters, inspire confidence in the formal errors. Distances to dwarfs cooler than ~6000 K appear to be unbiased, but those to hotter dwarfs tend to be too small by ~10% of the formal errors. Distances to giants tend to be too large by about the same amount. The median distance error in the whole sample of 216,000 stars is 28% and the error distribution is similar for both giants and dwarfs. Roughly half the stars in the RAVE survey are giants. The giant fraction is largest at low latitudes and in directions towards the Galactic Centre. Near the plane the metallicity distribution is remarkably narrow and centred on [M/H]-0.04 dex; with increasing |z| it broadens out and its median moves to [M/H] ~ -0.5. Mean age as a function of distance from the Galactic centre and distance |z| from the Galactic plane shows the anticipated increase in mean age with |z|.

We present results on the star-formation activity of an optically obscured region containing an embedded cluster (S255-IR) and molecular gas between two evolved HII regions S255 and S257. We have studied the complex using optical, near-infrared (NIR) imaging, optical spectroscopy and radio continnum mapping at 15 GHz, along with Spitzer-IRAC results. It is found that the main exciting sources of the evolved HII regions S255 and S257 and the compact HII regions associated with S255-IR are of O9.5 - B3 V nature, consistent with previous observations. Our NIR observations reveal 109 likely young stellar object (YSO) candidates in an area of ~ 4'.9 x 4'.9 centered on S255-IR, which include 69 new YSO candidates. Our observations increased the number of previously identified YSOs in this region by 32%. To see the global star formation, we constructed the V-I/V diagram for 51 optically identified IRAC YSOs in an area of ~ 13' x 13' centered on S255-IR. We suggest that these YSOs have an approximate age between 0.1 - 4 Myr, indicating a non-coeval star formation. Using spectral energy distribution models, we constrained physical properties and evolutionary status of 31 and 16 YSO candidates outside and inside the gas ridge, respectively. The models suggest that the sources associated within the gas ridge are of younger population (mean age ~ 1.2 Myr) than the sources outside the gas ridge (mean age ~ 2.5 Myr). The positions of the young sources inside the gas ridge at the interface of the HII regions S255 and S257, favor a site of induced star formation.

We present the third data release of the RAdial Velocity Experiment (RAVE) which is the first milestone of the RAVE project, releasing the full pilot survey. The catalog contains 83,072 radial velocity measurements for 77,461 stars in the southern celestial hemisphere, as well as stellar parameters for 39,833 stars. This paper describes the content of the new release, the new processing pipeline, as well as an updated calibration for the metallicity based upon the observation of additional standard stars. Spectra will be made available in a future release. The data release can be accessed via the RAVE webpage: http://www.rave-survey.org.

We present here a fast code for creating a synthetic survey of the Milky Way. Given one or more color-magnitude bounds, a survey size and geometry, the code returns a catalog of stars in accordance with a given model of the Milky Way. The model can be specified by a set of density distributions or as an N-body realization. We provide fast and efficient algorithms for sampling both types of models. As compared to earlier sampling schemes which generate stars at specified locations along a line of sight, our scheme can generate a continuous and smooth distribution of stars over any given volume. The code is quite general and flexible and can accept input in the form of a star formation rate, age metallicity relation, age velocity dispersion relation and analytic density distribution functions. Theoretical isochrones are then used to generate a catalog of stars and support is available for a wide range of photometric bands. As a concrete example we implement the Besancon Milky Way model for the disc. For the stellar halo we employ the simulated stellar halo N-body models of Bullock & Johnston (2005). In order to sample N-body models, we present a scheme that disperses the stars spawned by an N-body particle, in such a way that the phase space density of the spawned stars is consistent with that of the N-body particles. The code is ideally suited to generating synthetic data sets that mimic near future wide area surveys such as GAIA, LSST and HERMES. As an application we study the prospect of identifying structures in the stellar halo with a simulated GAIA survey. We plan to make the code publicly available at http://galaxia.sourceforge.net.

The ubiquity of substructure in the stellar halo has already been demonstrated by the SDSS and 2MASS and future surveys promise to explore the halo in ever more detail. This paper examines what can be learnt from current and future photometric-databases using group-finding techniques. We compare groups recovered from a sample of M-giants from 2MASS with those found in synthetic surveys of simulated $\Lambda$CDM stellar halos and demonstrate broad consistency. We also find that these recovered groups are likely to represent the majority of high-luminosity ($\log(L/L_{Sun})>6.67$) satellites accreted within the last 10 Gyr and on orbits with apocenters within 100 kpc. However the sensitivity of the M-giant survey to accretion events that were either ancient, from low-luminosity objects or those on radial orbits is limited because of the low number of stars, bias towards high-metallicity stars and the shallow depth. We examine the extent to which these limitations are addressed by current and future surveys, in particular catalogues of main-sequence turn-off (MSTO) stars from SDSS and LSST, and of RR-Lyrae stars from LSST or PanSTARRS. The MSTO surveys are more sensitive to low-luminosity events ($\log(L/L_{Sun}~5$ or less) than the 2MASS M-giant sample, while the RR-Lyrae surveys, with superior depth, are good at detecting events on highly eccentric orbits. When combined we expect these photometric surveys to provide a comprehensive picture of the last 10 Gyr of Galactic accretion. Events older than this are too phase mixed to be discovered and would require additional information such as velocity and chemical abundance.

We identify a new, nearby (0.5 < d < 10 kpc) stream in data from the RAdial Velocity Experiment (RAVE). As the majority of stars in the stream lie in the constellation of Aquarius we name it the Aquarius Stream. We identify 15 members of the stream lying between 30 < l < 75 and -70< b <-50, with heliocentric line-of-sight velocities V_los~-200 km/s. The members are outliers in the radial velocity distribution, and the overdensity is statistically significant when compared to mock samples created with both the Besan\c{c}on Galaxy model and newly-developed code Galaxia. The metallicity distribution function and isochrone fit in the log g - T_eff plane suggest the stream consists of a 10 Gyr old population with [m/H]~-1.0. We explore relations to other streams and substructures, finding the stream cannot be identified with known structures: it is a new, nearby substructure in the Galaxy's halo. Using a simple dynamical model of a dissolving satellite galaxy we account for the localization of the stream. We find that the stream is dynamically young and therefore likely the debris of a recently disrupted dwarf galaxy or globular cluster. The Aquarius stream is thus a specimen of ongoing hierarchical Galaxy formation, rare for being right in the solar suburb.

A density based hierarchical group-finding algorithm is used to identify stellar halo structures in a catalog of M-giants from the Two Micron All Sky Survey (2MASS). The intrinsic brightness of M-giant stars means that this catalog probes deep into the halo where substructures are expected to be abundant and easy to detect. Our analysis reveals 16 structures at high Galactic latitude (greater than 15 degree), of which 10 have been previously identified. Among the six new structures two could plausibly be due to masks applied to the data, one is associated with a strong extinction region and one is probably a part of the Monoceros ring. Another one originates at low latitudes, suggesting some contamination from disk stars, but also shows protrusions extending to high latitudes, implying that it could be a real feature in the stellar halo. The last remaining structure is free from the defects discussed above and hence is very likely a satellite remnant. Although the extinction in the direction of the structure is very low, the structure does match a low temperature feature in the dust maps. While this casts some doubt on its origin, the low temperature feature could plausibly be due to real dust in the structure itself. The angular position and distance of this structure encompass the Pisces overdensity traced by RR Lyraes in Stripe 82 of the Sloan Digital Sky Survey (SDSS). However, the 2MASS M-giants indicate that the structure is much more extended than what is visible with the SDSS, with the point of peak density lying just outside Stripe 82. The morphology of the structure is more like a cloud than a stream and reminiscent of that seen in simulations of satellites disrupting along highly eccentric orbits.

The chemical abundance patterns of the oldest stars in the Galaxy are expected to contain residual signatures of the first stars in the early universe. Numerous studies attempt to explain the intrinsic abundance scatter observed in some metal-poor populations in terms of chemical inhomogeneities dispersed throughout the early Galactic medium due to discrete enrichment events. Just how the complex data and models are to be interpreted with respect to "progenitor yields" remains an open question. Here we show that stochastic chemical evolution models to date have overlooked a crucial fact. Essentially all stars today are born in highly homogeneous star clusters and it is likely that this was also true at early times. When this ingredient is included, the overall scatter in the abundance plane [Fe/H] vs. [X/Fe] (C-space), where X is a nucleosynthetic element, can be much less than derived from earlier models. Moreover, for moderately flat cluster mass functions (gamma < 2), and/or for mass functions with a high mass cut-off (M_max > 10^5 M_sun), stars exhibit a high degree of clumping in C-space that can be identified even in relatively small data samples. Since stellar abundances can be modified by mass transfer in close binaries, clustered signatures are essential for deriving the yields of the first supernovae. We present a statistical test to determine whether a given set of observations exhibit such behaviour. Our initial work focusses on two dimensions in C-space, but we show that the clustering signal can be greatly enhanced by additional abundance axes. The proposed experiment will be challenging on existing 8-10m telescopes, but relatively straightforward for a multi-object echelle spectrograph mounted on a 25-40m telescope.

In this paper, we discuss the properties of CMDs, age, metallicity and radial velocities of eight massive LMC clusters using data taken from FORS2 multiobject spectrograph at the 8.2-meter VLT/UT1. The strong near-infrared Ca II triplet (CaT) lines of RGB stars obtained from the high S/N spectra are used to determine the metallicity and radial velocity of cluster members. We report for the first time spectroscopically determined metallicity values for four clusters based on the mean [Fe/H] value of ~10 cluster members each. We found two concentrations in the distribution of ages of the target clusters. Six have ages between 0.8-2.2 Gyr and the other two, NGC 1754 and NGC 1786, are very old. The metallicity of the six intermediate age clusters, with a mean age of 1.5 Gyr, is -0.49 with a scatter of only 0.04. This tight distribution suggests that a close encounter between the LMC and SMC may have caused not only the restart of cluster formation in the LMC but the generation of the central bar. The metallicity for the two old clusters is similar to that of the other old, metal-poor LMC clusters. We find that the LMC cluster system exhibits disk-like rotation with no clusters appearing to have halo kinematics and there is no evidence of a metallicity gradient in the LMC, in contrast with the stellar population of the MW and M33, where the metallicity decreases as galactocentric distance increases. The LMC's stellar bar may be the factor responsible for the dilution of any kind of gradient in the LMC.

It is generally agreed that interstellar dust grains consist of two main components, namely, silicates and graphites. Some models, like MRN model, assume these grains to be homogeneous spheres following a power law size distribution. This paper presents, in the framework of Mie theory, a parametrization of extinction spectrum curves of the silicates and the graphites separately in terms of frequency and the minimum and maximum of sizes in the distribution. Analytic expressions in ultraviolet and far-ultraviolet are presented for both types of grains.The values of maximum and minimum sizes for which these equations are valid have been identified. These equations can be useful in a number of situations involving silicate and graphite grains.

Using UBVRI Halpha CCD photometric observations and the archival NIR and X-ray data, we have carried out a multi-wavelength study of a young star cluster NGC 7419. An age of 22.5+/-3.0 Myr and a distance of 3230^{+330}_{-430} pc are derived for the cluster with a higher value of color excess ratio E(U-B)/E(B-V) than the normal one. There is an evidence for mass segregation in this dynamically relaxed cluster with mass function slope is in agreement with the Salpeter value. NIR and Halpha excess support the existence of a young (< 2 Myr) stellar population of Herbig Ae/Be stars (> 3.0 M_sun) indicating a second episode of star formation in the cluster region. Using XMM-Newton observations, we found several X-ray sources in the cluster region but none of the Herbig Ae/Be stars is detected in X-rays. We compare the distribution of upper limits for Herbig Ae/Be stars with the X-ray distribution functions of the T-Tauri and the Herbig Ae/Be stars from previous studies, and found that the X-ray emission level of these Herbig Ae/Be stars is not more than L_X ~5.2 x 10^{30} erg/s, which is not significantly higher than for the T-Tauri stars. Therefore, X-ray emission from Herbig Ae/Be stars could be the result of either unresolved companion stars or a process similar to T-Tauri stars. We report an extended X-ray emission from the cluster region NGC 7419, with a total L_X estimate of ~ 1.8 x 10^31 erg/s/arcmin^2. Investigation of dust and CO map of 1 degree region around the cluster indicates the presence of a foreground dust cloud which is most likely associated with star forming region Sh2-154. This cloud harbors uniformly distributed pre main sequence stars (0.1-2.0M_sun) and the star formation in this cloud depend mostly upon the primordial fragmentation.

This paper explores the mapping between the observable properties of a stellar halo in phase- and abundance-space and the parent galaxy's accretion history in terms of the characteristic epoch of accretion and mass and orbits of progenitor objects. The study utilizes a suite of eleven stellar halo models constructed within the context of a standard LCDM cosmology. The results demonstrate that coordinate-space studies are sensitive to the recent (0-8 Gyears ago) merger histories of galaxies (this timescale corresponds to the last few to tens of percent of mass accretion for a Milky-Way-type galaxy). Specifically, the {\it frequency, sky coverage} and {\it fraction of stars} in substructures in the stellar halo as a function of surface brightness are indicators of the importance of recent merging and of the luminosity function of infalling dwarfs. The {\it morphology} of features serves as a guide to the orbital distribution of those dwarfs. Constraints on the earlier merger history (> 8 Gyears ago) can be gleaned from the abundance patterns in halo stars: within our models, dramatic differences in the dominant epoch of accretion or luminosity function of progenitor objects leave clear signatures in the [alpha/Fe] and [Fe/H] distributions of the stellar halo - halos dominated by very early accretion have higher average [alpha/Fe], while those dominated by high luminosity satellites have higher [Fe/H]. This intuition can be applied to reconstruct much about the merger histories of nearby galaxies from current and future data sets.

We have used data from the Sloan Digital Sky Survey (SDSS) Data Release 5 to explore the overall structure and substructure of the stellar halo of the Milky Way using about 4 million color-selected main sequence turn-off stars. We fit oblate and triaxial broken power-law models to the data, and found a `best-fit' oblateness of the stellar halo 0.5<c/a<0.8, and halo stellar masses between Galactocentric radii of 1 and 40kpc of (3.7+/-1.2)x10^8 M_sun. The density profile of the stellar halo is approximately r^{-3}; it is possible that the power law slope is shallower inside 20kpc and steeper outside that radius. Yet, we found that all smooth and symmetric models were very poor fits to the distribution of stellar halo stars because the data exhibit a great deal of spatial substructure. We quantified deviations from a smooth oblate/triaxial model using the RMS of the data around the model profile on scales >~100pc, after accounting for the (known) contribution of Poisson uncertainties. The fractional RMS deviation of the actual stellar distribution from any smooth, parameterized halo model is >~40%: hence, the stellar halo is highly structured. We compared the observations with simulations of galactic stellar halos formed entirely from the accretion of satellites in a cosmological context by analysing the simulations in the same way as the data. While the masses, overall profiles, and degree of substructure in the simulated stellar halos show considerable scatter, the properties and degree of substructure in the Milky Way's halo match well the properties of a `typical' stellar halo built exclusively out of the debris from disrupted satellite galaxies. Our results therefore point towards a picture in which an important fraction of the Milky Way's stellar halo has been accreted from satellite galaxies.

We present a method to numerically estimate the densities of a discretely sampled data based on binary space partitioning tree. We start with a root node containing all the particles and then recursively divide each node into two nodes each containing roughly equal number of particles,until each of the nodes contains only one particle. The volume of such a leaf node provides an estimate of the local density. We implement an entropy-based node splitting criterion that results in a significant improvement in the estimation of densities compared to earlier work. The method is completely metric free and can be applied to arbitrary number of dimensions. We apply this method to determine the phase space densities of dark matter halos obtained from cosmological N-body simulations. We find that contrary to earlier studies, the volume distribution function $v(f)$ of phase space density $f$ does not have a constant slope but rather a small hump at high phase space densities. We demonstrate that a model in which a halo is made up by a superposition of Hernquist spheres is not capable in explaining the shape of $v(f)$ vs $f$ relation, whereas a model which takes into account the contribution of the main halo separately roughly reproduces the behavior as seen in simulations. The use of the presented method is not limited to calculation of phase space densities, but can be used as a general-purpose data-mining tool and due to its speed and accuracy it is ideally suited for analysis of large multidimensional data sets.

Using homogeneous CCD photometric data from the 105-cm Kiso Schmidt telescope covering a 50' x 50' field, we study the mass functions (MFs) of nine open clusters. The ages and Galactocentric distances of the target clusters vary from 16 - 2000 Myr and 9-10.8 kpc, respectively. The values of MF slopes vary from -1.1 to -2.1. The classical value derived by Salpeter (1955) for the slope of the IMF is \Gamma = -1.35. The MFs in the outer regions of the clusters are found to be steeper than in the inner regions, indicating the presence of mass segregation in the clusters.The MF slopes (in the outer region as well as the whole cluster) undergo an exponential decay with the evolutionary parameter \tau (= age/ relaxation time). It seems that the evaporation of low-mass members from outer regions of the clusters is not significant at larger Galactocentric distances. It is concluded that the initial mass function (IMF) in the anticentre direction of the Galaxy might have been steeper than the IMF in the opposite direction. A comparison of the observed CMDs of the clusters with synthetic CMDs gives a photometric binary content of ~40%.

We present $UBVI_c$ CCD photometry of the young open cluster Stock 8 with the aim to study the basic properties and star formation scenario in this region. The radius of the cluster is found to be $\sim 6^{\prime}$ ($\sim 3.6$ pc) and the reddening within the cluster region varies from $E(B-V)=0.40$ to 0.60 mag. The cluster is located at a distance of $2.05 \pm 0.10$ kpc. Using H$\alpha$ slitless spectroscopy and 2MASS NIR data we identified H$\alpha$ emission and NIR excess young stellar objects (YSOs), respectively. The colour-magnitude diagrams of these YSOs reveal that majority of these objects have ages between 1 to 5 Myr indicating a non-coeval star formation in the cluster. Massive stars in the cluster region reveal an average age of $\le$ 2 Myr. In the cluster region ($r \le 6^\prime$) the slope of the mass function (MF), $\Gamma$, in the mass range $\sim 1.0 \le M/M_\odot < 13.4$ can be represented by a power law having a slope of $-1.38\pm0.12$, which agrees well with Salpeter value (-1.35). In the mass range $0.3 \le M/M_\odot < 1.0$, the MF is also found to follow a power law with a shallower slope of $\Gamma = -0.58\pm 0.23$ indicating a break in the slope of the IMF at $\sim 1 M_\odot$. The slope of the $K$-band luminosity function for the cluster ($r \le 6^\prime$) is found to be $0.31\pm0.02$, which is smaller than the average value ($\sim$ 0.4) obtained for embedded star clusters. A significant number of YSOs are distributed along a Nebulous Stream towards the east side of the cluster. A small cluster is embedded in the Nebulous Stream. The YSOs lying in the Nebulous Stream and in the embedded cluster are found to be younger than the stars in the cluster Stock 8. It appears that star formation activity in the Nebulous Stream and embedded cluster may be independent from that of Stock 8.

We present $UBVI_C$ CCD photometry of the young open cluster Be 59 with the aim to study the star formation scenario in the cluster. The radial extent of the cluster is found to be $\sim$ 10 arcmin (2.9 pc). The interstellar extinction in the cluster region varies between $E(B-V) \simeq$ 1.4 to 1.8 mag. The ratio of total-to-selective extinction in the cluster region is estimated as $3.7\pm0.3$. The distance of the cluster is found to be $1.00\pm0.05$ kpc. Using near-infrared colours and slitless spectroscopy, we have identified young stellar objects (YSOs) in the open cluster Be 59 region. The ages of these YSOs range between $<1$ Myr to $\sim$ 2 Myr, whereas the mean age of the massive stars in the cluster region is found to be $\sim$ 2 Myr. There is evidence for second generation star formation outside the boundary of the cluster, which may be triggered by massive stars in the cluster. The slope of the initial mass function, $\Gamma$, in the mass range $2.5 < M/M_\odot \le 28$ is found to be $-1.01\pm0.11$ which is shallower than the Salpeter value (-1.35), whereas in the mass range $1.5 < M/M_\odot \le 2.5$ the slope is almost flat. The slope of the K-band luminosity function is estimated as $0.27\pm0.02$, which is smaller than the average value ($\sim$0.4) reported for young embedded clusters. Approximately 32% of H$\alpha$ emission stars of Be 59 exhibit NIR excess indicating that inner disks of the T-Tauri star (TTS) population have not dissipated. The MSX and IRAS-HIRES images around the cluster region are also used to study the emission from unidentified infrared bands and to estimate the spatial distribution of optical depth of warm and cold interstellar dust.

We present a comprehensive multi-wavelength study of the star-forming region NGC 1893 to explore the effects of massive stars on low-mass star formation. Using near-infrared colours, slitless spectroscopy and narrow-band $H\alpha$ photometry in the cluster region we have identified candidate young stellar objects (YSOs) distributed in a pattern from the cluster to one of the nearby nebulae Sim 129. The $V, (V-I)$ colour-magnitude diagram of the YSOs indicates that majority of these objects have ages between 1 to 5 Myr. The spread in the ages of the YSOs may indicate a non-coeval star formation in the cluster. The slope of the KLF for the cluster is estimated to be $0.34\pm0.07$, which agrees well with the average value ($\sim 0.4$) reported for young clusters. For the entire observed mass range $0.6 < M/M_\odot \le 17.7$ the value of the slope of the initial mass function, $`\Gamma$', comes out to be $-1.27\pm0.08$, which is in agreement with the Salpeter value of -1.35 in the solar neighborhood. However, the value of $`\Gamma$' for PMS phase stars (mass range $0.6 < M/M_\odot \le 2.0$) is found to be $-0.88\pm0.09$ which is shallower than the value ($-1.71\pm0.20$) obtained for MS stars having mass range $2.5 < M/M_\odot \le 17.7$ indicating a break in the slope of the mass function at $\sim 2 M_\odot$. Estimated $`\Gamma$' values indicate an effect of mass segregation for main-sequence stars, in the sense that massive stars are preferentially located towards the cluster center. The estimated dynamical evolution time is found to be greater than the age of the cluster, therefore the observed mass segregation in the cluster may be the imprint of the star formation process. There is evidence for triggered star formation in the region, which seems to govern initial morphology of the cluster.

We present colour-magnitude diagrams of open clusters, located in the range $112^\circ < l < 252^\circ$, manifesting stellar populations in the background of clusters. Some of the populations are found to be located beyond the Perseus arm and may be the part of Norma-Cygnus (outer) arm. The outer arm seems to be continued from $l\sim120^\circ$ to $l\sim235^\circ$. The background populations follow the downward warp of the Galactic plane around $l\sim240^\circ$.

GRB 050730 is a long duration high-redshift burst (z=3.967) discovered by Swift. The afterglow shows variability and is well monitored over a wide wavelength range. We present comprehensive temporal and spectral analysis of the afterglow of GRB 050730 including observations from the millimeter to X-rays. We use multi-wavelength afterglow data to understand the temporal and spectral decay properties with superimposed variability of this high redshift burst. Five telescopes were used to study the decaying afterglow of GRB 050730 in the B, V, r', R, i', I, J and K photometric pass bands. A spectral energy distribution was constructed at 2.9 hours post-burst in the K, J, I, R, V and B bands. X-ray data from the satellites Swift and XMM-Newton were used to study the afterglow evolution at higher energies. The early afterglow shows variability at early times and shows a steepening at ~0.1 days (8.6 ks) in the B, V, r', R, i', I, J and K passbands. The early afterglow light curve decayed with alpha_1 = -0.60+/-0.07 and alpha_2 = -1.71+/-0.06 based on R and I band data. A millimeter detection of the afterglow around 3 days after the burst shows an excess in comparison to predictions. The early X-ray light curve observed by Swift is complex and contains flares. At late times the X-ray light curve can be fit by a powerlaw alpha_x = -2.5+/-0.15 which is steeper than the optical light curve. A spectral energy distribution (SED) was constructed at \~2.9 hours after the burst. An electron energy index, p, of ~ 2.3 was calculated using the SED and the photon index from the X-ray afterglow spectra and indicates that the synchrotron cooling frequency nu_c is above observed frequencies.

In this paper we study the evolution of core and corona of nine open clusters using the projected radial density profiles derived from homogeneous CCD photometric data obtained through the 105-cm Kiso Schmidt telescope. The age and galactocentric distance of the target clusters varies from 16 Myr to 2000 Myr and 9 kpcto 10.8 kpc respectively. Barring Be 62, which is young open cluster, other clusters show a uniform reddening across the cluster region. The reddening in Be 62varies from $E(B-V)_{min}$= 0.70 mag to $E(B-V)_{max}$= 1.00 mag. The corona of six of the clusters in the present sample is found to be elongated, however on the basis of the present sample it is not possible to establish any correlation between the age and shape of the core. The elongated core in the case of young cluster Be 62 may reflect the initial conditions in the parental molecular cloud. The other results of the present study are (i) Core radius `$r_c$' and corona size $`r_{cn}$'/cluster radius $`r_{cl}$' are linearly correlated. (ii) The $r_c/r_{cn}/r_{cl}$ are linearly correlated with the number of stars in that region. (iii) In the age range 10-1000 Myr, the core and corona shrink with age. (iv) We find that in the galactocentric distance range 9 - 10 kpc, the core and corona/cluster extent of the clusters increase with the galactocentric distance.

(Abridged) We report results of a series of non radiative N-body/SPH simulations in a LCDM cosmology. We find that the spin of the baryonic component is on average larger than that of the dark matter (DM) component and we find this effect to be more pronounced at lower redshifts. A significant fraction f of gas has negative angular momentum and this fraction is found to increase with redshift. We describe a toy model in which the tangential velocities of particles are smeared by Gaussian random motions. This model is successful in explaining some of the angular momentum properties. We compare and contrast various techniques to determine the angular momentum distributions (AMDs). We show that broadening of velocity dispersions is unsuitable for making comparisons between gas and DM. We smooth the angular momentum of the particles over a fixed number of neighbors. We find that an analytical function based on gamma distribution can be used to describe a wide variety of profiles, with just one parameter \alpha. The distribution of the shape parameter $\alpha$ for both gas and DM follows roughly a log-normal distribution. The mean and standard deviation of log(\alpha) for gas is -0.04 and 0.11 respectively. About 90-95% of halos have \alpha<1.3, while exponential disks in NFW halos would require 1.3<\alpha<1.6. This implies that a typical halo in simulations has an excess of low angular momentum material as compared to that of observed exponential disks, a result which is consistent with the findings of earlier works. \alpha for gas is correlated with that of DM but they have a significant scatter <\alpha_Gas/\alpha_DM>=1.09 \pm 0.2. \alpha_Gas is also biased towards slightly higher values compared to \alpha_DM.

Pachmarhi Array of \v Cerenkov Telescopes (PACT) consists of a 5$\times$5 array of \v Cerenkov telescopes deployed over an area of 100 $m$ $\times$ 80 $m$, in the form of a rectangular matrix. The experiment is based on atmospheric \v Cerenkov technique using wavefront sampling technique. Each telescope consists of 7 parabolic mirrors mounted para-axially on an equatorial mount. At the focus of each mirror a fast phototube is mounted. In this experiment a large number of parameters have to be measured and recorded from each of the 175 phototubes in the shortest possible time. Further, the counting rates from each phototube as well as the analog sum of the 7 phototubes from each telescope (royal sum) need to be monitored at regular intervals during the run. In view of the complexity of the system, the entire array is divided into four smaller sectors each of which is handled by an independent field signal processing centre (FSPC) housed in a control room that collects, processes and records information from nearby six telescopes that belong to that sector. The distributed data acquisition system (DDAS) developed for the purpose consists of stand-alone sector data acquisition system (SDAS) in each of the four FSPC's and a master data acquisition system (MDAS). MDAS running in the master signal processing centre (MSPC) records data from each of the 25 telescopes. The data acquisition and monitoring PCs (SDAS and MDAS) are networked through LAN. The entire real time software for DDAS is developed in C under $linux$ environment. Most of the hardware in DDAS are designed and fabricated in-house. The design features and the performance of the entire system along with some other auxiliary systems to facilitate the entire observations will be presented.