We analyze the kinematics and spatial distribution of fundamental-mode 15,599 RR Lyrae (RRL) stars in the Galactic bulge by combining OGLE-IV photometric data and Gaia DR2 proper motions. We show that the longitudinal proper motion and the line-of-sight velocity give similar results for the rotation in the Galactic central regions. The angular velocity of bulge RRLs is found to be around $35$ km s$^{-1}$ kpc$^{-1}$, significantly smaller than that for the majority of bulge stars ($50-60$ km s$^{-1}$ kpc$^{-1}$); bulge RRLs have larger velocity dispersion (120$-$140 km s$^{-1}$) than younger stars. The dependence of the kinematics of the bulge RRLs on their metallicities is shown by their rotation curves and spatial distributions. Metal-poor RRLs ([Fe/H]<$-1$) show a smaller bar angle than metal-rich ones. We also find clues suggesting that RRLs in the bulge are not dominated by halo stars. These results might explain some previous conflicting results over bulge RRLs and help understand the chemodynamical evolution of the Galactic bulge.

The grand-design face-on spiral galaxy M51 is an excellent laboratory for studying magnetic fields in galaxies. We present new observations of M51 using the VLA at the frequency range of S-band (2-4GHz), to shed new light on the transition region between the disk and halo. We present images of the distributions of the total intensity, polarized intensity, degree of polarization, and rotation measure (RM). The RM distribution in S-band shows a fluctuating pattern without any apparent large-scale structure. We discuss a model of the depolarization of synchrotron radiation in a multi-layer magneto-ionic medium and compare the model predictions to the polarization data of M51 between 1-8GHz. Since the model predictions strongly differ within the wavelength range of the S-band, the new data are essential. The parameters of the model are adjusted to fit to the data of polarization fractions in a few selected regions. In three spiral arm regions, the turbulent field in the disk dominates with strengths between 18muG and 24muG, while the regular field strengths are 8-16muG. In one inter-arm region, the regular field strength of 18muG exceeds that of the turbulent field of 11muG. The regular field strengths in the halo are 3-5muG. The observed RMs in the disk-halo transition region are probably dominated by tangled regular fields, as predicted from models of evolving dynamos, and/or vertical fields, as predicted from numerical simulations of Parker instabilities or galactic winds. Both types of magnetic fields have frequent reversals on scales similar to or larger than the beam size (550pc) that contribute to an increase of the RM dispersion and to distortions of any large-scale pattern of the regular field. Our study devises new ways of analyzing and interpreting broadband multi-frequency polarization data that will be applicable to future data from, for example, the Square Kilometre Array.

The Square Kilometre Array (SKA) will answer fundamental questions about the origin, evolution, properties, and influence of magnetic fields throughout the Universe. Magnetic fields can illuminate and influence phenomena as diverse as star formation, galactic dynamics, fast radio bursts, active galactic nuclei, large-scale structure, and Dark Matter annihilation. Preparations for the SKA are swiftly continuing worldwide, and the community is making tremendous observational progress in the field of cosmic magnetism using data from a powerful international suite of SKA pathfinder and precursor telescopes. In this contribution, we revisit community plans for magnetism research using the SKA, in the light of these recent rapid developments. We focus in particular on the impact that new radio telescope instrumentation is generating, thus advancing our understanding of key SKA magnetism science areas, as well as the new techniques that are required for processing and interpreting the data. We discuss these recent developments in the context of the ultimate scientific goals for the SKA era.

Aims. The orbital structure of galaxies is strongly influenced by the accuracy of the force calculation during orbit integration. We explore the accuracy of force calculations for two expansion methods and determine which one is preferable for orbit integration. Methods. We specifically compare two methods, one was introduced by Hernquist \& Ostriker (HO), which uses a spherical coordinate system and was built specifically for the Hernquist model, and the other by Vasiliev \& Athanassoula (CylSP) has a cylindrical coordinate system. Our comparisons include the Dehnen profile, its triaxial extension (of which the Hernquist profile is a special case) and a multicomponent system including a bar and disk density distributions for both analytical models and $N$-body realizations.We specifically compare two methods, one was introduced by Hernquist \& Ostriker (HO), which uses a spherical coordinate system and was built specifically for the Hernquist model, and the other by Vasiliev \& Athanassoula (CylSP) has a cylindrical coordinate system. Our comparisons include the Dehnen profile, its triaxial extension (of which the Hernquist profile is a special case) and a multicomponent system including a bar and disk density distributions for both analytical models and $N$-body realizations. Results. For the generalized Dehnen density, the CylSP method is more accurate than the HO method for nearly all inner power-law indices and shapes at all radii. For $N$-body realizations of the Dehnen model, or snapshots of an $N$-body simulation, the CylSP method is more accurate than the HO method in the central region for the oblate, prolate, and triaxial Hernquist profiles if the particle number is more than $5\times 10^5$. For snapshots of the Hernquist models with spherical shape, the HO method is preferred. (Abridged)

We apply the vertical Jeans equation to the kinematics of Milky Way stars in the solar neighbourhood to measure the local dark matter density. More than 90,000 G- and K-type dwarf stars are selected from the cross-matched sample of LAMOST DR5 and Gaia DR2 for our analyses. The mass models applied consist of a single exponential stellar disc, a razor thin gas disc and a constant dark matter density. We first consider the simplified vertical Jeans equation which ignores the tilt term and assumes a flat rotation curve. Under a Gaussian prior on the total stellar surface density, the local dark matter density inferred from Markov Chain Monte Carlo simulations is $0.0133_{-0.0022}^{+0.0024}\ {\rm M}_{\odot}\,{\rm pc}^{-3}$. The local dark matter densities for subsamples in an azimuthal angle range of $-10^{\circ} < \phi < 5^{\circ}$ are consistent within their 1$\sigma$ errors. However, the northern and southern subsamples show a large discrepancy due to plateaux in the northern and southern vertical velocity dispersion profiles. These plateaux may be the cause of the different estimates of the dark matter density between the north and south. Taking the tilt term into account has little effect on the parameter estimations and does not explain the north and south asymmetry. Taking half of the difference of $\sigma_{z}$ profiles as unknown systematic errors, we then obtain consistent measurements for the northern and southern subsamples. We discuss the influence of the vertical data range, the scale height of the tracer population, the vertical distribution of stars and the sample size on the uncertainty of the determination of the local dark matter density.

Galaxy properties are known to correlate most tightly with the galaxy effective stellar velocity dispersion $\sigma_{\rm e}$. Here we look for {\em additional} trends at fixed $\sigma_{\rm e}$ using 1339 galaxies ($M_\ast \gtrsim 6\times10^9$ M$_\odot$) with different morphologies in the MaNGA (DR14) sample with integral-field spectroscopy data. We focus on the gradients ($\gamma_{\rm rms} \equiv \sigma(R_{\rm e}/4)/\sigma_{\rm e}$) of the stellar root-mean-square velocity ($V_{\rm rms} \equiv \sqrt{V^2 + \sigma^2}$), which we show traces the total mass density gradient $\gamma_{\rm tot}$ derived from dynamical models and, more weakly, the bulge fraction. We confirm that $\gamma_{\rm rms}$ increases with $\sigma_{\rm e}$, age and metallicity. We additionally find that these correlations still exist at fixed $\sigma_{\rm e}$, where galaxies with larger $\gamma_{\rm rms}$ are found to be older and more metal-rich. It means that mass density gradients contain information of the stellar population which is not fully accounted for by $\sigma_{\rm e}$. This result puts an extra constraint on our understanding of galaxy quenching. We compare our results with galaxies in the IllustrisTNG hydrodynamical simulations and find that, at fixed $\sigma_{\rm e}$, similar trends exist with age, the bulge fraction, and the total mass density slope but, unlike observations, no correlation with metallicity can be detected in the simulations.

We study the mass distribution and kinematics of the giant elliptical galaxy M87 (NGC 4486) using discrete chemo-dynamical, axisymmetric Jeans equation modelling. Our catalogue comprises 894 globular clusters (GCs) extending to a projected radius of $\sim 430$ kpc with line-of-sight velocities and colours, and Multi Unit Spectroscopic Explorer (MUSE) integral field unit data within the central $2.4$ kpc of the main galaxy. The gravitational potential for our models is a combination of a luminous matter potential with a varying mass-to-light ratio for the main galaxy, a supermassive black hole and a dark matter (DM) potential with a cusped or cored DM halo. The best-fitting models with either a cusped or a cored DM halo show no significant differences and both are acceptable. We obtain a total mass of $(2.16 \pm 0.38) \times 10^{13} M_{\odot}$ within $\sim$ 400 kpc. By including the stellar mass-to-light ratio gradient, the DM fraction increases from $\sim$ 26 percent (with no gradient) to $\sim$ 73 percent within $1\,R_e^{\rm maj}$ (major axis of half-light isophote, 14.2 kpc), and from $\sim$ 84 percent to $\sim$ 94 percent within $5\,R_e^{\rm maj}$ (71.2 kpc). Red GCs have moderate rotation with $V_{\rm max}/\sigma \sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/\sigma \sim$ 0.1. Red GCs have tangential velocity dispersion anisotropy, while blue GCs are consistent with being nearly isotropic. Our results suggest that red GCs are more likely to be born in-situ, while blue GCs are more likely to be accreted.

We develop an analytic mass model for lensing galaxies, based on a broken power-law (BPL) density profile, which is a power-law profile with a mass deficit or surplus in the central region. Under the assumption of an elliptically symmetric surface mass distribution, the deflection angle and magnification can be evaluated analytically for this new model. We compute the theoretical prediction for various quantities, including the volume and surface mass density profiles of the galaxies, and the aperture and luminosity-weighted line-of-sight velocity dispersions, and compare them to those measured from the Illustris simulation. We find an excellent agreement between our model prediction and the simulation, which validates our modeling. The high efficiency and accuracy of our model manifests itself as a promising tool for studying properties of galaxies with strong lensing.

We present population-orbit superposition models for external galaxies based on Schwarzschild's orbit-superposition method, by tagging the orbits with age and metallicity. The models fit the density distributions, as well as kinematic, age and metallicity maps from Integral Field Unit (IFU) spectroscopy observations. We validate the method and demonstrate its power by applying it to mock data, similar to those obtained by the Multi-Unit Spectroscopic Explorer (MUSE) IFU on the Very Large Telescope (VLT). These mock data are created from Auriga galaxy simulations, viewed at three different inclination angles ($\vartheta=40^o, 60^o, 80^o$). Constrained by MUSE-like mock data, our model can recover the galaxy's stellar orbit distribution projected in orbital circularity $\lambda_z$ vs. radius $r$, the intrinsic stellar population distribution in age $t$ vs. metallicity $Z$, and the correlation between orbits' circularity $\lambda_z$ and stellar age $t$. A physically motivated age-metallicity relation improves recovering the intrinsic stellar population distributions. We decompose galaxies into cold, warm and hot + counter-rotating components based on their orbit circularity distribution, and find that the surface density, mean velocity, velocity dispersion, age and metallicity maps of each component from our models well reproduce those from simulation, especially for projections close to edge-on. These galaxies exhibit strong global age vs. $\sigma_z$ relation, which is well recovered by our model. The method has the power to reveal the detailed build-up of stellar structures in galaxies, and offers a complement to local resolved, and high-redshift studies of galaxy evolution.

We report the first discovery of a thick-disk planet, LHS 1815b (TOI-704b, TIC 260004324), detected in the TESS survey. LHS 1815b transits a bright (V = 12.19 mag, K = 7.99 mag) and quiet M dwarf located $ 29.87\pm0.02 pc$ away with a mass of $0.502\pm0.015 M_{\odot}$ and a radius of $0.501\pm0.030 R_{\odot}$. We validate the planet by combining space and ground-based photometry, spectroscopy, and imaging. The planet has a radius of $1.088\pm 0.064 R_{\oplus}$ with a $3 \sigma$ mass upper-limit of $8.7 M_{\oplus}$. We analyze the galactic kinematics and orbit of the host star LHS1815 and find that it has a large probability ($P_{thick}/P_{thin} = 6482$) to be in the thick disk with a much higher expected maximal height ($Z_{max} = 1.8 kpc$) above the Galactic plane compared with other TESS planet host stars. Future studies of the interior structure and atmospheric properties of planets in such systems using for example the upcoming James Webb Space Telescope (JWST), can investigate the differences in formation efficiency and evolution for planetary systems between different Galactic components (thick and thin disks, and halo).

This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $\Delta\Sigma_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $\Delta\Sigma_1$ to the established bulge-type indicator $\Delta\langle\mu_e\rangle$ from Gadotti (2009) and show that classifying by $\Delta\Sigma_1$ agrees well with $\Delta\langle\mu_e\rangle$. $\Delta\Sigma_1$ requires no bulge-disk decomposition and can be measured on SDSS images out to $z = 0.07$. Bulge types using it are mapped onto twenty different structural and stellar-population properties for 12,000 SDSS central galaxies with masses 10.0 < log $M_*$/$M_{\odot}$ < 10.4. New trends emerge from this large sample. Structural parameters show fairly linear log-log relations vs. $\Delta\Sigma_1$ and $\Delta\langle\mu_e\rangle$ with only moderate scatter, while stellar-population parameters show a highly non-linear "elbow" in which specific star-formation rate remains roughly flat with increasing central density and then falls rapidly at the elbow, where galaxies begin to quench. P-bulges occupy the low-density end of the horizontal arm of the elbow and are universally star-forming, while C-bulges occupy the elbow and the vertical branch and exhibit a wide range of star-formation rates at fixed density. The non-linear relation between central density and star-formation rate has been seen before, but this mapping onto bulge class is new. The wide range of star-formation rates in C-bulges helps to explain why bulge classifications using different parameters have sometimes disagreed in the past. The elbow-shaped relation between density and stellar indices suggests that central structure and stellar-populations evolve at different rates as galaxies begin to quench.

We investigate the fundamental plane (FP) evolution of early-type galaxies in the IllustrisTNG-100 simulation (TNG100) from redshift $z=0$ to $z=2$. We find that a tight plane relation already exists as early as $z=2$. Its scatter stays as low as $\sim 0.08$ dex across this redshift range. Both slope parameters $b$ and $c$ (where $R \propto \sigma^b I^c$ with $R$, $\sigma$, and $I$ being the typical size, velocity dispersion, and surface brightness) of the plane evolve mildly since $z=2$, roughly consistent with observations. The FP residual $\rm Res$ ($\equiv\,a\,+\,b\log \sigma\,+\,c\log I\,-\,\log R$, where $a$ is the zero point of the FP) is found to strongly correlate with stellar age, indicating that stellar age can be used as a crucial fourth parameter of the FP. However, we find that $4c+b+2=\delta$, where $\delta \sim 0.8$ for FPs in TNG, rather than zero as is typically inferred from observations. This implies that a tight power-law relation between the dynamical mass-to-light ratio $M_{\rm dyn}/L$ and the dynamical mass $M_{\rm dyn}$ (where $M_{\rm dyn}\equiv 5\sigma^2R/G$, with $G$ being the gravitational constant) is not present in the TNG100 simulation. Recovering such a relation requires proper mixing between dark matter and baryons, as well as star formation occurring with correct efficiencies at the right mass scales. This represents a powerful constraint on the numerical models, which has to be satisfied in future hydrodynamical simulations.

