In an effort to probe the origin of surface brightness profile (SBP) breaks widely observed in nearby disk galaxies, we carry out a comparative study of stellar population profiles of 635 disk galaxies selected from the MaNGA spectroscopic survey. We classify our galaxies into single exponential (TI), down-bending (TII) and up-bending (TIII) SBP types, and derive their spin parameters and radial profiles of age/metallicity-sensitive spectral features. Most TII (TIII) galaxies have down-bending (up-bending) star formation rate (SFR) radial profiles, implying that abrupt radial changes of SFR intensities contribute to the formation of both TII and TIII breaks. Nevertheless, a comparison between our galaxies and simulations suggests that stellar migration plays a significant role in weakening down-bending $\Sigma_{\star}$ profile breaks. While there is a correlation between the break strengths of SBPs and age/metallicity-sensitive spectral features for TII galaxies, no such correlation is found for TIII galaxies, indicating that stellar migration may not play a major role in shaping TIII breaks, as is evidenced by a good correspondence between break strengths of $\Sigma_{\star}$ and surface brightness profiles of TIII galaxies. We do not find evidence for galaxy spin being a relevant parameter for forming different SBP types, nor do we find significant differences between the asymmetries of galaxies with different SBP types, suggesting that environmental disturbances or satellite accretion in the recent past do not significantly influence the break formation. By dividing our sample into early and late morphological types, we find that galaxies with different SBP types follow nearly the same tight stellar mass-$R_{25}$ relation, which makes the hypothesis that stellar migration alone can transform SBP types from TII to TI and then to TIII highly unlikely.

M51 and NGC 5195 is an interacting system that can be explored in great details with ground-based telescopes. The H II regions in M51 were observed using the 2.16 m telescope of the National Astronomical Observatories of the Chinese Academy of Sciences and the 6.5 m Multiple Mirror Telescope with spatial resolution of less than $\sim100$ pc. We obtain a total of 113 spectra across the galaxy and combine the literature data of Croxall et al. to derive a series of physical properties, including the gas-phase extinction, stellar population age, star formation rate (SFR) surface density, and oxygen abundance. The spatial distributions and radial profiles of these properties are investigated in order to study the characteristics of M51 and the clues to the formation and evolution of this galaxy. M51 presents a mild radial extinction gradient. The lower gas-phase extinction in the north spiral arms compared to the south arms are possibly caused by the past encounters with the companion galaxy of NGC 5195. A number of H II regions have the stellar age between 50 and 500 Myr, consistent with the recent interaction history by simulations in the literatures. The SFR surface density presents a mild radial gradient, which is ubiquitous in spiral galaxies. There is a negative metallicity gradient of $-0.08$ dex $R_{e}^{-1}$ in the disk region, which is also commonly found in many spiral galaxies. It is supported by the "inside-out" scenario of galaxy formation. We find a positive abundance gradient of 0.26 dex $R_{e}^{-1}$ in the inner region. There are possible reasons causing the positive gradient, including the freezing of the chemical enrichment due to the star-forming quenching in the bulge and the gas infall and dilution due to the pseudobulge growth and/or galactic interaction.

Using the star cluster catalogs from the Hubble Space Telescope program Legacy ExtraGalactic UV survey (LEGUS) and 8 $\mu$m images from the IRAC camera on the Spitzer Space Telescope for 5 galaxies within 5 Mpc, we investigate how the 8 $\mu$m dust luminosity correlates with the stellar age on the 30--50 pc scale of star forming regions. We construct a sample of 97 regions centered at local peaks of 8 $\mu$m emission, each containing one or more young star cluster candidates from the LEGUS catalogs. We find a tight anti-correlation with a Pearson correlation coefficient of $r=-0.84\pm0.05$ between the mass-normalized dust-only 8 $\mu$m luminosity and the age of stellar clusters younger than 1 Gyr; the 8 $\mu$m luminosity decreases with increasing age of the stellar population. Simple assumptions on a combination of stellar and dust emission models reproduce the observed trend. We also explore how the scatter of the observed trend depends on assumptions of stellar metallicity, PAH abundance, fraction of stellar light absorbed by dust, and instantaneous versus continuous star formation models. We find that variations in stellar metallicity have little effect on the scatter, while PAH abundance and the fraction of dust-absorbed light bracket the full range of the data. We also find that the trend is better explained by continuous star formation, rather than instantaneous burst models. We ascribe this result to the presence of multiple star clusters with different ages in many of the regions. Upper limits of the dust-only 8 $\mu$m emission as a function of age are provided.

The recent detections of the gravitational wave (GW) event GW170817 and its electromagnetic counterpart GRB 170817A produced by a binary neutron star (NS) merger is a new milestone of multimessenger astronomy. The time interval between these two signals has attracted widespread attention from physicists. In the braneworld scenario, GWs could propagate through the bulk while electromagnetic waves (EMWs) are bounded on the brane, i.e., our Universe. Therefore, the trajectories of GWs and EMWs may follow different pathes. If GWs and EMWs are originated simultaneously from the same source on the brane, they are expected to arrive at the observer successively. Consequently, the time delay between GW170817 and GRB 170817A may carry the information of the extra dimension. In this paper, we try to investigate the phenomenon in the context of a five-dimensional dS ($\text{dS}_5$) spacetime. We first study two special Universe models, i.e., de Sitter and Einstein-de Sitter models, and calculate the gravitation horizon radius in each case. For the real Universe, we then consider the $\Lambda$CDM model. Our results show that for the de Sitter model of the Universe, the $\text{dS}_5$ radius could not contribute to the time delay. With the data of the observation, we constrain the $\text{dS}_5$ radius to $\ell\gtrsim7.5\times10^{2}\,\text{Tpc}$ for the Einstein-de Sitter model and $\ell\gtrsim2.4\times10^{3}\,\text{Tpc}$ for the $\Lambda$CDM model. After considering the uncertainty in the source redshift and the time-lags given by different astrophysical processes of the binary NS merger, we find that our constraints are not sensitive to the redshift in the range of (0.005, 0.01) and the time-lag in the range of (-100s, 1.734s).

We introduce a novel methodology for establishing the presence of Standing Accretion Shock Instabilities (SASI) in the dynamics of a core collapse supernova from the observed neutrino event rate at water- or ice-based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test-statistics; it is also employed to assess the potential to estimate the frequency and the amplitude of the SASI modulations of the neutrino signal. The parameter estimation errors are consistent with the minimum possible errors as evaluated from the inverse of the Fisher information matrix, and close to the theoretical minimum for the SASI amplitude. Using results from a core-collapse simulation of a 15 solar-mass star by Kuroda $\it {et\, al.}$ (2017) as a test bed for the method, we find that SASI can be identified with high confidence for a distance to the supernova of up to $\sim 6$ kpc for IceCube and and up to $\sim 3$ kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will aid the investigation of a future galactic supernova.

An increasing number of young circumstellar disks show strikingly ordered (sub)millimeter polarization orientations along the minor axis, which is strong evidence for polarization due to scattering by ~0.1 mm sized grains. To test this mechanism further, we model the dust continuum and polarization data of HD 163296, one of the best observed disks with prominent rings and gaps, using the RADMC-3D radiative transfer code. We find that scattering by grains with a maximum size of 90$\mu$m can simultaneously reproduce the polarization observed at 0.87 mm (ALMA Band 7) and the unusually low spectral index of $\alpha$ ~ 1.5 between 0.87 and 1.25 mm (ALMA Band 6) in the optically thick inner disk as a result of more efficient scattering at a shorter wavelength. The relatively low spectral index of $\alpha$ ~ 2.5 inferred for the optically thin gaps is also reproduced by the same (relatively small) grains, as a result of telescope beam averaging of the gaps (with an intrinsic $\alpha$ ~ 4) and their adjacent optically thick rings (where $\alpha$ << 2). In this case, the long-standing tension between the grain sizes inferred from polarization and spectral index disappears because the relatively low $\alpha$ values are illusory and do not require large mm-sized grains. In addition, the polarization fraction has a unique pattern of azimuthal variation: higher along the major axis than the minor axis in the gaps but higher along the minor axis than the major axis in the rings. We find a rapidly declining polarization spectrum (with the fraction $p \propto \lambda^{-3}$ approximately) in the gaps, which becomes flattened or even inverted towards short wavelengths in the optically thick rings. These contrasting behaviors in the rings and gaps provide further tests of scattering-induced polarization that can be tested via multi-wavelength observations that resolve the disk substructure.