From targeted observations of ground-state OH masers towards 702 Multibeam (MMB) survey 6.7-GHz methanol masers, between Galactic longitudes 186$^{\circ}$ through the Galactic centre to 20$^{\circ}$, made as part of the `MAGMO' project, we present the physical and polarisation properties of the 1720-MHz OH maser transition, including the identification of Zeeman pairs. We present 10 new and 23 previously catalogued 1720-MHz OH maser sources detected towards star formation regions. In addition, we also detected 16 1720-MHz OH masers associated with supernova remnants and two sites of diffuse OH emission. Towards the 33 star formation masers, we identify 44 Zeeman pairs, implying magnetic field strengths ranging from $-$11.4 to $+$13.2 mG, and a median magnetic field strength of $|B_{LOS}|$ $\sim$ 6 mG. With limited statistics, we present the in-situ magnetic field orientation of the masers and the Galactic magnetic field distribution revealed by the 1720-MHz transition. We also examine the association statistics of 1720-MHz OH SFR masers with other ground-state OH masers, excited-state OH masers, class I and class II methanol masers and water masers, and compare maser positions with mid-infrared images of the parent star forming regions. Of the 33 1720-MHz star formation masers, ten are offset from their central exciting sources, and appear to be associated with outflow activity.

We explore the isothermal total density profiles of early-type galaxies (ETGs) in the IllustrisTNG simulation. For the selected 559 ETGs at $z = 0$ with stellar mass $10^{10.7}\mathrm{M}_{\odot} \leqslant M_{\ast} \leqslant 10^{11.9}\mathrm{M}_{\odot}$, the total power-law slope has a mean of $\langle\gamma^{\prime}\rangle = 2.011 \pm 0.007$ and a scatter of $\sigma_{\gamma^{\prime}} = 0.171$ over the radial range 0.4 to 4 times the stellar half mass radius. Several correlations between $\gamma^{\prime}$ and galactic properties including stellar mass, effective radius, stellar surface density, central velocity dispersion, central dark matter fraction and in-situ-formed stellar mass ratio are compared to observations and other simulations, revealing that IllustrisTNG reproduces many correlation trends, and in particular, $\gamma^{\prime}$ is almost constant with redshift below $z = 2$. Through analyzing IllustrisTNG model variations we show that black hole kinetic winds are crucial to lowering $\gamma^{\prime}$ and matching observed galaxy correlations. The effects of stellar winds on $\gamma^{\prime}$ are subdominant compared to AGN feedback, and differ due to the presence of AGN feedback from previous works. The density profiles of the ETG dark matter halos are well-described by steeper-than-NFW profiles, and they are steeper in the full physics (FP) run than their counterparts in the dark matter only (DMO) run. Their inner density slopes anti-correlates (remain constant) with the halo mass in the FP (DMO) run, and anti-correlates with the halo concentration parameter $c_{200}$ in both types of runs. The dark matter halos of low-mass ETGs are contracted whereas high-mass ETGs are expanded, suggesting that variations in the total density profile occur through the different halo responses to baryons.

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

At $q=1.81\pm 0.20 \times 10^{-5}$, KMT-2018-BLG-0029Lb has the lowest planet-host mass ratio $q$ of any microlensing planet to date by more than a factor of two. Hence, it is the first planet that probes below the apparent "pile-up" at $q=5$--10 $\times 10^{-5}$. The event was observed by {\it Spitzer}, yielding a microlens-parallax $\pi_{\rm E}$ measurement. Combined with a measurement of the Einstein radius $\theta_{\rm E}$ from finite-source effects during the caustic crossings, these measurements imply masses of the host $M_{\rm host}=1.14^{+0.10}_{-0.12}\, M_\odot$ and planet $M_{\rm planet} = 7.59^{+0.75}_{-0.69}\,M_\oplus$, system distance $D_L = 3.38^{+0.22}_{-0.26}\,\,{\rm kpc}$ and projected separation $a_\perp = 4.27^{+0.21}_{-0.23}\,{\rm au}$. The blended light, which is substantially brighter than the microlensed source, is plausibly due to the lens and could be observed at high resolution immediately.

In our earlier 2019 paper, we evaluated the reliability of Schwarzschild's orbit-superposition dynamical modelling method in estimating the internal mass distribution, intrinsic stellar shapes and orbit distributions of early-type galaxies (ETGs) taken from the Illustris cosmological simulation. We now apply the same techniques to galaxies taken from the integral-field survey Mapping Nearby Galaxies with APO (MaNGA), using a sample of 149 ETGs in the mass range of $10^{9.90}\sim10^{11.80} M_{\odot}$ and made up of 105 central and 44 satellite galaxies. We find that low-mass ETGs with $\log(M_*/M_{\odot})<11.1$ have an average dark matter fraction of $\sim0.2$ within one effective radius $R_{\rm e}$, tend to be oblate-like, and are dominated by rotation about their minor axis. High-mass ETGs with $\log(M_*/M_{\odot})>11.1$ have an average dark matter fraction of $\sim0.4$ within one effective radius $R_{\rm e}$, tend to be prolate-like, and are dominated by rotation about their major axis and by centrophilic orbits. The changes of internal structures within one $R_{\rm e}$ are dominated by the total stellar mass of the individual galaxies. We find no differences of internal structures between central and satellite ETGs for the same stellar masses. However, for similar stellar mass and colour distributions, we find that ETGs more prolate-like, or with more hot orbits, tend to have higher close neighbour counts at $r_p\sim40$ kpc.

We report the discovery of a super-Jovian planet in the microlensing event KMT-2016-BLG-1836, which was found by the Korea Microlensing Telescope Network's high-cadence observations (\Gamma ~ 4~{hr}^{-1}). The planet-host mass ratio q ~ 0.004. A Bayesian analysis indicates that the planetary system is composed of a super-Jovian M_{planet} = 2.2_{-1.1}^{+1.9} M_{J} planet orbiting an M or K dwarf M_{\rm host} = 0.49_{-0.25}^{+0.38} M_{Sun}, at a distance of D_{L} = 7.1_{-2.4}^{+0.8} kpc. The projected planet-host separation is 3.5^{+1.1}_{-0.9} AU, implying that the planet is located beyond the snowline of the host star. Future high-resolution images can potentially strongly constrain the lens brightness and thus the mass and distance of the planetary system. Without considering detailed detection efficiency, selection or publication biases, we find a potential "mass ratio desert" at -3.7 \lesssim \log q \lesssim -3.0 for the 31 published KMTNet planets.

We present the analysis of the binary-lens microlensing event OGLE-2013-BLG-0911. The best-fit solutions indicate the binary mass ratio of q~0.03 which differs from that reported in Shvartzvald+2016. The event suffers from the well-known close/wide degeneracy, resulting in two groups of solutions for the projected separation normalized by the Einstein radius of s~0.15 or s~7. The finite source and the parallax observations allow us to measure the lens physical parameters. The lens system is an M-dwarf orbited by a massive Jupiter companion at very close (M_{host}=0.30^{+0.08}_{-0.06} M_{Sun}, M_{comp}=10.1^{+2.9}_{-2.2} M_{Jup}, a_{exp}=0.40^{+0.05}_{-0.04} au) or wide (M_{host}=0.28^{+0.10}_{-0.08} M_{Sun}, M_{comp}=9.9^{+3.8}_{-3.5}M_{Jup}, a_{exp}=18.0^{+3.2}_{-3.2} au) separation. Although the mass ratio is slightly above the planet-brown dwarf (BD) mass-ratio boundary of q=0.03 which is generally used, the median physical mass of the companion is slightly below the planet-BD mass boundary of 13M_{Jup}. It is likely that the formation mechanisms for BDs and planets are different and the objects near the boundaries could have been formed by either mechanism. It is important to probe the distribution of such companions with masses of ~13M_{Jup} in order to statistically constrain the formation theories for both BDs and massive planets. In particular, the microlensing method is able to probe the distribution around low-mass M-dwarfs and even BDs which is challenging for other exoplanet detection methods.

We present deep rest-frame UV spectroscopic observations using the Gran Telescopio Canarias of six gravitationally lensed Lya emitters (LAEs) at $2.36<z<2.82$ selected from the BELLS GALLERY survey. By taking the magnifications into account, we show that LAEs can be as luminous as L(Lya) = 30x10$^{42}$ erg s-1 and M(UV) = -23 (AB) without invoking an AGN component, in contrast with previous findings. We measure Lya rest-frame equivalent widths, EW(Lya), ranging from 16\AA to 50\AA and Lya escape fractions, fesc(Lya), from 10% to 40%. Large EW(Lya) and fesc(Lya) are found predominantly in LAEs showing weak low-ionization ISM absorption (EW < 1\AA) and narrow Lya profiles (< 300 km s-1 FWHM) with their peak close (< 80 km s-1) to their systemic redshifts, suggestive of less scatter from low HI column densities that favours the escape of Lya photons. We infer stellar metallicities of Z/Zsun ~ 0.2 in almost all LAEs by comparing the P-Cygni profiles of the wind lines NV1240\AA and CIV1549\AA with those from stellar synthesis models. We also find a trend between M(UV) and the velocity offset of ISM absorption lines, such as the most luminous LAEs experience stronger outflows. The most luminous LAEs show star formation rates up to 180 Msun yr-1, yet they appear relatively blue ($\beta$(UV) ~ -1.8 to -2.0) showing evidence of little dust attenuation (E(B-V) = 0.10-0.14). These luminous LAEs may be particular cases of young starburst galaxies that have had no time to form large amounts of dust. If so, they are ideal laboratories to study the early phase of massive star formation, stellar and dust mass growth, and chemical enrichment histories of starburst galaxies at high-z.

In the standard cosmological model, dark matter drives the structure formation and constructs potential wells within which galaxies may form. The baryon fraction in dark halos can reach the universal value (15.7%) in massive clusters and decreases rapidly as the mass of the system decreases. The formation of dwarf galaxies is sensitive both to baryonic processes and the properties of dark matter owing to the shallow potential wells in which they form. In dwarf galaxies in the Local Group, dark matter dominates the mass content even within their optical-light half-radii (r_e ~ 1 kpc). However, recently it has been argued that not all dwarf galaxies are dominated by dark matter. Here we report 19 dwarf galaxies that could consist mainly of baryons up to radii well beyond r_e, at which point they are expected to be dominated by dark matter. Of these, 14 are isolated dwarf galaxies, free from the influence of nearby bright galaxies and high dense environments. This result provides observational evidence that could challenge the formation theory of low-mass galaxies within the framework of standard cosmology. Further observations, in particular deep imaging and spatially-resolved kinematics, are needed to constrain the baryon fraction better in such galaxies.

We report the lens mass and distance measurements of the nearby microlensing event TCP J05074264+2447555. We measure the microlens parallax vector ${\pi}_{\rm E}$ using Spitzer and ground-based light curves with constraints on the direction of lens-source relative proper motion derived from Very Large Telescope Interferometer (VLTI) GRAVITY observations. Combining this ${\pi}_{\rm E}$ determination with the angular Einstein radius $\theta_{\rm E}$ measured by VLTI GRAVITY observations, we find that the lens is a star with mass $M_{\rm L} = 0.495 \pm 0.063~M_{\odot}$ at a distance $D_{\rm L} = 429 \pm 21~{\rm pc}$. We find that the blended light basically all comes from the lens. The lens-source proper motion is $\mu_{\rm rel,hel} = 26.55 \pm 0.36~{\rm mas\,yr^{-1}}$, so with currently available adaptive-optics (AO) instruments, the lens and source can be resolved in 2021. This is the first microlensing event whose lens mass is unambiguously measured by interferometry + satellite parallax observations, which opens a new window for mass measurements of isolated objects such as stellar-mass black holes.

We study the evolutionary trend of the total density profile of early-type galaxies (ETGs) in IllustrisTNG. To this end, we trace ETGs from $z=0$ to $z=4$ and measure the power-law slope $\gamma^{\prime}$ of the total density profile for their main progenitors. We find that their $\gamma^{\prime}$ steepen on average during $z\sim4-2$, then becoming shallower until $z=1$, after which they remain almost constant, aside from a residual trend of becoming shallower towards $z=0$. We also compare to a statistical sample of ETGs at different redshifts, selected based on their luminosity profiles and stellar masses. Due to different selection effects, the average slopes of the statistical samples follow a modified evolutionary trend. They monotonically decrease since $z=3$, and after $z\approx 1$, they remain nearly invariant with a mild increase towards $z=0$. These evolutionary trends are mass-dependent for both samples, with low-mass galaxies having in general steeper slopes than their more massive counterparts. Galaxies that transitioned to ETGs more recently have steeper mean slopes as they tend to be smaller and more compact at any given redshift. By analyzing the impact of mergers and AGN feedback on the progenitors' evolution, we conjecture a multi-phase path leading to isothermality in ETGs: dissipation associated with rapid wet mergers tends to steepen $\gamma^{\prime}$ from $z=4$ to $z=2$, whereas subsequent AGN feedback (especially in the kinetic mode) makes $\gamma^{\prime}$ shallower again from $z=2$ to $z=1$. Afterwards, passive evolution from $z=1$ to $z=0$, mainly through gas-poor mergers, mildly decreases $\gamma^{\prime}$ and maintains the overall mass distribution close to isothermal.

In this paper we present a detailed analysis of the Faraday depth (FD) spectrum and its clean components obtained through the application of the commonly used technique of Faraday rotation measure synthesis to analyze spectro-polarimetric data. In order to directly compare the Faraday depth spectrum with physical properties of a magneto-ionic medium, we generated synthetic broad-bandwidth spectro-polarimetric observations from magnetohydrodynamic (MHD) simulations of a transonic, isothermal, compressible turbulent medium. We find that correlated magnetic field structures give rise to a combination of spiky, localized peaks at certain FD values, and broad structures in the FD spectrum. Although the majority of these spiky FD structures appear narrow, giving an impression of a Faraday thin medium, we show that they arise from strong synchrotron emissivity at that FD. Strong emissivity at a FD can arise because of both strong spatially-local polarized synchrotron emissivity at a FD or accumulation of weaker emissions along the distance through a medium that have Faraday depths within half the width of the rotation measure spread function. Such a complex Faraday depth spectrum is a natural consequence of MHD turbulence when the lines of sight pass through a few turbulent cells. This therefore complicates the convention of attributing narrow FD peaks to presence of a Faraday rotating medium along the line of sight. Our work shows that it is difficult to extract the FD along a line of sight from the Faraday depth spectrum using standard methods for a turbulent medium in which synchrotron emission and Faraday rotation occur simultaneously.

We apply gravity-based and density-based methods to identify clouds in numerical simulations of the star-forming, three-phase interstellar medium (ISM), and compare their properties and their global correlation with the star formation rate over time. The gravity-based method identifies bound objects, which have masses M ~ 10^3 - 10^4 M_solar at densities n_H ~ 100 cm^-3, and traditional virial parameters alpha_v ~ 0.5 - 5. For clouds defined by a density threshold n_H,min , the average virial parameter decreases, and the fraction of material that is genuinely bound increases, at higher n_H,min. Surprisingly, these clouds can be unbound even when alpha_v < 2, and high mass clouds (10^4 - 10^6 M_solar) are generally unbound. This suggests that the traditional alpha_v is at best an approximate measure of boundedness in the ISM. All clouds have internal turbulent motions increasing with size as sigma ~ 1 km/s(R/ pc)^1/2, similar to observed relations. Bound structures comprise a small fraction of the total simulation mass, with star formation efficiency per free-fall time epsilon_ff ~ 0.4. For n_H,min = 10 - 100 cm^-3, epsilon_ff ~ 0.03 - 0.3, increasing with density. Temporal correlation analysis between SFR(t) and aggregate mass M(n_H,min;t) at varying n_H,min shows that time delays to star formation are t_delay ~ t_ff(n_H,min). Correlation between SFR(t) and M(nH,min;t) systematically tightens at higher n_H,min. Considering moderate-density gas, selecting against high virial parameter clouds improves correlation with SFR, consistent with previous work. Even at high n_H,min, the temporal dispersion in (SFR-epsilon_ff M/t_ff )/<SFR> is ~ 50%, due to the large-amplitude variations and inherent stochasticity of the system.

We report the discovery and the analysis of the short (tE < 5 days) planetary microlensing event, OGLE-2015-BLG-1771. The event was discovered by the Optical Gravitational Lensing Experiment (OGLE), and the planetary anomaly (at I ~ 19) was captured by The Korea Microlensing Telescope Network (KMTNet). The event has three surviving planetary models that explain the observed light curves, with planet-host mass ratio q \~ 5.4 * 10^{-3}, 4.5 * 10^{-3} and 4.5 * 10^{-2}, respectively. The first model is the best-fit model, while the second model is disfavored by \Delta\chi^2 ~ 3. The last model is strongly disfavored by \Delta\chi^2 ~ 15 but not ruled out. A Bayesian analysis using a Galactic model indicates that the first two models are probably composed of a Saturn-mass planet orbiting a late M dwarf, while the third one could consist of a super-Jovian planet and a mid-mass brown dwarf. The source-lens relative proper motion is mu_rel ~ 9 mas/yr, so the source and lens could be resolved by current adaptive-optics (AO) instruments in 2021 if the lens is luminous.

The microlensing event OGLE-2015-BLG-0448 was observed by Spitzer and lay within the tidal radius of the globular cluster NGC 6558. The event had moderate magnification and was intensively observed, hence it had the potential to probe the distribution of planets in globular clusters. We measure the proper motion of NGC 6558 ($\mu_{\rm cl}$(N,E) = (+0.36+-0.10, +1.42+-0.10) mas/yr) as well as the source and show that the lens is not a cluster member. Even though this particular event does not probe the distribution of planets in globular clusters, other potential cluster lens events can be verified using our methodology. Additionally, we find that microlens parallax measured using OGLE photometry is consistent with the value found based on the light curve displacement between Earth and Spitzer.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

We present ALMA observations of the CO(1-0) line and 3-mm continuum emission in eight ultraluminous infrared (IR) quasi-stellar objects (QSOs) at z = 0.06-0.19. All eight IR QSO hosts are clearly resolved in their CO molecular gas emission with a median source size of 3.2 kpc, and seven out of eight sources are detected in 3-mm continuum, which is found to be more centrally concentrated with respect to molecular gas with sizes of 0.4-1.0 kpc. Our observations reveal a diversity of CO morphology and kinematics for the IR QSO systems which can be roughly classified into three categories, rotating gas disk with ordered velocity gradient, compact CO peak with disturbed velocity, and multiple CO distinct sources undergoing a merger between luminous QSO and a companion galaxy separated by a few kpc. The molecular gas in three of IR QSO hosts are found to be rotation-dominated with the ratio of the maximum rotation velocity to the local velocity dispersion of $V_{\rm rot}/\sigma=4-6$. Basic estimates of the dynamical masses within the CO-emitting regions give masses between $7.4\times10^9$ and $6.9\times10^{10}$ $M_\odot$. We find an increasing trend between BH mass accretion rate and star formation rate (SFR) over three orders of magnitude in far-IR luminosity/SFR, in line with the correlation between QSO bolometric luminosity and SF activity, indicative of a likely direct connection between AGN and SF activity over galaxy evolution timescales.

Planet population synthesis models predict an abundance of planets with semi-major axes between 1-10 au, yet they lie at the edge of the detection limits of most planet finding techniques. Discovering these planets and studying their distribution is critical to understanding the physical processes that drive planet formation. ROME/REA is a gravitational microlensing project whose main science driver is to discover exoplanets in the cold outer regions of planetary systems. To achieve this, it uses a novel approach combining a multi-band survey with reactive follow-up observations, exploiting the unique capabilities of the Las Cumbres Observatory (LCO) global network of robotic telescopes combined with a Target and Observation Manager (TOM) system. We present the main science objectives and a technical overview of the project, including initial results.

We report the discovery of an exoplanet in microlensing event OGLE-2015-BLG-1649. The planet/host-star mass ratio is $q =7.2 \times 10^{-3}$ and the projected separation normalized by the Einstein radius is $s = 0.9$. The upper limit of the lens flux is obtained from adaptive optics observations by IRCS/Subaru, which excludes the probability of a G-dwarf or more massive host star and helps to put a tighter constraint on the lens mass as well as commenting on the formation scenarios of giant planets orbiting low-mass stars. We conduct a Bayesian analysis including constraints on the lens flux to derive the probability distribution of the physical parameters of the lens system. We thereby find that the masses of the host star and planet are $M_{L} = 0.34 \pm 0.19 M_{\odot}$ and $M_{p} = 2.5^{+1.5}_{-1.4} M_{Jup}$, respectively. The distance to the system is $D_{L} = 4.23^{+1.51}_{-1.64}$kpc. The projected star-planet separation is $a_{\perp} = 2.07^{+0.65}_{-0.77}$AU. The lens-source relative proper motion of the event is quite high, at $\sim 7.1 \, {\rm mas/yr}$. Therefore, we may be able to determine the lens physical parameters uniquely or place much stronger constraints on them by measuring the color-dependent centroid shift and/or the image elongation with additional high resolution imaging already a few years from now.

The Radio Arc is one of the brightest systems of non-thermal filaments (NTFs) in the Galactic Center, located near several prominent HII regions (Sickle and Pistol) and the Quintuplet stellar cluster. We present observations of the Arc NTFs using the S-, C-, and X-bands of the Very Large Array interferometer. Our images of total intensity reveal large-scale helical features that surround the Arc NTFs, very narrow sub-filamentation, and compact sources along the NTFs. The distribution of polarized intensity is confined to a relatively small area along the NTFs. There are elongated polarized structures that appear to lack total intensity counterparts. We detect a range of rotation measure values from -1000 to -5800 rad m$\rm^{-2}$, likely caused by external Faraday rotation along the line of sight. After correcting for Faraday rotation, the intrinsic magnetic field orientation is found to generally trace the extent of the NTFs. However, the intrinsic magnetic field in several regions of the Arc NTFs shows an ordered pattern that is rotated with respect to the extent of the NTFs. We suggest this changing pattern may be caused by an additional magnetized source along the line of sight, so that we observe two field systems superposed in our observations. We suggest that the large scale helical segments near the Radio Arc could be components of such a source causing these changes in intrinsic magnetic field, and some variations in the polarization and rotation measure values along the NTFs.

Chondrules are silicate spheroids found in meteorites, serving as important fossil records of the early solar system. In order to form chondrules, chondrule precursors must be heated to temperatures much higher than the typical conditions in the current asteroid belt. One proposed mechanism for chondrule heating is the passage through bow shocks of highly eccentric planetesimals in the protoplanetary disk in the early solar system. However, it is difficult for planetesimals to gain and maintain such high eccentricities. In this paper, we present a new scenario in which planetesimals in the asteroid belt region are excited to high eccentricities by the Jovian sweeping secular resonance in a depleting disk, leading to efficient formation of chondrules. We study the orbital evolution of planetesimals in the disk using semi-analytic models and numerical simulations. We investigate the dependence of eccentricity excitation on the planetesimal's size as well as the physical environment, and calculate the probability for chondrule formation. We find that 50 - 2000 km planetesimals can obtain eccentricities larger than 0.6 and cause effective chondrule heating. Most chondrules form in high velocity shocks, in low density gas, and in the inner disk. The fraction of chondrule precursors which become chondrules is about 4 - 9 % between 1.5 - 3 AU. Our model implies that the disk depletion timescale is $\tau_\mathrm{dep}\approx 1~\mathrm{Myr}$, comparable to the age spread of chondrules; and that Jupiter formed before chondrules, no more than 0.7 Myr after the formation of the CAIs.

We use the Illustris simulation to study the distributions of ellipticities and kinematic misalignments of galactic projections, as well as the intrinsic shapes and rotation of the simulated galaxies. Our results for the projections of galaxies display clear trends of an overall increase of kinematic misalignment and a slight decrease of ellipticity for fast rotators with increasing masses, while revealing no clear distinction between slow rotators of different mass. It is also found that the number of very slow rotators with large ellipticities is much larger than found in observations. The intrinsic properties of the galaxies are then analysed. The results for the intrinsic shapes of the galaxies are mostly consistent with previous results inferred from observational data. The distributions of intrinsic misalignment of the galaxies suggest that some of the galaxies produced by Illustris have significant rotation around their medium axes. Further analysis reveals that most of these galaxies display signs of non-equilibrium. We then study the evolution of the intrinsic misalignments and shapes of three specific Illustris galaxies, which we consider as typical ones, along the main progenitor line of their merger trees, revealing how mergers influence the intrinsic shapes and kinematics: the spin axis in general stays close to the shortest axis, and tends to quickly relax to such an equilibrium state within a few dynamical times of the galaxy after major perturbations; triaxiality and intrinsic flatness in general decrease with time, however, sometimes increases occur that are clearly seen to correlate with major merger events.

We present the analysis of the event OGLE-2017-BLG-1186 from the 2017 Spitzer microlensing campaign. This is a remarkable microlensing event because its source is photometrically bright and variable, which makes it possible to perform an asteroseismic analysis using ground-based data. We find that the source star is an oscillating red giant with average timescale of $\sim 9$ d. The asteroseismic analysis also provides us source properties including the source angular size ($\sim 27~\mu{\rm as}$) and distance ($\sim 11.5$ kpc), which are essential for inferring the properties of the lens. When fitting the light curve, we test the feasibility of Gaussian Processes (GPs) in handling the correlated noise caused by the variable source. We find that the parameters from the GP model are generally more loosely constrained than those from the traditional $\chi^2$ minimization method. We note that this event is the first microlensing system for which asteroseismology and GPs have been used to account for the variable source. With both finite-source effect and microlens parallax measured, we find that the lens is likely a $\sim 0.045~M_{\odot}$ brown dwarf at distance $\sim 9.0$ kpc, or a $\sim 0.073~M_{\odot}$ ultracool dwarf at distance $\sim 9.8$ kpc. Combining the estimated lens properties with a Bayesian analysis using a Galactic model, we find a $\sim 35$ per cent probability for the lens to be a bulge object and $\sim 65$ per cent to be a background disc object.

We present the status of the Dark Energy Spectroscopic Instrument (DESI) and its plans and opportunities for the coming decade. DESI construction and its initial five years of operations are an approved experiment of the US Department of Energy and is summarized here as context for the Astro2020 panel. Beyond 2025, DESI will require new funding to continue operations. We expect that DESI will remain one of the world's best facilities for wide-field spectroscopy throughout the decade. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

We report a multiplanetary system found from the analysis of microlensing event OGLE-2018-BLG-1011, for which the light curve exhibits a double-bump anomaly around the peak. We find that the anomaly cannot be fully explained by the binary-lens or binary-source interpretations and its description requires the introduction of an additional lens component. The 3L1S (3 lens components and a single source) modeling yields three sets of solutions, in which one set of solutions indicates that the lens is a planetary system in a binary, while the other two sets imply that the lens is a multiplanetary system. By investigating the fits of the individual models to the detailed light curve structure, we find that the multiple-planet solution with planet-to-host mass ratios $\sim 9.5\times 10^{-3}$ and $\sim 15\times 10^{-3}$ are favored over the other solutions. From the Bayesian analysis, we find that the lens is composed of two planets with masses $1.8^{+3..4}_{-1.1}~M_{\rm J}$ and $2.8^{+5.1}_{-1.7}~M_{\rm J}$ around a host with a mass $0.18^{+0.33}_{-0.10}~M_\odot$ and located at a distance $7.1^{+1.1}_{-1.5}~{\rm kpc}$. The estimated distance indicates that the lens is the farthest system among the known multiplanetary systems. The projected planet-host separations are $a_{\perp,2}=1.8^{+2.1}_{-1.5}~{\rm au}$ ($0.8^{+0.9}_{-0.6}~{\rm au}$) and $a_{\perp,3}=0.8^{+0.9}_{-0.6}~{\rm au}$, where the values of $a_{\perp,2}$ in and out the parenthesis are the separations corresponding to the two degenerate solutions, indicating that both planets are located beyond the snow line of the host, as with the other four multiplanetary systems previously found by microlensing.

We present the analysis of stellar binary microlensing event OGLE-2015-BLG-0060 based on observations obtained from 13 different telescopes. Intensive coverage of the anomalous parts of the light curve was achieved by automated follow-up observations from the robotic telescopes of the Las Cumbres Observatory. We show that, for the first time, all main features of an anomalous microlensing event are well covered by follow-up data, allowing us to estimate the physical parameters of the lens. The strong detection of second-order effects in the event light curve necessitates the inclusion of longer-baseline survey data in order to constrain the parallax vector. We find that the event was most likely caused by a stellar binary-lens with masses $M_{\star1} = 0.87 \pm 0.12 M_{\odot}$ and $M_{\star2} = 0.77 \pm 0.11 M_{\odot}$. The distance to the lensing system is 6.41 $\pm 0.14$ kpc and the projected separation between the two components is 13.85 $\pm 0.16$ AU. Alternative interpretations are also considered.

Motivated by the recently discovered kinematic "Hubble sequence" shown by the stellar orbit-circularity distribution of 260 CALIFA galaxies, we make use of a comparable galaxy sample at z = 0 with a stellar mass range from 5E9 to 5E11 solar masses, selected from the IllustrisTNG simulation and study their stellar orbit compositions in relation to a number of other fundamental galaxy properties.We find that the TNG100 simulation broadly reproduces the observed fractions of different orbital components and their stellar mass dependencies. In particular, the mean mass dependencies of the luminosity fractions for the kinematically warm and hot orbits are well reproduced within model uncertainties of the observed galaxies. The simulation also largely reproduces the observed peak and trough features at a stellar mass of 1-2E10 solar masses, in the mean distributions of the cold- and hot-orbit fractions, respectively, indicating fewer cooler orbits and more hotter orbits in both more- and less-massive galaxies beyond such a mass range. Several marginal disagreements are seen between the simulation and observations: the average cold-orbit (counter-rotating) fractions of the simulated galaxies below (above) a stellar mass of 6E10 solar masses, are systematically higher than the observational data by < 10% (absolute orbital fraction); the simulation also seems to produce more scatter for the cold-orbit fraction and less so for the non-cold orbits at any given galaxy mass. Possible causes that stem from the adopted heating mechanisms are discussed.

Sharpless 2-27 (Sh2-27) is a nearby ${\rm H\,{\small II}}$ region excited by $\zeta$Oph. We present observations of polarized radio emission from 300 to 480$\,$MHz towards Sh2-27, made with the Parkes 64$\,$m Radio Telescope as part of the Global Magneto-Ionic Medium Survey. These observations have an angular resolution of $1.35^{\circ}$, and the data are uniquely sensitive to magneto-ionic structure on large angular scales. We demonstrate that background polarized emission towards Sh2-27 is totally depolarized in our observations, allowing us to investigate the foreground. We analyse the results of Faraday tomography, mapping the magnetised interstellar medium along the 165$\,$pc path to Sh2-27. The Faraday dispersion function in this direction has peaks at three Faraday depths. We consider both Faraday thick and thin models for this observation, finding that the thin model is preferred. We further model this as Faraday rotation of diffuse synchrotron emission in the Local Bubble and in two foreground neutral clouds. The Local Bubble extends for 80$\,$pc in this direction, and we find a Faraday depth of $-0.8 \pm 0.4\,$rad$\,$m$^{-2}$. This indicates a field directed away from the Sun with a strength of $-2.5\pm1.2\,\mu$G. The near and far neutral clouds are each about 30$\,$pc thick, and we find Faraday depths of $-6.6\pm0.6\,$rad$\,$m$^{-2}$ and $+13.7\pm0.8\,$rad$\,$m$^{-2}$, respectively. We estimate that the line-of-sight magnetic strengths in the near and far cloud are $B_{\parallel, \text{near}} \approx -15\,\mu\text{G}$ and $B_{\parallel, \text{far}} \approx +30\,\mu\text{G}$. Our results demonstrate that Faraday tomography can be used to investigate the magneto-ionic properties of foreground features in front of nearby ${\rm H\,{\small II}}$ regions.

We report the discovery of the planet OGLE-2018-BLG-0532Lb, with very obvious signatures in the light curve that lead to an estimate of the planet-host mass ratio $q=M_{\rm planet}/M_{\rm host}\simeq 1\times10^{-4}$. Although there are no obvious systematic residuals to this double-lens/single-source (2L1S) fit, we find that $\chi^2$ can be significantly improved by adding either a third lens (3L1S, $\Delta\chi^2=81$) or second source (2L2S, $\Delta\chi^2=65$) to the lens-source geometry. After thorough investigation, we conclude that we cannot decisively distinguish between these two scenarios and therefore focus on the robustly-detected planet. However, given the possible presence of a second planet, we investigate to what degree and with what probability such additional planets may affect seemingly single-planet light curves. Our best estimates for the properties of the lens star and the secure planet are: a host mass $M\sim 0.25\,M_\odot$, system distance $D_L\sim 1\,$kpc and planet mass $m_{p,1}= 8\,M_\oplus$ with projected separation $a_{1,\perp}=1.4\,$au. However, there is a relatively bright $I=18.6$ (and also relatively blue) star projected within $<50\,$mas of the lens, and if future high-resolution images show that this is coincident with the lens, then it is possible that it is the lens, in which case, the lens would be both more massive and more distant than the best-estimated values above.

The NRAO VLA Sky Survey (NVSS) Rotation Measure (RM) catalogue is invaluable for the study of cosmic magnetism. However, the RM values reported in it can be affected by n$\pi$-ambiguity, resulting in deviations of the reported RM from the true values by multiples of +-652.9 rad m-2. We therefore set off to observationally constrain the fraction of sources in the RM catalogue affected by this ambiguity. New broadband spectro-polarimetric observations were performed with the Karl G. Jansky Very Large Array (VLA) at 1--2 GHz, with 23 n$\pi$-ambiguity candidates selected by their peculiarly high |RM| values. We identified nine sources with erroneous RM values due to n$\pi$-ambiguity and 11 with reliable RM values. In addition, we found two sources to be unpolarised and one source to be inconsistent with neither n$\pi$-ambiguity nor reliable RM cases. By comparing the statistical distributions of the above two main classes, we devised a measure of how much a source's RM deviates from that of its neighbours: $\Delta/\sigma$, which we found to be a good diagnostic of n$\pi$-ambiguity. With this, we estimate that there are at least 50 sources affected by n$\pi$-ambiguity among the 37,543 sources in the catalogue. Finally, we explored the Faraday complexities of our sources revealed by our broadband observations.

The NRAO VLA Sky Survey (NVSS) Rotation Measure (RM) catalogue has enabled numerous studies in cosmic magnetism, and will continue being a unique dataset complementing future polarisation surveys. Robust comparisons with these new surveys will however require further understandings in the systematic effects present in the NVSS RM catalogue. In this paper, we make careful comparisons between our new on-axis broadband observations with the Karl G. Jansky Very Large Array and the NVSS RM results for 23 sources. We found that two unpolarised sources were reported as polarised at about 0.5% level in the RM catalogue, and noted significant differences between our newly derived RM values and the catalogue values for the remaining 21 sources. These discrepancies are attributed to off-axis instrumental polarisation in the NVSS RM catalogue. By adopting the 0.5% above as the typical off-axis instrumental polarisation amplitude, we quantified its effect on the reported RMs with a simulation, and found that on average the RM uncertainties in the catalogue have to be increased by $\approx$ 10% to account for the off-axis instrumental polarisation effect. This effect is more substantial for sources with lower fractional polarisation, and is a function of the source's true RM. Moreover, the distribution of the resulting RM uncertainty is highly non-Gaussian. With the extra RM uncertainty incorporated, we found that the RM values from the two observations for most (18 out of 21) of our polarised targets can be reconciled. The remaining three are interpreted as showing hints of time variabilities in RM.

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

We evaluate the capabilities of Schwarzschild's orbit-superposition method by applying it to galaxies from the large scale, high resolution Illustris simulation. Nine early-type galaxies with a range of triaxiality are selected, and we create mock integral field unit data for five line-of-sight projections of each galaxy. Each of the 45 mock data sets is taken as an independent observed galaxy. Using van den Bosch's 2008 triaxial Schwarzschild implementation, we assess model estimates of various galaxy properties, covering mass profiles, intrinsic shapes, stellar orbit distributions and velocity anisotropies. Total mass within $\overline{R_{\rm e}}$ is recovered well with average deviations within $\pm15$ percent. Stellar mass is underestimated by $\sim24$ percent and dark matter overestimated by $\sim38$ percent (assuming an NFW dark matter profile and allowing for degeneracy between stellar mass and dark matter mass). Using a gNFW profile, these values improve to $\sim13$ percent for stellar mass and $\sim18$ percent for dark matter. Axis ratio estimates show a moderate bias of $\Delta (b/a)=0.07$ and $\Delta (c/a)=0.14$ ($a\ge b\ge c$). Distributions of the orbit circularities $\lambda_z$ and $\lambda_x$, representing rotation about the minor and major axes, are well reconstructed. Separating orbits into thermal categories, our models match the average fractions of these categories to within $10$ percent. Velocity anisotropy is well estimated with values matching in the inner regions but becoming slightly radially biased in the outer regions. Overall, the galaxy property estimates we obtained using Schwarzschild modelling are not implausible and are representative of the simulated galaxies we modelled.

We report the mass and distance measurements of two single-lens events from the 2017 Spitzer microlensing campaign. The ground-based observations yield the detection of finite-source effects, and the microlens parallaxes are derived from the joint analysis of ground-based observations and Spitzer observations. We find that the lens of OGLE-2017-BLG-1254 is a $0.60 \pm 0.03 M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.11~\text{kpc}$, where $D_{\rm LS}$ is the distance between the lens and the source. The second event, OGLE-2017-BLG-1161, is subject to the known satellite parallax degeneracy, and thus is either a $0.51^{+0.12}_{-0.10} M_{\odot}$ star with $D_{\rm LS} = 0.40 \pm 0.12~\text{kpc}$ or a $0.38^{+0.13}_{-0.12} M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.19~\text{kpc}$. Both of the lenses are therefore isolated stars in the Galactic bulge. By comparing the mass and distance distributions of the eight published Spitzer finite-source events with the expectations from a Galactic model, we find that the Spitzer sample is in agreement with the probability of finite-source effects occurrence in single lens events.

We derive the mass weighted total density slopes within the effective (half-light) radius, $\gamma'$, for more than 2000 nearby galaxies from the SDSS-IV MaNGA survey using Jeans-anisotropic-models applied to IFU observations. Our galaxies span a wide range of the stellar mass ($10^9$ $M_{\rm \odot}< M_* < 10^{12}$ M$_{\odot}$) and the velocity dispersion (30 km/s $< \sigma_v <$ 300 km/s). We find that for galaxies with velocity dispersion $\sigma_v>100$ km/s, the density slope has a mean value $\langle \gamma^{\prime} \rangle = 2.24$ and a dispersion $\sigma_{\gamma}=0.22$, almost independent of velocity dispersion. A clear turn over in the $\gamma'-\sigma_v$ relation is present at $\sigma\sim 100$ km/s, below which the density slope decreases rapidly with $\sigma_v$. Our analysis shows that a large fraction of dwarf galaxies (below $M_* = 10^{10}$ M$_{\odot}$) have total density slopes shallower than 1, which implies that they may reside in cold dark matter halos with shallow density slopes. We compare our results with that of galaxies in hydrodynamical simulations of EAGLE, Illustris and IllustrisTNG projects, and find all simulations predict shallower density slopes for massive galaxies with high $\sigma_v$. Finally, we explore the dependence of $\gamma'$ on the positions of galaxies in halos, namely centrals vs. satellites, and find that for the same velocity dispersion, the amplitude of $\gamma'$ is higher for satellite galaxies by about 0.1.

The Wide Field Infrared Survey Telescope (WFIRST) will monitor $\sim 2$ deg$^2$ toward the Galactic bulge in a wide ($\sim 1-2~\mu$m) W149 filter at 15-minute cadence with exposure times of $\sim$50s for 6 seasons of 72 days each, for a total $\sim$41,000 exposures taken over $\sim$432 days, spread over the 5-year prime mission. This will be one of the deepest exposures of the sky ever taken, reaching a photon-noise photometric precision of 0.01 mag per exposure and collecting a total of $\sim 10^9$ photons over the course of the survey for a W149$_{\rm AB}\sim 21$ star. Of order $4 \times 10^7$ stars will be monitored with W149$_{\rm AB}$<21, and 10$^8$ stars with W145$_{\rm AB}$<23. The WFIRST microlensing survey will detect $\sim$54,000 microlensing events, of which roughly 1% ($\sim$500) will be due to isolated black holes, and $\sim$3% ($\sim$1600) will be due to isolated neutron stars. It will be sensitive to (effectively) isolated compact objects with masses as low as the mass of Pluto, thereby enabling a measurement of the compact object mass function over 10 orders of magnitude. Assuming photon-noise limited precision, it will detect $\sim 10^5$ transiting planets with sizes as small as $\sim 2~R_\oplus$, perform asteroseismology of $\sim 10^6$ giant stars, measure the proper motions to $\sim 0.3\%$ and parallaxes to $\sim 10\%$ for the $\sim 6 \times 10^6$ disk and bulge stars in the survey area, and directly detect $\sim 5 \times 10^3$ Trans-Neptunian objects (TNOs) with diameters down to $\sim 10$ km, as well as detect $\sim 10^3$ occulations of stars by TNOs during the survey. All of this science will completely serendipitous, i.e., it will not require modifications of the WFIRST optimal microlensing survey design. Allowing for some minor deviation from the optimal design, such as monitoring the Galactic center, would enable an even broader range of transformational science.

We report observations of the binary microlensing event OGLE-2018-BLG-0022, provided by the ROME/REA Survey, which indicate that the lens is a low-mass binary star consisting of M3 (0.375+/-0.020 Msun) and M7 (0.098+/-0.005 Msun) components. The lens is unusually close, at 0.998+/-0.047 kpc, compared with the majority of microlensing events, and despite its intrinsically low luminosity, it is likely that AO observations in the near future will be able to provide an independent confirmation of the lens masses.

Free-floating planets are the remnants of violent dynamical rearrangements of planetary systems. It is possible that even our own solar system ejected a large planet early in its evolution. WFIRST will have the ability to detect free-floating planets over a wide range of masses, but it will not be able to directly measure their masses. Microlensing parallax observations can be used to measure the masses of isolated objects, including free-floating planets, by observing their microlensing events from two locations. The intra-L2 separation between WFIRST and Euclid is large enough to enable microlensing parallax measurements, especially given the exquisite photometric precision that both spacecraft are capable of over wide fields. In this white paper we describe how a modest investment of observing time could yield hundreds of parallax measurements for WFIRST's bound and free-floating planets. We also describe how a short observing campaign of precursor observations by Euclid can improve WFIRST's bound planet and host star mass measurements.

[abridged] WFIRST is uniquely capable of finding planets with masses as small as Mars at separations comparable to Jupiter, i.e., beyond the current ice lines of their stars. These planets fall between the close-in planets found by Kepler and the wide separation gas giants seen by direct imaging and ice giants inferred from ALMA observations. Furthermore, the smallest planets WFIRST can detect are smaller than the planets probed by RV and Gaia at comparable separations. Interpreting planet populations to infer the underlying formation and evolutionary processes requires combining results from multiple detection methods to measure the full variation of planets as a function of planet size, orbital separation, and host star mass. Microlensing is the only way to find planets from 0.5 to 5M_E at 1 to 5au. The case for a microlensing survey from space has not changed in the past 20 yrs: space allows wide-field diffraction-limited observations that resolve main-sequence stars in the bulge, which allows the detection and characterization of the smallest signals including those from planets with masses at least as small as Mars. What has changed is that ground-based (GB) microlensing is reaching its limits, underscoring the scientific necessity for a space-based survey. GB microlensing has found a break in the mass-ratio distribution at about a Neptune, implying that these are the most common microlensing planet and that planets smaller than this are rare. However, GB microlensing reaches its detection limits only slightly below the observed break. WFIRST will measure the shape of the mass-ratio function below the break by finding numerous smaller planets: 500 Neptunes, 600 gas giants, 200 Earths, and planets as small as 0.1M_E. Because it will also measure host masses and distances, WFIRST will also track the behavior of the planet distribution as a function of separation and host star mass.

Microlensing is the only method that can detect and measure mass of wide orbit, low mass, solar system analog exoplanets. Mass measurements of such planets would yield massive science on planet formation, exoplanet demographics, free floating planets, planet frequencies towards the galaxy. High res follow-up observations of past microlens targets provide a mass measurement of microlens planets and hosts at an uncertainty of <20%. This will be primary method for mass measurement with WFIRST. We advocate for the fact that high resolution observations with AO, HST and JWST(in future) remain necessary in coming decade to develop the methods, to determine the field and filter selection, understand the systematics and to develop a robust pipeline to release high quality data products from WFIRST microlensing survey such that the astronomy community can promptly engage in the science. We also support future high res obs with US ELTs with advanced Laser AO systems in context of enhancing the science return of WFIRST microlensing survey. We endorse the 2018 Exoplanet Science Strategy report published by the National Academy. This white paper extends and complements the material presented therein. In particular, this white paper supports the recommendation of the National Academy Exoplanet Science Strategy report that: NASA should launch WFIRST to conduct its microlensing survey of distant planets and to demonstrate the technique of coronagraphic spectroscopy on exoplanet targets. This white paper also supports to the finding from that report which states "A number of activities, including precursor and concurrent observations using ground- and space-based facilities, would optimize the scientific yield of the WFIRST microlensing survey."

We analyze the stellar initial mass functions (IMF) of a large sample of early type galaxies (ETGs) provided by MaNGA. The large number of IFU spectra of individual galaxies provide high signal-to-noise composite spectra that are essential for constraining IMF and to investigate possible radial gradients of the IMF within individual galaxies. The large sample of ETGs also make it possible to study how the IMF shape depends on various properties of galaxies. We adopt a novel approach to study IMF variations in ETGs, use Bayesian inferences based on full spectrum fitting. The Bayesian method provides a statistically rigorous way to explore potential degeneracy in spectrum fitting, and to distinguish different IMF models with Bayesian evidence. We find that the IMF slope depends systematically on galaxy velocity dispersion, in that galaxies of higher velocity dispersion prefer a more bottom-heavy IMF, but the dependence is almost entirely due to the change of metallicity, $Z$, with velocity dispersion. The IMF shape also depends on stellar age, $A$, but the dependence is completely degenerate with that on metallicity through a combination $AZ^{-1.42}$. Using independent age and metallicity estimates we find that the IMF variation is produced by metallicity instead of age. The IMF near the centers of massive ETGs appears more bottom-heavy than that in the outer parts, while a weak opposite trend is seen for low-mass ETGs. Uncertainties produced by star formation history, dust extinction, $\alpha$-element abundance enhancement and noise in the spectra are tested.

We use a sample of ~3000 galaxies from the MaNGA MPL-7 internal data release to study the alpha abundance distribution within low-redshift galaxies. We use the Lick index ratio Mgb/<Fe> as an alpha abundance indicator to study relationships between the alpha abundance distribution and galaxy properties such as effective stellar velocity dispersion within 0.3 effective radii (sigma_*), galaxy environment, and dark matter halo formation time (z_f). We find that (1) all galaxies show a tight correlation between Mgb/<Fe> and sigma_*; (2) `old' (H_beta < 3) low-sigma_* galaxies in high local density environment and inner regions within galaxy groups are enhanced in Mgb/<Fe>, while `young' (H_beta>3) galaxies and high-mass galaxies show no or less environmental dependence; (3) `old' galaxies with high-z_f show enhanced Mgb/<Fe> over low- and medium-z_f; (4) Mgb/<Fe> gradients are close to zero and show dependence on sigma_* but no obvious dependence on the environment or z_f. Our study indicates that stellar velocity dispersion or galaxy mass is the main parameter driving the Mgb/<Fe> enhancement, although environments appear to have modest effects, particularly for low- and medium-mass galaxies.

Three basic ingredients are required to generate a simple stellar population (SSP) library, i.e., an initial mass function (IMF), a stellar evolution model/isochrones, and an empirical/theoretical stellar spectral library. However, there are still some uncertainties to the determination and understanding of these ingredients. We perform the spectral fitting to test the relative parameter offsets between these uncertainties using two different stellar population models, two different empirical stellar libraries, two different isochrones, and the Salpeter and Chabrier IMFs. Based on these setups, we select five SSP libraries generated with the Galaxev/STELIB and Vazdekis/MILES models, and apply them to the pPXF full-spectrum fitting of both MaNGA and mock spectra. We find that: 1) Compared to the Galaxev/STELIB model, spectral fitting qualities with the Vazdekis/MILES model have significant improvements for those metal-rich (especially over-solar) spectra, which cause better reduced $\chi^2$ distributions and more precisely fitted absorption lines. This might due to the lack of metal rich stars in the empirical STELIB library, or code improvement of the Vazdekis model. 2) When applying the Vazdekis/MILES model for spectral fitting, the IMF variation will lead to not only a systematic offset in $M_*/L_r$, but also offsets in age and metallicity, and these offsets increase with increasing stellar population ages. However, the IMF-variation caused metallicity offsets disappear in the case of Galaxev/STELIB based libraries. 3) The Padova2000 model provides a better match to the MaNGA galaxy spectra at [M/H]$_L<-1.0$, while the BaSTI model match the local galaxy spectra better at [M/H]$_L>-1.0$. Current tests suggest that spectral fitting with the Vazdekis/MILES+BaSTI combination would be a better choice for local galaxies.

With about 12 000 OB type stars selected from the LAMOST and Gaia survey, we study their 3 dimensional velocity distribution over the range of galactocentric radius from 6 to 15 kpc in the Galactic disk plane. A clear ripple pattern in the radial velocity ($V_R$) map is shown. The median $V_R$ reaches $-8$ km s$^{-1}$ at $R\sim9$ kpc, then increases to $\sim0$ km s$^{-1}$ at $R\sim12$ kpc, and later declines to below $-10$ km s$^{-1}$ beyond $R\sim13$ kpc. The median azimuthal velocity ($V_\phi$) map shows a similar pattern but has roughly $1/4$ phase difference with the radial velocity. Although the ripple of negative $V_R$ at $\sim9$ kpc extends to about 40$^\circ$ in the azimuth angle, it does not align with either the Local or the Perseus spiral arm. Moreover, the farther ripple beyond 13 kpc does not match the Outer spiral arm either. This indicates that the non-axisymmetric kinematic features are not induced by perturbations of known spiral structures. The central rotating bar can not lead to such patterns in the outer disk either. External perturbation of a dwarf galaxy or a dark matter sub-halo can induce such patterns but requires more evidence from both observations and simulations. The $V_\phi$ map in the $Z$--$V_Z$ plane of the OB stars is also investigated. Despite asymmetry to some degree, no spiral pattern is found. This is reasonable since most of the OB stars have ages much younger than 100 Myrs, which is smaller than one orbital period around the Galactic center.

We report the discovery of a low mass-ratio planet $(q = 1.3\times10^{-4})$, i.e., 2.5 times higher than the Neptune/Sun ratio. The planetary system was discovered from the analysis of the KMT-2017-BLG-0165 microlensing event, which has an obvious short-term deviation from the underlying light curve produced by the host of the planet. Although the fit improvement with the microlens parallax effect is relatively low, one component of the parallax vector is strongly constrained from the light curve, making it possible to narrow down the uncertainties of the lens physical properties. A Bayesian analysis yields that the planet has a super-Neptune mass $(M_{2}=34_{-12}^{+15}~M_{\oplus})$ orbiting a Sun-like star $(M_{1}=0.76_{-0.27}^{+0.34}~M_{\odot})$ located at $4.5~{\rm kpc}$. The blended light is consistent with these host properties. The projected planet-host separation is $a_{\bot}={3.45_{-0.95}^{+0.98}}~{\rm AU}$, implying that the planet is located outside the snowline of the host, i.e., $a_{sl}\sim2.1~{\rm AU}$. KMT-2017-BLG-0165Lb is the sixteenth microlensing planet with mass ratio $q<3\times10^{-4}$. Using the fifteen of these planets with unambiguous mass-ratio measurements, we apply a likelihood analysis to investigate the form of the mass-ratio function in this regime. If we adopt a broken power law for the form of this function, then the break is at $q_{\rm br}\simeq0.55\times10^{-4}$, which is much lower than previously estimated. Moreover, the change of the power law slope, $\zeta>3.3$ is quite severe. Alternatively, the distribution is also suggestive of a "pile-up" of planets at Neptune-like mass ratios, below which there is a dramatic drop in frequency.

The Smith Cloud is a high velocity cloud (HVC) with an orbit suggesting it has made at least one passage through the Milky Way disk. A magnetic field found around this cloud has been thought to provide extra stability as it passes through the Galactic halo. We use the Karl G. Jansky Very Large Array to measure Faraday rotation measures (RMs) towards 1105 extragalactic background point sources behind and next to the Smith Cloud to constrain the detailed geometry and strength of its magnetic field. The RM pattern across the cloud gives a detailed morphology of the magnetic field structure which indicates a field draped over the ionized gas and compressed at the head of the cloud. We constrain the peak line-of-sight magnetic field strength to > +5 uG and create a model of the magnetic field to demonstrate that a draped configuration can qualitatively explain the morphology of the observed RMs.

Faraday rotation occurs along every line of sight in the Galaxy; Rotation Measure (RM) synthesis allows a three-dimensional representation of the interstellar magnetic field. This study uses data from the Global Magneto-Ionic Medium Survey, a combination of single-antenna spectro-polarimetric studies, including northern sky data from the DRAO 26-m Telescope (1270-1750 MHz) and southern sky data from the Parkes 64-m Telescope (300-480 MHz). From the synthesized Faraday spectral cubes we compute the zeroth, first, and second moments to find the total polarized emission, mean and RM-width of the polarized emission. From DRAO first moments we find a weak vertical field directed from Galactic North to South, but Parkes data reveal fields directed towards the Sun at high latitudes in both hemispheres: the two surveys clearly sample different volumes. DRAO second moments show feature widths in Faraday spectra increasing with decreasing positive latitudes, implying that longer lines of sight encounter more Faraday rotating medium, but this is not seen at negative latitudes. Parkes data show the opposite: at positive latitudes the second moment decreases with decreasing latitude, but not at negative latitudes. Comparing first moments with RMs of pulsars and extragalactic sources and a study of depolarization together confirm that the DRAO survey samples to larger distances than the Parkes data. Emission regions in the DRAO survey are typically 700 to 1000 pc away, slightly beyond the scale-height of the magneto-ionic medium; emission detected in the Parkes survey is entirely within the magneto-ionic disk, less than 500 pc away.

We explore orbit properties of 35 prolate-triaxial galaxies selected from the Illustris cosmological hydrodynamic simulation. We present a detailed study of their orbit families, and also analyse relations between the relative abundance of the orbit families and the spin parameter, triaxiality, the ratio of the angular momentum and the baryon fraction. We find that box orbits dominate the orbit structure for most prolate-triaxial galaxies, especially in the central region. The fraction of irregular orbits in the prolate-triaxial galaxies is small, and it increases with galaxy radius. Both the x-tube and z-tube orbits are important in prolate-triaxial galaxies, especially in the outer regions. The fraction of box orbits for prolate-triaxial galaxies ($0.7<T<1$) decreases with the triaxiality of the stars, while the fraction of x-tube orbits increases with the triaxiality for a given axis ratio. The fraction of box orbits increases and the fractions of x-tube and z-tube orbits weakly decrease with increasing baryon fraction. These results help to understand the structure of prolate-triaxial galaxies and provide cross-checks for constructing dynamical models of prolate-triaxial galaxies by using the Schwarzschild or the Made-to-Measure methods.

Cosmological simulations predict that an intergalactic magnetic field (IGMF) pervades the large scale structure (LSS) of the Universe. Measuring the IGMF is important to determine its origin (i.e. primordial or otherwise). Using data from the LOFAR Two Metre Sky Survey (LoTSS), we present the Faraday rotation measure (RM) and depolarisation properties of the giant radio galaxy J1235+5317, at a redshift of $z = 0.34$ and 3.38 Mpc in size. We find a mean RM difference between the lobes of $2.5\pm0.1$ rad/m$^2$ , in addition to small scale RM variations of ~0.1 rad/m$^2$ . From a catalogue of LSS filaments based on optical spectroscopic observations in the local universe, we find an excess of filaments intersecting the line of sight to only one of the lobes. Associating the entire RM difference to these LSS filaments leads to a gas density-weighted IGMF strength of ~0.3 {\mu}G. However, direct comparison with cosmological simulations of the RM contribution from LSS filaments gives a low probability (~5%) for an RM contribution as large as 2.5 rad/m$^2$ , for the case of IGMF strengths of 10 to 50 nG. It is likely that variations in the RM from the Milky Way (on 11' scales) contribute significantly to the mean RM difference, and a denser RM grid is required to better constrain this contribution. In general, this work demonstrates the potential of the LOFAR telescope to probe the weak signature of the IGMF. Future studies, with thousands of sources with high accuracy RMs from LoTSS, will enable more stringent constraints on the nature of the IGMF.

We present an analysis of the microlensing event OGLE-2015-BLG-0232. This event is challenging to characterize for two reasons. First, the light curve is not well sampled during the caustic crossing due to the proximity of the full Moon impacting the photometry quality. Moreover, the source brightness is difficult to estimate because this event is blended with a nearby K dwarf star. We found that the light curve deviations are likely due to a close brown dwarf companion (i.e., s = 0.55 and q = 0.06), but the exact nature of the lens is still unknown. We finally discuss the potential of follow-up observations to estimate the lens mass and distance in the future.

We report the discovery of a giant planet in the KMT-2016-BLG-1397 microlensing event, which was found by The Korea Microlensing Telescope Network (KMTNet) alone. The time scale of this event is t_E = 40.0 +- 0.5 days and the mass ratio between the lens star and its companion is q = 0.016 +- 0.002. The planetary perturbation in the light curve is a smooth bump, resulting in the classical binary-lens/binary-source (2L1S/1L2S) degeneracy. We measure the V - I color of the (putative) two sources in the 1L2S model, and then effectively rule out the binary source solution. The finite-source effect is marginally detected. Combined with the limits on the blend flux and the probability distribution of the source size normalized by the Einstein radius rho, a Bayesian analysis yields the lens mass M_L = 0.45+0.33-0.28 M_sun, at distance of D_L = 6.60+1.10-1.30 kpc. Thus the companion is a super-Jupiter of a mass m_p = 7.0+5.2-4.3 M_J , at a projected separation r = 5.1+1.5-1.7 AU, indicating that the planet is well beyond the snow line of the host star.

We extend Schwarzschild's dynamical modelling method to model absorption line strength data as well as the more usual luminosity and kinematic data. Our approach draws on earlier published work by the first author with the Syer and Tremaine made-to-measure (M2M) dynamical modelling method and uses similar ideas to create a chemo-Schwarzschild method. We apply our extended Schwarzschild method to the same four early type galaxies (NGC 1248, NGC 3838, NGC 4452, NGC 4551) as the chemo-M2M work, and are able to recover successfully the 2D absorption line strength for the three lines we model (Hbeta, Fe5015, Mgb). We believe that this is the first time Schwarzschild's method has been used in this way. The techniques developed can be applied to modelling other aspects of galaxies, for example age and metallicity data coming from stellar population modelling, and are not specific to absorption line strength data.

An excellent laboratory for studying large scale magnetic fields is the grand de- sign face-on spiral galaxy M51. Due to wavelength-dependent Faraday depolarization, linearly polarized synchrotron emission at different radio frequencies gives a picture of the galaxy at different depths: Observations at L-band (1-2 GHz) probe the halo region while at C- and X- band (4-8 GHz) the linearly polarized emission probe the disk region of M51. We present new observations of M51 using the Karl G. Jansky Very Large Array (VLA) at S-band (2-4 GHz), where previously no polarization observations existed, to shed new light on the transition region between the disk and the halo. We discuss a model of the depolarization of synchrotron radiation in a multilayer magneto-ionic medium and compare the model predictions to the multi-frequency polarization data of M51 between 1-8GHz. The new S-band data are essential to distinguish between different models. Our study shows that the initial model parameters, i.e. the total reg- ular and turbulent magnetic field strengths in the disk and halo of M51, need to be adjusted to successfully fit the models to the data.

The unprecedented sensitivity, angular resolution and broad bandwidth coverage of Square Kilometre Array (SKA) radio polarimetric observations will allow us to address many long-standing mysteries in cosmic magnetism science. I will highlight the unique capabilities of the SKA to map the warm hot intergalactic medium, reveal detailed 3-dimensional structures of magnetic fields in local galaxies and trace the redshift evolution of galactic magnetic fields.

With rotation measure (RM) towards 37,543 polarised sources, the Taylor et al. (2009) RM catalogue has been widely exploited in studies of the foreground magneto-ionic media. However, due to limitations imposed by observations in survey mode in the narrowband era, the listed RM values are inevitably affected by various systematic effects. With new Karl G. Jansky Very Large Array (VLA) broadband spectro-polarimetric observations at L-band, we set off to observationally examine the robustness of the Taylor catalogue. This would facilitate combinations and comparisons of it with results from current and future polarisation surveys such as Polarization Sky Survey of the Universe's Magnetism (POSSUM), VLA Sky Survey (VLASS), and the eventual Square Kilometre Array (SKA). Our on-axis pointed observations, in conjunction with simulations, allowed us to estimate the impact of off-axis polarisation leakage on the measured RM values. This demonstrates the importance to properly calibrate for the off-axis leakage terms in future all-sky polarisation surveys, in order to obtain high fidelity polarisation information from sources down to low fractional polarisation.

To constrain the large-scale magnetic field strengths in cosmologically distant galax- ies, we derive the probability distribution function of Faraday rotation measure (RM) when random lines of sight pass through a sample of disc galaxies, with axisymmetric large-scale magnetic fields. We find that the width of the RM distribution of the galaxy sample is directly related to the mean large-scale field strength of the galaxy population, provided the dispersion within the sample is lower than the mean value. In the absence of additional constraints on parameters describing the magneto-ionic medium of the intervening galaxies, and in the situation where RMs produced in the intervening galaxies have already been statistically isolated from other RM contributions along the lines of sight, our simple model of the magneto-ionic medium in disc galaxies suggests that the mean large-scale magnetic field of the population can be measured to within ~ 50% accuracy.

The MaNGA project has obtained IFU data for several thousand nearby galaxies, including barred galaxies. With the two dimensional spectral and kinematic information provided by IFUs, we can measure the pattern speed of a barred galaxy, which determines the bar dynamics. We apply the non-parametric method proposed by Tremaine \& Weinberg to estimate the bar pattern speed for 53 barred galaxies, making this the largest sample studied so far in this way. Our sample is selected from the MaNGA first public data release as part of SDSS Data Release 13 according mainly to the axis ratio and position angle difference between the bar and disc, while kinematic data is from the later SDSS Data Release 14. We have used both the photometric position angle from the photometric image and the kinematic position angle from the stellar velocity map to derive the pattern speed. Combining three independent bar length measurements and the circular velocity from Jeans Anisotropic modelling (JAM), we also determine the dimensionless ratio $\cal{R}$ of the corotation radius to the bar length. We find that the galaxy's position angle is the main uncertainty in determining the bar pattern speed. The kinematic position angle leads to fewer ultrafast bars than the photometric position angle, and this could be due to the method of measuring the kinematic position angle. We study the dependence of $\cal{R}$ values on galaxy properties such as the dark matter fraction from JAM modelling and the stellar age and metallicity from stellar population synthesis (SPS). A positive correlation between the bar length and bar strength is found: the longer the bar, the stronger the bar. However, no other significant correlations are found. This may result from errors in deriving the $\cal{R}$ values or from the complex formation and slowdown processes of galactic bars.

Gravitational waves emitted from compact binary coalescence can be subject to wave diffraction if they are gravitationally lensed by an intervening mass clump whose Schwarzschild timescale matches the wave period. Waves in the ground-based frequency band $f\sim 10$--$10^3\,$Hz are sensitive to clumps with masses $M_E \sim 10^2$--$10^3\,M_\odot$ enclosed within the impact parameter. These can be the central parts of low mass $M_L \sim 10^3$--$10^6\,M_\odot$ dark matter halos, which are predicted in Cold Dark Matter scenarios but are challenging to observe. Neglecting finely-tuned impact parameters, we focus on lenses aligned generally on the Einstein scale for which multiple lensed images may not form in the case of an extended lens. In this case, diffraction induces amplitude and phase modulations whose sizes $\sim 10\%$--$20\%$ are small enough so that standard matched filtering with unlensed waveforms do not degrade, but are still detectable for events with high signal-to-noise ratio. We develop and test an agnostic detection method based on dynamic programming, which does not require a detailed model of the lensed waveforms. For pseudo-Jaffe lenses aligned up to the Einstein radius, we demonstrate that a pair of fully upgraded aLIGO/Virgo detectors can extract diffraction imprints from binary black hole mergers out to $z_s \sim 0.2$--$0.3$. The prospect will improve dramatically for a third-generation detector for which binary black hole mergers out to $z_s \sim 2$--$4$ will all become valuable sources.

Data from the newly-commissioned \textit{Transiting Exoplanet Survey Satellite} (TESS) has revealed a "hot Earth" around LHS 3844, an M dwarf located 15 pc away. The planet has a radius of $1.32\pm 0.02$ $R_\oplus$ and orbits the star every 11 hours. Although the existence of an atmosphere around such a strongly irradiated planet is questionable, the star is bright enough ($I=11.9$, $K=9.1$) for this possibility to be investigated with transit and occultation spectroscopy. The star's brightness and the planet's short period will also facilitate the measurement of the planet's mass through Doppler spectroscopy.

We propose a new strategy of finding strongly-lensed supernovae (SNe) by monitoring known galaxy-scale strong-lens systems. Strongly lensed SNe are potentially powerful tools for the study of cosmology, galaxy evolution, and stellar populations, but they are extremely rare. By targeting known strongly lensed starforming galaxies, our strategy significantly boosts the detection efficiency for lensed SNe compared to a blind search. As a reference sample, we compile the 128 galaxy-galaxy strong-lens systems from the Sloan Lens ACS Survey (SLACS), the SLACS for the Masses Survey, and the Baryon Oscillation Spectroscopic Survey Emission-Line Lens Survey. Within this sample, we estimate the rates of strongly-lensed Type Ia SN (SNIa) and core-collapse SN (CCSN) to be $1.23 \pm 0.12$ and $10.4 \pm 1.1$ events per year, respectively. The lensed SN images are expected to be widely separated with a median separation of 2 arcsec. Assuming a conservative fiducial lensing magnification factor of 5 for the most highly magnified SN image, we forecast that a monitoring program with a single-visit depth of 24.7 mag (5$\sigma$ point source, $r$ band) and a cadence of 5 days can detect 0.49 strongly-lensed SNIa event and 2.1 strongly-lensed CCSN events per year within this sample. Our proposed targeted-search strategy is particularly useful for prompt and efficient identifications and follow-up observations of strongly-lensed SN candidates. It also allows telescopes with small field of views and limited time to efficiently discover strongly-lensed SNe with a pencil-beam scanning strategy.

The Beijing-Arizona Sky Survey (BASS) is a new wide-field legacy imaging survey in the northern Galactic cap using the 2.3m Bok telescope. The survey will cover about 5400 deg$^2$ in the $g$ and $r$ bands, and the expected 5$\sigma$ depths (corrected for the Galactic extinction) in the two bands are 24.0 and 23.4 mag, respectively. BASS started observations in January 2015, and has completed about 41% of the whole area as of July 2016. The first data release contains both calibrated images and photometric catalogs obtained in 2015 and 2016. The depths of single-epoch images in the two bands are 23.4 and 22.9 mag, and the full depths of three epochs are about 24.1 and 23.5 mag, respectively.

The Wide Field InfraRed Survey Telescope (WFIRST) is the next NASA astrophysics flagship mission, to follow the James Webb Space Telescope (JWST). The WFIRST mission was chosen as the top-priority large space mission of the 2010 astronomy and astrophysics decadal survey in order to achieve three primary goals: to study dark energy via a wide-field imaging survey, to study exoplanets via a microlensing survey, and to enable a guest observer program. Here we assess the ability of the several WFIRST designs to achieve the goal of the microlensing survey to discover a large sample of cold, low-mass exoplanets with semimajor axes beyond roughly one AU, which are largely impossible to detect with any other technique. We present the results of a suite of simulations that span the full range of the proposed WFIRST architectures, from the original design envisioned by the decadal survey, to the current design, which utilizes a 2.4-m telescope donated to NASA. By studying such a broad range of architectures, we are able to determine the impact of design trades on the expected yields of detected exoplanets. In estimating the yields we take particular care to ensure that our assumed Galactic model predicts microlensing event rates that match observations, consider the impact that inaccuracies in the Galactic model might have on the yields, and ensure that numerical errors in lightcurve computations do not bias the yields for the smallest mass exoplanets. For the nominal baseline WFIRST design and a fiducial planet mass function, we predict that a total of ${\sim}1400$ bound exoplanets with mass greater than ${\sim}0.1~M_{\oplus}$ should be detected, including ${\sim}200$ with mass ${\lesssim}3~M_{\oplus}$. WFIRST should have sensitivity to planets with mass down to ${\sim}0.02~M_{\oplus}$, or roughly the mass of Ganymede.

By inverting the distributions of galaxies' apparent ellipticities and misalignment angles (measured around the projected half-light radius $R_{\rm e}$) between their photometric and kinematic axes, we study the intrinsic shape distribution of 189 slow rotator early-type galaxies with stellar masses $2\times 10^{11} M_{\odot}<M_\ast<2\times 10^{12} M_{\odot}$, extracted from a sample of about 2200 galaxies with integral-field stellar kinematics from the DR14 of the SDSS-IV MaNGA IFU survey. Thanks to the large sample of slow rotators, Graham+18 showed that there is clear structure in the misalignment angle distribution, with two peaks at both $0^{\circ}$ and $90^{\circ}$ misalignment (characteristic of oblate and prolate rotation respectively). Here we invert the observed distribution from Graham+18. The large sample allows us to go beyond the known fact that slow rotators are weakly triaxial and to place useful constraints on their intrinsic triaxiality distribution (around $1R_{\rm e}$) for the first time. The shape inversion is generally non-unique. However, we find that, for a wide set of model assumptions, the observed distribution clearly requires a dominant triaxial-oblate population. For some of our models, the data suggest a hint for a minor triaxial-prolate population, but a dominant prolate population is ruled out.

Using mock spectra based on Vazdekis/MILES library fitted within the wavelength region 3600-7350\AA, we analyze the bias and scatter on the resulting physical parameters induced by the choice of fitting algorithms and observational uncertainties, but avoid effects of those model uncertainties. We consider two full-spectrum fitting codes: pPXF and STARLIGHT, in fitting for stellar population age, metallicity, mass-to-light ratio, and dust extinction. With pPXF we find that both the bias in the population parameters and the scatter in the recovered logarithmic values follows the expected trend. The bias increases for younger ages and systematically makes recovered ages older, $M_*/L_r$ larger and metallicities lower than the true values. For reference, at S/N=30, and for the worst case ($t=10^8$yr), the bias is 0.06 dex in $M_*/L_r$, 0.03 dex in both age and [M/H]. There is no significant dependence on either E(B-V) or the shape of the error spectrum. Moreover, the results are consistent for both our 1-SSP and 2-SSP tests. With the STARLIGHT algorithm, we find trends similar to pPXF, when the input E(B-V)<0.2 mag. However, with larger input E(B-V), the biases of the output parameter do not converge to zero even at the highest S/N and are strongly affected by the shape of the error spectra. This effect is particularly dramatic for youngest age, for which all population parameters can be strongly different from the input values, with significantly underestimated dust extinction and [M/H], and larger ages and $M_*/L_r$. Results degrade when moving from our 1-SSP to the 2-SSP tests. The STARLIGHT convergence to the true values can be improved by increasing Markov Chains and annealing loops to the "slow mode". For the same input spectrum, pPXF is about two order of magnitudes faster than STARLIGHT's "default mode" and about three order of magnitude faster than STARLIGHT's "slow mode".

The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefront controlled technology program, and two targeted science programs: a program to study dark energy, and a statistical census of exoplanets with a microlensing survey, which uses nearly one quarter of WFIRST's observing time in the current design reference mission. The New Worlds, New Horizons (decadal survey) midterm assessment summarizes the science case for the WFIRST exoplanet microlensing survey with this statement: "WFIRST's microlensing census of planets beyond 1 AU will perfectly complement Kepler's census of compact systems, and WFIRST will also be able to detect free-floating planets unbound from their parent stars\rlap."

K2 Campaign 9 (K2C9) was the first space-based microlensing parallax survey capable of measuring microlensing parallaxes of free-floating planet candidate microlensing events. Simultaneous to K2C9 observations we conducted the K2C9 Canada-France-Hawaii Telescope Multi-Color Microlensing Survey (K2C9-CFHT MCMS) in order to measure the colors of microlensing source stars to improve the accuracy of K2C9's parallax measurements. We describe the difference imaging photometry analysis of the K2C9-CFHT MCMS observations, and present the project's first data release. This includes instrumental difference flux lightcurves of 217 microlensing events identified by other microlensing surveys, reference image photometry calibrated to PanSTARRS data release 1 photometry, and tools to convert between instrumental and calibrated flux scales. We derive accurate analytic transformations between the PanSTARRS bandpasses and the Kepler bandpass, as well as angular diameter-color relations in the PanSTARRS bandpasses. To demonstrate the use of our data set, we analyze ground-based and K2 data of a short timescale microlensing event, OGLE-2016-BLG-0795. We find the event has a timescale $t_{\rm E}=4.5 \pm 0.1$~days and microlens parallax $\pi_{\rm E}=0.12 \pm 0.03$ or $0.97 \pm 0.04$, subject to the standard satellite parallax degeneracy. We argue that the smaller value of the parallax is more likely, which implies that the lens is likely a stellar-mass object in the Galactic bulge as opposed to a super-Jupiter mass object in the Galactic disk.

Microlensing can access planet populations that no other method can probe: cold wide-orbit planets beyond the snow line, planets in both the Galactic bulge and disk, and free floating planets (FFPs). The demographics of each population will provide unique constraints on planet formation. Over the past 5 years, U.S. microlensing campaigns with Spitzer and UKIRT have provided a powerful complement to international ground-based microlensing surveys, with major breakthroughs in parallax measurements and probing new regions of the Galaxy. The scientific vitality of these projects has also promoted the development of the U.S. microlensing community. In the 2020s, the U.S. can continue to play a major role in ground-based microlensing by leveraging U.S. assets to complement ongoing ground-based international surveys. LSST and UKIRT microlensing surveys would probe vast regions of the Galaxy, where planets form under drastically different conditions. Moreover, while ground-based surveys will measure the planet mass-ratio function beyond the snow line, adaptive optics (AO) observations with ELTs would turn all of these mass ratios into masses and also distinguish between very wide-orbit planets and genuine FFPs. To the extent possible, cooperation of U.S. scientists with international surveys should also be encouraged and supported.

Deriving the Faraday rotation measure (RM) of quasar absorption line systems, which are tracers of high-redshift galaxies intervening background quasars, is a powerful tool for probing magnetic fields in distant galaxies. Statistically comparing the RM distributions of two quasar samples, with and without absorption line systems, allows one to infer magnetic field properties of the intervening galaxy population. Here, we have derived the analytical form of the probability distribution function (PDF) of RM produced by a single galaxy with an axisymmetric large-scale magnetic field. We then further determine the PDF of RM for one random sight line traversing each galaxy in a population with a large-scale magnetic field prescription. We find that the resulting PDF of RM is dominated by a Lorentzian with a width that is directly related to the mean axisymmetric large-scale field strength $\langle B_0 \rangle$ of the galaxy population if the dispersion of $B_0$ within the population is smaller than $\langle B_0 \rangle$. Provided that RMs produced by the intervening galaxies have been successfully isolated from other RM contributions along the line of sight, our simple model suggests that $\langle B_0 \rangle$ in galaxies probed by quasar absorption line systems can be measured within $\approx50$ per cent accuracy without additional constraints on the magneto-ionic medium properties of the galaxies. Finally, we discuss quasar sample selection criteria that are crucial to reliably interpret observations, and argue that within the limitations of the current database of absorption line systems, high-metallicity damped Lyman-$\alpha$ absorbers are best suited to study galactic dynamo action in distant disc galaxies.

In this paper, we study the strong gravitational lensing of gravitational waves (GWs) from a statistical perspective, with particular focus on the high frequency GWs from stellar binary black hole coalescences. These are most promising targets for ground-based detectors such as Advanced Laser Interferometer Gravitational Wave Observatory (aLIGO) and the proposed Einstein Telescope (ET) and can be safely treated under the geometrical optics limit for GW propagation. We perform a thorough calculation of the lensing rate, by taking account of effects caused by the ellipticity of lensing galaxies, lens environments, and magnification bias. We find that in certain GW source rate scenarios, we should be able to observe strongly lensed GW events once per year ($\sim1~\text{yr}^{-1}$) in the aLIGO survey at its design sensitivity; for the proposed ET survey, the rate could be as high as $\sim80~\text{yr}^{-1}$. These results depend on the estimate of GW source abundance, and hence can be correspondingly modified with an improvement in our understanding of the merger rate of stellar binary black holes. We also compute the fraction of four-image lens systems in each survey, predicting it to be $\sim30$ per cent for the aLIGO survey and $\sim6$ per cent for the ET survey. Finally, we evaluate the possibility of missing some images due to the finite survey duration, by presenting the probability distribution of lensing time delays. We predict that this selection bias will be insignificant in future GW surveys, as most of the lens systems ($\sim90$ per cent) will have time delays less than $\sim1$ month, which will be far shorter than survey durations.

We analyze the binary gravitational microlensing event OGLE-2017-BLG-1130 (mass ratio q~0.45), the first published case in which the binary anomaly was only detected by the Spitzer Space Telescope. This event provides strong evidence that some binary signals can be missed by observations from the ground alone but detected by Spitzer. We therefore invert the normal procedure, first finding the lens parameters by fitting the space-based data and then measuring the microlensing parallax using ground-based observations. We also show that the normal four-fold space-based degeneracy in the single-lens case can become a weak eight-fold degeneracy in binary-lens events. Although this degeneracy is resolved in event OGLE-2017-BLG-1130, it might persist in other events.

We measure $\lambda_{R_e}$, a proxy for galaxy specific stellar angular momentum within one effective radius, and the ellipticity, $\epsilon$, for about 2300 galaxies of all morphological types observed with integral field spectroscopy as part of the MaNGA survey, the largest such sample to date. We use the $(\lambda_{R_e}, \epsilon)$ diagram to separate early-type galaxies into fast and slow rotators. We also visually classify each galaxy according to its optical morphology and two-dimensional stellar velocity field. Comparing these classifications to quantitative $\lambda_{R_e}$ measurements reveals tight relationships between angular momentum and galaxy structure. In order to account for atmospheric seeing, we use realistic models of galaxy kinematics to derive a general approximate analytic correction for $\lambda_{R_e}$. Thanks to the size of the sample and the large number of massive galaxies, we unambiguously detect a clear bimodality in the $(\lambda_{R_e}, \epsilon)$ diagram which may result from fundamental differences in galaxy assembly history. There is a sharp secondary density peak inside the region of the diagram with low $\lambda_{R_e}$ and $\epsilon < 0.4$, previously suggested as the definition for slow rotators. Most of these galaxies are visually classified as non-regular rotators and have high velocity dispersion. The intrinsic bimodality must be stronger, as it tends to be smoothed by noise and inclination. The large sample of slow rotators allows us for the first time to unveil a secondary peak at +/-90 degrees in their distribution of the misalignments between the photometric and kinematic position angles. We confirm that genuine slow rotators start appearing above a stellar mass of 2\times10^{11} M_{\odot}$ where a significant number of high-mass fast rotators also exist.

We perform full spectrum fitting stellar population analysis and Jeans Anisotropic modelling (JAM) of the stellar kinematics for about 2000 early-type galaxies (ETGs) and spiral galaxies from the MaNGA DR14 sample. Galaxies with different morphologies are found to be located on a remarkably tight mass plane which is close to the prediction of the virial theorem, extending previous results for ETGs. By examining an inclined projection (`the mass-size' plane), we find that spiral and early-type galaxies occupy different regions on the plane, and their stellar population properties (i.e. age, metallicity and stellar mass-to-light ratio) vary systematically along roughly the direction of velocity dispersion, which is a proxy for the bulge fraction. Galaxies with higher velocity dispersions have typically older ages, larger stellar mass-to-light ratios and are more metal rich, which indicates that galaxies increase their bulge fractions as their stellar populations age and become enriched chemically. The age and stellar mass-to-light ratio gradients for low-mass galaxies in our sample tend to be positive ($\rm centre<outer$), while the gradients for most massive galaxies are negative. The metallicity gradients show a clear peak around velocity dispersion $\log_{10} \sigma_{\rm e}\approx 2.0$, which corresponds to the critical mass $\sim 3\times 10^{10}M_{\odot}$ of the break in the mass-size relation. Spiral galaxies with large mass and size have the steepest gradients, while the most massive ETGs, especially above the critical mass $M_{\rm crit}\ge 2\times 10^{11} M_{\odot}$, where slow rotator ETGs start dominating, have much flatter gradients. This may be due to differences in their evolution histories, e.g. mergers.

Cosmic ray pressure gradients transfer energy and momentum to extraplanar gas in disk galaxies, potentially driving significant mass loss as galactic winds. This may be particularly important for launching high-velocity outflows of "cool" (T < 10^4 K) gas. We study cosmic-ray driven disk winds using a simplified semi-analytic model assuming streamlines follow the large-scale gravitational potential gradient. We consider scaled Milky Way-like potentials including a disk, bulge, and halo with a range of halo velocities V_H = 50-300 km/s, and streamline footpoints with radii in the disk R_0=1-16 kpc at height 1 kpc. Our solutions cover a wide range of footpoint gas velocity u_0, magnetic-to-cosmic-ray pressure ratio, gas-to-cosmic-ray pressure ratio, and angular momentum. Cosmic ray streaming at the Alfv\'en speed enables the effective sound speed C_eff to increase from the footpoint to a critical point where C_eff,c = u_c ~ V_H; this differs from thermal winds in which C_eff decreases outward. The critical point is typically at a height of 1-6 kpc from the disk, increasing with V_H, and the asymptotic wind velocity exceeds the escape speed of the halo. Mass loss rates are insensitive to the footpoint values of the magnetic field and angular momentum. In addition to numerical parameter space exploration, we develop and compare to analytic scaling relations. We show that winds have mass loss rates per unit area up to ~ Pi_0 V_H^-5/3 u_0^2/3 where Pi_0 is the footpoint cosmic ray pressure and u_0 is set by the upwelling of galactic fountains. The predicted wind mass-loss rate exceeds the star formation rate for V_H < 200 km/s and u_0 = 50 km/s, a typical fountain velocity.