We study the relations between gas-phase metallicity ($Z$), local stellar mass surface density ($\Sigma_*$), and the local star formation surface density ($\Sigma_{\rm SFR}$) in a sample of 1120 star-forming galaxies from the MaNGA survey. At fixed $\Sigma_{*}$ the local metallicity increases as decreasing of $\Sigma_{\rm SFR}$ or vice versa for metallicity calibrators of N2 and O3N2. Alternatively, at fixed $\Sigma_{\rm SFR}$ metallicity increases as increasing of $\Sigma_{*}$, but at high mass region, the trend is flatter. However, the dependence of metallicity on $\Sigma_{\rm SFR}$ is nearly disappeared for N2O2 and N2S2 calibrators. We investigate the local metallicity against $\Sigma_{\rm SFR}$ with different metallicity calibrators, and find negative/positive correlations depending on the choice of the calibrator. We demonstrate that the O32 ratio (or ionization parameter) is probably dependent on star formation rate at fixed local stellar mass surface density. Additional, the shape of $\Sigma_*$ -- $Z$ -- $\Sigma_{\rm SFR}$ (FMR) depends on metallicity calibrator and stellar mass range. Since the large discrepancy between the empirical fitting-based (N2, O3N2) to electronic temperature metallicity and the photoionization model-dependent (N2O2, N2S2) metallicity calibrations, we conclude that the selection of metallicity calibration affects the existence of FMR on $\Sigma_{\rm SFR}$.

We present a study on the radial profiles of D4000,luminosity-weighted stellar ages $\tau_L$,and luminosity-weighted stellar metallicities $[Z/H]_L$ of 3654 nearby galaxies($0.01<z<0.15$)using the IFU spectroscopic data from the MaNGA survey available in the SDSS DR15,in an effort to explore the connection between median stellar population radial gradients($\nabla$D4000,$\nabla\tau_L,\nabla[Z/H]_L$)out to~$1.5R_e$ and various galaxy properties,including stellar mass($M_\star$),specific star formation rate(sSFR),morphologies,and local environment. We find that $M_\star$ is the single most predictive physical property for$\nabla$D4000 and$\nabla[Z/H]_L$. The most predictive properties for $\nabla\tau_L$ are sSFR,and to a lesser degree,$M_\star$. The environmental parameters,including local galaxy overdensities and central-satellite division,have virtually no correlation with stellar population radial profiles for the whole sample,but the $\nabla$D4000 of star-forming satellite galaxies with$M_\star\lesssim 10^{10}M_\odot$exhibit a significant positive correlation with galaxy overdensities. Galaxies with lower sSFR have on average steeper negative stellar population gradients,and this sSFR dependence is stronger for more massive star-forming galaxies. The negative correlation between the median stellar population gradients and$M_\star$ are best described largely as segmented relationships, whereby median gradients of galaxies with$\log M_\star\lesssim 10$(with the exact value depending on sSFR)have much weaker mass dependence than galaxies with higher$M_\star$. While the dependence of the radial gradients of ages and metallicities on T-Types and central stellar mass surface densities are generally not significant,galaxies with later T-Types or lower central mass densities tend to have significantly lower D4000,younger$\tau_L$ and lower$[Z/H]_L$ across the radial ranges probed in this study.

We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object very similar to mildly active Solar System comets with an out-gassing rate of $\sim$10$^{27}$ mol/sec. The photometry, taken in filters spanning the visible and NIR range shows a gradual brightening trend of $\sim0.03$ mags/day since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The lightcurve from recent and pre-discovery data reveals a brightness trend suggesting the recent onset of significant H$_2$O sublimation with the comet being active with super volatiles such as CO at heliocentric distances $>$6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 04, we estimate a diameter of 2I's nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/'Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size-distribution resulting in a slope of $\sim$3.4$\pm$1.2, similar to Solar System comets and bodies produced from collisional equilibrium.

We present the pilot study component of the Fluorescent Lyman-Alpha Structures in High-z Environments (FLASHES) Survey; the largest integral-field spectroscopy survey to date of the circumgalactic medium at $z=2.3-3.1$. We observed 48 quasar fields between 2015 and 2018 with the Palomar Cosmic Web Imager (Matuszewski et al. 2010). Extended HI Lyman-$\mathrm{\alpha}$ emission is discovered around 42/48 of the observed quasars, ranging in projected, flux-weighted radius from 21-71 proper kiloparsecs (pkpc), with 26 nebulae exceeding $100\mathrm{~pkpc}$ in effective diameter. The circularly averaged surface brightness radial profile peaks at a maximum of $\mathrm{1\times 10^{-17}~erg~s^{-1}~cm^{-2}~arcsec^{-2}}$ ($2\times10^{-15}~\mathrm{erg~s^{-1}~cm^{-2}~arcsec^{-2}}$ adjusted for cosmological dimming) and luminosities range from $1.9\times10^{43}~\mathrm{erg~s^{-1}}$ to $-14.1\times10^{43}~\mathrm{erg~s^{-1}}$. The emission appears to have a highly eccentric morphology and a maximum covering factor of $50\%$ ($60\%$ for giant nebulae). On average, the nebular spectra are red-shifted with respect to both the systemic redshift and Ly$\alpha$ peak of the quasar spectrum. The integrated spectra of the nebulae mostly have single or double-peaked line shapes with global dispersions ranging from $167~\mathrm{km~s^{-1}}$ to $690~\mathrm{km~s^{-1}}$, though the individual (Gaussian) components of lines with complex shapes mostly appear to have dispersions $\leq 400$ $\mathrm{km~s^{-1}}$, and the flux-weighted velocity centroids of the lines vary by thousands of $ \mathrm{km~s^{-1}}$ with respect to the systemic QSO redshifts. Finally, the root-mean-square velocities of the nebulae are found to be consistent with gravitational motions expected in dark matter halos of mass $\mathrm{M_h \simeq10^{12.5} M_\odot}$. We compare these results to existing surveys at both higher and lower redshift.

We compile a sample of about 157,000 spaxels from the Mapping Nearby Galaxies at the Apache Point Observatory survey to derive the average dust attenuation curve for subgalactic star-forming regions of local star-forming galaxies (SFGs) in the optical wavelength, following the method of \cite{Calzetti1994}. We obtain a $D_n(4000)$-independent average attenuation curve for spaxels with $1.1\leq D_n(4000)<1.3$, which is similar to the one derived from either local starbursts or normal SFGs. We examine whether and how the shape of the average attenuation curve changes with several local and global physical properties. For spaxels with $1.2\leq D_n(4000)<1.3$, we find no dependence on either local or global physical properties for the shape of the average attenuation curve. However, for spaxels with younger stellar population ($1.1\leq D_n(4000)<1.2$), shallower average attenuation curves are found for star-forming regions with smaller stellar mass surface density, smaller star formation rate surface density, or those residing in the outer region of galaxies. These results emphasize the risk of using one single attenuation curve to correct the dust reddening for all types of star-forming regions, especially for those with fairly young stellar population.

Polarized continuum emission from aligned grains in disks around young stellar objects can be used to probe the magnetic field, radiation anisotropy, or drift between dust and gas, depending on whether the non-spherical grains are aligned magnetically, radiatively or mechanically. We show that it can also be used to probe another key disk property -- the temperature gradient -- along sight lines that are optically thick, independent of the grain alignment mechanism. We first illustrate the technique analytically using a simple 1D slab model, which yields an approximate formula that relates the polarization fraction to the temperature gradient with respect to the optical depth tau at the tau=1 surface. The formula is then validated using models of stellar irradiated disks with and without accretion heating. The promises and challenges of the technique are illustrated with a number of Class 0 and I disks with ALMA dust polarization data, including NGC 1333 IRAS4A1, IRAS 16293B, BHB 07-11, L1527, HH 212 and HH 111. We find, in particular, that the sight lines passing through the near-side of a highly inclined disk trace different temperature gradient directions than those through the far-side, which can lead to a polarization orientation on the near-side that is orthogonal to that on the far-side, and that the HH 111 disk may be such a case. Our technique for probing the disk temperature gradient through dust polarization can complement other methods, particularly those using molecular lines.

The Beijing-Arizona Sky Survey (BASS) is a wide and deep imaging survey to cover a 5400 deg$^2$ area in the Northern Galactic Cap with the 2.3m Bok telescope using two filters ($g$ and $r$ bands). The Mosaic $z$-band Legacy Survey (MzLS) covers the same area in $z$ band with the 4m Mayall telescope. These two surveys will be used for spectroscopic targeting of the Dark Energy Spectroscopic Instrument (DESI). The BASS survey observations were completed in 2019 March. This paper describes the third data release (DR3) of BASS, which contains the photometric data from all BASS and MzLS observations between 2015 January and 2019 March. The median astrometric precision relative to {\it Gaia} positions is about 17 mas and the median photometric offset relative to the PanSTARRS1 photometry is within 5 mmag. The median $5\sigma$ AB magnitude depths for point sources are 24.2, 23.6, and 23.0 mag for $g$, $r$, and $z$ bands, respectively. The photometric depth within the survey area is highly homogeneous, with the difference between the 20\% and 80\% depth less than 0.3 mag. The DR3 data, including raw data, calibrated single-epoch images, single-epoch photometric catalogs, stacked images, and co-added photometric catalogs, are publicly accessible at \url{http://batc.bao.ac.cn/BASS/doku.php?id=datarelease:home}.

With the first observations of debris disks as well as proposed planets around white dwarfs, the question of how rocky planets around such stellar remnants can be characterized and probed for signs of life becomes tangible. White dwarfs are similar in size to Earth and have relatively stable environments for billions of years after initial cooling, making them intriguing targets for exoplanet searches and terrestrial planet atmospheric characterization. Their small size and the resulting large planet transit signal allows observations with next generation telescopes to probe the atmosphere of such rocky planets, if they exist. We model high-resolution transmission spectra for planets orbiting white dwarfs from as they cool from 6,000-4,000 K, for i) planets receiving equivalent irradiation to modern Earth, and ii) planets orbiting at the distance around a cooling white dwarf which allows for the longest continuous time in the habitable zone. All high-resolution transmission spectra will be publicly available online upon publication of this study and can be used as a tool to prepare and interpret upcoming observations with JWST, the Extremely Large Telescopes as well as mission concepts like Origins, HabEx, and LUVOIR.