Galaxy formation entails the hierarchical assembly of mass, along with the condensation of baryons and the ensuing, self-regulating star formation. The stars form a collisionless system whose orbit distribution retains dynamical memory that can constrain a galaxy's formation history. The ordered-rotation dominated orbits with near maximum circularity $\lambda_z \simeq1$ and the random-motion dominated orbits with low circularity $\lambda_z \simeq0$ are called kinematically cold and kinematically hot, respectively. The fraction of stars on `cold' orbits, compared to the fraction of stars on `hot' orbits, speaks directly to the quiescence or violence of the galaxies' formation histories. Here we present such orbit distributions, derived from stellar kinematic maps via orbit-based modelling for a well defined, large sample of 300 nearby galaxies. The sample, drawn from the CALIFA survey, includes the main morphological galaxy types and spans the total stellar mass range from $10^{8.7}$ to $10^{11.9}$ solar masses. Our analysis derives the orbit-circularity distribution as a function of galaxy mass, $p(\lambda_z~|~M_\star)$, and its volume-averaged total distribution, $p(\lambda_z)$. We find that across most of the considered mass range and across morphological types, there are more stars on `warm' orbits defined as $0.25\le \lambda_z \le 0.8$ than on either `cold' or `hot' orbits. This orbit-based "Hubble diagram" provides a benchmark for galaxy formation simulations in a cosmological context.

Polarised radio emission from PSR J1745-2900 has already been used to investigate the strength of the magnetic field in the Galactic Centre, close to Sagittarius A*. Here we report how persistent radio emission from this magnetar, for over four years since its discovery, has revealed large changes in the observed Faraday rotation measure, by up to 3500 rad m$^{-2}$ (a five per cent fractional change). From simultaneous analysis of the dispersion measure, we determine that these fluctuations are dominated by variations in the projected magnetic field, rather than the integrated free electron density, along the changing line of sight to the rapidly moving magnetar. From a structure function analysis of rotation measure variations, and a recent epoch of rapid change of rotation measure, we determine a minimum scale of magnetic fluctuations of size ~ 2 au at the Galactic Centre distance, inferring PSR J1745-2900 is just ~ 0.1 pc behind an additional scattering screen.

We study 379 central and 159 satellite early-type galaxies with two-dimensional kinematics from the integral-field survey Mapping Nearby Galaxies at APO (MaNGA) to determine how their angular momentum content depends on stellar and halo mass. Using the Yang et. al. (2007) group catalog, we identify central and satellite galaxies in groups with halo masses in the range 10^12.5 h^-1 M_sun < M_200b < 10^15 h^-1 M_sun. As in previous work, we see a sharp dependence on stellar mass, in the sense that ~ 70% of galaxies with stellar mass M_* > 10^11 h^-2 M_sun tend to have very little rotation, while nearly all galaxies at lower mass show some net rotation. The ~ 30% of high-mass galaxies that have significant rotation do not stand out in other galaxy properties except for a higher incidence of ionized gas emission. Our data are consistent with recent simulation results suggesting that major merging and gas accretion have more impact on the rotational support of lower-mass galaxies. When carefully matching the stellar mass distributions, we find no residual differences in angular momentum content between satellite and central galaxies at the 20\% level. Similarly, at fixed mass, galaxies have consistent rotation properties across a wide range of halo mass. However, we find that errors in classification of centrals and satellites with group finders systematically lowers differences between satellite and central galaxies at a level that is comparable to current measurement uncertainties. To improve constraints, the impact of group finding methods will have to be forward modeled via mock catalogs.

We report the discovery of OGLE-2016-BLG-1190Lb, which is likely to be the first Spitzer microlensing planet in the Galactic bulge/bar, an assignation that can be confirmed by two epochs of high-resolution imaging of the combined source-lens baseline object. The planet's mass M_p= 13.4+-0.9 M_J places it right at the deuterium burning limit, i.e., the conventional boundary between "planets" and "brown dwarfs". Its existence raises the question of whether such objects are really "planets" (formed within the disks of their hosts) or "failed stars" (low mass objects formed by gas fragmentation). This question may ultimately be addressed by comparing disk and bulge/bar planets, which is a goal of the Spitzer microlens program. The host is a G dwarf M_host = 0.89+-0.07 M_sun and the planet has a semi-major axis a~2.0 AU. We use Kepler K2 Campaign 9 microlensing data to break the lens-mass degeneracy that generically impacts parallax solutions from Earth-Spitzer observations alone, which is the first successful application of this approach. The microlensing data, derived primarily from near-continuous, ultra-dense survey observations from OGLE, MOA, and three KMTNet telescopes, contain more orbital information than for any previous microlensing planet, but not quite enough to accurately specify the full orbit. However, these data do permit the first rigorous test of microlensing orbital-motion measurements, which are typically derived from data taken over <1% of an orbital period.

As radio polarised emission from astrophysical objects traverse through foreground magnetised plasma, the physical conditions along the lines of sight are encrypted in the form of Rotation Measure (RM). We performed broadband spectro-polarimetric observations of high Rotation Measure (|RM| >~ 300 rad m-2) sources away from the Galactic plane (|b| > 10 deg) selected from the NVSS RM catalogue. The main goals are to verify the NVSS RM values, which could be susceptible to n{\pi}-ambiguity, as well as to identify the origin of the extreme RM values. We show that 40 % of our sample suffer from n{\pi}-ambiguity in the NVSS RM catalogue. There are also hints of RM variabilities over ~20 years epoch for most of our sources, as revealed by comparing the RM values of the two studies in the same frequency ranges after correcting for n{\pi}-ambiguity. At last, we demonstrate the possibility of applying QU-fitting to study the ambient media of AGNs.

We present the analysis of OGLE-2016-BLG-0613, for which the lensing light curve appears to be that of a typical binary-lens event with two caustic spikes but with a discontinuous feature on the trough between the spikes. We find that the discontinuous feature was produced by a planetary companion to the binary lens. We find 4 degenerate triple-lens solution classes, each composed of a pair of solutions according to the well-known wide/close planetary degeneracy. One of these solution classes is excluded due to its relatively poor fit. For the remaining three pairs of solutions, the most-likely primary mass is about $M_1\sim 0.7\,M_\odot$ while the planet is a super-Jupiter. In all cases the system lies in the Galactic disk, about half-way toward the Galactic bulge. However, in one of these three solution classes, the secondary of the binary system is a low-mass brown dwarf, with relative mass ratios (1 : 0.03 : 0.003), while in the two others the masses of the binary components are comparable. These two possibilities can be distinguished in about 2024 when the measured lens-source relative proper motion will permit separate resolution of the lens and source.

Context. The ability to automatically select scientifically-important transient events from an alert stream of many such events, and to conduct follow-up observations in response, will become increasingly important in astronomy. With wide-angle time domain surveys pushing to fainter limiting magnitudes, the capability to follow-up on transient alerts far exceeds our follow-up telescope resources, and effective target prioritization becomes essential. The RoboNet-II microlensing program is a pathfinder project which has developed an automated target selection process (RoboTAP) for gravitational microlensing events which are observed in real-time using the Las Cumbres Observatory telescope network. Aims. Follow-up telescopes typically have a much smaller field-of-view compared to surveys, therefore the most promising microlens- ing events must be automatically selected at any given time from an annual sample exceeding 2000 events. The main challenge is to select between events with a high planet detection sensitivity, aiming at the detection of many planets and characterizing planetary anomalies. Methods. Our target selection algorithm is a hybrid system based on estimates of the planet detection zones around a microlens. It follows automatic anomaly alerts and respects the expected survey coverage of specific events. Results. We introduce the RoboTAP algorithm, whose purpose is to select and prioritize microlensing events with high sensitivity to planetary companions. In this work, we determine the planet sensitivity of the RoboNet follow-up program and provide a working example of how a broker can be designed for a real-life transient science program conducting follow-up observations in response to alerts, exploring the issues that will confront similar programs being developed for the Large Synoptic Survey Telescope (LSST) and other time domain surveys.

A cold high-velocity (HV, $\sim$ 200 km/s) peak was first reported in several Galactic bulge fields based on the APOGEE commissioning observations. Both the existence and the nature of the high-velocity peak are still under debate. Here we revisit this feature with the latest APOGEE DR13 data. We find that most of the low latitude bulge fields display a skewed Gaussian distribution with a HV shoulder. However, only 3 out of 53 fields show distinct high-velocity peaks around 200 km/s. The velocity distribution can be well described by Gauss-Hermite polynomials, except the three fields showing clear HV peaks. We find that the correlation between the skewness parameter ($h_{3}$) and the mean velocity ($\bar{v}$), instead of a distinctive HV peak, is a strong indicator of the bar. It was recently suggested that the HV peak is composed of preferentially young stars. We choose three fields showing clear HV peaks to test this hypothesis using the metallicity, [$\alpha$/M] and [C/N] as age proxies. We find that both young and old stars show HV features. The similarity between the chemical abundances of stars in the HV peaks and the main component indicates that they are not systematically different in terms of chemical abundance or age. In contrast, there are clear differences in chemical space between stars in the Sagittarius dwarf and the bulge stars. The strong HV peaks off-plane are still to be explained properly, and could be different in nature.

Extended gaps in the debris disks of both Vega and Fomalhaut have been observed. These structures have been attributed to tidal perturbations by multiple super-Jupiter gas giant planets. Within the current observational limits, however, no such massive planets have been detected. Here we propose a less stringent `lone-planet' scenario to account for the observed structure with a single eccentric gas giant and suggest that clearing of these wide gaps is induced by its sweeping secular resonance. During the depletion of the disk gas, the planet's secular resonance propagates inward and clears a wide gap over an extended region of the disk. Although some residual intermediate-size planetesimals may remain in the gap, their surface density is too low to either produce super-Earths or lead to sufficiently frequent disruptive collisions to generate any observable dusty signatures. The main advantage of this lone-planet sweeping-secular-resonance model over the previous multiple gas giant tidal truncation scenario is the relaxed requirement on the number of gas giants. The observationally inferred upper mass limit can also be satisfied provided the hypothetical planet has a significant eccentricity. A significant fraction of solar or more massive stars bear gas giant planets with significant eccentricities. If these planets acquired their present-day kinematic properties prior to the depletion of their natal disks, their sweeping secular resonance would effectively impede the retention of neighboring planets and planetesimals over a wide range of orbital semi-major axes.

We study galaxy shapes in the Illustris cosmological hydrodynamic simulation. We find that massive galaxies have a higher probability of being prolate. For galaxies with stellar mass larger than $10^{11}\rm M_{\odot}$, 35 out of total 839 galaxies are prolate. For 21 galaxies with stellar mass larger than $10^{12}\rm M_{\odot}$, 9 are prolate, 4 are triaxial while the others are close to being oblate. There are almost no prolate galaxies with stellar mass smaller than $3\times10^{11}\rm M_{\odot}$. We check the merger history of the prolate galaxies, and find that they are formed by major dry mergers. All the prolate galaxies have at least one such merger, with most having mass ratios between $1:1$ and $1:3$. The gas fraction (gas mass to total baryon mass) of the progenitors is 0-3 percent for nearly all these mergers, except for one whose second progenitor contains $\sim 15\%$ gas mass, while its main progenitor still contains less than $5\%$. For the 35 massive prolate galaxies that we find, 18 of them have minor axis rotation, and their angular momenta mostly come from the spin angular momenta of the progenitors (usually that of the main progenitor). We analyse the merger orbits of these prolate galaxies and find that most of them experienced a nearly radial merger orbit. Oblate galaxies with major dry mergers can have either radial or circular merger orbits. We further discuss various properties of these prolate galaxies, such as spin parameter $\lambda_{\rm R}$, spherical anisotropy parameter $\beta$, dark matter fraction, as well as inner density slopes for the stellar, dark matter and total mass distributions.

We present a morphological study of the 17 lensed Lyman-$\alpha$ emitter (LAE) galaxies of the Baryon Oscillation Spectroscopic Survey Emission-Line Lens Survey (BELLS) for the GALaxy-Ly$\alpha$ EmitteR sYstems (BELLS GALLERY) sample. This analysis combines the magnification effect of strong galaxy-galaxy lensing with the high resolution of the Hubble Space Telescope ($HST$) to achieve a physical resolution of $\sim$80 pc for this $2<z<3$ LAE sample, allowing a detailed characterization of the LAE rest-frame ultraviolet continuum surface brightness profiles and substructure. We use lens-model reconstructions of the LAEs to identify and model individual clumps, which we subsequently use to constrain the parameters of a generative statistical model of the LAE population. Since the BELLS GALLERY sample is selected primarily on the basis of Lyman-$\alpha$ emission, the LAEs that we study here are likely to be directly comparable to those selected in wide-field narrow-band LAE surveys, in contrast with the lensed LAEs identified in cluster lensing fields. We find an LAE clumpiness fraction of approximately 88%, significantly higher than found in previous (non-lensing) studies. We find a well-resolved characteristic clump half-light radii of $\sim$350 pc, a scale comparable to the largest H II regions seen in the local universe. This statistical characterization of LAE surface-brightness profiles will be incorporated into future lensing analyses using the BELLS GALLERY sample to constrain the incidence of dark-matter substructure in the foreground lensing galaxies.

Magnetic fields play a pivotal role in the physics of interstellar medium in galaxies, but there are few observational constraints on how they evolve across cosmic time. Spatially resolved synchrotron polarization maps at radio wavelengths reveal well-ordered large-scale magnetic fields in nearby galaxies that are believed to grow from a seed field via a dynamo effect. To directly test and characterize this theory requires magnetic field strength and geometry measurements in cosmologically distant galaxies, which are challenging to obtain due to the limited sensitivity and angular resolution of current radio telescopes. Here, we report the cleanest measurements yet of magnetic fields in a galaxy beyond the local volume, free of the systematics traditional techniques would encounter. By exploiting the scenario where the polarized radio emission from a background source is gravitationally lensed by a foreground galaxy at z = 0.439 using broadband radio polarization data, we detected coherent $\mu$G magnetic fields in the lensing disk galaxy as seen 4.6 Gyrs ago, with similar strength and geometry to local volume galaxies. This is the highest redshift galaxy whose observed coherent magnetic field property is compatible with a mean-field dynamo origin.

We present the analysis of the binary gravitational microlensing event MOA-2015-BLG-020. The event has a fairly long timescale (about 63 days) and thus the light curve deviates significantly from the lensing model that is based on the rectilinear lens-source relative motion. This enables us to measure the microlensing parallax through the annual parallax effect. The microlensing parallax parameters constrained by the ground-based data are confirmed by the Spitzer observations through the satellite parallax method. By additionally measuring the angular Einstein radius from the analysis of the resolved caustic crossing, the physical parameters of the lens are determined. It is found that the binary lens is composed of two dwarf stars with masses $M_1 = 0.606 \pm 0.028M_\odot$ and $M_2 = 0.125 \pm 0.006M_\odot$ in the Galactic disk. Assuming the source star is at the same distance as the bulge red clump stars, we find the lens is at a distance $D_L = 2.44 \pm 0.10 kpc$. In the end, we provide a summary and short discussion of all published microlensing events in which the annual parallax effect is confirmed by other independent observations.

PHL 6625 is a luminous quasi-stellar object (QSO) at z = 0.3954 located behind the nearby galaxy NGC 247 (z = 0.0005). Hubble Space Telescope (HST) observations revealed an arc structure associated with it. We report on spectroscopic observations with the Very Large Telescope (VLT) and multiwavelength observations from the radio to the X-ray band for the system, suggesting that PHL 6625 and the arc are a close pair of merging galaxies, instead of a strong gravitational lens system. The QSO host galaxy is estimated to be (4-28) x 10^10 M_sun, and the mass of the companion galaxy of is estimated to be M_* = (6.8 +/- 2.4) x 10^9 M_sun, suggesting that this is a minor merger system. The QSO displays typical broad emission lines, from which a black hole mass of about (2-5) x 10^8 M_sun and an Eddington ratio of about 0.01-0.05 can be inferred. The system represents an interesting and rare case where a QSO is associated with an ongoing minor merger, analogous to Arp 142.

The rest-frame UV-optical (i.e., $NUV-B$) color is sensitive to both low-level recent star formation (specific star formation rate - sSFR) and dust. In this Letter, we extend our previous work on the origins of $NUV-B$ color gradients in star-forming galaxies (SFGs) at $z\sim1$ to those at $z\sim2$. We use a sample of 1335 large (semi-major axis radius $R_{\rm SMA}>0.''18$) SFGs with extended UV emission out to $2R_{\rm SMA}$ in the mass range $M_{\ast} = 10^{9}-10^{11}M_{\odot}$ at $1.5<z<2.8$ in the CANDELS/GOODS-S and UDS fields. We show that these SFGs generally have negative $NUV-B$ color gradients (redder centres), and their color gradients strongly increase with galaxy mass. We also show that the global rest-frame $FUV-NUV$ color is approximately linear with $A_{\rm V}$, which is derived by modeling the observed integrated FUV to NIR spectral energy distributions of the galaxies. Applying this integrated calibration to our spatially-resolved data, we find a negative dust gradient (more dust extinguished in the centers), which steadily becomes steeper with galaxy mass. We further find that the $NUV-B$ color gradients become nearly zero after correcting for dust gradients regardless of galaxy mass. This indicates that the sSFR gradients are negligible and dust reddening is likely the principal cause of negative UV-optical color gradients in these SFGs. Our findings support that the buildup of the stellar mass in SFGs at the Cosmic Noon is self-similar inside $2R_{\rm SMA}$.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We report on the first results from a large-scale observing campaign aiming to use astrometric microlensing to detect and place limits on the mass of single objects, including stellar remnants. We used the Hubble Space Telescope to monitor stars near the Galactic Center for 3 years, and we measured the brightness and positions of $\sim$2 million stars at each observing epoch. In addition to this, we monitored the same pointings using the VIMOS imager on the Very Large Telescope. The stars we monitored include several bright microlensing events observed from the ground by the OGLE collaboration. In this paper, we present the analysis of our photometric and astrometric measurements for 6 of these events, and derive mass constraints for the lens in each of these. Although these constraints are limited by the photometric precision of ground-based data, and our ability to determine the lens distance, we were able to constrain the size of the Einstein ring radius thanks to our precise astrometric measurements, the first routine measurements of this type from a large-scale observing program. This demonstrates the power of astrometric microlensing as a tool to constrain the masses of stars, stellar remnants, and, in the future, of extrasolar planets, using precise ground- and space-based observations.

With a sample of 48,161 K giant stars selected from the LAMOST DR 2 catalogue, we construct torus models in a large volume extending, for the first time, from the solar vicinity to a Galactocentric distance of $\sim 20$ kpc, reaching the outskirts of the Galactic disc. We show that the kinematics of the K giant stars match conventional models, e.g. as created by Binney in 2012, in the Solar vicinity. However such two-disc models fail if they are extended to the outer regions, even if an additional disc component is utilised. If we loosen constraints in the Sun's vicinity, we find that an effective thick disc model could explain the anti-centre of the MW. The LAMOST data imply that the sizes of the Galactic discs are much larger, and that the outer disc is much thicker, than previously thought, or alternatively that the outer structure is not a conventional disc at all. However, the velocity dispersion $\sigma_{0z}$ of the kinematically thick disc in the best-fitting model is about 80 km s$^{-1}$ and has a scale parameter $R_{\sigma}$ for an exponential distribution function of $\sim 19$ kpc. Such a height $\sigma_{0z}$ is strongly rejected by current measurements in the solar neighbourhood, and thus a model beyond quasi-thermal, two or three thin or thick discs is required.

This paper uses radial colour profiles to infer the distributions of dust, gas and star formation in z=0.4-1.4 star-forming main sequence galaxies. We start with the standard UVJ-based method to estimate dust extinction and specific star formation rate (sSFR). By replacing J with I band, a new calibration method suitable for use with ACS+WFC3 data is created (i.e. UVI diagram). Using a multi-wavelength multi-aperture photometry catalogue based on CANDELS, UVI colour profiles of 1328 galaxies are stacked in stellar mass and redshift bins. The resulting colour gradients, covering a radial range of 0.2--2.0 effective radii, increase strongly with galaxy mass and with global $A_V$. Colour gradient directions are nearly parallel to the Calzetti extinction vector, indicating that dust plays a more important role than stellar population variations. With our calibration, the resulting $A_V$ profiles fall much more slowly than stellar mass profiles over the measured radial range. sSFR gradients are nearly flat without central quenching signatures, except for $M_*>10^{10.5} M_{\odot}$, where central declines of 20--25 per cent are observed. Both sets of profiles agree well with previous radial sSFR and (continuum) $A_V$ measurements. They are also consistent with the sSFR profiles and, if assuming a radially constant gas-to-dust ratio, gas profiles in recent hydrodynamic models. We finally discuss the striking findings that SFR scales with stellar mass density in the inner parts of galaxies, and that dust content is high in the outer parts despite low stellar-mass surface densities there.

We report the discovery of a microlensing planet --- MOA-2016-BLG-227Lb --- with a large planet/host mass ratio of $q \simeq 9 \times 10^{-3}$. This event was located near the $K2$ Campaign 9 field that was observed by a large number of telescopes. As a result, the event was in the microlensing survey area of a number of these telescopes, and this enabled good coverage of the planetary light curve signal. High angular resolution adaptive optics images from the Keck telescope reveal excess flux at the position of the source above the flux of the source star, as indicated by the light curve model. This excess flux could be due to the lens star, but it could also be due to a companion to the source or lens star, or even an unrelated star. We consider all these possibilities in a Bayesian analysis in the context of a standard Galactic model. Our analysis indicates that it is unlikely that a large fraction of the excess flux comes from the lens, unless solar type stars are much more likely to host planets of this mass ratio than lower mass stars. We recommend that a method similar to the one developed in this paper be used for other events with high angular resolution follow-up observations when the follow-up observations are insufficient to measure the lens-source relative proper motion.

We perform Jeans anisotropic modeling (JAM) on elliptical and spiral galaxies from the MaNGA DR13 sample. By comparing the stellar mass-to-light ratios estimated from stellar population synthesis (SPS) and from JAM, we find a similar systematic variation of the initial mass function (IMF) as in the earlier $\rm ATLAS^{3D}$ results. Early type galaxies (elliptical and lenticular) with lower velocity dispersions within one effective radius are consistent with a Chabrier-like IMF while galaxies with higher velocity dispersions are consistent with a more bottom heavy IMF such as the Salpeter IMF. Spiral galaxies have similar systematic IMF variations, but with slightly different slopes and larger scatters, due to the uncertainties caused by higher gas fractions and extinctions for these galaxies. Furthermore, we examine the effects of stellar mass-to-light ratio gradients on our JAM modeling, and find that the trends from our results becomes stronger after considering the gradients.

The Beijing-Arizona Sky Survey (BASS) is a wide-field two-band photometric survey of the Northern Galactic Cap using the 90Prime imager on the 2.3 m Bok telescope at Kitt Peak. It is a four-year collaboration between the National Astronomical Observatory of China and Steward Observatory, the University of Arizona, serving as one of the three imaging surveys to provide photometric input catalogs for target selection of the Dark Energy Spectroscopic Instrument (DESI) project. BASS will take up to 240 dark/grey nights to cover an area of about 5400 deg$^2$ in the $g$ and $r$ bands. The 5$\sigma$ limiting AB magnitudes for point sources in the two bands, corrected for the Galactic extinction, are 24.0 and 23.4 mag, respectively. BASS, together with other DESI imaging surveys, will provide unique science opportunities that cover a wide range of topics in both Galactic and extragalactic astronomy.

The Fan Region is one of the dominant features in the polarized radio sky, long thought to be a local (distance < 500 pc) synchrotron feature. We present 1.3-1.8 GHz polarized radio continuum observations of the region from the Global Magneto-Ionic Medium Survey (GMIMS) and compare them to maps of Halpha and polarized radio continuum intensity from 0.408-353 GHz. The high-frequency (> 1 GHz) and low-frequency (< 600 MHz) emission have different morphologies, suggesting a different physical origin. Portions of the 1.5 GHz Fan Region emission are depolarized by about 30% by ionized gas structures in the Perseus Arm, indicating that this fraction of the emission originates >2 kpc away. We argue for the same conclusion based on the high polarization fraction at 1.5 GHz (about 40%). The Fan Region is offset with respect to the Galactic plane, covering -5{\deg} < b < +10{\deg}; we attribute this offset to the warp in the outer Galaxy. We discuss origins of the polarized emission, including the spiral Galactic magnetic field. This idea is a plausible contributing factor although no model to date readily reproduces all of the observations. We conclude that models of the Galactic magnetic field should account for the > 1 GHz emission from the Fan Region as a Galactic-scale, not purely local, feature.

We report the discovery of a very bright (r = 20.16), highly magnified, and yet intrinsically very luminous Ly{\alpha} emitter (LAE) at z = 2.82. This system comprises four images in the observer plane with a maximum separation of ~ 6" and it is lensed by a z = 0.55 massive early-type galaxy. It was initially identified in the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey for GALaxy-Ly{\alpha} EmitteR sYstems (BELLS GALLERY) survey, and follow-up imaging and spectroscopic observations using the Gran Telescopio Canarias (GTC) and William Herschel Telescope (WHT) confirmed the lensing nature of this system. A lens model using a singular isothermal ellipsoid in an external shear field reproduces quite well the main features of the system, yielding an Einstein radius of 2.95" +/- 0.10", and a total magnification factor for the LAE of 8.8 +/- 0.4. This LAE is one of the brightest and most luminous galaxy-galaxy strong lenses known. We present initial imaging and spectroscopy showing the basic physical and morphological properties of this lensed system.

We present \textsl{Hubble Space Telescope} (\textsl{HST}) F606W-band imaging observations of 21 galaxy-Ly$\alpha$ emitter lens candidates in the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) for GALaxy-Ly$\alpha$ EmitteR sYstems (BELLS GALLERY) survey. 17 systems are confirmed to be definite lenses with unambiguous evidence of multiple imaging. The lenses are primarily massive early-type galaxies (ETGs) at redshifts of approximately $0.55$, while the lensed sources are Ly$\alpha$ emitters (LAEs) at redshifts from 2 to 3. Although the \textsl{HST} imaging data are well fit by smooth lens models consisting of singular isothermal ellipsoids in an external shear field, a thorough exploration of dark substructures in the lens galaxies is required. The Einstein radii of the BELLS GALLERY lenses are on average $60\%$ larger than those of the BELLS lenses because of the much higher source redshifts which will allow a detailed investigation of the radius evolution of the mass profile in ETGs. With the aid of the average $\sim 13 \times$ lensing magnification, the LAEs are resolved to comprise individual star-forming knots of a wide range of properties with characteristic sizes from less than 100 pc to several kpc, rest-frame far UV apparent AB magnitudes from 29.6 to 24.2, and typical projected separations of 500 pc to 2 kpc.

We present a study on the stellar age and metallicity distributions for 1105 galaxies using the STARLIGHT software on MaNGA integral field spectra. We derive age and metallicity gradients by fitting straight lines to the radial profiles, and explore their correlations with total stellar mass M*, NUV-r colour and environments, as identified by both the large scale structure (LSS) type and the local density. We find that the mean age and metallicity gradients are close to zero but slightly negative, which is consistent with the inside-out formation scenario. Within our sample, we find that both the age and metallicity gradients show weak or no correlation with either the LSS type or local density environment. In addition, we also study the environmental dependence of age and metallicity values at the effective radii. The age and metallicity values are highly correlated with M* and NUV-r and are also dependent on LSS type as well as local density. Low-mass galaxies tend to be younger and have lower metallicity in low-density environments while high-mass galaxies are less affected by environment.

Decomposing an arbitrary electron energy distribution into sums of Maxwellian and power law components is an efficient method to calculate synchrotron emission and absorption. We use this method to study the effect of non-thermal electrons on submm images and spectra of the Galactic center black hole, Sgr A*. We assume a spatially uniform functional form for the electron distribution function and use a semi-analytic radiatively inefficient accretion flow and a 2D general relativistic MHD snapshot as example models of the underlying accretion flow structure. We develop simple analytic models which allow us to generalize from the numerical examples. A high energy electron component containing a small fraction (few per cent) of the total internal energy (e.g. a "power law tail") can produce a diffuse halo of emission, which modifies the observed image size and structure. A population of hot electrons with a larger energy fraction (e.g. resulting from a diffusion in electron energy space) can dominate the emission, so that the observed images and spectra are well approximated by considering only a single thermal component for a suitable choice of the electron temperature. We discuss the implications of these results for estimating accretion flow or black hole parameters from images and spectra, and for the identification of the black hole "shadow" in future mm-VLBI data. In particular, the location of the first minimum in visibility profiles does not necessarily correspond to the shadow size as sometimes assumed.

We present a study of the magnetic field properties of NGC 4038/9 (the `Antennae' galaxies), the closest example of a late stage merger of two spiral galaxies. Wideband polarimetric observations were performed using the Karl G. Jansky Very Large Array between 2 and 4 GHz. Rotation measure synthesis and Faraday depolarization analysis was performed to probe the magnetic field strength and structure at spatial resolution of $\sim1$ kpc. Highly polarized emission from the southern tidal tail is detected with intrinsic fractional polarization close to the theoretical maximum ($0.62\pm0.18$), estimated by fitting the Faraday depolarization with a volume that is both synchrotron emitting and Faraday rotating containing random magnetic fields. Magnetic fields are well aligned along the tidal tail and the Faraday depths shows large-scale smooth variations preserving its sign. This suggests the field in the plane of the sky to be regular up to $\sim20$ kpc, which is the largest detected regular field structure on galactic scales. The equipartition field strength of $\sim8.5~\mu$G of the regular field in the tidal tail is reached within a few 100 Myr, likely generated by stretching of the galactic disc field by a factor of 4--9 during the tidal interaction. The regular field strength is greater than the turbulent fields in the tidal tail. Our study comprehensively demonstrates, although the magnetic fields within the merging bodies are dominated by strong turbulent magnetic fields of $\sim20~\mu$G in strength, tidal interactions can produce large-scale regular field structure in the outskirts.

The dynamics and evolution of any galactic structure are strongly influenced by the properties of the orbits that constitute it. In this paper, we compare two orbit classification schemes, one by Laskar (NAFF) , and the other by Carpintero and Aguilar (CA), by applying both of them to orbits obtained by following individual particles in a numerical simulation of a barred galaxy. We find that, at least for our case and some provisos, the main frequencies calculated by the two methods are in good agreement: for $80\%$ of the orbits the difference between the results of the two methods is less than $5\%$ for all three main frequencies. However, it is difficult to evaluate the amount of regular or chaotic bar orbits in a given system. The fraction of regular orbits obtained by the NAFF method strongly depends on the critical frequency drift parameter, while in the CA method the number of fundamental frequencies strongly depends on the frequency difference parameter $L_{\rm r}$ and the maximum integer used for searching the linear independence of the fundamental frequencies. We also find that, for a given particle, in general the projection of its motion along the bar minor axis is more regular than the other two projections, while the projection along the intermediate axis is the least regular.

The Thirty Meter Telescope (TMT) first light instrument IRIS (Infrared Imaging Spectrograph) will complete its preliminary design phase in 2016. The IRIS instrument design includes a near-infrared (0.85 - 2.4 micron) integral field spectrograph (IFS) and imager that are able to conduct simultaneous diffraction-limited observations behind the advanced adaptive optics system NFIRAOS. The IRIS science cases have continued to be developed and new science studies have been investigated to aid in technical performance and design requirements. In this development phase, the IRIS science team has paid particular attention to the selection of filters, gratings, sensitivities of the entire system, and science cases that will benefit from the parallel mode of the IFS and imaging camera. We present new science cases for IRIS using the latest end-to-end data simulator on the following topics: Solar System bodies, the Galactic center, active galactic nuclei (AGN), and distant gravitationally-lensed galaxies. We then briefly discuss the necessity of an advanced data management system and data reduction pipeline.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

We introduce the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) for GALaxy-Ly$\alpha$ EmitteR sYstems (BELLS GALLERY) Survey, which is a Hubble Space Telescope program to image a sample of galaxy-scale strong gravitational lens candidate systems with high-redshift Ly$\alpha$ emitters (LAEs) as the background sources. The goal of the BELLS GALLERY Survey is to illuminate dark substructures in galaxy-scale halos by exploiting the small-scale clumpiness of rest-frame far-UV emission in lensed LAEs, and to thereby constrain the slope and normalization of the substructure-mass function. In this paper, we describe in detail the spectroscopic strong-lens selection technique, which is based on methods adopted in the previous Sloan Lens ACS (SLACS) Survey, BELLS, and SLACS for the Masses Survey. We present the BELLS GALLERY sample of the 21 highest-quality galaxy--LAE candidates selected from $\approx 1.4 \times 10^6$ galaxy spectra in the BOSS of the Sloan Digital Sky Survey III. These systems consist of massive galaxies at redshifts of approximately 0.5 strongly lensing LAEs at redshifts from 2--3. The compact nature of LAEs makes them an ideal probe of dark substructures, with a substructure-mass sensitivity that is unprecedented in other optical strong-lens samples. The magnification effect from lensing will also reveal the structure of LAEs below 100 pc scales, providing a detailed look at the sites of the most concentrated unobscured star formation in the universe. The source code used for candidate selection is available for download as a part of this release.

We present a combined analysis of the observations of the gravitational microlensing event OGLE-2015-BLG-0479 taken both from the ground and by the {\it Spitzer Space Telescope}. The light curves seen from the ground and from space exhibit a time offset of $\sim 13$ days between the caustic spikes, indicating that the relative lens-source positions seen from the two places are displaced by parallax effects. From modeling the light curves, we measure the space-based microlens parallax. Combined with the angular Einstein radius measured by analyzing the caustic crossings, we determine the mass and distance of the lens. We find that the lens is a binary composed of two G-type stars with masses $\sim 1.0\ M_\odot$ and $\sim 0.9\ M_\odot$ located at a distance $\sim 3$ kpc. In addition, we are able to constrain the complete orbital parameters of the lens thanks to the precise measurement of the microlens parallax derived from the joint analysis. In contrast to the binary event OGLE-2014-BLG-1050, which was also observed by {\it Spitzer}, we find that the interpretation of OGLE-2015-BLG-0479 does not suffer from the degeneracy between $(\pm,\pm)$ and $(\pm,\mp)$ solutions, confirming that the four-fold parallax degeneracy in single-lens events collapses into the two-fold degeneracy for the general case of binary-lens events. The location of the blend in the color-magnitude diagram is consistent with the lens properties, suggesting that the blend is the lens itself. The blend is bright enough for spectroscopy and thus this possibility can be checked from future follow-up observations.

Simultaneous observations of microlensing events from multiple locations allow for the breaking of degeneracies between the physical properties of the lensing system, specifically by exploring different regions of the lens plane and by directly measuring the "microlens parallax". We report the discovery of a 30-55$M_J$ brown dwarf orbiting a K dwarf in microlensing event OGLE-2015-BLG-1319. The system is located at a distance of $\sim$5 kpc toward the Galactic bulge. The event was observed by several ground-based groups as well as by $Spitzer$ and $Swift$, allowing the measurement of the physical properties. However, the event is still subject to an 8-fold degeneracy, in particular the well-known close-wide degeneracy, and thus the projected separation between the two lens components is either $\sim$0.25 AU or $\sim$45 AU. This is the first microlensing event observed by $Swift$, with the UVOT camera. We study the region of microlensing parameter space to which $Swift$ is sensitive, finding that while for this event $Swift$ could not measure the microlens parallax with respect to ground-based observations, it can be important for other events. Specifically, for detecting nearby brown dwarfs and free-floating planets in high magnification events.

We calculate the microlensing event rate and typical time-scales for the free-floating planet (FFP) population that is predicted by the core accretion theory of planet formation. The event rate is found to be ~$1.8\times 10^{-3}$ of that for the stellar population. While the stellar microlensing event time-scale peaks at around 20 days, the median time-scale for FFP events (~0.1 day) is much shorter. Our values for the event rate and the median time-scale are significantly smaller than those required to explain the \cite{Sum+11} result, by factors of ~13 and ~16, respectively. The inclusion of planets at wide separations does not change the results significantly. This discrepancy may be too significant for standard versions of both the core accretion theory and the gravitational instability model to explain satisfactorily. Therefore, either a modification to the planet formation theory is required, or other explanations to the excess of short-time-scale microlensing events are needed. Our predictions can be tested by ongoing microlensing experiment such as KMTNet, and by future satellite missions such as WFIRST and Euclid.

The rest-frame UV-optical (i.e., NUV-B) color index is sensitive to the low-level recent star formation and dust extinction, but it is insensitive to the metallicity. In this Letter, we have measured the rest-frame NUV-B color gradients in ~1400 large ($\rm r_e>0.18^{\prime\prime}$), nearly face-on (b/a>0.5) main-sequence star-forming galaxies (SFGs) between redshift 0.5 and 1.5 in the CANDELS/GOODS-S and UDS fields. With this sample, we study the origin of UV-optical color gradients in the SFGs at z~1 and discuss their link with the buildup of stellar mass. We find that the more massive, centrally compact, and more dust extinguished SFGs tend to have statistically more negative raw color gradients (redder centers) than the less massive, centrally diffuse, and less dusty SFGs. After correcting for dust reddening based on optical-SED fitting, the color gradients in the low-mass ($M_{\ast} <10^{10}M_{\odot}$) SFGs generally become quite flat, while most of the high-mass ($M_{\ast} > 10^{10.5}M_{\odot}$) SFGs still retain shallow negative color gradients. These findings imply that dust reddening is likely the principal cause of negative color gradients in the low-mass SFGs, while both increased central dust reddening and buildup of compact old bulges are likely the origins of negative color gradients in the high-mass SFGs. These findings also imply that at these redshifts the low-mass SFGs buildup their stellar masses in a self-similar way, while the high-mass SFGs grow inside out.

With the aim of exploring the fast evolutionary path from the blue cloud of star-forming galaxies to the red sequence of quiescent galaxies in the local universe, we select a local advanced merging infrared luminous and ultraluminous galaxy (adv-merger (U)LIRGs) sample and perform careful dust extinction corrections to investigate their positions in the SFR-$M_{\ast}$, u-r and NUV-r color-mass diagrams. The sample consists of 89 (U)LIRGs at the late merger stage, obtained from cross-correlating the IRAS Point Source Catalog Redshift Survey and 1 Jy ULIRGs samples with the Sloan Digital Sky Survey DR7 database. Our results show that $74\%\pm 5\%$ of adv-merger (U)LIRGs are localized above the $1\, \sigma$ line of the local star-forming galaxy main sequence. We also find that all adv-merger (U)LIRGs are more massive than and as blue as the blue cloud galaxies after corrections of Galactic and internal dust extinctions, with $95\%\pm 2\%$ and $81\%\pm 4\%$ of them outside the blue cloud on the u-r and NUV-r color-mass diagrams, respectively. These results, combined with the short timescale for exhausting the molecular gas reservoir in adv-merger (U)LIRGs ($3\times 10^{7}$ to $3\times 10^{8}$ years), imply that the adv-merger (U)LIRGs are likely at the starting point of the fast evolutionary track previously proposed by several groups. While the number density of adv-merger (U)LIRGs is only $\sim 0.1\%$ of the blue cloud star-forming galaxies in the local universe, this evolutionary track may play a more important role at high redshift.

We report on the mass and distance measurements of two single-lens events from the 2015 \emph{Spitzer} microlensing campaign. With both finite-source effect and microlens parallax measurements, we find that the lens of OGLE-2015-BLG-1268 is very likely a brown dwarf. Assuming that the source star lies behind the same amount of dust as the Bulge red clump, we find the lens is a $45\pm7$ $M_{\rm J}$ brown dwarf at $5.9\pm1.0$ kpc. The lens of of the second event, OGLE-2015-BLG-0763, is a $0.50\pm0.04$ $M_\odot$ star at $6.9\pm1.0$ kpc. We show that the probability to definitively measure the mass of isolated microlenses is dramatically increased once simultaneous ground- and space-based observations are conducted.

$K2$'s Campaign 9 ($K2$C9) will conduct a $\sim$3.7 deg$^{2}$ survey toward the Galactic bulge from 7/April through 1/July of 2016 that will leverage the spatial separation between $K2$ and the Earth to facilitate measurement of the microlens parallax $\pi_{\rm E}$ for $\gtrsim$127 microlensing events. These will include several that are planetary in nature as well as many short-timescale microlensing events, which are potentially indicative of free-floating planets (FFPs). These satellite parallax measurements will in turn allow for the direct measurement of the masses of and distances to the lensing systems. In this white paper we provide an overview of the $K2$C9 space- and ground-based microlensing survey. Specifically, we detail the demographic questions that can be addressed by this program, including the frequency of FFPs and the Galactic distribution of exoplanets, the observational parameters of $K2$C9, and the array of resources dedicated to concurrent observations. Finally, we outline the avenues through which the larger community can become involved, and generally encourage participation in $K2$C9, which constitutes an important pathfinding mission and community exercise in anticipation of $WFIRST$.

We report Faraday rotation measurements of 11 extragalactic radio sources with lines of sight through the Rosette Nebula, a prominent HII region associated with the star cluster NGC 2244. It is also a prototypical example of a "stellar bubble" produced by the winds of the stars in NGC 2244. The goal of these measurements is to better determine the strength and structure of the magnetic field in the nebula. We calculate the rotation measure (RM) through two methods, a least-squares fit to $\chi$( $\lambda^2$) and Rotation Measure Synthesis. In conjunction with our results from Savage et al. (2013), we find an excess RM due to the shell of the nebula of +40 to +1200 rad m$^{-2}$ above a background RM of +147 rad m$^{-2}$. We discuss two forms of a simple shell model intended to reproduce the magnitude of the observed RM as a function of distance from the center of the Rosette Nebula. The models represent different physical situations for the magnetic field within the shell of the nebula. The first assumes that there is an increase in the magnetic field strength and plasma density at the outer radius of the HII region, such as would be produced by a strong magnetohydrodynamic shock wave. The second model assumes that any increase in the RM is due solely to an increase in the density, and the Galactic magnetic field is unaffected in the shell. We employ a Bayesian analysis to distinguish between the two forms of the model.

We report the detection and mass measurement of a binary lens OGLE-2015-BLG-1285La,b, with the more massive component having $M_1>1.35\,M_\odot$ (80% probability). A main-sequence star in this mass range is ruled out by limits on blue light, meaning that a primary in this mass range must be a neutron star or black hole. The system has a projected separation $r_\perp= 6.1\pm 0.4\,{\rm AU}$ and lies in the Galactic bulge. These measurements are based on the "microlens parallax" effect, i.e., comparing the microlensing light curve as seen from $Spitzer$, which lay at $1.25\,{\rm AU}$ projected from Earth, to the light curves from four ground-based surveys, three in the optical and one in the near infrared. Future adaptive optics imaging of the companion by 30m class telescopes will yield a much more accurate measurement of the primary mass. This discovery both opens the path and defines the challenges to detecting and characterizing black holes and neutron stars in wide binaries, with either dark or luminous companions. In particular, we discuss lessons that can be applied to future $Spitzer$ and $Kepler$ K2 microlensing parallax observations.

We assess the effectiveness of the Jeans-Anisotropic-MGE (JAM) technique with a state-of-the-art cosmological hydrodynamic simulation, the Illustris project. We perform JAM modelling on 1413 simulated galaxies with stellar mass M^* > 10^{10}M_{sun}, and construct an axisymmetric dynamical model for each galaxy. Combined with a Markov Chain Monte Carlo (MCMC) simulation, we recover the projected root-mean-square velocity (V_rms) field of the stellar component, and investigate constraints on the stellar mass-to-light ratio, M^*/L, and the fraction of dark matter f_{DM} within 2.5 effective radii (R_e). We find that the enclosed total mass within 2.5 R_e is well constrained to within 10%. However, there is a degeneracy between the dark matter and stellar components with correspondingly larger individual errors. The 1 sigma scatter in the recovered M^*/L is 30-40% of the true value. The accuracy of the recovery of M^*/L depends on the triaxial shape of a galaxy. There is no significant bias for oblate galaxies, while for prolate galaxies the JAM-recovered stellar mass is on average 18% higher than the input values. We also find that higher image resolutions alleviate the dark matter and stellar mass degeneracy and yield systematically better parameter recovery.

Measurement of the local dark matter density plays an important role in both Galactic dynamics and dark matter direct detection experiments. However, the estimated values from previous works are far from agreeing with each other. In this work, we provide a well-defined observed sample with 1427 G \& K type main-sequence stars from the LAMOST spectroscopic survey, taking into account selection effects, volume completeness, and the stellar populations. We apply a vertical Jeans equation method containing a single exponential stellar disk, a razor thin gas disk, and a constant dark matter density distribution to the sample, and obtain a total surface mass density of $\rm {78.7 ^{+3.9}_{-4.7}\ M_{\odot}\ pc^{-2}}$ up to 1 kpc and a local dark matter density of $0.0159^{+0.0047}_{-0.0057}\,\rm M_{\odot}\,\rm pc^{-3}$. We find that the sampling density (i.e. number of stars per unit volume) of the spectroscopic data contributes to about two-thirds of the uncertainty in the estimated values. We discuss the effect of the tilt term in the Jeans equation and find it has little impact on our measurement. Other issues, such as a non-equilibrium component due to perturbations and contamination by the thick disk population, are also discussed.

We report the detection of a Cold Neptune m_planet=21+/-2MEarth orbiting a 0.38MSol M dwarf lying 2.5-3.3 kpc toward the Galactic center as part of a campaign combining ground-based and Spitzer observations to measure the Galactic distribution of planets. This is the first time that the complex real-time protocols described by Yee et al. (2015), which aim to maximize planet sensitivity while maintaining sample integrity, have been carried out in practice. Multiple survey and follow-up teams successfully combined their efforts within the framework of these protocols to detect this planet. This is the second planet in the Spitzer Galactic distribution sample. Both are in the near-to-mid disk and clearly not in the Galactic bulge.

Time delays in gravitational lenses can be used to determine the Hubble constant and the lens potential. In future surveys, many gravitational lenses can be discovered, and their time delays and image positions can in principle be measured. Using an elliptical power-law potential, we show that combinations of image positions and time delays for quadruple lenses yield simple analytical expressions that are connected with observable quantities. These relations can be used to obtain the approximate axis ratio q, the Einstein radius and the slope. We apply this method to RX J1131-1231, and show that our analytical results match the full numerical determinations approximately. Our approach can quickly determine rough values of lens parameters, which can then be used as initial guesses for further refinement through numerical modelling and may be useful for automated lens search in large surveys.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

We know that magnetic fields are pervasive across all scales in the Universe and over all of cosmic time and yet our understanding of many of the properties of magnetic fields is still limited. We do not yet know when, where or how the first magnetic fields in the Universe were formed, nor do we fully understand their role in fundamental processes such as galaxy formation or cosmic ray acceleration or how they influence the evolution of astrophysical objects. The greatest challenge to addressing these issues has been a lack of deep, broad bandwidth polarimetric data over large areas of the sky. The Square Kilometre Array will radically improve this situation via an all-sky polarisation survey that delivers both high quality polarisation imaging in combination with observations of 7-14 million extragalactic rotation measures. Here we summarise how this survey will improve our understanding of a range of astrophysical phenomena on scales from individual Galactic objects to the cosmic web.

We present an analysis of the diffuse soft X-ray emission from the nuclear region of M51 combining both XMM-Newton RGS and Chandra data. Most of the RGS spectrum of M51 can be fitted with a thermal model with a temperature of $\sim0.5$ keV except for the OVII triplet, which is forbidden-line dominated. The Fe L-shell lines peak around the southern cloud, where the OVIII and NVII Lya lines also peak. In contrast, the peak of the OVII forbidden line is about 10$"$ offset from that of the other lines, indicating that it is from a spatially distinct component. The spatial distribution of the OVII triplet mapped by the Chandra data shows that most of the OVII triplet flux is located at faint regions near edges, instead of the southern cloud where other lines peak. This distribution of the OVII triplet is inconsistent with the photoionization model. Other mechanisms that could produce the anomalous OVII triplet, including a recombining plasma and charge exchange X-ray emission, are discussed.

We explore the kinematics (both the radial velocity and the proper motion) of the vertical X-shaped feature in the Milky Way with an N-body bar/bulge model. From the solar perspective, the distance distribution of particles is double-peaked in fields passing through the X-shape. The separation and amplitude ratio between the two peaks qualitatively match the observed trends towards the Galactic bulge. We confirm clear signatures of cylindrical rotation in the pattern of mean radial velocity across the bar/bulge region. We also find possible imprints of coherent orbital motion inside the bar structure in the radial velocity distribution along l=0 degree, where the near and far sides of the bar/bulge show excesses of approaching and receding particles. The coherent orbital motion is also reflected in the slight displacement of the zero-velocity-line in the mean radial velocity, and the displacement of the maximum/minimum in the mean longitudinal proper motion across the bulge region. We find some degree of anisotropy in the stellar velocity within the X-shape, but the underlying orbital family of the X-shape cannot be clearly distinguished. Two potential applications of the X-shape in previous literature are tested, i.e., bulge rotation and Galactic center measurements. We find that the proper motion difference between the two sides of the X-shape can be used to estimate the mean azimuthal streaming motion of the bulge, but not the pattern speed of the bar. We also demonstrate that the Galactic center can be located with the X-shape, but the accuracy depends on the fitting scheme, the number of fields, and their latitudinal coverage.

Discrepancies between the observed and model-predicted radio flux ratios are seen in a number of quadruply-lensed quasars. The most favored interpretation of these anomalies is that CDM substructures present in lensing galaxies perturb the lens potentials and alter image magnifications and thus flux ratios. So far no consensus has emerged regarding whether or not the predicted CDM substructure abundance fully accounts for the lensing flux anomaly observations. Accurate modeling relies on a realistic lens sample in terms of both the lens environment and internal structures and substructures. In this paper we construct samples of generalised and specific lens potentials, to which we add (rescaled) subhalo populations from the galaxy-scale Aquarius and the cluster-scale Phoenix simulation suites. We further investigate the lensing effects from subhalos of masses several orders of magnitude below the simulation resolution limit. The resulting flux ratio distributions are compared to the currently best available sample of radio lenses. The observed anomalies in B0128+437, B0712+472 and B1555+375 are more likely to be caused by propagation effects or oversimplified lens modeling, signs of which are already seen in the data. Among the quadruple systems that have closely located image triplets/pairs, the anomalous flux ratios of MG0414+0534 can be reproduced by adding CDM subhalos to its macroscopic lens potential, with a probability of 5%-20%; for B0712+472, B1422+231, B1555+375 and B2045+265, these probabilities are only of a few percent. We hence find that CDM substructures are unlikely to be the whole reason for radio flux anomalies. We discuss other possible effects that might also be at work.

Faraday rotation of polarised background sources is a unique probe of astrophysical magnetic fields in a diverse range of foreground objects. However, to understand the properties of the polarised sources themselves and of depolarising phenomena along the line of sight, we need to complement Faraday rotation data with polarisation observations over very broad bandwidths. Just as it is impossible to properly image a complex source with limited u-v coverage, we can only meaningfully understand the magneto-ionic properties of polarised sources if we have excellent coverage in $\lambda^2$-space. We here propose a set of broadband polarisation surveys with the Square Kilometre Array, which will provide a singular set of scientific insights on the ways in which galaxies and their environments have evolved over cosmic time.

Magnetic fields are an important ingredient of the interstellar medium (ISM). Besides their importance for star formation, they govern the transport of cosmic rays, relevant to the launch and regulation of galactic outflows and winds, which in turn are pivotal in shaping the structure of halo magnetic fields. Mapping the small-scale structure of interstellar magnetic fields in many nearby galaxies is crucial to understand the interaction between gas and magnetic fields, in particular how gas flows are affected. Elucidation of the magnetic role in, e.g., triggering star formation, forming and stabilising spiral arms, driving outflows, gas heating by reconnection and magnetising the intergalactic medium has the potential to revolutionise our physical picture of the ISM and galaxy evolution in general. Radio polarisation observations in the very nearest galaxies at high frequencies (>= 3 GHz) and with high spatial resolution (<= 5") hold the key here. The galaxy survey with SKA1 that we propose will also be a major step to understand the galactic dynamo, which is important for models of galaxy evolution and for astrophysical magnetohydrodynamics in general. Field amplification by turbulent gas motions, which is crucial for efficient dynamo action, has been investigated so far only in simulations, while compelling evidence of turbulent fields from observations is still lacking.

Magnetic fields play an important role in shaping the structure and evolution of the interstellar medium (ISM) of galaxies, but the details of this relationship remain unclear. With SKA1, the 3D structure of galactic magnetic fields and its connection to star formation will be revealed. A highly sensitive probe of the internal structure of the magnetoionized ISM is the partial depolarization of synchrotron radiation from inside the volume. Different configurations of magnetic field and ionized gas within the resolution element of the telescope lead to frequency-dependent changes in the observed degree of polarization. The results of spectro-polarimetric observations are tied to physical structure in the ISM through comparison with detailed modeling, supplemented with the use of new analysis techniques that are being actively developed and studied within the community such as Rotation Measure Synthesis. The SKA will enable this field to come into its own and begin the study of the detailed structure of the magnetized ISM in a sample of nearby galaxies, thanks to its extraordinary wideband capabilities coupled with the combination of excellent surface brightness sensitivity and angular resolution.

Magnetic fields in the Milky Way are present on a wide variety of sizes and strengths, influencing many processes in the Galactic ecosystem such as star formation, gas dynamics, jets, and evolution of supernova remnants or pulsar wind nebulae. Observation methods are complex and indirect; the most used of these are a grid of rotation measures of unresolved polarized extragalactic sources, and broadband polarimetry of diffuse emission. Current studies of magnetic fields in the Milky Way reveal a global spiral magnetic field with a significant turbulent component; the limited sample of magnetic field measurements in discrete objects such as supernova remnants and HII regions shows a wide variety in field configurations; a few detections of magnetic fields in Young Stellar Object jets have been published; and the magnetic field structure in the Galactic Center is still under debate. The SKA will unravel the 3D structure and configurations of magnetic fields in the Milky Way on sub-parsec to galaxy scales, including field structure in the Galactic Center. The global configuration of the Milky Way disk magnetic field, probed through pulsar RMs, will resolve controversy about reversals in the Galactic plane. Characteristics of interstellar turbulence can be determined from the grid of background RMs. We expect to learn to understand magnetic field structures in protostellar jets, supernova remnants, and other discrete sources, due to the vast increase in sample sizes possible with the SKA. This knowledge of magnetic fields in the Milky Way will not only be crucial in understanding of the evolution and interaction of Galactic structures, but will also help to define and remove Galactic foregrounds for a multitude of extragalactic and cosmological studies.

We present a study of the magneto-ionic medium in the Whirlpool galaxy (M51) using new wide-band multi-configuration polarization data at L band (1-2 GHz) obtained at the Karl G. Jansky Very Large Array. By fitting the observed diffuse complex polarization $Q$+$iU$ as a function of wavelength directly to various depolarization models, we find that polarized emission from M51 at 1-2 GHz originates from the top of the synchrotron disk and then experiences Faraday rotation in the near-side thermal halo of the galaxy. Thus, the scale height of the thermal gas must exceed that of the synchrotron emitting gas at L band. The observed Faraday depth distribution at L band is consistent with a halo field that comprises of a plane-parallel bisymmetric component and a vertical component which produces a Faraday rotation of $\sim$ $-$9 rad m$^{-2}$. The derived rotation measure structure functions indicate a characteristic scale of rotation measure fluctuations of less than 560 pc in the disk and approximately 1 kpc in the halo. The outer scale of turbulence of 1 kpc found in the halo of M51 is consistent with superbubbles and the Parker instability being the main energy injection mechanisms in galactic halos.

We use about 15,000 F/G nearby dwarf stars selected from the LAMOST pilot survey to map the U-V velocity distribution in the solar neighbourhood. An extreme deconvolution algorithm is applied to reconstruct an empirical multi-Gaussian model. In addition to the well known substructures, e.g., Sirius, Coma Berenices, Hyades-Pleiades over-densities, several new substructures are unveiled. A ripple-like structure from (U, V) = (-120, -5) to (103, -32)km/s is clearly seen in the U-V distribution. This structure seems associated with resonance induced by the Galactic bar, since it is extended in U while having a small dispersion in V at the same time. A ridge structure between (U, V) = (-60, 40) and (-15, 15) km/s is also found. Although similar substructures have been seen in the Hipparcos data, their origin is still unclear. Another compact over-density is seen at (U, V) = (-102, -24). With this large data sample, we find that the substructure located at V~70 km/s and the Arcturus group are essentially parallel in V, which may indicate that they originate from an unrelaxed disk component perturbed by the rotating bar.

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

Lensing flux-ratio anomalies are most likely caused by gravitational lensing by small-scale dark matter structures. These anomalies offer the prospect of testing a fundamental prediction of the cold dark matter (CDM) cosmological model: the existence of numerous substructures that are too small to host visible galaxies. In two previous studies we found that the number of subhalos in the six high-resolution simulations of CDM galactic halos of the Aquarius project is not sufficient to account for the observed frequency of flux ratio anomalies seen in selected quasars from the CLASS survey. These studies were limited by the small number of halos used, their narrow range of masses (1-2E12 solar masses) and the small range of lens ellipticities considered. We address these shortcomings by investigating the lensing properties of a large sample of halos with a wide range of masses in two sets of high resolution simulations of cosmological volumes and comparing them to a currently best available sample of radio quasars. We find that, as expected, substructures do not change the flux-ratio probability distribution of image pairs and triples with large separations, but they have a significant effect on the distribution at small separations. For such systems, CDM substructures can account for a substantial fraction of the observed flux-ratio anomalies. For large close-pair separation systems, the discrepancies existing between the observed flux ratios and predictions from smooth halo models are attributed to simplifications inherent in these models which do not take account of fine details in the lens mass distributions.

The Exoplanet Euclid Legacy Survey (ExELS) proposes to determine the frequency of cold exoplanets down to Earth mass from host separations of ~1 AU out to the free-floating regime by detecting microlensing events in Galactic Bulge. We show that ExELS can also detect large numbers of hot, transiting exoplanets in the same population. The combined microlensing+transit survey would allow the first self-consistent estimate of the relative frequencies of hot and cold sub-stellar companions, reducing biases in comparing "near-field" radial velocity and transiting exoplanets with "far-field" microlensing exoplanets. The age of the Bulge and its spread in metallicity further allows ExELS to better constrain both the variation of companion frequency with metallicity and statistically explore the strength of star-planet tides. We conservatively estimate that ExELS will detect ~4100 sub-stellar objects, with sensitivity typically reaching down to Neptune-mass planets. Of these, ~600 will be detectable in both Euclid's VIS (optical) channel and NISP H-band imager, with ~90% of detections being hot Jupiters. Likely scenarios predict a range of 2900-7000 for VIS and 400-1600 for H-band. Twice as many can be expected in VIS if the cadence can be increased to match the 20-minute H-band cadence. The separation of planets from brown dwarfs via Doppler boosting or ellipsoidal variability will be possible in a handful of cases. Radial velocity confirmation should be possible in some cases, using 30-metre-class telescopes. We expect secondary eclipses, and reflection and emission from planets to be detectable in up to ~100 systems in both VIS and NISP-H. Transits of ~500 planetary-radius companions will be characterised with two-colour photometry and ~40 with four-colour photometry (VIS,YJH), and the albedo of (and emission from) a large sample of hot Jupiters in the H-band can be explored statistically.

The China Dark Matter Experiment reports results on light WIMP dark matter searches at the China Jinping Underground Laboratory with a germanium detector array with a total mass of 20 g. The physics threshold achieved is 177 eVee ("ee" represents electron equivalent energy) at 50% signal efficiency. With 0.784 kg-days of data, exclusion region on spin-independent coupling with the nucleon is derived, improving over our earlier bounds at WIMP mass less than 4.6 GeV.

We study the dynamics of the giant elliptical galaxy M87 from the central to the outermost regions with the made-to-measure (M2M) method. We use a new catalogue of 922 globular cluster line-of- sight velocities extending to a projected radius of 180 kpc (equivalent to 25 M87 effective radii), and SAURON integral field unit data within the central 2.4 kpc. 263 globular clusters, mainly located beyond 40 kpc, are newly observed by the Next Generation Virgo Survey (NGVS). For the M2M modelling, the gravitational potential is taken as a combination of a luminous matter potential with a constant stellar mass-to-light ratio and a dark matter potential modelled as a logarithmic potential. Our best dynamical model returns a stellar mass-to-light ratio in the I band of M/LI = 6.0(+ -0.3) M_sun/L_sun with a dark matter potential scale velocity of 591(+ -50) km/s and scale radius of 42(+ -10) kpc. We determine the total mass of M87 within 180 kpc to be (1.5 + - 0.2) 10^13 M_sun. The mass within 40 kpc is smaller than previous estimates determined using globular cluster kinematics that did not extend beyond 45 kpc. With our new globular cluster velocities at much larger radii, we see that globular clusters around 40 kpc show an anomalously large velocity dispersion which affected previous results. The mass we derived is in good agreement with that inferred from ROSAT X-ray observation out to 180 kpc. Within 30 kpc our mass is also consistent with that inferred from Chandra and XMM-Newton X-ray observations, while within 120 kpc it is about 20% smaller. The model velocity dispersion anisotropy beta parameter for the globular clusters in M87 is small, varying from -0.2 at the centre to 0.2 at 40 kpc, and gradually decreasing to zero at 120 kpc.

IRIS (InfraRed Imaging Spectrograph) is a first light near-infrared diffraction limited imager and integral field spectrograph being designed for the future Thirty Meter Telescope (TMT). IRIS is optimized to perform astronomical studies across a significant fraction of cosmic time, from our Solar System to distant newly formed galaxies (Barton et al. [1]). We present a selection of the innovative science cases that are unique to IRIS in the era of upcoming space and ground-based telescopes. We focus on integral field spectroscopy of directly imaged exoplanet atmospheres, probing fundamental physics in the Galactic Center, measuring 10^4 to 10^10 Msun supermassive black hole masses, resolved spectroscopy of young star-forming galaxies (1 < z < 5) and first light galaxies (6 < z < 12), and resolved spectroscopy of strong gravitational lensed sources to measure dark matter substructure. For each of these science cases we use the IRIS simulator (Wright et al. [2], Do et al. [3]) to explore IRIS capabilities. To highlight the unique IRIS capabilities, we also update the point and resolved source sensitivities for the integral field spectrograph (IFS) in all five broadband filters (Z, Y, J, H, K) for the finest spatial scale of 0.004" per spaxel. We briefly discuss future development plans for the data reduction pipeline and quicklook software for the IRIS instrument suite.

We undertake the first study of two-planet microlensing models recovered from simulations of microlensing events generated by realistic multi-planet systems in which 292 planetary events including 16 two-planet events were detected from 6690 simulated light curves. We find that when two planets are recovered, their parameters are usually close to those of the two planets in the system most responsible for the perturbations. However, in one of the 16 examples, the apparent mass of both detected planets was more than doubled by the unmodeled influence of a third, massive planet. This fraction is larger than, but statistically consistent with, the roughly 1.5% rate of serious mass errors due to unmodeled planetary companions for the 274 cases from the same simulation in which a single planet is recovered. We conjecture that an analogous effect due to unmodeled stellar companions may occur more frequently. For seven out of 23 cases in which two planets in the system would have been detected separately, only one planet was recovered because the perturbations due to the two planets had similar forms. This is a small fraction (7/274) of all recovered single-planet models, but almost a third of all events that might plausibly have led to two-planet models. Still, in these cases, the recovered planet tends to have parameters similar to one of the two real planets most responsible for the anomaly.

We conduct the first microlensing simulation in the context of planet formation model. The planet population is taken from the Ida & Lin core accretion model for $0.3M_\odot$ stars. With $6690$ microlensing events, we find for a simplified Korea Microlensing Telescopes Network (KMTNet) the fraction of planetary events is $2.9\%$ , out of which $5.5\%$ show multiple-planet signatures. The number of super-Earths, super-Neptunes and super-Jupiters detected are expected to be almost equal. Our simulation shows that high-magnification events and massive planets are favored by planet detections, which is consistent with previous expectation. However, we notice that extremely high-magnification events are less sensitive to planets, which is possibly because the 10 min sampling of KMTNet is not intensive enough to capture the subtle anomalies that occur near the peak. This suggests that while KMTNet observations can be systematically analyzed without reference to any follow-up data, follow-up observations will be essential in extracting the full science potential of very high-magnification events. The uniformly high-cadence observations expected for KMTNet also result in $\sim 55\%$ of all detected planets being non-caustic-crossing, and more low-mass planets even down to Mars-mass being detected via planetary caustics. We also find that the distributions of orbital inclinations and planet mass ratios in multiple-planet events agree with the intrinsic distributions.

Recently the commissioning APOGEE observations of the Galactic bulge reported that a significant fraction of stars ($\sim10%$) are in a cold ($\sigma_{\rm V} \approx 30$ km/s) high velocity peak (Galactocentric radial velocity $\approx 200$ km/s). These stars are speculated to reflect the stellar orbits in the Galactic bar. In this study, we use two $N$-body models of the Milky Way-like disk galaxy with different bar strengths to critically examine this possibility. The general trends of the Galactocentric radial velocity distribution in observations and simulations are similar, but neither our models nor the BRAVA data reveal a statistically significant cold high velocity peak. A Monte Carlo test further suggests that it is possible for a spurious high velocity peak to appear if there are only a limited number of stars observed. Thus, the reported cold high velocity peak, even if it is real, is unlikely due to stars on the bar-supporting orbits. Our models do predict an excess of stars with high radial velocity, but not in a distinct peak. In the distance--velocity diagram, the high velocity particles in different fields exist at a similar distance $\sim8.5 \pm 1$ kpc away from the Sun. This result may be explained with geometric intersections between the line-of-sight and the particle orbits; high velocity stars naturally exist approximately at the tangent point, without constituting a distinct peak. We further demonstrate that even without the presence of a bar structure, particle motions in an axisymmetric disk can also exhibit an excess of high velocity stars.

Sub-parsec binary massive black holes (BBHs) are long anticipated to exist in many QSOs but remain observationally elusive. In this paper, we propose a novel method to probe sub-parsec BBHs through microlensing of lensed QSOs. If a QSO hosts a sub-parsec BBH in its center, it is expected that the BBH is surrounded by a circum-binary disk, each component of the BBH is surrounded by a small accretion disk, and a gap is opened by the secondary component in between the circum-binary disk and the two small disks. Assuming such a BBH structure, we generate mock microlensing light curves for some QSO systems that host BBHs with typical physical parameters. We show that microlensing light curves of a BBH QSO system at the infrared-optical-UV bands can be significantly different from those of corresponding QSO system with a single massive black hole (MBH), mainly because of the existence of the gap and the rotation of the BBH (and its associated small disks) around the center of mass. We estimate the half-light radii of the emission region at different wavelengths from mock light curves and find that the obtained half-light radius vs. wavelength relations of BBH QSO systems can be much flatter than those of single MBH QSO systems at a wavelength range determined by the BBH parameters, such as the total mass, mass ratio, separation, accretion rates, etc. The difference is primarily due to the existence of the gap. Such unique features on the light curves and half-light radius-wavelength relations of BBH QSO systems can be used to select and probe sub-parsec BBHs in a large number of lensed QSOs to be discovered by current and future surveys, including the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), the Large Synoptic Survey telescope (LSST) and Euclid.

At a distance of 8 kpc, the center of our Galaxy is the nearest galactic nucleus, and has been the subject of numerous key projects undertaken by great observatories such as Chandra, Spitzer, and Herschel. However, there are still no surveys of molecular gas properties in the Galactic center with less than 30" (1 pc) resolution. There is also no sensitive polarization survey of this region, despite numerous nonthermal magnetic features apparently unique to the central 300 parsecs. In this paper, we outline the potential the VLASS has to fill this gap. We assess multiple considerations in observing the Galactic center, and recommend a C-band survey with 10 micro-Jy continuum RMS and sensitive to molecular gas with densities greater than 10^4 cm^{-3}, covering 17 square degrees in both DnC and CnB configurations ( resolution ~5"), totaling 750 hours of observing time. Ultimately, we wish to note that the upgraded VLA is not just optimized for fast continuum surveys, but has a powerful correlator capable of simultaneously observing continuum emission and dozens of molecular and recombination lines. This is an enormous strength that should be fully exploited and highlighted by the VLASS, and which is ideally suited for surveying the center of our Galaxy.

A VLA Sky Survey of the extragalactic sky at S band (2-4 GHz) with polarization information can uniquely probe the magneto-ionic medium in a wide range of astrophysical environments over cosmic time. For a shallow all-sky survey, we expect to detect over 4 million sources in total intensity $>$ 0.45 mJy beam$^{-1}$ and over 2.2$\times$10$^5$ sources in polarized intensity. With these new observations, we expect to discover new classes of polarized radio sources in very turbulent astrophysical environments and those with extreme values of Faraday depth. Moreover, by determining reliable Faraday depths and by modeling depolarization effects, we can derive properties of the magneto-ionic medium associated with AGNs, absorption line systems and galaxies, addressing the following unresolved questions: (1) What is the covering fraction, the degree of turbulence and the origin of absorption line systems? (2) What is the thermal content in AGNs and radio galaxies? (3) How do AGNs and galaxies evolve over cosmic time? (4) What causes the increase in percentage polarization with decreasing flux densities at the low flux density end of the polarized source count? (5) What is the growth rate of large-scale magnetic fields in galaxies?

We study microlensing light curves by a triple lens, in particular, by a primary star plus two planets. A four-fold degeneracy is confirmed in the light curves, similar to the close and wide degeneracy found in a double lens. Furthermore, we derive a set of equations for triple-lens in the external shear approximation. By using these external shear equations, we identify two kinds of continuous degeneracies which may confused double and triple lenses, i.e. the continuous external shear degeneracy among triple-lens systems and the double-triple lens degeneracy. These degeneracies are particularly important in high magnification events, and thus some caution needs to be applied when one infers the fraction of stars hosting multiple planets from microlensing. We study the dependence of the degeneracies on the lensing parameters (e.g., source trajectory) and give recipes about how the degeneracies should be explored with real data.

We study gravitational lensing when plasma surrounds the lens. An extra deflection angle is induced by the plasma in addition to the deflection generated by gravity. An inhomogeneous plasma distribution generates a greater effect than a homogeneous one, and may cause significant effects to be detected in low frequency radio observations (a few hundred MHz). In particular, the lensed image positions will be different between optical and radio observations. The change of position due to a plasma can reach a few tens of milli-arcsec, which is readily detectable. One can use the position difference in different frequencies to estimate the density of plasma in the lens. The magnification ratios between multiple images are mainly determined by other properties of the lens, and are only weakly affected by the plasma. More importantly, we find that the strong lensing time delay will be affected by the plasma. Estimation of the Hubble constant from the time delay in low radio frequency observation may be slightly biased due to plasma in the lens. Unfortunately, the ionosphere of the Earth strongly affects low frequency radio observations. Thus our ability to detect the effect depends on how well we are able to calibrate out the ionosphere.

We study three-dimensional microlensing where two lenses are located at different distances along the line of sight. We formulate the lens equation in complex notations and recover several previous results. There are in total either 4 or 6 images, with an equal number of images with positive and negative parities. We find that the sum of signed magnifications for six image configurations is unity. Furthermore, we show that the light curves can be qualitatively different from those for binary lensing in a single plane. In particular, the magnifications between a `U'-shaped caustic crossing can be close to unity, rather than having a minimum magnification of 3 as in the single plane binary lensing. There is only a small probability three-dimensional microlensing events will be seen in microlensing toward the Galactic centre. It is more likely they will be first seen in cosmological microlensing.

We study the spatial distribution of the Fe 6.4 and 6.7 keV lines in the nuclear region of M82 using the Chandra archival data with a total exposure time of 500 ks. The deep exposure provides a significant detection of the Fe 6.4 keV line. Both the Fe 6.4 and 6.7 keV lines are diffuse emissions with similar spatial extent, but their morphology do not exactly follow each other. Assuming a thermal collisional-ionization-equilibrium model, the fitted temperatures are around 5-6 keV and the Fe abundances are about 0.4-0.6 solar value. We also report the spectrum of a point source, which shows a strong Fe 6.7 keV line and is likely a supernova remnant or a superbubble. The fitted Fe abundance of the point source is 1.7 solar value. It implies that part of the iron may be depleted from the X-ray emitting gases as the predicted abundance is about 5 solar value assuming complete mixing. If this is a representative case of the Fe enrichment, a mild mass-loading of a factor of 3 will make the Fe abundance of the point source in agreement with that of the hot gas, which then implies that most of the hard X-ray continuum (2-8 keV) of M82 has a thermal origin. In addition, the Fe 6.4 keV line is consistent with the fluorescence emission irradiated by the hard photons from nuclear point sources.

We report the first detection of magnetic fields associated with the Smith High Velocity Cloud. We use a catalog of Faraday rotation measures towards extragalactic radio sources behind the Smith Cloud, new HI observations from the Green Bank Telescope, and a spectroscopic map of H{\alpha} from the Wisconsin H-Alpha Mapper Northern Sky Survey. There are enhancements in rotation measure of approximately 100 rad m^(-2) which are generally well correlated with decelerated H{\alpha} emission. We estimate a lower limit on the line-of-sight component of the field of approximately 8 {\mu}G along a decelerated filament; this is a lower limit due to our assumptions about the geometry. No RM excess is evident in sightlines dominated by HI or H{\alpha} at the velocity of the Smith Cloud. The smooth H{\alpha} morphology of the emission at the Smith Cloud velocity suggests photoionization by the Galactic ionizing radiation field as the dominant ionization mechanism, while the filamentary morphology and high (approximately 1 Rayleigh) H{\alpha} intensity of the lower-velocity magnetized ionized gas suggests an ionization process associated with shocks due to interaction with the Galactic interstellar medium. The presence of the magnetic field may contribute to the survival of high velocity clouds like the Smith Cloud as they move from the Galactic halo to the disk. We expect these data to provide a test for magnetohydrodynamic simulations of infalling gas.

Using the potential from N-body simulations, we construct the Galactic bar models with the Schwarzschild method. By varying the pattern speed and the position angle of the bar, we find that the best-fit bar model has pattern speed $\Omega_{\rm p}=40\ \rm{km\ s^{-1}\ kpc^{-1}}$, and bar angle $\theta_{\rm bar}=45^{\circ}$. $N$-body simulations show that the best-fit model is stable for more than 1.5 Gyrs. Combined with the results in Wang et al. (2012), we find that the bar angle and/or the pattern speed are not well constrained by BRAVA data in our Schwarzschild models. The proper motions predicted from our model are slightly larger than those observed in four fields. In the future, more kinematic data from the ground and space-based observations will enable us to refine our model of the Milky Way bar.

The China Dark matter Experiment collaboration reports the first experimental limit on WIMP dark matter from 14.6 kg-day of data taken with a 994 g p-type point-contact germanium detector at the China Jinping underground Laboratory where the rock overburden is more than 2400 m. The energy threshold achieved was 400 eVee. According to the 14.6 kg-day live data, we placed the limit of N= 1.75 * 10^{-40} cm^{2} at 90% confidence level on the spin-independent cross-section at WIMP mass of 7 GeV before differentiating bulk signals from the surface backgrounds.

We present a study of the luminosity density distribution of the Galactic bar using number counts of red clump giants (RCGs) from the OGLE-III survey. The data were recently published by Nataf et al. (2013) for 9019 fields towards the bulge and have $2.94\times 10^6$ RC stars over a viewing area of $90.25 \,\textrm{deg}^2$. The data include the number counts, mean distance modulus ($\mu$), dispersion in $\mu$ and full error matrix, from which we fit the data with several tri-axial parametric models. We use the Markov Chain Monte Carlo (MCMC) method to explore the parameter space and find that the best-fit model is the $E_3$ model, with the distance to the GC is 8.13 kpc, the ratio of semi-major and semi-minor bar axis scale lengths in the Galactic plane $x_{0},y_{0}$, and vertical bar scale length $z_0$, is $x_0:y_0:z_0 \approx 1.00:0.43:0.40$ (close to being prolate). The scale length of the stellar density profile along the bar's major axis is $\sim$ 0.67 kpc and has an angle of $29.4^\circ$, slightly larger than the value obtained from a similar study based on OGLE-II data. The number of estimated RC stars within the field of view is $2.78 \times 10^6$, which is systematically lower than the observed value. We subtract the smooth parametric model from the observed counts and find that the residuals are consistent with the presence of an X-shaped structure in the Galactic centre, the excess to the estimated mass content is $\sim 5.8%$. We estimate the total mass of the bar is $\sim 1.8 \times 10^{10} M_\odot$. Our results can be used as a key ingredient to construct new density models of the Milky Way and will have implications on the predictions of the optical depth to gravitational microlensing and the patterns of hydrodynamical gas flow in the Milky Way.

We perform photometric analyses for a bright early-type galaxy (ETG) sample with 2949 galaxies ($M_{\rm r}<-22.5$ mag) in the redshift range of 0.05 to 0.15, drawn from the SDSS DR7 with morphological classification from Galaxy Zoo 1. We measure the Petrosian and isophotal magnitudes, as well as the corresponding half-light radius for each galaxy. We find that for brightest galaxies ($M_{\rm r}<-23$ mag), our Petrosian magnitudes, and isophotal magnitudes to 25 ${\rm mag/arcsec^2}$ and 1\% of the sky brightness are on average 0.16 mag, 0.20 mag, and 0.26 mag brighter than the SDSS Petrosian values, respectively. In the first case the underestimations are caused by overestimations in the sky background by the SDSS PHOTO algorithm, while the latter two are also due to deeper photometry. Similarly, the typical half-light radii ($r_{50}$) measured by the SDSS algorithm are smaller than our measurements. As a result, the bright-end of the $r$-band luminosity function is found to decline more slowly than previous works. Our measured luminosity densities at the bright end are more than one order of magnitude higher than those of Blanton et al. (2003), and the stellar mass densities at $M_{\ast}\sim 5\times10^{11} M_{\odot}$ and $M_{\ast}\sim 10^{12} M_{\odot}$ are a few tenths and a factor of few higher than those of Bernardi et al. (2010). These results may significantly alleviate the tension in the assembly of massive galaxies between observations and predictions of the hierarchical structure formation model.

The CDEX Collaboration has been established for direct detection of light dark matter particles, using ultra-low energy threshold p-type point-contact germanium detectors, in China JinPing underground Laboratory (CJPL). The first 1 kg point-contact germanium detector with a sub-keV energy threshold has been tested in a passive shielding system located in CJPL. The outputs from both the point-contact p+ electrode and the outside n+ electrode make it possible to scan the lower energy range of less than 1 keV and at the same time to detect the higher energy range up to 3 MeV. The outputs from both p+ and n+ electrode may also provide a more powerful method for signal discrimination for dark matter experiment. Some key parameters, including energy resolution, dead time, decay times of internal X-rays, and system stability, have been tested and measured. The results show that the 1 kg point-contact germanium detector, together with its shielding system and electronics, can run smoothly with good performances. This detector system will be deployed for dark matter search experiments.

{\it Spitzer} images of Tycho's supernova remnant in the mid-infrared reveal limb-brightened emission from the entire periphery of the shell and faint filamentary structures in the interior. As with other young remnants, this emission is produced by dust grains, warmed to $\sim 100$ K in the post-shock environment by collisions with energetic electrons and ions. The ratio of the 70 to 24 $\mu$m fluxes is a diagnostic of the dust temperature, which in turn is a sensitive function of the plasma density. We find significant variations in the 70/24 flux ratio around the periphery of Tycho's forward shock, implying order-of-magnitude variations in density. While some of these are likely localized interactions with dense clumps of the interstellar medium, we find an overall gradient in the ambient density surrounding Tycho, with densities 3-10 times higher in the NE than in the SW. This large density gradient is qualitatively consistent with the variations in the proper motion of the shock observed in radio and X-ray studies. Overall, the mean ISM density around Tycho is quite low ($\sim 0.1-0.2$ cm$^{-3}$), consistent with the lack of thermal X-ray emission observed at the forward shock. We perform two-dimensional hydrodynamic simulations of a Type Ia SN expanding into a density gradient in the ISM, and find that the overall round shape of the remnant is still easily acheivable, even for explosions into significant gradients. However, this leads to an offset of the center of the explosion from the geometric center of the remnant of up to 20%, although lower values of 10% are preferred. The best match with hydrodynamical simulations is achieved if Tycho is located at a large (3-4 kpc) distance in a medium with a mean preshock density of $\sim 0.2$ cm$^{-3}$. Such preshock densities are obtained for highly ($\ga 50$%) porous ISM grains.

We have made a serendipitous discovery of a massive cD galaxy at z=1.096 in a candidate rich cluster in the HUDF area of GOODS-South. This brightest cluster galaxy is the most distant cD galaxy confirmed to date. Ultra-deep HST/WFC3 images reveal an extended envelope starting from ~10 kpc and reaching ~70 kpc in radius along the semi-major axis. The spectral energy distributions indicate that both its inner component and outer envelope are composed of an old, passively-evolving stellar population. The cD galaxy lies on the same mass-size relation as the bulk of quiescent galaxies at similar redshifts. The cD galaxy has a higher stellar mass surface density but a similar velocity dispersion to those of more-massive, nearby cDs. If the cD galaxy is one of the progenitors of today's more massive cDs, its size and stellar mass have had to increase on average by factors of $3.4\pm1.1$ and $3.3\pm1.3$ over the past ~8 Gyrs, respectively. Such increases in size and stellar mass without being accompanied by significant increases in velocity dispersion are consistent with evolutionary scenarios driven by both major and minor dry mergers. If such cD envelopes originate from dry mergers, our discovery of even one example proves that some BCGs entered the dry merger phase at epochs earlier than z=1. Our data match theoretical models which predict that the continuance of dry mergers at z<1 can result in structures similar to those of massive cD galaxies seen today. Moreover, our discovery is a surprise given that the extreme depth of the HUDF is essential to reveal such an extended cD envelope at z>1 and, yet, the HUDF covers only a minuscule region of sky. Adding that cDs are rare, Our serendipitous discovery hints that such cDs may be more common than expected. [Abridged]

Euclid is the second M-class mission of the ESA Cosmic Vision programme, with the principal science goal of studying dark energy. Euclid is also expected to undertake additional Legacy Science programmes. One proposal is the Exoplanet Euclid Legacy Survey (ExELS) which will be the first survey able to measure the abundance of exoplanets down to Earth mass for host separations from ~1AU out to the free-floating (unbound) regime. The cold and free-floating exoplanet regimes represent a crucial discovery space for testing planet formation theories. ExELS will use the gravitational microlensing technique and will detect 1000 microlensing events per month over 1.6 deg^2 of the Galactic bulge. We assess how many of these events will have detectable planetary signatures using a detailed multi-wavelength microlensing simulator (MABuLS) which incorporates the Besancon Galactic model with 3D extinction. MABuLS is the first theoretical simulation of microlensing to treat the effects of point spread function (PSF) blending self-consistently with the underlying Galactic model. We use MABuLS, together with current numerical models for the Euclid PSFs, to explore a number of designs and de-scope options for ExELS, including the exoplanet yield as a function of filter choice and slewing time, and the effect of systematic photometry errors. Using conservative extrapolations of current empirical exoplanet mass functions determined from ground-based microlensing and radial velocity surveys, ExELS can expect to detect a few hundred cold exoplanets around mainly G, K and M-type stellar hosts, including ~45 Earth-mass planets and ~6 Mars-mass planets for an observing programme totalling 10 months. ExELS will be capable of measuring the cold exoplanet mass function down to Earth mass or below, with orbital separations from ~1AU to the free-floating regime. (Abridged)

We study the lensing magnification effect on background galaxies. Differential magnification due to different magnifications of different source regions of a galaxy will change the lensed composite spectra. The derived properties of the background galaxies are therefore biased. For simplicity, we model galaxies as a superposition of an axis-symmetric bulge and a face-on disk in order to study the differential magnification effect on the composite spectra. We find that some properties derived from the spectra (e.g., velocity dispersion, star formation rate and metallicity) are modified. Depending on the relative positions of the source and the lens, the inferred results can be either over- or under-estimates of the true values. In general, for an extended source at strong lensing regions with high magnifications, the inferred physical parameters (e.g. metallicity) can be strongly biased. Therefore detailed lens modelling is necessary to obtain the true properties of the lensed galaxies.

We present a Faraday rotation measure (RM) study of the diffuse, polarized, radio emission from the giant lobes of the nearest radio galaxy, Centaurus A. After removal of the smooth Galactic foreground RM component, using an ensemble of background source RMs located outside the giant lobes, we are left with a residual RM signal associated with the giant lobes. We find the most likely origin of this residual RM is from thermal material mixed throughout the relativistic lobe plasma. The alternative possibility of a thin-skin/boundary layer of magnetoionic material swept up by the expansion of the lobes is highly unlikely since it requires, at least, an order of magnitude enhancement of the swept up gas over the expected intragroup density on these scales. Strong depolarisation observed from 2.3 to 0.96 GHz also supports the presence of a significant amount of thermal gas within the lobes; although depolarisation solely due to RM fluctuations in a foreground Faraday screen on scales smaller than the beam cannot be ruled out. Considering the internal Faraday rotation scenario, we find a thermal gas number density of ~10^{-4} cm^{-3} implying a total gas mass of ~10^{10} M_solar within the lobes. The thermal pressure associated with this gas (with temperature kT ~ 0.5 keV, obtained from recent X-ray results) is approximately equal to the non-thermal pressure, indicating that over the volume of the lobes, there is approximate equipartition between the thermal gas, radio-emitting electrons and magnetic field (and potentially any relativistic protons present).

We study the effects on the number counts of sub-millimetre galaxies due to gravitational lensing. We explore the effects on the magnification cross section due to halo density profiles, ellipticity and cosmological parameter (the power-spectrum normalisation $\sigma_8$). We show that the ellipticity does not strongly affect the magnification cross section in gravitational lensing while the halo radial profiles do. Since the baryonic cooling effect is stronger in galaxies than clusters, galactic haloes are more concentrated. In light of this, a new scenario of two halo population model is explored where galaxies are modeled as a singular isothermal sphere profile and clusters as a Navarro, Frenk and White (NFW) profile. We find the transition mass between the two has modest effects on the lensing probability. The cosmological parameter $\sigma_8$ alters the abundance of haloes and therefore affects our results. Compared with other methods, our model is simpler and more realistic. The conclusions of previous works is confirm that gravitational lensing is a natural explanation for the number count excess at the bright end.

We demonstrate how the Syer & Tremaine made-to-measure method of stellar dynamical modelling can be adapted to model a rotating galactic bar. We validate our made-to-measure changes using observations constructed from the existing Shen et al. (2010) N-body model of the Milky Way bar, together with kinematic observations of the Milky Way bulge and bar taken by the Bulge Radial Velocity Assay (BRAVA). Our results for a combined determination of the bar angle and bar pattern speed (~30 degrees and ~40 km/s/kpc) are consistent with those determined by the N-body model. Whilst the made-to-measure techniques we have developed are applied using a particular N-body model and observational data set, they are in fact general and could be applied to other Milky Way modelling scenarios utilising different N-body models and data sets. Additionally, we use the exercise as a vehicle for illustrating how N-body and made-to-measure methods might be combined into a more effective method.

The masses of stars including stellar remnants are almost exclusively known from binary systems. In this work, we study gravitational microlensing of faint background galaxies by isolated neutron stars (pulsars). We show that the resulting surface brightness distortions can be used to determine the masses of neutron star. Due to different evolutionary histories, isolated neutron stars may have different masses from those in binary systems, and thus provide unique insight into their equation of states under extreme physical conditions. We search for existing pulsar catalogs and find one promising pair of a nearby pulsar and a background galaxy. This method will become more practical for the next generation optical and radio surveys and telescopes.

We use Schwarzschild's orbit-superposition technique to construct self-consistent models of the Galactic bar. Using $\chi^2$ minimisation, we find that the best-fit Galactic bar model has a pattern speed $\Omega_{\rm p}=60 \rm{km s^{-1} kpc^{-1}}$, disk mass $\rm{M_{\rm d}=1.0\times10^{11}M_{\odot}}$ and bar angle $\theta_{\rm bar}=20^{\circ}$ for an adopted bar mass $\rm{M_{\rm bar}=2\times10^{10}M_{\odot}}$. The model can reproduce not only the three-dimensional and projected density distributions but also velocity and velocity dispersion data from the BRAVA survey. We also predict the proper motions in the range $l=[-12^{\circ},12^{\circ}]$, $b=[-10^{\circ},10^{\circ}]$, which appear to be higher than observations in the longitudinal direction. The model is stable within a timescale of 0.5 Gyr, but appears to deviate from steady-state on longer timescales. Our model can be further tested by future observations such as those from GAIA.

We present a study of the magnetic field of the Large Magellanic Cloud (LMC), carried out using diffuse polarized synchrotron emission data at 1.4 GHz acquired at the Parkes Radio Telescope and the Australia Telescope Compact Array. The observed diffuse polarized emission is likely to originate above the LMC disk on the near side of the galaxy. Consistent negative rotation measures (RMs) derived from the diffuse emission indicate that the line-of-sight magnetic field in the LMC's near-side halo is directed coherently away from us. In combination with RMs of extragalactic sources that lie behind the galaxy, we show that the LMC's large scale magnetic field is likely to be of quadrupolar geometry, consistent with the prediction of dynamo theory. On smaller scales, we identify two brightly polarized filaments southeast of the LMC, associated with neutral hydrogen arms. The filaments' magnetic field potentially aligns with the direction towards the Small Magellanic Cloud. We suggest that tidal interactions between the Small and the Large Magellanic Clouds in the past 10^9 years is likely to have shaped the magnetic field in these filaments.

We present measurements of the velocity dispersion profile (VDP) for galaxy groups in the final data release of the Sloan Digital Sky Survey (SDSS). For groups of given mass we estimate the redshift-space cross-correlation function (CCF) with respect to a reference galaxy sample, xi(r_p, pi), the projected CCF, w_p(r_p), and the real-space CCF, xi(r). The VDP is then extracted from the redshift distortion in xi(r_p, pi), by comparing xi(r_p, pi) with xi(r). We find that the velocity dispersion (VD) within virial radius (R_200) shows a roughly flat profile, with a slight increase at radii below ~0.3 R_200 for high mass systems. The average VD within the virial radius, sigma_v, is a strongly increasing function of central galaxy mass. We apply the same methodology to N-body simulations with the concordance Lambda cold dark matter cosmology but different values of the density fluctuation parameter sigma_8, and we compare the results to the SDSS results. We show that the sigma_v-M_* relation from the data provides stringent constraints on both sigma_8 and sigma_ms, the dispersion in log M_* of central galaxies at fixed halo mass. Our best-fitting model suggests sigma_8 = 0.86 +/- 0.03 and sigma_ms = 0.16 +/- 0.03. The slightly higher value of sigma_8 compared to the WMAP7 result might be due to a smaller matter density parameter assumed in our simulations. Our VD measurements also provide a direct measure of the dark matter halo mass for central galaxies of different luminosities and masses, in good agreement with the results obtained by Mandelbaum et al. (2006) from stacking the gravitational lensing signals of the SDSS galaxies.

Since the first discovery of microlensing events nearly two decades ago, gravitational microlensing has accumulated tens of TBytes of data and developed into a powerful astrophysical technique with diverse applications. The review starts with a theoretical overview of the field and then proceeds to discuss the scientific highlights. (1) Microlensing observations toward the Magellanic Clouds rule out the Milky Way halo being dominated by MAssive Compact Halo Objects (MACHOs). This confirms most dark matter is non-baryonic, consistent with other observations. (2) Microlensing has discovered about 20 extrasolar planets (16 published), including the first two Jupiter-Saturn like systems and the only "cold Neptunes" yet detected. They probe a different part of the parameter space and will likely provide the most stringent test of core accretion theory of planet formation. (3) Microlensing provides a unique way to measure the mass of isolated stars, including brown dwarfs to normal stars. Half a dozen or so stellar mass black hole candidates have also been proposed. (4) High-resolution, target-of-opportunity spectra of highly-magnified dwarf stars provide intriguing "age" determinations which may either hint at enhanced helium enrichment or unusual bulge formation theories. (5) Microlensing also measured limb-darkening profiles for close to ten giant stars, which challenges stellar atmosphere models. (6) Data from surveys also provide strong constraints on the geometry and kinematics of the Milky Way bar (through proper motions); the latter indicates predictions from current models appear to be too anisotropic compared with observations. The future of microlensing is bright given the new capabilities of current surveys and forthcoming new telescope networks from the ground and from space. Some open issues in the field are identified and briefly discussed.

We present a study of the Milky Way disk and halo magnetic field, determined from observations of Faraday rotation measure (RM) towards 641 polarized extragalactic radio sources in the Galactic longitude range 100-117 degs, within 30 degs of the Galactic plane. For |b| < 15 degs, we observe a symmetric RM distribution about the Galactic plane. This is consistent with a disk field in the Perseus arm of even parity across the Galactic mid-plane. In the range 15<|b|<30 degs, we find median rotation measures of -15+/-4 rad/m^2 and -62+/-5 rad/m^2 in the northern and southern Galactic hemispheres, respectively. If the RM distribution is a signature of the large-scale field parallel to the Galactic plane, this suggests that the halo magnetic field toward the outer Galaxy does not reverse direction across the mid-plane. The variation of RM as a function of Galactic latitude in this longitude range is such that RMs become more negative at larger |b|. This is consistent with an azimuthal magnetic field of strength 2 microGauss (7 microGauss) at a height 0.8-2 kpc above (below) the Galactic plane between the local and the Perseus spiral arm. We propose that the Milky Way could possess spiral-like halo magnetic fields similar to those observed in M51.

We explore the dynamical signatures imprinted by baryons on dark matter haloes during the formation process using the OverWhelmingly Large Simulations (OWLS), a set of state-of-the-art high-resolution cosmological hydrodynamical simulations. We present a detailed study of the effects of the implemented feedback prescriptions on the orbits of dark matter particles, stellar particles and subhaloes, analysing runs with no feedback, with stellar feedback and with feedback from supermassive black holes. We focus on the central regions (0.25 r_{200}) of haloes with virial masses ~ 6 x 10^{13} (~ 7 x 10^{11}) Msun/h at z = 0(2). We also investigate how the orbital content (relative fractions of the different orbital types) of these haloes depends on several key parameters such as their mass, redshift and dynamical state. The results of spectral analyses of the orbital content of these simulations are compared, and the change in fraction of box, tube and irregular orbits is quantified. Box orbits are found to dominate the orbital structure of dark matter haloes in cosmological simulations. There is a strong anticorrelation between the fraction of box orbits and the central baryon fraction. While radiative cooling acts to reduce the fraction of box orbits, strong feedback implementations result in a similar orbital distribution to that of the dark matter only case. The orbital content described by the stellar particles is found to be remarkably similar to that drawn from the orbits of dark matter particles, suggesting that either they have forgotten their dynamical history, or that subhaloes bringing in stars are not biased significantly with respect to the main distribution. The orbital content of the subhaloes is in broad agreement with that seen in the outer regions of the particle distributions.

We present a spatial analysis of the soft X-ray and H{\alpha} emissions from the outflow of the starburst galaxy M82. We find that the two emissions are tightly correlated on various scales. The O VII triplet of M82, as resolved by X-ray grating observations of XMM-Newton, is dominated by the forbidden line, inconsistent with the thermal prediction. The O VII triplet also shows some spatial variations. We discuss three possible explanations for the observed O VII triplet, including the charge exchange at interfaces between the hot outflow and neutral cool gas, a collisional non-equilibrium-ionization recombining plasma, and resonance scattering.

We identify a total of 120 early-type Brightest Cluster Galaxies (BCGs) at 0.1<z<0.4 in two recent large cluster catalogues selected from the Sloan Digital Sky Survey (SDSS). They are selected with strong emission lines in their optical spectra, with both H{\alpha} and [O II]{\lambda}3727 line emission, which indicates significant ongoing star formation. They constitute about ~ 0.5% of the largest, optically-selected, low-redshift BCG sample, and the fraction is a strong function of cluster richness. Their star formation history can be well described by a recent minor and short starburst superimposed on an old stellar component, with the recent episode of star formation contributing on average only less than 1 percent of the total stellar mass. We show that the more massive star-forming BCGs in richer clusters tend to have higher star formation rate (SFR) and specific SFR (SFR per unit galaxy stellar mass). We also compare their statistical properties with a control sample selected from X-ray luminous clusters, and show that the fraction of star-forming BCGs in X-ray luminous clusters is almost one order of magnitude larger than that in optically-selected clusters. BCGs with star formation in cooling flow clusters usually have very flat optical spectra and show the most active star formation, which may be connected with cooling flows.

The microlensing event OGLE-2008-BLG-510 is characterised by an evident asymmetric shape of the peak, promptly detected by the ARTEMiS system in real time. The skewness of the light curve appears to be compatible both with binary-lens and binary-source models, including the possibility that the lens system consists of an M dwarf orbited by a brown dwarf. The detection of this microlensing anomaly and our analysis demonstrates that: 1) automated real-time detection of weak microlensing anomalies with immediate feedback is feasible, efficient, and sensitive, 2) rather common weak features intrinsically come with ambiguities that are not easily resolved from photometric light curves, 3) a modelling approach that finds all features of parameter space rather than just the `favourite model' is required, and 4) the data quality is most crucial, where systematics can be confused with real features, in particular small higher-order effects such as orbital motion signatures. It moreover becomes apparent that events with weak signatures are a silver mine for statistical studies, although not easy to exploit. Clues about the apparent paucity of both brown-dwarf companions and binary-source microlensing events might hide here.

We aim to summarize the current state of knowledge regarding Galactic Faraday rotation in an all-sky map of the Galactic Faraday depth. For this we have assembled the most extensive catalog of Faraday rotation data of compact extragalactic polarized radio sources to date. In the map making procedure we use a recently developed algorithm that reconstructs the map and the power spectrum of a statistically isotropic and homogeneous field while taking into account uncertainties in the noise statistics. This procedure is able to identify some rotation angles that are offset by an integer multiple of pi. The resulting map can be seen as an improved version of earlier such maps and is made publicly available, along with a map of its uncertainty. For the angular power spectrum we find a power law behavior with a power law index of -2.14 for a Faraday sky where an overall variance profile as a function of Galactic latitude has been removed, in agreement with earlier work. We show that this is in accordance with a 3D Fourier power spectrum P(k) proportional to k^-2.14 of the underlying field n_e times B_r under simplifying geometrical and statistical assumptions.

We report CO detections in 17 out of 19 infrared ultraluminous QSO (IR QSO) hosts observed with the IRAM 30m telescope. The cold molecular gas reservoir in these objects is in a range of 0.2--2.1$\times 10^{10}M_\odot$ (adopting a CO-to-${\rm H_2}$ conversion factor $\alpha_{\rm CO}=0.8 M_\odot {\rm (K km s^{-1} pc^2)^{-1}}$). We find that the molecular gas properties of IR QSOs, such as the molecular gas mass, star formation efficiency ($L_{\rm FIR}/L^\prime_{\rm CO}$) and the CO (1-0) line widths, are indistinguishable from those of local ultraluminous infrared galaxies (ULIRGs). A comparison of low- and high-redshift CO detected QSOs reveals a tight correlation between L$_{\rm FIR}$ and $L^\prime_{\rm CO(1-0)}$ for all QSOs. This suggests that, similar to ULIRGs, the far-infrared emissions of all QSOs are mainly from dust heated by star formation rather than by active galactic nuclei (AGNs), confirming similar findings from mid-infrared spectroscopic observations by {\it Spitzer}. A correlation between the AGN-associated bolometric luminosities and the CO line luminosities suggests that star formation and AGNs draw from the same reservoir of gas and there is a link between star formation on $\sim$ kpc scale and the central black hole accretion process on much smaller scales.

The flux-ratio anomalies observed in multiply-lensed quasar images are most plausibly explained as the result of perturbing structures superposed on the underlying smooth matter distribution of the primary lens. The cold dark matter cosmological model predicts that a large number of substructures should survive inside larger halos but, surprisingly, this population alone has been shown to be insufficient to explain the observed distribution of the flux ratios of quasar's multiple images. Other halos (and their own subhalos) projected along the line of sight to the primary lens have been considered as additional source of perturbation. In this work, we use ray tracing through the Millennium II simulation to investigate the importance of projection effects due to halos and subhalos of mass m>1E8 Msun/h and extend our analysis to lower masses, m>1E6 Msun/h, using Monte-Carlo halo distributions. We find that the magnitude of the violation depends strongly on the density profile and concentration of the intervening halos, but clustering plays only a minor role. For a typical lensing geometry (lens at redshift 0.6 and source at redshift 2), background haloes (behind the main lens) are more likely to cause a violation than foreground halos. We conclude that line-of-sight structures can be as important as intrinsic substructures in causing flux-ratio anomalies. The combined effect of perturbing structures within the lens and along the line of sight in the LCDM universe results in a cusp-violation probability of 20-30%. This alleviates the discrepancy between models and current data, but a larger observational sample is required for a stronger test of the theory.

We study moderate gravitational lensing where a background galaxy is magnified substantially, but not multiply imaged, by an intervening galaxy. We focus on the case where both the lens and source are elliptical galaxies. The signatures of moderate lensing include isophotal distortions and systematic shifts in the fundamental plane and Kormendy relation, which can potentially be used to statistically determine the galaxy mass profiles. These effects are illustrated using Monte Carlo simulations of galaxy pairs where the foreground galaxy is modelled as a singular isothermal sphere model and observational parameters appropriate for the Large Synoptic Survey Telescope (LSST). The range in radius probed by moderate lensing will be larger than that by strong lensing, and is in the interesting regime where the density slope may be changing.

We take a sample of 24 elliptical and lenticular galaxies previously analysed by the SAURON project using three-integral dynamical models created with Schwarzschild's method, and re-analyse them using the made-to-measure (M2M) method of dynamical modelling. We obtain good agreement between the two methods in determining the dynamical mass-to-light (M/L) ratios for the galaxies with over 80% of ratios differing by < 10% and over 95% differing by < 20%. We show that (M/L)_M2M is approximately equal to (M/L)_Sch. For the global velocity dispersion anisotropy parameter delta, we find similar values but with fewer of the made-to-measure models tangentially anisotropic by comparison with their SAURON Schwarzschild counterparts. Our investigation is the largest comparative application of the made-to-measure method to date.

We constrain properties of cluster haloes by performing likelihood analysis using lensing shear and flexion data. We test our analysis using two mock cluster haloes: an isothermal ellipsoid (SIE) model and a more realistic elliptical Navarro-Frenk-White (eNFW) model. For both haloes, we find that flexion is more sensitive to the halo ellipticity than shear. The introduction of flexion information significantly improves the constraints on halo ellipticity, orientation and mass. We also point out that there is a degeneracy between the mass and the ellipticity of SIE models in the lensing signal.

Properties of hot gas outflows from galaxies are generally measured from associated X-ray line emission assuming that it represents atomic transitions in thermally excited hot gas. X-ray line emission, however, can also arise from the charge exchange between highly ionised ions and neutral species. The K\alpha\ triplet of He-like ions can be used as a powerful diagnostic, because the charge exchange X-ray emission (CXE) favours the inter-combination and forbidden lines, while the thermal emission favours the resonance line. We analyse the OVII triplet of a sample of nine nearby star-forming galaxies observed by the XMM-Newton reflection grating spectrometers. For most galaxies, the forbidden lines are comparable to or stronger than the resonance lines, which is in contrast to the thermal prediction. For NGC 253, M51, M83, M61, NGC 4631, and the Antennae galaxy, the observed line ratios are consistent with the ratio of the CXE; for M94 and NGC 2903, the observed ratios indicate multiple origins; for M82, different regions show different line ratios, also indicating multiple origins. We discuss other possible mechanisms that can produce a relatively strong forbidden line, such as a collisional non-equilibrium-ionization recombining/ionizing plasma, which are not favoured. These results suggest that the CXE may be a common process and contribute a significant fraction of the soft X-ray line emission for galaxies with massive star formation.

Most galaxies host central supermassive black holes. As two galaxies merge, the black holes also merge. The final single black hole may suffer a kick due to asymmetric gravitational radiation and may not be at the centre of the galaxy; off-centre black holes may also be produced by other means such as sustained acceleration due to asymmetric jet power. We model the main galaxy as a singular isothermal sphere and the black hole as an off-centre point lens, and study the critical curves and caustics using complex notation. We identify the critical parameters that govern the transitions in the topology of critical curves, caustics and pseudo-caustics, and find the number of images can be two, three, four and five. We show examples of image configurations, including cases where three highly de-magnified images are found close to the centre. The perturbation on the image magnification due to the black hole scales linearly with its mass in the off-centre case, and quadratically when the black hole is at the centre. Such images are difficult to observe unless high-contrast and high-resolution imaging facilities (e.g., the Square Kilometer Array in the radio) become available.

Recent observations indicate that many if not all galaxies host massive central black holes (BHs). In this paper we explore the influence of supermassive binary black holes (SMBBHs) on their actions as gravitational lenses. When lenses are modelled as singular isothermal ellipsoids, binary black holes change the critical curves and caustics differently as a function of distance. Each black hole can in principle create at least one additional image, which, if observed, provides evidence of black holes. By studying how SMBBHs affect the cumulative distribution of magnification for images created by black holes, we find that the cross section for at least one such additional image to have a magnification larger than $10^{-5}$ is comparable to the cross section for producing multiple-images in singular isothermal lenses. Such additional images may be detectable with high-resolution and large dynamic range maps of multiply-imaged systems from future facilities, such as the SKA. The probability of detecting at least one image (two images) with magnification above $10^{-3}$ is $\sim 0.2 \fBH$ ($\sim 0.05 \fBH$) in a multiply-imaged lens system, where $\fBH$ is the fraction of galaxies housing binary black holes. We also study the effects of SMBBHs on the core images when galaxies have shallower central density profiles (modelled as non-singular isothermal ellipsoids). We find that the cross section of the usually faint core images is further suppressed by SMBBHs. Thus their presence should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses.

Microlensing is most sensitive to binary lenses with relatively large orbital separations, and as such, typical binary microlensing events show little or no orbital motion during the event. However, despite the strength of binary microlensing features falling off rapidly as the lens separation decreases, we show that it is possible to detect repeating features in the lightcurve of binary microlenses that complete several orbits during the microlensing event. We investigate the lightcurve features of such Rapidly Rotating Lens (RRL) events. We derive analytical limits on the range of parameters where these effects are detectable, and confirm these numerically. Using a population synthesis Galactic model we estimate the RRL event rate for a ground-based and space-based microlensing survey to be 0.32fb and 7.8fb events per year respectively, assuming year-round monitoring and where fb is the binary fraction. We detail how RRL event parameters can be quickly estimated from their lightcurves, and suggest a method to model RRL events using timing measurements of lightcurve features. Modelling RRL lightcurves will yield the lens orbital period and possibly measurements of all orbital elements including the inclination and eccentricity. Measurement of the period from the lightcurve allows a mass-distance relation to be defined, which when combined with a measurement of microlens parallax or finite source effects, can yield a mass measurement to a two-fold degeneracy. With sub-percent accuracy photometry it is possible to detect planetary companions, but the likelihood of this is very small.

The cold dark matter theory predicts triaxial dark matter haloes. The radial distribution of halo ellipticity depends on baryonic processes and the nature of dark matter particles (collisionless or collisional). Here we show that we can use lensing flexion ratios to measure the halo ellipticity as a function of radius. We introduce a weight function and study the relationship between the first and second order statistics of flexion ratios, both of which can be used to reduce the bias in the estimate of ellipticity. we perform numerical tests for our method, and demonstrate that it can reduce the bias and determine the halo ellipticity as a function of radius. We also point out that the minimum mean flexion ratio can be used to trace the centres of galaxy clusters.

Starburst galaxies are primary feedback sources of mechanical energy and metals, which are generally measured from associated X-ray emission lines assuming that they are from the thermal emission of the outflowing hot gas. Such line emission, however, can also arise from the charge exchange X-ray emission (CXE) between highly ionized ions and neutral species. To understand the feedback of energy and metals, it is crucial to determine the origin of the X-ray emission lines and to distinguish the contributions from the CXE and the thermal emission. The origin of the lines can be diagnosed by the K{\alpha} triplets of He-like ions, because the CXE favors the inter-combination and forbidden lines, while the thermal emission favors the resonance line. We analyze the triplets of O VII, Ne IX, and Mg XI observed in the XMM- Newton reflection grating spectra of the starburst galaxy M82. The flux contribution of the CXE is 90%, 50%, and 30% to the O VII, Ne IX, and Mg XI triplet, respectively. Averaged over all the three triplets, the contribution of the CXE is \sim 50% of the total observed triplet flux. To correctly understand the hot outflow of starburst galaxies, it is necessary to include the CXE. Based on the measured CXE contributions to the O VII, Ne IX, and Mg XI triplets, we estimate the relative abundances of O, Ne, and Mg of the outflow and find they are similar to the solar ratios.

The observed number counts of high-redshift galaxy candidates have been used to build up a statistical description of star-forming activity at redshift z >~ 7, when galaxies reionized the Universe. Standard models predict that a high incidence of gravitational lensing will probably distort measurements of flux and number of these earliest galaxies. The raw probability of this happening has been estimated to be ~ 0.5 percent, but can be larger owing to observational biases. Here we report that gravitational lensing is likely to dominate the observed properties of galaxies with redshifts of z >~ 12, when the instrumental limiting magnitude is expected to be brighter than the characteristic magnitude of the galaxy sample. The number counts could be modified by an order of magnitude, with most galaxies being part of multiply imaged systems, located less than 1 arcsec from brighter foreground galaxies at z ~ 2. This lens-induced association of high-redshift and foreground galaxies has perhaps already been observed among a sample of galaxy candidates identified at z ~ 10.6. Future surveys will need to be designed to account for a significant gravitational lensing bias in high-redshift galaxy samples.

We have determined 194 Faraday rotation measures (RMs) of polarized extragalactic radio sources using new, multi-channel polarization observations at frequencies around 1.4 GHz from the Very Large Array (VLA) in the Galactic plane at $17^\circ \leq l \leq 63^\circ$ and $205^\circ \leq l \leq 253^\circ$. This catalog fills in gaps in the RM coverage of the Galactic plane between the Canadian Galactic Plane Survey and Southern Galactic Plane Survey. Using this catalog we have tested the validity of recently-proposed axisymmetric and bisymmetric models of the large-scale (or regular) Galactic magnetic field, and found that of the existing models we tested, an axisymmetric spiral model with reversals occurring in rings (as opposed to along spiral arms) best matched our observations. Building on this, we have performed our own modeling, using RMs from both extragalactic sources and pulsars. By developing independent models for the magnetic field in the outer and inner Galaxy, we conclude that in the inner Galaxy, the magnetic field closely follows the spiral arms, while in the outer Galaxy, the field is consistent with being purely azimuthal.Furthermore, the models contain no reversals in the outer Galaxy, and together seem to suggest the existence of a single reversed region that spirals out from the Galactic center.

A standard binary microlensing event lightcurve allows just two parameters of the lensing system to be measured: the mass ratio of the companion to its host, and the projected separation of the components in units of the Einstein radius. However, other exotic effects can provide more information about the lensing system. Orbital motion in the lens is one such effect, which if detected, can be used to constrain the physical properties of the lens. To determine the fraction of binary lens lightcurves affected by orbital motion (the detection efficiency) we simulate lightcurves of orbiting binary star and star-planet (planetary) lenses and simulate the continuous, high-cadence photometric monitoring that will be conducted by the next generation of microlensing surveys that are beginning to enter operation. The effect of orbital motion is measured by fitting simulated lightcurve data with standard static binary microlensing models; lightcurves that are poorly fit by these models are considered to be detections of orbital motion. We correct for systematic false positive detections by also fitting the lightcurves of static binary lenses. For a continuous monitoring survey without intensive follow-up of high magnification events, we find the orbital motion detection efficiency for planetary events with caustic crossings to be 0.061+-0.010, consistent with observational results, and 0.0130+-0.0055 for events without caustic crossings (smooth events). Similarly for stellar binaries, the orbital motion detection efficiency is 0.098+-0.011 for events with caustic crossings and is 0.048+-0.006 for smooth events. These result in combined (caustic crossing and smooth) orbital motion detection efficiencies of 0.029+-0.005 for planetary lenses and 0.070+-0.006 for stellar binary lenses. We also investigate how various microlensing parameters affect the orbital motion detectability. [Abridged]

Lensing flux-ratio anomalies have been frequently observed and taken as evidence for the presence of abundant dark matter substructures in lensing galaxies, as predicted by the cold dark matter (CDM) model of cosmogony. In previous work, we examined the cusp-caustic relations of the multiple images of background quasars lensed by galaxy-scale dark matter haloes, using a suite of high-resolution N-body simulations (the Aquarius simulations). In this work, we extend our previous calculations to incorporate both the baryonic and diffuse dark components in lensing haloes. We include in each lensing simulation: (1) a satellite galaxy population derived from a semi-analytic model applied to the Aquarius haloes, (2) an empirical Milky-Way globular cluster population and (3) satellite streams (diffuse dark component) identified in the simulations. Accounting for these extra components, we confirm our earlier conclusion that the abundance of intrinsic substructures (dark or bright, bound or diffuse) is not sufficient to explain the observed frequency of cusp-caustic violations in the CLASS survey. We conclude that the observed effect could be the result of the small number statistics of CLASS, or intergalactic haloes along the line of sight acting as additional sources of lensing flux anomalies. Another possibility is that this discrepancy signals a failure of the CDM model.

We present a study of the vertical magnetic field of the Milky Way towards the Galactic poles, determined from observations of Faraday rotation toward more than 1000 polarized extragalactic radio sources at Galactic latitudes |b| > 77 degs, using the Westerbork Radio Synthesis Telescope and the Australia Telescope Compact Array. We find median rotation measures (RMs) of 0.0 +/- 0.5 rad/m^2 and +6.3 +/- 0.7 rad/m^2 toward the north and south Galactic poles, respectively, demonstrating that there is no coherent vertical magnetic field in the Milky Way at the Sun's position. If this is a global property of the Milky Way's magnetism, then the lack of symmetry across the disk rules out pure dipole or quadrupole geometries for the Galactic magnetic field. The angular fluctuations in RM seen in our data show no preferred scale within the range ~ 0.1 to 25 degs. The observed standard deviation in RM of ~ 9 rad/m^2 then implies an upper limit of ~1microGauss on the strength of the random magnetic field in the warm ionized medium at high Galactic latitudes.

We re-derive the made-to-measure method of Syer and Tremaine 1996, for modelling stellar systems and individual galaxies, and demonstrate how extensions to the made-to-measure method may be implemented and used. We illustrate the enhanced made-to-measure method by determining the mass-to-light ratio of a galaxy modelled as a Plummer sphere. From the standard galactic observables of surface brightness and line-of-sight velocity dispersion together with the h_4 Gauss-Hermite coefficient of the line-of-sight velocity distribution, we successfully recover the true mass-to-light ratio of our toy galaxy. Using kinematic data from Kleyna et al 2002, we then estimate the mass-to-light ratio of the dwarf spheroidal galaxy Draco achieving a V-band value of 539 \pm 136 M_{\odot} / L_{\odot}. We describe the main aspects of creating a made-to-measure galaxy model and show how the key modelling parameters may be determined.

There is a remarkable synergy between requirements for Dark Energy probes by cosmic shear measurements and planet hunting by microlensing. Employing weak and strong gravitational lensing to trace and detect the distribution of matter on cosmic and Galactic scales, but as well as to the very small scales of exoplanets is a unique meeting point from cosmology to exoplanets. It will use gravity as the tool to explore the full range of masses not accessible by any other means. EUCLID is a 1.2m telescope with optical and IR wide field imagers and slitless spectroscopy, proposed to ESA Cosmic Vision to probe for Dark Energy, Baryonic acoustic oscillation, galaxy evolution, and an exoplanet hunt via microlensing. A 3 months microlensing program will already efficiently probe for planets down to the mass of Mars at the snow line, for free floating terrestrial or gaseous planets and habitable super Earth. A 12+ months survey would give a census on habitable Earth planets around solar like stars. This is the perfect complement to the statistics that will be provided by the KEPLER satellite, and these missions combined will provide a full census of extrasolar planets from hot, warm, habitable, frozen to free floating.

The incidence of sub-galactic level substructures is an important quantity, as it is a generic prediction of high-resolution Cold Dark Matter (CDM) models which is susceptible to observational test. Confrontation of theory with observations is currently in an uncertain state. In particular, gravitational lens systems appear to show evidence for flux ratio anomalies, which are expected from CDM substructures although not necessarily in the same range of radius as observed. However, the current small samples of lenses suggest that the lens galaxies in these systems are unusually often accompanied by luminous galaxies. Here we investigate a large sample of unlensed elliptical galaxies from the COSMOS survey, and determine the fraction of objects with satellites, in excess of background counts, as a function of satellite brightness and separation from the primary object. We find that the incidence of luminous satellites within 20 kpc is typically a few tenths of one percent for satellites of a few tenths of the primary flux, comparable to what is observed for the wider but shallower SDSS survey. Although the environments of lenses in the SLACS survey are compatible with this observation, the CLASS radio survey lenses are significantly in excess of this.

We use a high-resolution $N$-body simulation to investigate the influence of background galaxy properties, including redshift, size, shape and clustering, on the efficiency of forming giant arcs by gravitational lensing of rich galaxy clusters. Two large sets of ray-tracing simulations are carried out for 10 massive clusters at two redshifts, i.e. $z_{\rm l} \sim 0.2$ and 0.3. The virial mass ($M_{\rm vir}$) of the simulated lens clusters at $z\sim0.2$ ranges from $6.8\times10^{14} h^{-1} {M_{\odot}}$ to $1.1\times 10^{15} h^{-1} M_{\odot}$. The information of background galaxies brighter than 25 magnitude in the $I$-band is taken from Cosmological Evolution Survey (COSMOS) imaging data. Around $1.7\times 10^5$ strong lensing realizations with these images as background galaxies have been performed for each set. We find that the efficiency for forming giant arcs for $z_{\rm l}=0.2$ clusters is broadly consistent with observations. The efficiency of producing giant arcs by rich clusters is weakly dependent on the source size and clustering. Our principal finding is that a small proportion ($\sim 1/3$) of galaxies with elongated shapes (e.g. ellipticity $\epsilon=1-b/a>0.5$) can boost the number of giant arcs substantially. Compared with recent studies where a uniform ellipticity distribution from 0 to 0.5 is used for the sources, the adoption of directly observed shape distribution increases the number of giant arcs by a factor of $\sim2$. Our results indicate that it is necessary to account for source information and survey parameters (such as point-spread-function, seeing) to make correct predictions of giant arcs and further to constrain the cosmological parameters.(abridged)

We use galaxy catalogues constructed by combining high-resolution N-body simulations with semi-analytic models of galaxy formation to study the properties of Long Gamma-Ray Burst (LGRB) host galaxies. We assume that LGRBs originate from the death of massive young stars and analyse how results are affected by different metallicity constraints on the progenitor stars. As expected, the host sample with no metallicity restriction on the progenitor stars provides a perfect tracer of the cosmic star formation history. When LGRBs are required to be generated by low-metallicity stars, they trace a decreasing fraction of the cosmic star formation rate at lower redshift, as a consequence of the global increase in metallicity. We study the properties of host galaxies up to high redshift (~9), finding that they typically have low-metallicity (Z<0.5 Z_sun) and that they are small (M<10^9 M_sun), bluer and younger than the average galaxy population, in agreement with observational data. They are also less clustered than typical L_* galaxies in the Universe, and their descendents are massive, red and reside in groups of galaxies with halo mass between 10^{13} M_sun to 10^{14} M_sun.

We use high-resolution Aquarius simulations of Milky Way-sized haloes in the LCDM cosmology to study the effects of dark matter substructures on gravitational lensing. Each halo is resolved with ~ 10^8 particles (at a mass resolution ~ 10^3-4 M_sun/h) within its virial radius. Subhaloes with masses larger than 10^5 M_sun/h are well resolved, an improvement of at least two orders of magnitude over previous lensing studies. We incorporate a baryonic component modelled as a Hernquist profile and account for the response of the dark matter via adiabatic contraction. We focus on the "anomalous" flux ratio problem, in particular on the violation of the cusp-caustic relation due to substructures. We find that subhaloes with masses less than ~ 10^8 M_sun/h play an important role in causing flux anomalies; such low mass subhaloes have been unresolved in previous studies. There is large scatter in the predicted flux ratios between different haloes and between different projections of the same halo. In some cases, the frequency of predicted anomalous flux ratios is comparable to that observed for the radio lenses, although in most cases it is not. The probability for the simulations to reproduce the observed violations of the cusp lenses is about 0.001. We therefore conclude that the amount of substructure in the central regions of the Aquarius haloes is insufficient to explain the observed frequency of violations of the cusp-caustic relation. These conclusions are based purely on our dark matter simulations which ignore the effect of baryons on subhalo survivability.

We study the clustering properties of barred galaxies using data from the Sloan Digital Sky Survey (SDSS). We compute projected redshift-space two-point cross-correlation functions w_p(r_p) for a sample of nearly 1000 galaxies for which we have performed detailed structural decompositions using the methods described in Gadotti (2009). The sample includes 286 barred galaxies. The clustering of barred and unbarred galaxies of similar stellar mass is indistinguishable over all the scales probed (~20kpc - 30Mpc). This result also holds even if the sample is restricted to bars with bluer g-i colours (and hence younger ages). Our result also does not change if we split our sample of barred galaxies according to bar-to-total luminosity ratio, bar boxyness, effective surface brightness, length, or the shape of the surface density profile within the bar. There is a hint that red, elliptical bars are more strongly clustered than red and less elliptical bars, on scales >~ 1 Mpc, although the statistical significance is not high. We conclude that there is no significant evidence that bars are a product of mergers or interactions. We tentatively interpret the stronger clustering of the more elliptical bars as evidence that they are located in older galaxies, which reside in more massive haloes.

We search for ongoing major dry-mergers in a well selected sample of local Brightest Cluster Galaxies (BCGs) from the C4 cluster catalogue. 18 out of 515 early-type BCGs with redshift between 0.03 and 0.12 are found to be in major dry-mergers, which are selected as pairs (or triples) with $r$-band magnitude difference $\dm<1.5$ and projected separation $\rp<30$ kpc, and showing signatures of interaction in the form of significant asymmetry in residual images. We find that the fraction of BCGs in major dry-mergers increases with the richness of the clusters, consistent with the fact that richer clusters usually have more massive (or luminous) BCGs. We estimate that present-day early-type BCGs may have experienced on average $\sim 0.6 (\tmerge/0.3\Gyr)^{-1}$ major dry-mergers and through this process increases their luminosity (mass) by $15% (\tmerge/0.3\Gyr)^{-1} (\fmass/0.5)$ on average since $z=0.7$, where $\tmerge$ is the merging timescale and $\fmass$ is the mean mass fraction of companion galaxies added to the central ones. We also find that major dry-mergers do not seem to elevate radio activities in BCGs. Our study shows that major dry-mergers involving BCGs in clusters of galaxies are not rare in the local Universe, and they are an important channel for the formation and evolution of BCGs.

Based on the Sloan Digital Sky Survey DR6 (SDSS) and Millennium Simulation (MS) we investigate the alignment between galaxies and large-scale structure. For this purpose we develop two new statistical tools, namely the alignment correlation function and the cos(2theta)-statistic. The former is a two-dimensional extension of the traditional two-point correlation function and the latter is related to the ellipticity correlation function used for cosmic shear measurements. Both are based on the cross correlation between a sample of galaxies with orientations and a reference sample which represents the large-scale structure. Applied to the SDSS galaxy catalog the alignment correlation function reveals an overabundance of reference galaxies along the major axes of red, luminous (L > L*) galaxies out to projected separations of 60 Mpc/h. No alignment signal is detected for blue galaxies. The cos(2\theta)-statistic yields very similar results. Starting from a MS semi-analytic galaxy catalog we assign an orientation to each red, luminous and central galaxy, based on the central region of the host halo. Alternatively, we use the orientation of the host halo itself. We find a mean projected misalignment between a halo and its central region of ~25 deg. Agreement with the SDSS results is good if the central orientations are used. Using the halo orientations overestimates the observed alignment by more than a factor of 2. The large volume of the MS allows to generate two-dimensional maps of the alignment correlation function which show the reference galaxy distribution to be flattened parallel to the orientations of red luminous galaxies with axis ratios of ~0.5 and ~0.75 for halo and central orientations,respectively. These ratios are almost independent of scale out to 60 Mpc/h.

Microlensing events are usually selected among single-peaked non-repeating light curves in order to avoid confusion with variable stars. However, a microlensing event may exhibit a second microlensing brightening episode when the source or/and the lens is a binary system. A careful analysis of these repeating events provides an independent way to study the statistics of wide binary stars and to detect extrasolar planets. Previous theoretical studies predicted that 0.5 - 2 % of events should repeat due to wide binary lenses. We present a systematic search for such events in about 4000 light curves of microlensing candidates detected by the Optical Gravitational Lensing Experiment (OGLE) towards the Galactic Bulge from 1992 to 2007. The search reveals a total of 19 repeating candidates, with 6 clearly due to a wide binary lens. As a by-product we find that 64 events (~2% of the total OGLE-III sample) have been miss-classified as microlensing; these miss-classified events are mostly nova or other types of eruptive stars. The number and importance of repeating events will increase considerably when the next-generation wide-field microlensing experiments become fully operational in the future.

The basic concepts of gravitational microlensing are introduced. We start with the lens equation, and then derive the image positions and magnifications. The statistical quantities of optical depth and event rate are then described. We finish with a summary and a list of challenges and open questions. A problem set is given for students to practice.

Substructures, expected in cold dark matter haloes, have been proposed to explain the anomalous flux ratios in gravitational lenses. About 25% of lenses in the Cosmic Lens All-Sky Survey (CLASS) appear to have luminous satellites within ~ 5 kpc/h of the main lensing galaxies, which are usually at redshift z ~ 0.2-1. In this work we use the Millennium Simulation combined with galaxy catalogues from semi-analytical techniques to study the predicted frequency of such satellites in simulated haloes. The fraction of haloes that host bright satellites within the (projected) central regions is similar for red and blue hosts and is found to increase as a function of host halo mass and redshift. Specifically, at z = 1, about 11% of galaxy-sized haloes (with masses between 10^{12} M_sun/h and 10^{13} M_sun/h) host bright satellite galaxies within a projected radius of 5 kpc/h. This fraction increases to about 17% (25%) if we consider bright (all) satellites of only group-sized haloes (with masses between 10^{13} M_sun/h and 10^{14} M_sun/h). These results are roughly consistent with the fraction (~ 25%) of CLASS lensing galaxies observed to host luminous satellites. At z = 0, only ~ 3% of galaxy-sized haloes host bright satellite galaxies. The fraction rises to ~ 6%, (10%) if we consider bright (all) satellites of only group-sized haloes at z = 0. However, most of the satellites found in the inner regions are `orphan' galaxies where the dark matter haloes have been completely stripped. Thus the agreement crucially depends on the true survival rate of these `orphan' galaxies. We also discuss the effects of numerical resolution and cosmologies on our results.