The state-of-the-art geometry models of stars/dust suggest that dust attenuation toward nebular regions ($A_{V,gas}$) is always larger than that of stellar regions ($A_{V,star}$). Utilizing the newly released integral field spectroscopic data from the MaNGA survey, we investigate whether and how the $A_{V,star}/A_{V,gas}$ ratio varies from subgalactic to galactic scales. On a subgalactic scale, we report a stronger correlation between $A_{V,star}$ and $A_{V,gas}$ for more active HII regions. The local $A_{V,star}/A_{V,gas}$ is found to have moderate nonlinear correlations with three tracers of diffuse ionized gas (DIG), as well as indicators of gas-phase metallicity and ionization. The DIG regions tend to have larger $A_{V,star}/A_{V,gas}$ compared to classic HII regions excited by young OB stars. Metal-poor regions with a higher ionized level suffer much less nebular attenuation and thus have larger $A_{V,star}/A_{V,gas}$ ratios. A low-$A_{V,gas}$ and high-$A_{V,star}/A_{V,gas}$ sequence, which can be resolved into DIG-dominated and metal-poor regions, on the three BPT diagrams is found. Based on these observations, we suggest that besides the geometry of stars/dust, local physical conditions such as metallicity and ionized level also play an important role in determining the $A_{V,star}/A_{V,gas}$. On a galactic scale, the global $A_{V,star}/A_{V,gas}$ ratio has strong correlations with stellar mass ($M_*$), moderate correlations with SFR and metallicity, and weak correlations with inclination and specific SFR. Galaxies with larger $M_*$ and higher SFR that are more metal-rich tend to have smaller $A_{V,star}/A_{V,gas}$ ratios. Such correlations form a decreasing trend of $A_{V,star}/A_{V,gas}$ along the star-forming main sequence and mass-metallicity relation. The dust growth process accompanied by galaxy growth might be one plausible explanation for our observations.

We discuss the potential of detecting thermal neutrinos from matter-rich binary mergers, via a decades-long multi-messenger campaign involving a Mt-scale water Cherenkov neutrino detector and one or more next generation gravitational wave detectors, capable of observing mergers up to redshift $z\sim 2$. The search of neutrinos in time-coincidence with gravitational wave detections will allow to identify single neutrinos from individual mergers above the background, and to study their distributions in energy, redshift and type (double neutron-star or neutron-star-black hole merger) of the candidate sources. We find that, for merger rates consistent with current LIGO-Virgo constraints, and for a $100~{\rm Mt\cdot yr}$ exposure, between ${\mathcal O(10^{-1})}$ and ${\mathcal O(10)}$ neutrino events are expected. For extreme cases of mergers with more than $10^{52}$ ergs emitted in $\mathrel{{\bar \nu}_e}$, the number of events can be as large as $\sim 100$, with sensitivity to mergers up to redshift $z\sim 0.5$ or so. Such scenarios can already be tested with a $10~{\rm Mt\cdot yr}$ exposure, resulting in constraints on the post-merger evolution of the systems being considered.

The next generation of ground- and space-based telescopes will be able to observe rocky Earth-like planets in the near future, transiting their host star. We explore how the transmission spectrum of Earth changed through its geological history. These transmission spectra provide a template for how to characterize an Earth-like exoplanet - from a young prebiotic world to a modern Earth. They also allow us to explore at what point in its evolution a distant observer could identify life on our Pale Blue Dot and other worlds like it. We chose atmosphere models representative of five geological epochs of Earth's history, corresponding to a prebiotic high CO2-world 3.9 billion years ago (Ga), an anoxic world around 3.5 Ga, and 3 epochs through the rise of oxygen from 0.2 percent to present atmospheric levels of 21 percent. Our transmission spectra show atmospheric spectral features, which would show a remote observer that Earth had a biosphere since about 2 billion years ago. These high-resolution transmission spectral database of Earth through geological time from the VIS to the IR is available online and can be used as a tool to optimize our observation strategy, train retrieval methods, and interpret upcoming observations with JWST, the Extremely Large Telescopes and future mission concepts like Origins, HabEx, and LUOVIR.

The closest stars that harbor potentially habitable planets are cool M-stars. Upcoming ground- and space-based telescopes will be able to search the atmosphere of such planets for a range of chemicals. To facilitate this search and to inform upcoming observations, we model the high-resolution reflection spectra of two of the closest potentially habitable exoplanets for a range of terrestrial atmospheres and surface pressures for active and inactive phases of their host stars for both oxic and anoxic conditions: Proxima b, the closest potentially habitable exoplanet, and Trappist-1e, one of 3 Earth-size planets orbiting in the Habitable Zone of Trappist-1. We find that atmospheric spectral features, including biosignatures like O2 in combination with a reduced gas like CH4 for oxic atmospheres, as well as climate indicators like CO2 and H2O for all atmospheres, show absorption features in the spectra of Proxima b and Trappist-1e models. However for some features like oxygen, high-resolution observations will be critical to identify them in a planet's reflected flux. Thus these two planets will be among the best targets for upcoming observations of potential Earth-like planets in reflected light with planned Extremely Large Telescopes.

The spectra of emission-line galaxies (ELGs) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digit Sky Survey (SDSS) are used to study the mass-metallicity relation (MZR) at $z\sim0.8$. The selected sample contains about 180,000 massive star-forming galaxies with $0.6 < z < 1.05$ and $9 < {\rm log}(M_{\star}/M_{\odot}) < 12$. The spectra are stacked in bins of different parameters including redshift, stellar mass, star formation rate (SFR), specific star formation rate (sSFR), half-light radius, mass density, and optical color. The average MZR at $z\sim0.83$ has a downward evolution in the MZR from local to high-redshift universe, which is consistent with previous works. At a specified stellar mass, galaxies with higher SFR/sSFR and larger half-light radius have systematically lower metallicity. This behavior is reversed for galaxies with larger mass density and optical color. Among the above physical parameters, the MZR has the most significant dependency on SFR. Our galaxy sample at $0.6<z<1.05$ approximately follows the fundamental metallicity relation (FMR) in the local universe, although the sample inhomogeneity and incompleteness might have effect on our MZR and FMR.

The GJ 357 system harbors 3 planets orbiting a bright, nearby M2.5V star at 9.44pc. The innermost planet GJ 357 b (TOI-562.01) is a hot transiting Earth-size planet with Earth-like density, which receives about 12 times the irradiation Earth receives from the Sun, and was detected using data from TESS. Radial velocities discovered two more planets in the system at 9.12 (GJ 357 c) and 55.6 days (GJ 357 d), with minimum masses of 3.59+/-0.50 and 6.1+/-1 Earth masses, and an irradiation of 4.4 and 0.38 Earths irradiation, respectively. GJ 357 d receives slightly less stellar irradiation than Mars does in our own Solar System, which puts it in the Habitable Zone for its host star. GJ 357 d could not have been detected with TESS and whether it transits remains an open question. Here we model under what conditions GJ 357 d could sustain surface habitability and present planetary models as well as synthetic transmission, reflection and emission spectra for a range of models for GJ 357 d from water worlds to Earth-like models. With Earth-analog outgassing rates, GJ 357 d would be a frozen rocky world, however with an increased CO2 level, as would be expected if a geological cycle regulates carbon dioxide concentration like on Earth, the planet models show temperate surface conditions. If we can detect a transit of GJ 357 d, it would become the closest transiting, potentially habitable planet in the solar neighborhood. Even if GJ 357 d does not transit, the brightness of its star makes this planet in the Habitable Zone of a close-by M star a prime target for observations with Extremely Large telescopes as well as future space missions.

On 27 Apr 2015, when 67P/C-G was at 1.76 au from the Sun and moving towards perihelion, the OSIRIS and VIRTIS-M instruments on Rosetta observed the evolving dust and gas coma during a complete rotation of the comet. We aim to characterize the dust, H2O and CO2 gas spatial distribution in the inner coma. To do this we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the one calculated using a thermo-physical model. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O) and 4200 nm (CO2). We measured the average signal in a circular annulus, to study spatial variation around the comet, and in a sector of the annulus, to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet centre. The spatial correlation between dust and water, both coming from the sun-lit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0deg sub-solar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production, with respect to the value calculated using a simple thermo-physical model, is observed when the head lobe and regions of the Southern hemisphere with strong seasonal variations are illuminated. A drastic decrease in dust production, when the water production (both measured and from the model) displays a maximum, happens when typical Northern consolidated regions are illuminated and the Southern hemisphere regions with strong seasonal variations are instead in shadow.