We present a new joint analysis of pulsar dispersion measures and diffuse H-alpha emission in the Milky Way, which we use to derive the density, pressure and filling factor of the thick disk component of the warm ionised medium (WIM) as a function of height above the Galactic disk. By excluding sightlines at low Galactic latitude that are contaminated by HII regions and spiral arms, we find that the exponential scale-height of free electrons in the diffuse WIM is 1830 (+120, -250) pc, a factor of two larger than has been derived in previous studies. The corresponding inconsistent scale heights for dispersion measure and emission measure imply that the vertical profiles of mass and pressure in the WIM are decoupled, and that the filling factor of WIM clouds is a geometric response to the competing environmental influences of thermal and non-thermal processes. Extrapolating the properties of the thick-disk WIM to mid-plane, we infer a volume-averaged electron density 0.014 +- 0.001 cm^-3, produced by clouds of typical electron density 0.34 +- 0.06 cm^-3 with a volume filling factor 0.04 +- 0.01. As one moves off the plane, the filling factor increases to a maximum of ~30% at a height of approximately 1-1.5 kpc, before then declining to accommodate the increasing presence of hot, coronal gas. Since models for the WIM with a ~1 kpc scale-height have been widely used to estimate distances to radio pulsars, our revised parameters suggest that the distances to many high-latitude pulsars have been substantially underestimated.