Many studies have shown that RR Lyrae variable stars (RRL) are powerful stellar tracers of Galactic halo structure and satellite galaxies. The Dark Energy Survey (DES), with its deep and wide coverage (g ~ 23.5 mag) in a single exposure; over 5000 deg$^{2}$) provides a rich opportunity to search for substructures out to the edge of the Milky Way halo. However, the sparse and unevenly sampled multiband light curves from the DES wide-field survey (median 4 observations in each of grizY over the first three years) pose a challenge for traditional techniques used to detect RRL. We present an empirically motivated and computationally efficient template fitting method to identify these variable stars using three years of DES data. When tested on DES light curves of previously classified objects in SDSS stripe 82, our algorithm recovers 89% of RRL periods to within 1% of their true value with 85% purity and 76% completeness. Using this method, we identify 5783 RRL candidates, ~31% of which are previously undiscovered. This method will be useful for identifying RRL in other sparse multiband data sets.

Main-belt asteroid (6478) Gault unexpectedly sprouted two tails in late 2018 and early 2019, identifying it as a new active asteroid. Here we present observations obtained by the 1.2-m Zwicky Transient Facility survey telescope that provide detailed time-series coverage of the onset and evolution of Gault's activity. Gault exhibited two brightening events, with the first one starting on 2018 Oct. 18$\pm5$ days and a second one starting on 2018 Dec. 24$\pm1$ days. The amounts of mass released are $2\times10^7$ kg and $1\times10^6$ kg, respectively. Based on photometric measurements, each event persisted for about a month. Gault's color has not changed appreciably over time, with a pre-outburst color of $g_\mathrm{PS1}-r_\mathrm{PS1}=0.50\pm0.04$ and $g_\mathrm{PS1}-r_\mathrm{PS1}=0.46\pm0.04$ during the two outbursts. Simulations of dust dynamics shows that the ejecta consists of dust grains of up to 10 $\mu$m in size that are ejected at low velocities below $1~\mathrm{m~s^{-1}}$ regardless of particle sizes. This is consistent with non-sublimation-driven ejection events. The size distribution of the dust exhibits a broken power-law, with particles at 10--20 $\mu$m following a power-law of $-2.5$ to $-3.0$, while larger particles follow a steeper slope of $-4.0$. The derived properties can be explained by either rotational excitation of the nucleus or a merger of a near-contact binary, with the latter scenario to be statistically more likely.

Atmospheric transmission from Dome A, Antarctica, presents new possibilities in the field of terahertz astronomy, where space telescopes have been the only observational tools until now. Using atmospheric transmission measurements from Dome A with a Fourier transform spectrometer, transmission spectra and long-term stabilities have been analyzed at 1.461 THz, 3.393 THz, 5.786 THz and 7.1 THz, which show that important atmospheric windows for terahertz astronomy open for a reasonable length of time in the winter season. With large aperture terahertz telescopes and interferometers at Dome A, high angular resolution terahertz observations are foreseen of atomic fine-structure lines from ionized gas and a water ice feature from protoplanetary disks.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

We present observations of supernova (SN) 2017ens, discovered by the ATLAS survey and identified as a hot blue object through the GREAT program. The redshift z=0.1086 implies a peak brightness of M_g=-21.1 mag, placing the object within the regime of superluminous supernovae. We observe a dramatic spectral evolution, from initially being blue and featureless, to later developing features similar to those of the broadlined Type Ic SN 1998bw, and finally showing ~2000 km s^-1 wide H-alpha and H-beta emission. Relatively narrow Balmer emission (reminiscent of a SN IIn) is present at all times. We also detect coronal lines, indicative of a dense circumstellar medium. We constrain the progenitor wind velocity to ~50-60 km s^-1 based on P-Cygni profiles, which is far slower than those present in Wolf-Rayet stars. This may suggest that the progenitor passed through a luminous blue variable phase, or that the wind is instead from a binary companion red supergiant star. At late times we see the ~2000 km s^-1 wide H-alpha emission persisting at high luminosity (~3x10^40 erg s^-1) for at least 100 day, perhaps indicative of additional mass loss at high velocities that could have been ejected by a pulsational pair instability.

The stellar mass-metallicity relation ($M_* - Z$, MZR) indicates that the metallicities of galaxies increase with increasing stellar masses. The fundamental metallicity relation (FMR) suggests that the galaxies with higher star formation rates (SFRs) tend to have lower metallicities for a given stellar mass. To examine whether the MZR and FMR still hold at poorer metallicities and higher redshifts, we compile a sample of 35 star-forming galaxies (SFGs) at $0.6<z<0.9$ using the public spectral database ($\rm v5\_10\_0$) of emission-line galaxies from the extended Baryon Oscillation Spectroscopic Survey (eBOSS). These galaxies are identified for their significant auroral $\rm [OIII]\lambda4363$ emission line ($\rm S/N \geq 3$). With the electronic temperature metallicity calibration, we find nine SFGs are extremely metal-poor galaxies with $\rm 12 + log(O/H) \leq 7.69 \ (1/10 \ Z_\odot)$. The metallicity of the most metal-deficient galaxy is $7.35\pm 0.09$ (about 1/20 $Z_{\odot}$). Compared with the SFGs with normal metallicities in local and high redshift universe, our metal-poor SFGs have more than ten times higher SFRs at a fixed stellar mass. We create a new mass -- SFR relation for these metal-poor galaxies at $0.6<z<0.9$. Due to the higher SFRs and younger stellar ages, our metal-poor SFGs deviate from the MZR and FMR in the local universe toward lower metallicities, confirming the existence of FMR, as well as the cosmic evolution of MZR and FMR with redshift.

The metallicity and its relationship with other galactic properties is a fundamental probe of the evolution of galaxies. In this work, we select about 750,000 star-forming spatial pixels from 1122 blue galaxies in the MaNGA survey to investigate the global stellar mass - local stellar mass surface density - gas-phase metallicity ($M_*$ - $\Sigma_*$ - $Z$ ) relation. At a fixed $M_*$, the metallicity increases steeply with increasing $\Sigma_*$. Similarly, at a fixed $\Sigma_*$, the metallicity increases strongly with increasing $M_*$ at low mass end, while this trend becomes less obvious at high mass end. We find the metallicity to be more strongly correlated to $\Sigma_*$ than to $M_*$. Furthermore, we construct a tight (0.07 dex scatter) $M_*$ - $\Sigma_*$ - $Z$ relation, which reduces the scatter in the $\Sigma_*$ - $Z$ relation by about 30$\%$ for galaxies with $7.8 < {\rm log}(M_*/M_\odot) < 11.0$, while the reduction of scatter is much weaker for high-mass galaxies. This result suggests that, especially for low-mass galaxies, the $M_*$ - $\Sigma_*$ - $Z$ relation is largely more fundamental than the $M_*$ - $Z$ and $\Sigma_*$ - $Z$ relations, meaning that both $M_*$ and $\Sigma_*$ play important roles in shaping the local metallicity. We also find that the local metallicity is probably independent on the local star formation rate surface density at a fixed $M_*$ and $\Sigma_*$. Our results are consistent with the scenario that the local metallicities in galaxies are shaped by the combination of the local stars formed in the history and the metal loss caused by galactic winds.

The phase function of the dust coma of comet 67P has been determined from Rosetta/OSIRIS images \citep{Bertini17}. This function show a deep minimum at phase angles near 100$^\circ$, and a strong backscattering enhancement. These two properties cannot be reproduced by regular models of cometary dust, most of them based on wavelength-sized and randomly-oriented aggregate particles. We show, however, that an ensamble of oriented elongated particles of a wide variety of aspect ratios, with radii $r \gtrsim$10 $\mu$m, and whose long axes are perpendicular to the direction of the solar radiation, are capable of reproducing the observed phase function. These particles must be absorbing, with an imaginary part of the refractive index of about 0.1 to match the expected geometric albedo, and with porosity in the 60-70\% range.

We report observations of the reactivations of main-belt comets 238P/Read and 288P/(300163) 2006 VW139, that also track the evolution of each object's activity over several months in 2016 and 2017. We additionally identify and analyze archival SDSS data showing 288P to be active in 2000, meaning that both 238P and 288P have now each been confirmed to be active near perihelion on three separate occasions. From data obtained of 288P from 2012-2015 when it appeared inactive, we find best-fit R-band H,G phase function parameters of H_R=16.80+/-0.12 mag and G_R=0.18+/-0.11, corresponding to effective component radii of r_c=0.80+/-0.04 km, assuming a binary system with equally-sized components. Fitting linear functions to ejected dust masses inferred for 238P and 288P soon after their observed reactivations in 2016, we find an initial average net dust production rate of 0.7+/-0.3 kg/s and a best-fit start date of 2016 March 11 (when the object was at a true anomaly of -63 deg) for 238P, and an initial average net dust production rate of 5.6+/-0.7 kg/s and a best-fit start date of 2016 August 5 (when the object was at a true anomaly of -27 deg) for 288P. Applying similar analyses to archival data, we find similar start points for previous active episodes for both objects, suggesting that minimal mantle growth or ice recession occurred between the active episodes in question. Some changes in dust production rates between active episodes are detected, however. More detailed dust modeling is suggested to further clarify the process of activity evolution in main-belt comets.

The near-Earth asteroid 3200 Phaethon (1983 TB) is an attractive object not only from a scientific viewpoint but also because of JAXA's DESTINY+ target. The rotational lightcurve and spin properties were investigated based on the data obtained in the ground-based observation campaign of Phaethon. We aim to refine the lightcurves and shape model of Phaethon using all available lightcurve datasets obtained via optical observation, as well as our time-series observation data from the 2017 apparition. Using eight 1-2-m telescopes and an optical imager, we acquired the optical lightcurves and derived the spin parameters of Phaethon. We applied the lightcurve inversion method and SAGE algorithm to deduce the convex and non-convex shape model and pole orientations. We analysed the optical lightcurve of Phaethon and derived a synodic and a sidereal rotational period of 3.6039 h, with an axis ratio of a/b = 1.07. The ecliptic longitude (lambda) and latitude (beta) of the pole orientation were determined as (308, -52) and (322, -40) via two independent methods. A non-convex model from the SAGE method, which exhibits a concavity feature, is also presented.

The Rosetta space probe accompanied comet 67P/Churyumov-Gerasimenko for more than two years, obtaining an unprecedented amount of unique data of the comet nucleus and inner coma. This work focuses identifying the source regions of faint jets and outbursts and on studying the spectrophotometric properties of some outbursts. We use observations acquired with the OSIRIS/NAC camera during July-October 2015, that is, close to perihelion. More than 200 jets of different intensities were identified directly on the nucleus. Some of the more intense outbursts appear spectrally bluer than the comet dark terrain in the vivible-to-near-infrared region. We attribute this spectral behavior to icy grains mixed with the ejected dust. Some of the jets have an extremely short lifetime. They appear on the cometary surface during the color sequence observations, and vanish in less than some few minutes after reaching their peak. We also report a resolved dust plume observed in May 2016 at a resolution of 55 cm/pixel, which allowed us to estimate an optical depth of $\sim$0.65 and an ejected mass of $\sim$ 2200 kg. We present the results on the location, duration, and colors of active sources on the nucleus of 67P from the medium-resolution (i.e., 6-10 m/pixel) images acquired close to perihelion passage. The observed jets are mainly located close to boundaries between different morphological regions. Jets depart not only from cliffs, but also from smooth and dust-covered areas, from fractures, pits, or cavities that cast shadows and favor the recondensation of volatiles. This study shows that faint jets or outbursts continuously contribute to the cometary activity close to perihelion passage, and that these events are triggered by illumination conditions. Faint jets or outbursts are not associated with a particular terrain type or morphology.

We analyze sparsely-sampled near-infrared (JHKs) light curves of a sample of 1781 Mira variable candidates in M33, originally discovered using I-band time-series observations. We extend our single-band semi-parametric Gaussian process modeling of Mira light curves to a multi-band version and obtain improved period determinations. We use our previous results on near-infrared properties of candidate Miras in the LMC to classify the majority of the M33 sample into Oxygen- or Carbon-rich subsets. We derive Period-Luminosity relations for O-rich Miras and determine a distance modulus for M33 of 24.80 +- 0.06 mag.

We present a new mechanism for deep crustal heating in accreting neutron stars. Charged pions ($\pi^+$) are produced in nuclear collisions on the neutron star surface during active accretion and upon decay they provide a flux of neutrinos into the neutron star crust. For massive and/or compact neutron stars, neutrinos deposit $\approx 1\textrm{--} 2 \, \mathrm{MeV}$ of heat per accreted nucleon into the inner crust. The strength of neutrino heating is comparable to the previously known sources of deep crustal heating, such as from pycnonuclear fusion reactions, and is relevant for studies of cooling neutron stars. We model the thermal evolution of a transient neutron star in a low-mass X-ray binary, and in the particular case of the neutron star MXB~1659-29 we show that additional deep crustal heating requires a higher thermal conductivity for the neutron star inner crust. A better knowledge of pion production cross sections near threshold would improve the accuracy of our predictions.

We investigate the spatially-resolved star formation main sequence in star-forming galaxies (SFGs) using Integral Field Spectroscopic (IFS) observations from the Mapping Nearby Galaxies at the Apache Point Observatory (MaNGA) survey. We demonstrate that the correlation between the stellar mass surface density ($\Sigma_*$) and star formation rate surface density ($\Sigma_{\mathrm{SFR}}$) holds down to sub-galactic scale, leading to the Sub-Galactic Main Sequence (SGMS). By dividing galaxies into two populations based on their recent mass assembly modes, we find the resolved main sequence in galaxies with 'outside-in' mode is steeper than that in galaxies with 'inside-out' mode. This is also confirmed on a galaxy-by-galaxy level, where we find the distributions of SGMS slopes for individual galaxies are clearly separated for the two populations. When normalizing and stacking the SGMS of individual galaxies on one panel for the two populations, we find the inner regions of galaxies with 'inside-out' mode statistically exhibit a suppression in star formation, with a less significant trend in the outer regions of galaxies with 'outside-in' mode. In contrast, the inner regions of galaxies with 'outside-in' mode and the outer regions of galaxies with 'inside-out' mode follow a slightly sub-linear scaling relation with a slope $\sim$0.9, which is in good agreement with previous findings, suggesting that they are experiencing a universal regulation without influences of additional physical processes.

By using the Hectospec 6.5 m Multiple Mirror Telescope (MMT) and the 2.16 m telescope of National Astronomical Observatories, Chinese Academy of Sciences (NAOC), we obtained 188 high signal-to-noise ratio (S/N) spectra of HII regions in the nearby galaxy M101, which are the largest spectroscopic sample of HII regions for this galaxy so far. These spectra cover a wide range of regions on M101, which enables us to analyze two dimensional distributions of its physical properties. The physical parameters are derived from emission lines or stellar continuum, including stellar population age, electron temperature, oxygen abundance and etc. The oxygen abundances are derived using two empirical methods based on O3N2 and R$_{23}$ indicators, as well as the direct Te method when OIII$\lambda4363$ is available. By applying the harmonic decomposition analysis to the velocity field, we obtained line-of-sight rotation velocity of 71 km s$^{-1}$ and a position angle of 36 degree. The stellar age profile shows an old stellar population in galaxy center and a relative young stellar population in outer regions, suggesting an old bulge and a young disk. Oxygen abundance profile exhibits a clear break at $\sim$18 kpc, with a gradient of $-$0.0364 dex kpc$^{-1}$ in the inner region and $-$0.00686 dex kpc$^{-1}$ in the outer region. Our results agree with the "inside-out" disk growth scenario of M101.

Core-collapse supernova explosions are driven by a central engine that converts a small fraction of the gravitational binding energy released during core collapse to outgoing kinetic energy. The suspected mode for this energy conversion is the neutrino mechanism, where a fraction of the neutrinos emitted from the newly formed protoneutron star are absorbed by and heat the matter behind the supernova shock. Accurate neutrino-matter interaction terms are crucial for simulating these explosions. In this proceedings for IAUS 331, SN 1987A, 30 years later, we explore several corrections to the neutrino-nucleon scattering opacity and demonstrate the effect on the dynamics of the core-collapse supernova central engine via two dimensional neutrino-radiation-hydrodynamics simulations. Our results reveal that the explosion properties are sensitive to corrections to the neutral-current scattering cross section at the 10-20% level, but only for densities at or above $\sim 10^{12}$ g cm$^{-3}$

Core collapse supernova simulations can be sensitive to neutrino interactions near the neutrinosphere. This is the surface of last scattering. We model the neutrinosphere region as a warm unitary gas of neutrons. A unitary gas is a low density system of particles with large scattering lengths. We calculate modifications to neutrino scattering cross sections because of the universal spin and density correlations of a unitary gas. These correlations can be studied in laboratory cold atom experiments. We find significant reductions in cross sections, compared to free space interactions, even at relatively low densities. These reductions could reduce the delay time from core bounce to successful explosion in multidimensional supernova simulations.

On 3 July 2016, several instruments on board ESA's Rosetta spacecraft detected signs of an outburst event on comet 67P, at a heliocentric distance of 3.32 AU from the sun, outbound from perihelion. We here report on the inferred properties of the ejected dust and the surface change at the site of the outburst. The activity coincided with the local sunrise and continued over a time interval of 14 - 68 minutes. It left a 10m-sized icy patch on the surface. The ejected material comprised refractory grains of several hundred microns in size, and sub-micron-sized water ice grains. The high dust mass production rate is incompatible with the free sublimation of crystalline water ice under solar illumination as the only acceleration process. Additional energy stored near the surface must have increased the gas density. We suggest a pressurized sub-surface gas reservoir, or the crystallization of amorphous water ice as possible causes.