Thirteen exo-planets have been discovered using the gravitational microlensing technique (out of which 7 have been published). These planets already demonstrate that super-Earths (with mass up to ~10 Earth masses) beyond the snow line are common and multiple planet systems are not rare. In this White Paper we introduce the basic concepts of the gravitational microlensing technique, summarise the current mode of discovery and outline future steps towards a complete census of planets including Earth-mass planets. In the near-term (over the next 5 years) we advocate a strategy of automated follow-up with existing and upgraded telescopes which will significantly increase the current planet detection efficiency. In the medium 5-10 year term, we envision an international network of wide-field 2m class telescopes to discover Earth-mass and free-floating exo-planets. In the long (10-15 year) term, we strongly advocate a space microlensing telescope which, when combined with Kepler, will provide a complete census of planets down to Earth mass at almost all separations. Such a survey could be undertaken as a science programme on Euclid, a dark energy probe with a wide-field imager which has been proposed to ESA's Cosmic Vision Programme.

We analyse mid-infrared (MIR) spectroscopic properties for 19 ultra-luminous infrared quasars (IR QSOs) in the local universe based on the spectra from the Infrared Spectrograph on board the Spitzer Space Telescope. The MIR properties of IR QSOs are compared with those of optically-selected Palomar-Green QSOs (PG QSOs) and ultra-luminous infrared galaxies (ULIRGs). The average MIR spectral features from ~ 5 to 30um, including the spectral slopes, 6.2um PAH emission strengths and [NeII] 12.81um luminosities of IR QSOs, differ from those of PG QSOs. In contrast, IR QSOs and ULIRGs have comparable PAH and [NeII] luminosities. These results are consistent with IR QSOs being at a transitional stage from ULIRGs to classical QSOs. We also find that the colour index alpha(30, 15) is a good indicator of the relative contribution of starbursts to AGNs for all QSOs. Correlations between the [NeII] 12.81um and PAH 6.2um luminosities and those between the [NeII], PAH with 60um luminosities for ULIRGs and IR QSOs indicate that both [NeII] and PAH luminosities are approximate star formation rate indicators for IR QSOs and starburst-dominated galaxies; the scatters are, however, quite large (~ 0.7 to 0.8 dex). Finally the correlation between the EW(PAH 6.2um) and outflow velocities suggests that star formation activities are suppressed by feedback from AGNs and/or supernovae.