We present a statistical analysis of the distribution of large scale topographic features on comet 67P/Churyumov-Gerasimenko. We observe that the cumulative cliff height distribution across the surface follows a power law with a slope equal to -1.69 +- 0.02. When this distribution is studied independently for each region, we find a good correlation between the slope of the power law and the orbital erosion rate of the surface. For instance, the northern hemisphere topography is dominated by structures on the 100~m scale while the southern hemisphere topography, illuminated at perihelion, is dominated by 10~m scale terrain features. Our study suggest that the current size of a cliff is controlled not only by material cohesion but by the dominant erosional process in each region. This observation can be generalized to other comets, where we argue that primitive nuclei are characterized by the presence of large cliffs with a cumulative height power index equal to or above -1.5, while older, eroded cometary surfaces have a power index equal to or below -2.3. In effect, our model shows that a measure of the topography provides a quantitative assessment of a comet's erosional history, i.e. its evolutionary age.

The Southern hemisphere of the 67P/Churyumov-Gerasimenko comet has become visible from Rosetta only since March 2015. It was illuminated during the perihelion passage and therefore it contains the regions that experienced the strongest heating and erosion rate, thus exposing the subsurface most pristine material. In this work we investigate, thanks to the OSIRIS images, the geomorphology, the spectrophotometry and some transient events of two Southern hemisphere regions: Anhur and part of Bes. Bes is dominated by outcropping consolidated terrain covered with fine particle deposits, while Anhur appears strongly eroded with elongated canyon-like structures, scarp retreats, different kinds of deposits, and degraded sequences of strata indicating a pervasive layering. We discovered a new 140 m long and 10 m high scarp formed in the Anhur/Bes boundary during/after the perihelion passage, close to the area where exposed CO$_2$ and H$_2$O ices were previously detected. Several jets have been observed originating from these regions, including the strong perihelion outburst, an active pit, and a faint optically thick dust plume. We identify several areas with a relatively bluer slope (i.e. a lower spectral slope value) than their surroundings, indicating a surface composition enriched with some water ice. These spectrally bluer areas are observed especially in talus and gravitational accumulation deposits where freshly exposed material had fallen from nearby scarps and cliffs. The investigated regions become spectrally redder beyond 2 au outbound when the dust mantle became thicker, masking the underlying ice-rich layers.

Images of the nucleus and the coma (gas and dust) of comet 67P/Churyumov- Gerasimenko have been acquired by the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras since March 2014 using both the Wide Angle Camera (WAC) and the Narrow Angle Camera (NAC). We use images from the NAC camera to study a bright outburst observed in the southern hemisphere on July 29, 2015. The high spatial resolution of the NAC is needed to localize the source point of the outburst on the surface of the nucleus. The heliocentric distance is 1.25 au and the spacecraft-comet distance is 186 km. Aiming to better understand the physics that led to the outgassing, we used the Direct Simulation Monte Carlo (DSMC) method to study the gas flow close to the nucleus and the dust trajectories. The goal is to understand the mechanisms producing the outburst. We reproduce the opening angle of the outburst in the model and constrain the outgassing ratio between the outburst source and the local region. The outburst is in fact a combination of both gas and dust, in which the active surface is approximately 10 times more active than the average rate found in the surrounding areas. We need a number of dust particles 7.83 $\times$ 10$^{11}$ - 6.90 $\times$ 10$^{15}$ (radius 1.97 - 185 {\mu}m), which corresponds to a mass of dust 220 - 21 $\times$ 10$^{3}$kg.