We present a study of the magnetic field of the Small Magellanic Cloud (SMC), carried out using radio Faraday rotation and optical starlight polarization data. Consistent negative rotation measures (RMs) across the SMC indicate that the line-of-sight magnetic field is directed uniformly away from us with a strength 0.19 +/- 0.06 microGauss. Applying the Chandrasekhar-Fermi method to starlight polarization data yields an ordered magnetic field in the plane of the sky of strength 1.6 +/- 0.4 microGauss oriented at a position angle 4 +/- 12 degs, measured counter-clockwise from the great circle on the sky joining the SMC to the Large Magellanic Cloud (LMC). We construct a three-dimensional magnetic field model of the SMC, under the assumption that the RMs and starlight polarization probe the same underlying large-scale field. The vector defining the overall orientation of the SMC magnetic field shows a potential alignment with the vector joining the center of the SMC to the center of the LMC, suggesting the possibility of a "pan-Magellanic'' magnetic field. A cosmic-ray driven dynamo is the most viable explanation of the observed field geometry, but has difficulties accounting for the observed uni-directional field lines. A study of Faraday rotation through the Magellanic Bridge is needed to further test the pan-Magellanic field hypothesis.

The detection of fast quasi-periodic variability from accreting black holes and neutron stars has been used to constrain their masses, radii, and spins. If the observed oscillations are linear modes in the accretion disks, then bounds can be placed on the properties of the central objects by assuming that these modes are locally sub-Keplerian. If, on the other hand, the observed oscillations correspond to non-linear resonances between disk modes, then the properties of the central objects can be measured by assuming that the resonant modes are excited at the same radial annulus in the disk. In this paper, we use numerical simulations of vertically integrated, axisymmetric hydrodynamic accretion disks to provide examples of situations in which the assumptions implicit in both methods are not satisfied. We then discuss our results for the robustness of the mass and spin measurements of compact objects from variability studies.

Based on a cosmological N-body simulation we analyze spatial and kinematic alignments of satellite halos within six times the virial radius of group size host halos (Rvir). We measure three different types of spatial alignment: halo alignment between the orientation of the group central substructure (GCS) and the distribution of its satellites, radial alignment between the orientation of a satellite and the direction towards its GCS, and direct alignment between the orientation of the GCS and that of its satellites. In analogy we use the directions of satellite velocities and probe three further types of alignment: the radial velocity alignment between the satellite velocity and connecting line between satellite and GCS, the halo velocity alignment between the orientation of the GCS and satellite velocities and the auto velocity alignment between the satellites orientations and their velocities. We find that satellites are preferentially located along the major axis of the GCS within at least 6 Rvir (the range probed here). Furthermore, satellites preferentially point towards the GCS. The most pronounced signal is detected on small scales but a detectable signal extends out to 6 Rvir. The direct alignment signal is weaker, however a systematic trend is visible at distances < 2 Rvir. All velocity alignments are highly significant on small scales. Our results suggest that the halo alignment reflects the filamentary large scale structure which extends far beyond the virial radii of the groups. In contrast, the main contribution to the radial alignment arises from the adjustment of the satellite orientations in the group tidal field. The projected data reveal good agreement with recent results derived from large galaxy surveys. (abridged)

The photometry data base of the second phase of the OGLE microlensing experiment, OGLE-II, is a rich source of information about the kinematics and structure of the Galaxy. In this work we use the OGLE-II proper motion catalogue to identify candidate stars which have high proper motions. 521 stars with proper motion mu > 50 mas/yr in the OGLE-II proper motion catalogue (Sumi 2004) were cross-identified with stars in the MACHO high proper motion catalogue, and the DENIS and 2MASS infra-red photometry catalogues. Photometric distances were computed for stars with colours consistent with G/K and M type stars. 6 stars were newly identified as possible nearby (< 50 pc) M dwarfs.

We investigate the photometric properties of the early type Brightest Cluster Galaxies (BCGs) using a carefully selected sample of 85 BCGs from the C4 cluster catalogue with redshift less than 0.1. We perform accurate background subtractions, and surface photometry for these BCGs to 25 $\magsec$ in the Sloan $r$-band. By quantitatively analysing the gradient of the Petrosian profiles of BCGs, we find that a large fraction of BCGs have extended stellar envelopes in their outskirts; more luminous BCGs tend to have more extended stellar halos that are likely connected with mergers. A comparison sample of elliptical galaxies was chosen with similar apparent magnitude and redshift ranges, for which the same photometric analysis procedure is applied. We find that BCGs have steeper size-luminosity ($R \propto L^\alpha$) and Faber-Jackson ($L \propto \sigma^\beta$) relations than the bulk of early type galaxies. Furthermore, the power-law indices ($\alpha$ and $\beta$) in these relations increase as the isophotal limits become deeper. For isophotal limits from 22 to 25 $\magsec$, BCGs are usually larger than the bulk of early type galaxies, and a large fraction ($\sim 49%$) of BCGs have disky isophotal shapes. The differences in the scaling relations are consistent with a scenario where the dynamical structure and formation route of BCGs may be different from the bulk of early type galaxies, in particular dry (dissipationless) mergers may play a more important role in their formation; we highlight several possible dry merger candidates in our sample.

We present a direct detection of the gravitational lens that caused the microlensing event MACHO-95-BLG-37. This is the first fully resolved microlensing system involving a source in the Galactic bulge, and the second such system in general. The lens and source are clearly resolved in images taken with the High Resolution Channel of the Advanced Camera for Surveys on board the Hubble Space Telescope (HST) ~9 years after the microlensing event. The presently available data are not sufficient for the final, unambiguous identification of the gravitational lens and the microlensed source. While the light curve models combined with the high resolution photometry for individual objects indicate that the source is red and the lens is blue, the color-magnitude diagram for the line of sight and the observed proper motions strongly support the opposite case. The first scenario points to a metal-poor lens with mass M = ~0.6 M_Sun at the distance D_l = ~4 kpc. In the second scenario the lens could be a main-sequence star with M = 0.8 - 0.9 M_Sun about half-way to the Galactic bulge or in the foreground disk, depending on the extinction.

Gravitational microlensing surveys target very dense stellar fields in the local group. As a consequence the microlensed source stars are often blended with nearby unresolved stars. The presence of `blending' is a cause of major uncertainty when determining the lensing properties of events towards the Galactic centre. After demonstrating empirical cases of blending we utilize Monte Carlo simulations to probe the effects of blending. We generate artificial microlensing events using an HST luminosity function convolved to typical ground-based seeing, adopting a range of values for the stellar density and seeing. We find that a significant fraction of bright events are blended, contrary to the oft-quoted assumption that bright events should be free from blending. We probe the effect that this erroneous assumption has on both the observed event timescale distribution and the optical depth, using realistic detection criteria relevent to the different surveys. Importantly, under this assumption the latter quantity appears to be reasonably unaffected across our adopted values for seeing and density. The timescale distribution is however biased towards smaller values, even for the least dense fields. The dominant source of blending is from lensing of faint source stars, rather than lensing of bright source stars blended with nearby fainter stars. We also explore other issues, such as the centroid motion of blended events and the phenomena of `negative' blending. Furthermore, we breifly note that blending can affect the determination of the centre of the red clump giant region from an observed luminosity function. This has implications for a variety of studies, e.g. mapping extinction towards the bulge and attempts to constrain the parameters of the Galactic bar through red clump giant number counts. (Abridged)

We use high-resolution N-body numerical simulations to study the number of predicted large-separation multiply-imaged systems produced by clusters of galaxies in the SDSS photometric and spectroscopic quasar samples. We incorporate the condensation of baryons at the centre of clusters by (artificially) adding a brightest central galaxy (BCG) as a truncated isothermal sphere. We make predictions in two flat cosmological models: a LCDM model with a matter density $\Omega_0=0.3$, and $\sigma_8=0.9$ (LCDM0), and a model favoured by the WMAP three-year data with $\Omega_0=0.238$, and $\sigma_8=0.74$ (WMAP3). We found that the predicted multiply-imaged quasars with separation >10" is about 6.2 and 2.6 for the SDSS photometric (with an effective area 8000 deg$^2$) and spectroscopic (with an effective area 5000 deg$^2$) quasar samples respectively in the LCDM0 model; the predicted numbers of large-separation lensed quasars agree well with the observations. These numbers are reduced by a factor of 7 or more in the WMAP3 model, and are consistent with data at < 8% level. The predicted cluster lens redshift peaks around redshift 0.5, and 90% are between 0.3 and 1. We find that the BCG creates a central circular region, comparable to the Einstein ring of the BCG, where the central image disappears in the usual three-image and five-image configurations. If we include four image systems as an extreme case of five-image systems (with an infinitely demagnified central image), we find that 68% of the central images are fainter by a factor of 100 than the brightest image, and about 80% are within 1.5"of the BCG.

(abridged) The discovery of OGLE 2005-BLG-390Lb, the first cool rocky/icy exoplanet, impressively demonstrated the sensitivity of the microlensing technique to extra-solar planets below 10 M_earth. A planet of 1 M_earth in the same spot would have provided a detectable deviation with an amplitude of ~ 3 % and a duration of ~ 12 h. An early detection of a deviation could trigger higher-cadence sampling which would have allowed the discovery of an Earth-mass planet in this case. Here, we describe the implementation of an automated anomaly detector, embedded into the eSTAR system, that profits from immediate feedback provided by the robotic telescopes that form the RoboNet-1.0 network. It went into operation for the 2007 microlensing observing season. As part of our discussion about an optimal strategy for planet detection, we shed some new light on whether concentrating on highly-magnified events is promising and planets in the 'resonant' angular separation equal to the angular Einstein radius are revealed most easily. Given that sub-Neptune mass planets can be considered being common around the host stars probed by microlensing (preferentially M- and K-dwarfs), the higher number of events that can be monitored with a network of 2m telescopes and the increased detection efficiency for planets below 5 M_earth arising from an optimized strategy gives a common effort of current microlensing campaigns a fair chance to detect an Earth-mass planet (from the ground) ahead of the COROT or Kepler missions. The detection limit of gravitational microlensing extends even below 0.1 M_earth, but such planets are not very likely to be detected from current campaigns. However, these will be within the reach of high-cadence monitoring with a network of wide-field telescopes or a space-based telescope.

Using a large galaxy group catalogue based on the Sloan Digital Sky Survey Data Release 4 we measure three different types of intrinsic galaxy alignment within groups: halo alignment between the orientation of the brightest group galaxies (BGG) and the distribution of its satellite galaxies, radial alignment between the orientation of a satellite galaxy and the direction towards its BGG, and direct alignment between the orientation of the BGG and that of its satellites. In agreement with previous studies we find that satellite galaxies are preferentially located along the major axis. In addition, on scales r < 0.7 Rvir we find that red satellites are preferentially aligned radially with the direction to the BGG. The orientations of blue satellites, however, are perfectly consistent with being isotropic. Finally, on scales r < 0.1 \Rvir, we find a weak but significant indication for direct alignment between satellites and BGGs. We briefly discuss the implications for weak lensing measurements.

We study the growth of black holes and stellar population in spheroids at high redshift using several (sub)mm-loud QSO samples. Applying the same criteria established in an earlier work, we find that, similar to IR QSOs at low redshift, the far-infrared emission of these (sub)mm-loud QSOs mainly originates from dust heated by starbursts. By combining low-z IR QSOs and high-z (sub)mm-loud QSOs, we find a trend that the star formation rate ($\Mstardot$) increases with the accretion rate ($\Mdot$). We compare the values of $\Mstardot/\Mdot$ for submm emitting galaxies (SMGs), far-infrared ultraluminous/hyperluminous QSOs and typical QSOs, and construct a likely evolution scenario for these objects. The (sub)mm-loud QSO transition phase has both high $\Mdot$ and $\Mstardot$ and hence is important for establishing the correlation between the masses of black holes and spheroids.

Recent studies have shown that the distribution of satellite galaxies is preferentially aligned with the major axis of their central galaxy. The strength of this alignment has been found to depend strongly on the colours of the satellite and central galaxies, and only weakly on the mass of the halo in which the galaxies reside. In this paper we study whether these alignment signals, and their dependence on galaxy and halo properties, can be reproduced in a hierarchical structure formation model of a $\Lambda$CDM concordance cosmology. To that extent we use a large $N$-body simulation which we populate with galaxies following a semi-analytical model for galaxy formation. We find that if the orientation of the central galaxy is perfectly aligned with that of its dark matter halo, then the predicted central-satellite alignment signal is much stronger than observed. If, however, the minor axis of a central galaxy is perfectly aligned with the angular momentum vector of its dark matter halo, we can accurately reproduce the observed alignment strength as function of halo mass and galaxy color. Although this suggests that the orientation of central galaxies is governed by the angular momentum of their dark matter haloes, we emphasize that any other scenario in which the minor axes of central galaxy and halo are misaligned by $\sim 40^{\circ}$ (on average) will match the data equally well. Finally, we show that dependence of the alignment strength on the color of the central galaxy is most likely an artefact due to interlopers in the group catalogue. The dependence on the color of the satellite galaxies, on the other hand, is real and owes to the fact that red satellites are associated with subhaloes that were more massive at their time of accretion.

Red clump giants in the Galactic bulge are approximate standard candles and hence they can be used as distance indicators. We compute the proper motion dispersions of RCG stars in the Galactic bulge using the proper motion catalogue from the second phase of the Optical Gravitational Microlensing Experiment (OGLE-II, Sumi et al. 2004) for 45 fields. The proper motion dispersions are measured to a few per cent accuracy due to the large number of stars in the fields. The observational sample is comprised of 577736 stars. These observed data are compared to a state-of-the-art particle simulation of the Galactic bulge region. The predictions are in rough agreement with observations, but appear to be too anisotropic in the velocity ellipsoid. We note that there is significant field-to-field variation in the observed proper motion dispersions. This could either be a real feature, or due to some unknown systematic effect.

Red clump giant stars can be used as distance indicators to trace the mass distribution of the Galactic bar. We use RCG stars from 44 bulge fields from the OGLE-II microlensing collaboration database to constrain analytic tri-axial models for the Galactic bar. We find the bar major axis is oriented at an angle of 24 - 27 degrees to the Sun-Galactic centre line-of-sight. The ratio of semi-major and semi-minor bar axis scale lengths in the Galactic plane x_0, y_0, and vertical bar scale length z_0, is x_0 : y_0 : z_0 = 10 : 3.5 : 2.6, suggesting a slightly more prolate bar structure than the working model of Gerhard (2002) which gives the scale length ratios as x_0 : y_0 : z_0 = 10 : 4 : 3 .

Semi-analytical models of galaxy formation can be used to predict the evolution of the number density of early-type galaxies as a function of the circular velocity at the virial radius, v_{c,vir}. Gravitational lensing probability and separation distribution on the other hand are sensitive to the velocity dispersion (or circular velocity) at about the effective radius. We adopt the singular isothermal ellipsoid (SIE) lens model to estimate the velocity dispersion at the effective radius. We use radio lenses from the Cosmic Lens All-Sky Survey and the PMN-NVSS Extragalactic Lens Survey to study how the velocity dispersions, \sigma_SIE, are related to v_{c,vir}. When we include both the lensing probability and separation distribution as our lensing constraints, we find \sigma_SIE /(200\kms) = [(1.17_{-0.26}^{+0.40}) v_{c,vir}/ (200\kms)]^{0.22^{+0.05}_{-0.04}} for 200\kms \la \sigma_SIE \la 260\kms; at \sigma_SIE = 200\kms, the ratio \sqrt{2} \sigma_SIE / v_{c,vir} is about 1.65^{+0.57}_{-0.37} (68% CL) but decreases to 0.65_{-0.12}^{+0.15} (68% CL) for \sigma_SIE = 260\kms. These results are consistent with those of Seljak (2002) obtained from galaxy-galaxy weak lensing for galaxies of around L_*. However, our results clearly suggest that the ratio must vary significantly as \sigma_SIE is varied and are marginally discrepant with the Seljak results at \sigma_SIE = 260\kms. The scaling \sigma_SIE \propto v_{c,vir}^{0.22\pm 0.05} is broadly consistent with those from galaxy occupation statistics studies and the most recent galaxy-galaxy weak lensing study. (Abridged)

Using a set of high-resolution N-body/SPH cosmological simulations with identical initial conditions but run with different numerical setups, we investigate the influence of baryonic matter on the mass distribution of dark halos when radiative cooling is NOT included. We compare the concentration parameters of about 400 massive halos with virial mass from $10^{13}$ \Msun to $7.1 \times 10^{14}$ \Msun. We find that the concentration parameters for the total mass and dark matter distributions in non radiative simulations are on average larger by ~3% and 10% than those in a pure dark matter simulation. Our results indicate that the total mass density profile is little affected by a hot gas component in the simulations. After carefully excluding the effects of resolutions and spurious two-body heating between dark matter and gas particles, we conclude that the increase of the dark matter concentration parameters is due to interactions between baryons and dark matter. We demonstrate this with the aid of idealized simulations of two-body mergers. The results of individual halos simulated with different mass resolutions show that the gas profiles of densities, temperature and entropy are subjects of mass resolution of SPH particles. In particular, we find that in the inner parts of halos, as the SPH resolution increases the gas density becomes higher but both the entropy and temperature decrease.

We propose a new smoothing method for obtaining surface densities from discrete particle positions from numerical simulations. This is an essential step for many applications in gravitational lensing. This method is based on the ``scatter'' interpretation of the discrete density field in the Smoothed Particle Hydrodynamics. We use Monte Carlo simulations of uniform density fields and one isothermal ellipsoid to empirically derive the noise properties, and best smoothing parameters (such as the number of nearest neighbors used). A cluster from high-resolution simulations is then used to assess the reality of high-order singularities such as swallowtails and butterflies in caustics, which are important for the interpretation of substructures in gravitational lenses. We also compare our method with the Delaunay tesselation field estimator using the galaxy studied by Bradac et al. (2004), and find good agreements. We show that higher order singularities are not only connected with bound subhaloes but also with the satellite streams. However, the presence of high-order singularities are sensitive to not only the fluctuation amplitude of the surface density, but also the detailed form of the underlying smooth lensing potential (such as ellipticity and external shear).

We use high-resolution $N$-body simulations to investigate the optical depth of giant arcs with length-to-width ratio larger than 7.5 and 10 in the `standard' $\LCDM$ model with $\sigma_8=0.9$ and $\Omega_{\rm m,0}=0.3$ and a model based on three-year Wilkinson Microwave Anisotropy Probe (WMAP) data. We find that, in dark-matter only simulations, the lensing probability in the three-year WMAP model (with $\sigma_8=0.74$ and $\Omega_{\rm m,0}=0.238)$ decreases by a factor of $\sim 6$ compared with that in the `standard' $\LCDM$ model. The effects of baryonic cooling, star formation and feedbacks are uncertain, but we argue that baryons will only increase the the lensing cross-section by a moderate factor, $\sim 2$. We conclude that the low central value of $\sigma_8$ and $\Omega_{\rm m,0}$ preferred by the WMAP three-year data may be too low to be compatible with observations if conventional assumptions of the background source population are correct.

We present the first systematic search for microlensing events with variability in their baselines using data from the third phase of the Optical Gravitational Lensing Experiment (OGLE-III). A total of 137 candidates (88 new) was discovered toward the Galactic bulge. Among these, 21 have periodic oscillations in their baselines, 111 are irregular variables and 5 are potential long period detached eclipsing binaries. This is about 10% of the total number of constant baseline events. They are hence quite common and can be regarded as a new type of exotic events, which allow the determination of extra parameters of the events. We show that microlensing of variable stars may allow us to break the degeneracy between the blending parameter and magnification. We note that in some cases variability hidden in the baseline due to strong blending may be revealed in highly magnified events and resemble other exotic microlensing behavior, including planetary deviation. A new system (VEWS) for detecting ongoing variable baseline microlensing events is presented.

We revisit the shapes of isophotes for elliptical (E) and lenticular (S0) galaxies by studying 847 nearby early-type galaxies selected from the Sloan Digital Sky Survey Data Release 4 with velocity dispersions above 200 km/s. The IRAF task {\tt ellipse} was used to derive the deviations of the isophotes from pure ellipses (Fourier coefficients a_3/a and a_4/a), position angles and ellipticities as a function of radius. We show the statistical distributions of the a_3/a and a_4/a parameters as a function of velocity dispersion, magnitude, and colour. The a_4/a parameter is correlated with the ellipticity and absolute Petrosian magnitude of galaxies. No significant correlation was found between the a_4/a parameter with colour and velocity dispersion. A cross-correlation between the FIRST survey and the SDSS data reveals a lack of powerful radio emitters in discy E/S0s, as previously found by Bender et al. We also find that boxy E/S0s favor denser environments while discy E/S0s favor more isolated environments. The median values of changes in the ellipticity and position angle between one and one and a half Petrosian half-light radii in the isophotes are about -0.023 and 1.61 degree respectively. The average change in the position angle is much larger, about 4.12 degree, due to an extended tail. The change in ellipticity is weakly correlated with the ellipticity itself, with an increasing ellipticity for galaxies with higher ellipticity as the radius increases. The isophote parameters for the 847 galaxies are available online.

We study the morphology and star formation properties of 159 local luminous infrared galaxy (LIRG) using multi-color images from Data Release 2 (DR2) of the Sloan Digital Sky Survey (SDSS). The LIRGs are selected from a cross-correlation analysis between the IRAS survey and SDSS. They are all brighter than 15.9 mag in the r-band and below redshift ~ 0.1, and so can be reliably classified morphologically. We find that the fractions of interacting/merging and spiral galaxies are ~ 48% and ~ 40% respectively. Our results complement and confirm the decline (increase) in the fraction of spiral (interacting/merging) galaxies from z ~1 to z ~ 0.1, as found by Melbourne, Koo & Le Floc'h (2005). About 75% of spiral galaxies in the local LIRGs are barred, indicating that bars may play an important role in triggering star formation rates > 20 M_{sun}/yr in the local universe. Compared with high redshift LIRGs, local LIRGs have lower specific star formation rates, smaller cold gas fractions and a narrower range of stellar masses. Local LIRGs appear to be either merging galaxies forming intermediate mass ellipticals or spiral galaxies undergoing high star formation activities regulated by bars.

We present a proper motion mini-survey of 35 fields in the vicinity of Baade window, (l, b) = (1 deg, -4 deg), sampling roughly a 5 x 2.5 deg region of the Galactic bar. Our second epoch observations collected with the ACS/HRC instrument on board the Hubble Space Telescope were combined with the archival WFPC2/PC images. The resulting time baselines are in the range of 4 - 8 years. Precise proper motions of 15,863 stars were determined in the reference frame defined by the mean motion of stars with magnitudes between I_F814W = 16.5 - 21.5 along the line of sight. We clearly detect small gradients in proper motion dispersions (sigma_l, sigma_b) ~ (3.0, 2.5) mas/yr, and in the amount of anisotropy (sigma_l/sigma_b ~ 1.2). Both the longitude dispersion sigma_l and its ratio to the vertical motion sigma_b increase toward the Galactic plane. The decline of the anisotropy ratio sigma_l/sigma_b toward the minor axis of the bulge is mostly due to increasing sigma_b. We also find, for the first time, a significant negative covariance term in the transverse velocity field sigma_lb/(sigma_l*sigma_b) ~ -0.10. Our results extend by a factor of ~15 the number of the Galactic bar fields with good proper motion dispersions.

We use galaxy groups selected from the Sloan Digital Sky Survey to examine the alignment between the orientation of the central galaxy (defined as the brightest group member) and the distribution of satellite galaxies. By construction, we therefore only address the alignment on scales smaller than the halo virial radius. We find a highly significant alignment of satellites with the major axis of their central galaxy. This is in qualitative agreement with the recent study of Brainerd (2005), but inconsistent with several previous studies who detected a preferential minor axis alignment. The alignment strength in our sample is strongest between red central galaxies and red satellites. On the contrary, the satellite distribution in systems with a blue central galaxy is consistent with isotropic. We also find that the alignment strength is stronger in more massive haloes and at smaller projected radii from the central galaxy. In addition, there is a weak indication that fainter (relative to the central galaxy) satellites are more strongly aligned. We present a detailed comparison with previous studies, and discuss the implications of our findings for galaxy formation.