We present a summary of the campaign of remote observations that supported the European Space Agency's Rosetta mission. Telescopes across the globe (and in space) followed comet 67P/Churyumov-Gerasimenko from before Rosetta's arrival until nearly the end of mission in September 2016. These provided essential data for mission planning, large-scale context information for the coma and tails beyond the spacecraft, and a way to directly compare 67P with other comets. The observations revealed 67P to be a relatively `well behaved' comet, typical of Jupiter family comets and with activity patterns that repeat from orbit-to-orbit. Comparison between this large collection of telescopic observations and the in situ results from Rosetta will allow us to better understand comet coma chemistry and structure. This work is just beginning as the mission ends -- in this paper we present a summary of the ground-based observations and early results, and point to many questions that will be addressed in future studies.

(144977) 2005 EC127 is an V-/A-type inner-main-belt asteroid with a diameter of 0.6 +- 0.1 km. Asteroids of this size are believed to have rubble-pile structure, and, therefore, cannot have a rotation period shorter than 2.2 hours. However, our measurements show that asteroid 2005 EC127 completes one rotation in 1.65 +- 0.01 hours with a peak-to-peak light-curve variation of ~0.5 mag. Therefore, this asteroid is identified as a large super-fast rotator. Either a rubble-pile asteroid with a bulk density of ~6 g cm^-3 or an asteroid with an internal cohesion of 47 +- 30 Pa can explain 2005 EC127. However, the scenario of high bulk density is very unlikely for asteroids. To date, only six large super-fast rotators, including 2005 EC127, have been reported, and this number is very small when compared with the much more numerous fast rotators. We also note that none of the six reported large SFRs are classified as C-type asteroids.

The spectra of 413 star-forming (or HII) regions in M33 (NGC 598) were observed by using the multifiber spectrograph of Hectospec at the 6.5-m Multiple Mirror Telescope (MMT). By using this homogeneous spectra sample, we measured the intensities of emission lines and some physical parameters, such as electron temperatures, electron densities, and metallicities. Oxygen abundances were derived via the direct method (when available) and two empirical strong-line methods, namely, O3N2 and N2. In the high-metallicity end, oxygen abundances derived from O3N2 calibration were higher than those derived from N2 index, indicating an inconsistency between O3N2 and N2 calibrations. We presented a detailed analysis of the spatial distribution of gas-phase oxygen abundances in M33 and confirmed the existence of the axisymmetric global metallicity distribution widely assumed in literature. Local variations were also observed and subsequently associated with spiral structures to provide evidence of radial migration driven by arms. Our O/H gradient fitted out to 1.1 $R_{25}$ resulted in slopes of $-0.17\pm0.03$, $-0.19\pm0.01$, and $-0.16\pm0.17$ dex $R_{25}^{-1}$ utilizing abundances from O3N2, N2 diagnostics, and direct method, respectively.

We present a sample of 48 metal-poor galaxies at $ z < 0.14$ selected from 92,510 galaxies in the LAMOST survey. These galaxies are identified for their detection of the auroral emission line \oiii$\lambda$4363 above $3\sigma$ level, which allows a direct measurement of the electron temperature and the oxygen abundance. The emission line fluxes are corrected for internal dust extinction using Balmer decrement method. With electron temperature derived from \oiii$\lambda\lambda4959,5007/\oiii\lambda4363$ and electron density from $\sii\lambda6731/\sii\lambda6717$, we obtain the oxygen abundances in our sample which range from $\zoh= 7.63$ (0.09 $\Zsun$) to $8.46$ (0.6 $\Zsun$). We find an extremely metal-poor galaxy with $\zoh=7.63 \pm 0.01$. With multiband photometric data from FUV to NIR and $\ha$ measurements, we also determine the stellar masses and star formation rates, based on the spectral energy distribution fitting and $\ha$ luminosity, respectively. We find that our galaxies have low and intermediate stellar masses with $\rm 6.39 \le log(M/M_{\sun})\le 9.27$, and high star formation rates (SFRs) with $\rm -2.18 \le log(SFR/M_{\sun} yr^{-1}) \le 1.95$. We also find that the metallicities of our galaxies are consistent with the local $T_e$-based mass-metallicity relation, while the scatter is about 0.28 dex. Additionally, assuming the coefficient of $\rm \alpha =0.66$, we find most of our galaxies follow the local mass-metallicity-SFR relation, while a scatter about 0.24 dex exists, suggesting the mass-metallicity relation is weakly dependent on SFR for those metal-poor galaxies.

Nuclear pasta, with nucleons arranged into tubes, sheets, or other complex shapes, is expected in core collapse supernovae (SNe) at just below nuclear density. We calculate the additional opacity from neutrino-pasta coherent scattering using molecular dynamics simulations. We approximately include this opacity in simulations of SNe. We find that pasta slows neutrino diffusion and greatly increases the neutrino signal at late times of 10 or more seconds after stellar core collapse. This signal, for a galactic SN, should be clearly visible in large detectors such as Super-Kamiokande.

We generalize our virial approach to study the neutral current neutrino response of nuclear matter at low densities. In the long-wavelength limit, the virial expansion makes model-independent predictions for neutrino-nucleon scattering rates and the density S_V and spin S_A responses. We find S_A is significantly reduced from one even at low densities. We provide a simple fit S_A^f(n,T,Y_p) of the axial response as a function of density n, temperature T and proton fraction Y_p. This fit reproduces our model independent virial results at low densities and reproduces the Burrows and Sawyer random phase approximation (RPA) results at high densities. Our fit can be incorporated into supernova simulations in a straight forward manner. Preliminary one dimensional supernova simulations suggest that the reduction in the axial response may enhance neutrino heating rates in the gain region during the accretion phase of a core-collapse supernovae.

Template-based extrapolations from only one photometric band can be a cost-effective method to estimate the total infrared (IR) luminosities ($L_{\mathrm{IR}}$) of galaxies. By utilizing multi-wavelength data that covers across 0.35--500\,$\mathrm{\mu m}$ in GOODS-North and GOODS-South fields, we investigate the accuracy of this monochromatic extrapolated $L_{\mathrm{IR}}$ based on three IR spectral energy distribution (SED) templates (\citealt[CE01]{Chary2001}; \citealt[DH02]{Dale2002}; \citealt[W08]{Wuyts2008a}) out to $z\sim 3.5$. We find that the CE01 template provides the best estimate of $L_{\mathrm{IR}}$ in {\it Herschel}/PACS bands, while the DH02 template performs best in {\it Herschel}/SPIRE bands. To estimate $L_{\mathrm{IR}}$, we suggest that extrapolations from the available longest wavelength PACS band based on the CE01 template can be a good estimator. Moreover, if PACS measurement is unavailable, extrapolations from SPIRE observations but based on the \cite{Dale2002} template can also provide a statistically unbiased estimate for galaxies at $z\lesssim 2$. The emission of rest-frame 10--100\,$\mathrm{\mu m}$ range of IR SED can be well described by all the three templates, but only the DH02 template shows nearly unbiased estimate of the emission of the rest-frame submillimeter part.

During its two years mission around comet 67P/Churyumov-Gerasimenko, ESA's Rosetta spacecraft had the unique opportunity to follow closely a comet in the most active part of its orbit. Many studies have presented the typical features associated to the activity of the nucleus, such as localized dust and gas jets. Here we report on series of more energetic transient events observed during the three months surrounding the comet's perihelion passage in August 2015. We detected and characterized 34 outbursts with the Rosetta cameras, one every 2.4 nucleus rotation. We identified 3 main dust plume morphologies associated to these events: a narrow jet, a broad fan, and more complex plumes featuring both previous types together. These plumes are comparable in scale and temporal variation to what has been observed on other comets. We present a map of the outbursts source locations, and discuss the associated topography. We find that the spatial distribution sources on the nucleus correlates well with morphological region boundaries, especially in areas marked by steep scarps or cliffs. Outbursts occur either in the early morning or shortly after the local noon, indicating two potential processes: Morning outbursts may be triggered by thermal stresses linked to the rapid change of temperature, afternoon events are most likely related to the diurnal or seasonal heat wave reaching volatiles buried under the first surface layer. In addition, we propose that some events can be the result of a completely different mechanism, in which most of the dust is released upon the collapse of a cliff.

The terahertz and far-infrared (FIR) band, from approximately 0.3 THz to 15 THz (1 mm to 20 micron), is important for astrophysics as the thermal radiation of much of the universe peaks at these wavelengths and many spectral lines that trace the cycle of interstellar matter also lie within this band. However, water vapor renders the terrestrial atmosphere opaque to this frequency band over nearly all of the Earth's surface. Early radiometric measurements below 1 THz at Dome A, the highest point of the cold and dry Antarctic ice sheet, suggest that it may offer the best possible access for ground-based astronomical observations in the terahertz and FIR band. To address uncertainty in radiative transfer modelling, we carried out measurements of atmospheric radiation from Dome A spanning the entire water vapor pure rotation band from 20 micron to 350 micron wavelength by a Fourier transform spectrometer. Our measurements expose atmospheric windows having significant transmission throughout this band. Furthermore, by combining our broadband spectra with auxiliary data on the atmospheric state over Dome A, we set new constraints on the spectral absorption of water vapor at upper tropospheric temperatures important for accurately modeling the terrestrial climate. In particular, we find that current spectral models significantly underestimate the H2O continuum absorption.

Beginning in March 2014, the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras began capturing images of the nucleus and coma (gas and dust) of comet 67P/Churyumov-Gerasimenko using both the wide angle camera (WAC) and the narrow angle camera (NAC). The many observations taken since July of 2014 have been used to study the morphology, location, and temporal variation of the comet's dust jets. We analyzed the dust monitoring observations shortly after the southern vernal equinox on May 30 and 31, 2015 with the WAC at the heliocentric distance Rh = 1.53 AU, where it is possible to observe that the jet rotates with the nucleus. We found that the decline of brightness as a function of the distance of the jet is much steeper than the background coma, which is a first indication of sublimation. We adapted a model of sublimation of icy aggregates and studied the effect as a function of the physical properties of the aggregates (composition and size). The major finding of this article was that through the sublimation of the aggregates of dirty grains (radius a between 5 microm and 50 microm) we were able to completely reproduce the radial brightness profile of a jet beyond 4 km from the nucleus. To reproduce the data we needed to inject a number of aggregates between 8.5 x $10^{13}$ and 8.5 x $10^{10}$ for a = 5 microm and 50 microm respectively, or an initial mass of $H_2O$ ice around 22kg.

We present observations of SN 2015bn (= PS15ae = CSS141223-113342+004332 = MLS150211-113342+004333), a Type I superluminous supernova (SLSN) at redshift $z=0.1136$. As well as being one of the closest SLSNe I yet discovered, it is intrinsically brighter ($M_U\approx-23.1$) and in a fainter galaxy ($M_B\approx-16.0$) than other SLSNe at $z\sim0.1$. We used this opportunity to collect the most extensive dataset for any SLSN I to date, including densely-sampled spectroscopy and photometry, from the UV to the NIR, spanning $-$50 to +250 days from optical maximum. SN 2015bn fades slowly, but exhibits surprising undulations in the light curve on a timescale of 30-50 days, especially in the UV. The spectrum shows extraordinarily slow evolution except for a rapid transformation between +7 and +20-30 days. No narrow emission lines from slow-moving material are observed at any phase. We derive physical properties including the bolometric luminosity, and find slow velocity evolution and non-monotonic temperature and radial evolution. A deep radio limit rules out a healthy off-axis gamma-ray burst, and places constraints on the pre-explosion mass loss. The data can be consistently explained by a $\gtrsim10\,{\rm M}_\odot$ stripped progenitor exploding with $\sim 10^{51}\,$erg kinetic energy, forming a magnetar with a spin-down timescale of $\sim20$ days (thus avoiding a gamma-ray burst) that reheats the ejecta and drives ionization fronts. The most likely alternative scenario -- interaction with $\sim20\,{\rm M}_\odot$ of dense, inhomogeneous circumstellar material -- can be tested with continuing radio follow-up.

The Wide Angle Camera of the OSIRIS instrument on board the Rosetta spacecraft is equipped with several narrowband filters that are centered on the emission lines and bands of various fragment species. These are used to determine the evolution of the production and spatial distribution of the gas in the inner coma of comet 67P with time and heliocentric distance, here between 2.6 - 1.3 AU pre-perihelion. Our observations indicate that the emission observed in the OH, OI, CN, NH, and NH2 filters is mostly produced by dissociative electron impact excitation of different parent species. We conclude that CO2 rather than H2O is a significant source of the [OI] 630 nm emission. A strong plume-like feature observed in the in CN and [OI] filters is present throughout our observations. This plume is not present in OH emission and indicates a local enhancement of the CO2/H2O ratio by as much as a factor of 3. We observed a sudden decrease in intensity levels after March 2015, which we attribute to decreased electron temperatures in the first kilometers above the nucleus surface.

This paper reviews gravitational wave sources and their detection. One of the most exciting potential sources of gravitational waves are coalescing binary black hole systems. They can occur on all mass scales and be formed in numerous ways, many of which are not understood. They are generally invisible in electromagnetic waves, and they provide opportunities for deep investigation of Einstein's general theory of relativity. Sect. 1 of this paper considers ways that binary black holes can be created in the universe, and includes the prediction that binary black hole coalescence events are likely to be the first gravitational wave sources to be detected. The next parts of this paper address the detection of chirp waveforms from coalescence events in noisy data. Such analysis is computationally intensive. Sect. 2 reviews a new and powerful method of signal detection based on the GPU-implemented summed parallel infinite impulse response filters. Such filters are intrinsically real time alorithms, that can be used to rapidly detect and localise signals. Sect. 3 of the paper reviews the use of GPU processors for rapid searching for gravitational wave bursts that can arise from black hole births and coalescences. In sect. 4 the use of GPU processors to enable fast efficient statistical significance testing of gravitational wave event candidates is reviewed. Sect. 5 of this paper addresses the method of multimessenger astronomy where the discovery of electromagnetic counterparts of gravitational wave events can be used to identify sources, understand their nature and obtain much greater science outcomes from each identified event.

[Background] The "pasta" phase of nuclear matter may play an important role in the structure and evolution of neutron stars. Recent works suggest nuclear pasta has a high resistivity which could be explained by the presence of long lived topological defects. The defects act as impurities that decrease thermal and electrical conductivity of the pasta. [Purpose] To quantify how topological defects affect transport properties of nuclear pasta and estimate this effect using an impurity parameter $Q_{\text{imp}}$. [Methods] Contrast molecular dynamics simulations of up to 409\,600 nucleons arranged in parallel nuclear pasta slabs (perfect pasta) with simulations of pasta slabs connected by topological defects (impure pasta). From these simulations compare the viscosity and heat conductivity of perfect and impure pasta to obtain an effective impurity parameter $Q_{\text{imp}}$ due to the presence of defects. [Results] Both the viscosity and thermal conductivity calculated for both perfect and impure pasta are anisotropic, peaking along directions perpendicular to the slabs and reaching a minimum close to zero parallel to them. In our 409\,600 nucleon simulation topological defects connecting slabs of pasta reduce both the thermal conductivity and viscosity on average by about 37\%. We estimate an effective impurity parameter due to the defects of order $Q_{\text{imp}}\sim30$. [Conclusions] Topological defects in the pasta phase of nuclear matter have an effect similar to impurities in a crystal lattice. The irregularities introduced by the defects reduce the thermal and electrical conductivities and the viscosity of the system. This effect can be parameterized by a large impurity parameter $Q_{\text{imp}}\sim30$.

We analyze the physical properties and dynamical origin of a curved jet of comet 67P/Churyumov-Gerasimenko that was observed repeatedly in several nucleus rotations starting on May 30 and persisting until early August, 2015. We simulated the motion of dust grains ejected from the nucleus surface under the influence of the gravity and viscous drag effect of the expanding gas flow from the rotating nucleus. The formation of the curved jet is a combination of the size of the dust particles (~0.1-1 mm) and the location of the source region near the nucleus equator. This enhances the spiral feature of the collimated dust stream after the dust is accelerated to a terminal speed on the order of m/s.

Dust jets, i.e. fuzzy collimated streams of cometary material arising from the nucleus, have been observed in-situ on all comets since the Giotto mission flew by comet 1P/Halley in 1986. Yet their formation mechanism remains unknown. Several solutions have been proposed, from localized physical mechanisms on the surface/sub-surface (see review in Belton (2010)) to purely dynamical processes involving the focusing of gas flows by the local topography (Crifo et al. 2002). While the latter seems to be responsible for the larger features, high resolution imagery has shown that broad streams are composed of many smaller features (a few meters wide) that connect directly to the nucleus surface. We monitored these jets at high resolution and over several months to understand what are the physical processes driving their formation, and how this affects the surface. Using many images of the same areas with different viewing angles, we performed a 3-dimensional reconstruction of collimated jets, and linked them precisely to their sources on the nucleus. Results.We show here observational evidence that the Northern hemisphere jets of comet 67P arise from areas with sharp topographic changes and describe the physical processes involved. We propose a model in which active cliffs are the main source of jet-like features, and therefore the regions eroding the fastest on comets. We suggest that this is a common mechanism taking place on all comets.

In this paper we investigate the opportunities provided by the James Webb Space Telescope (JWST) for significant scientific advances in the study of solar system bodies and rings using stellar occultations. The strengths and weaknesses of the stellar occultation technique are evaluated in light of JWST's unique capabilities. We identify several possible JWST occultation events by minor bodies and rings, and evaluate their potential scientific value. These predictions depend critically on accurate a priori knowledge of the orbit of JWST near the Sun-Earth Lagrange-point 2 (L2). We also explore the possibility of serendipitous stellar occultations by very small minor bodies as a by-product of other JWST observing programs. Finally, to optimize the potential scientific return of stellar occultation observations, we identify several characteristics of JWST's orbit and instrumentation that should be taken into account during JWST's development.

The ultra-long period Cepheids (ULPCs) are classical Cepheids with pulsation periods exceeding $\approx 80$ days. The intrinsic brightness of ULPCs are ~1 to ~3 mag brighter than their shorter period counterparts. This makes them attractive in future distance scale work to derive distances beyond the limit set by the shorter period Cepheids. We have initiated a program to search for ULPCs in M31, using the single-band data taken from the Palomar Transient Factory, and identified eight possible candidates. In this work, we presented the VI-band follow-up observations of these eight candidates. Based on our VI-band light curves of these candidates and their locations in the color-magnitude diagram and the Period-Wesenheit diagram, we verify two candidates as being truly ULPCs. The six other candidates are most likely other kinds of long-period variables. With the two confirmed M31 ULPCs, we tested the applicability of ULPCs in distance scale work by deriving the distance modulus of M31. It was found to be $\mu_{M31,ULPC}=24.30\pm0.76$ mag. The large error in the derived distance modulus, together with the large intrinsic dispersion of the Period-Wesenheit (PW) relation and the small number of ULPCs in a given host galaxy, means that the question of the suitability of ULPCs as standard candles is still open. Further work is needed to enlarge the sample of calibrating ULPCs and reduce the intrinsic dispersion of the PW relation before re-considering ULPCs as suitable distance indicators.

Low-inclination Near-Earth Asteroid (297274) 1996 SK, which is also classified as a Potentially Hazardous Asteroid, has a highly eccentric orbit. It was studied by multi-wavelength photometry within the framework of an NEA color survey at the Lulin Observatory. We report here the finding of large color variation across the surface of (297274) 1996 SK within one asteroidal rotation period of $4.656\pm0.122$ hours and classify it as a S-type asteroid according to its average colors of $B-V=0.767\pm0.033$, $V-R=0.482\pm0.021$, $V-I=0.801\pm0.025$ and the corresponding relative reflectance spectrum. It might be indicative of differential space weathering effect or compositional inhomogeneity of the surface materials.

We reported the monitoring results on spectrophotometry, photometry and imaging of comet 103P/Hartley 2 obtained at Lulin (1m), Calar Alto (2.2m) and Beijing Astronomical (2.16m) Observatory from April to December 2010. We found that a dust feature at sunward direction was detected starting from the end of September until the beginning of December (our last observation from the Lulin and Calar Alto observatory). Two distinct sunward jet features in the processed images were observed on October 11 and after October 29 until November 2. In parallel, the CN images reveal two asymmetrical jet features which are nearly perpendicular to the Sun-nucleus direction and this asymmetrical features implies that the comet was in a nearly side-on view in late-October and early-November. Additional to the jet features, the average result of the C2-to-CN production rate ratio ranges from 0.7 to 1.5 which places 103P/Hartley 2 as being of typical cometary chemistry. We found that the $r_h$ dependence for the dust production rate, Afrho (5,000 km), is -3.75+-0.45 before perihelion and is -3.44 +-1.20 during post-perihelion period. We detected the higher dust reddening is around the optocenter and getting bluer outward along the sunward jet feature and concluded that the former one, higher dust reddening, could be associated with strong jet activity and the latter one, the lowering of the reddening, might imply the optical properties changed or could be associated with outburst. The average dust color did not appear to vary significantly as the comet passed through perihelion.

The dynamics of collisionless galaxy can be described by the Vlasov-Poisson system. By the Jean's theorem, all the spherically symmetric steady galaxy models are given by a distribution of {\Phi}(E,L), where E is the particle energy and L the angular momentum. In a celebrated Doremus-Feix-Baumann Theorem, the galaxy model {\Phi}(E,L) is stable if the distribution {\Phi} is monotonically decreasing with respect to the particle energy E. On the other hand, the stability of {\Phi}(E,L) remains largely open otherwise. Based on a recent abstract instability criterion of Guo-Lin, we constuct examples of unstable galaxy models of f(E,L) and f(E) in which f fails to be monotone in E.

To determine the stability and instability of a given steady galaxy configuration is one of the fundamental problems in the Vlasov theory for galaxy dynamics. In this article, we study the stability of isotropic spherical symmetric galaxy models $f_{0}(E)$, for which the distribution function $f_{0}$ depends on the particle energy $E$ only. In the first part of the article, we derive the first sufficient criterion for linear instability of $f_{0}(E):$ $f_{0}(E)$ is linearly unstable if the second-order operator \[ A_{0}\equiv-\Delta+4\pi\int f_{0}^{\prime}(E)\{I-\mathcal{P}\}dv \] has a negative direction, where $\mathcal{P}$ is the projection onto the function space $\{g(E,L)\},$ $L$ being the angular momentum [see the explicit formula (\ref{A0-radial})]. In the second part of the article, we prove that for the important King model, the corresponding $A_{0}$ is positive definite. Such a positivity leads to the nonlinear stability of the King model under all spherically symmetric perturbations.

The 1-m telescope at Lulin Observatory and the 0.76-m Katzman Automatic Imaging Telescope at Lick Observatory were used to observe the optical afterglow of the short-duration (1.2--1.5 s) gamma-ray burst (GRB) 040924. This object has a soft high-energy spectrum, thus making it an exceptional case, perhaps actually belonging to the short-duration tail of the long-duration GRBs. Our data, combined with other reported measurements, show that the early R-band light curve can be described by two power laws with index alpha = -0.7 (at t = 16-50 min) and alpha = -1.06 (at later times). The rather small difference in the spectral indices can be more easily explained by an afterglow model invoking a cooling break rather than a jet break.

We study in this paper the possible influence of the dark energy on the detection of the dark matter particles. In models of dark energy described by a dynamical scalar field such as the Quintessence, its interaction with the dark matter will cause the dark matter particles such as the neutralino vary as a function of space and time. Given a specific model of the Quintessence and its interaction in this paper we calculate numerically the corrections to the neutralino masses and the induced spectrum of the neutrinos from the annihilation of the neutralinos pairs in the core of the Sun. This study gives rise to a possibility of probing for dark energy in the experiments of detecting the dark matter particles.