A TeV spectral break in the total flux of cosmic-ray electrons and positrons (CREs) at which the spectral power index softens from $\sim 3$ to $\sim 4$ has been observed by H.E.S.S. and recently confirmed by DAMPE with a high significance of $6.6~\sigma$. Such an observation is apparently inconsistent with the data from other experiments such as Fermi-LAT, AMS-02 and CALET. We perform a global analysis to the latest CRE data including Fermi-LAT, AMS-02, CALET, DAMPE and H.E.S.S. with energy scale uncertainties taken into account to improve the consistency between the data sets. The fit result strongly favors the existence of the break at $\sim 1$ TeV with an even higher statistical significance of $13.3~\sigma$. In view of the tentative CRE break, we revisit a number of models of nearby sources, such as a single generic Pulsar Wind Nebula (PWN), known multiple PWNe from the ATNF catalog, and their combinations with either an additional Dark Matter (DM) component or a Supernova Remnant (SNR). We show that the CRE break at $\sim 1$ TeV, together with the known CR positron excess points towards the possibility that the nearby sources should be highly charge asymmetric. Among the models under consideration, the one with a PWN plus SNR is most favored by the current data. The favoured distance and age of the PWN and SNR sources are both within $0.6$ kpc and around $10^{5}$ yr respectively. Possible candidate sources include PSR J0954-5430, Vela and Monogem ring, etc. We find that for the models under consideration, the additional DM component is either unnecessary, or predicts too much photons in tension with the H.E.S.S. data of $\gamma$-rays from the direction of the Galactic Center. We also show that the current measurement of the anisotropies in the arrival direction of the CRE can be useful in determining the property of the sources.

A relationship between dust-reprocessed light from recent star formation and the amount of star-forming gas in a galaxy produces a correlation between WISE 12$\mu$m emission and CO line emission. Here we explore this correlation on kiloparsec scales with CO(1-0) maps from EDGE-CALIFA matched in resolution to WISE 12$\mu$m images. We find strong CO-12$\mu$m correlations within each galaxy (median Pearson $r = 0.85$) and we show that the scatter in the global CO-12$\mu$m correlation is largely driven by differences from galaxy to galaxy. The correlation is stronger than that between star formation rate and H$_2$ surface densities over the same set of pixels (median $r=0.71$). We explore multi-variable regression to predict $\Sigma(\mathrm{H_2})$ surface density using the WISE 12$\mu$m data combined with global and resolved galaxy properties, and provide the fit parameters for the best estimators. We find that $\Sigma(\mathrm{H_2})$ estimators that include $\Sigma(\mathrm{12\mu m})$ are able to predict $\Sigma(\mathrm{H_2})$ with $>10$% better accuracy than estimators that include resolved optical properties ($\Sigma(\mathrm{SFR}),\Sigma(\mathrm{M_*}),A_V$ and $12+\log\mathrm{O/H}$) instead of $\Sigma(\mathrm{12\mu m})$. The best single-property estimator is $\log \frac{\Sigma(\mathrm{H_2})}{\mathrm{M_\odot\>pc^{-2}}} = (0.49 \pm 0.01) + (0.71 \pm 0.01)\log \frac{\Sigma(\mathrm{12\mu m})}{\mathrm{L_\odot\>pc^{-2}}}$, with an average predictive accuracy of 0.19 dex per pixel, and intrinsic scatter of 0.17 dex. This correlation can be used to efficiently estimate H$_2$ surface densities down to at least $1 \> M_\odot \> \mathrm{pc^{-2}}$ on small spatial scales within nearby galaxies. This correlation may prove useful to probe even lower gas densities with the better mid-infrared sensitivities expected from the James Webb Space Telescope.

Upcoming large astronomical surveys are expected to capture an unprecedented number of strong gravitational lensing systems in the Universe. Deep learning is emerging as a promising practical tool in detection and quantification of these galaxy-scale image distortions. However, absence of large quantities of representative data from current astronomical surveys requires development of robust forward modeling of synthetic lensing images. Using a realistic and unbiased sample of the strong lenses created by using state-of-the-art extragalactic catalogs, we train a modular deep learning pipeline for uncertainty-quantified detection and modeling with intermediate image processing components for denoising and deblending the lensing systems. We demonstrate a higher degree of interpretability and controlled systematics due to domain-specific task modules that are trained with different stages of synthetic image generation. For lens detection and modeling, we obtain semantically meaningful latent spaces that separate classes and provide uncertainty estimates that explain the misclassified images and provide uncertainty bounds on the lens parameters. In addition, we obtain an improved performance---lens detection (classification) improved from 82% with the baseline to 94%, while the lens modeling (regression) accuracy improved by 25% over the baseline model.

Measurements of The Hubble-Lemaitre constant from early- and local-universe observations show a significant discrepancy. In an attempt to understand the origin of this mismatch, independent techniques to measure $H_0$ are required. One such technique, strong lensing time delays, is set to become a leading contender amongst the myriad methods due to forthcoming large strong lens samples. It is therefore critical to understand the systematic effects inherent in this method. In this paper, we quantify the influence of additional structures along the line-of-sight by adopting realistic lightcones derived from the \textit{CosmoDC2} semi-analytical extra-galactic catalogue. Using multiple lens plane ray-tracing to create a set of simulated strong lensing systems, we have investigated the impact of line-of-sight structures on time-delay measurements and in turn, on the inferred value of $H_0$. We have also tested the reliability of existing procedures for correcting for line-of-sight effects. We find that if the integrated contribution of the of line-of-sight structures is close to a uniform mass sheet, the bias in $H_0$ can be adequately corrected by including a constant external convergence $\kappa_{\rm ext}$ in the lens model. However, for realistic line-of-sight structures comprising many galaxies at different redshifts, this simple correction over-estimates the bias by a factor of approximately three. We therefore conclude that lens modelling must incorporate multiple lens planes to account for line-of-sight structures for accurate and precise inference of $H_0$.

Wolf-Rayet (WR) galaxies are a rare population of galaxies that host living high-mass stars during their WR phase (i.e. WR stars) and are thus expected to provide interesting constraints on the stellar Initial Mass Function, massive star formation, stellar evolution models, etc. Spatially resolved spectroscopy should in principle provide a more efficient way of identifying WR galaxies than single-fiber surveys of galactic centers such as SDSS-I & II, as WR stars should be more preferentially found in discs. Using IFU data from the ongoing SDSS-IV MaNGA survey, we have performed a thorough search for WR galaxies. We first identify H II regions in each datacube and carry out full spectral fitting to the stacked spectra. We then visually inspect the residual spectrum of each H II region and identify WR regions that present a significant "blue bump" at 4600-4750 A. The resulting WR catalog includes 267 WR regions of ~500pc (radius) sizes, distributed in 90 galaxies from the current sample of MaNGA (MaNGA Product Launch 7). We find WR regions are exclusively found in galaxies that show bluest colors and highest star formation rates for their mass. Most WR galaxies have late-type morphologies and show relatively large asymmetry in their images, implying that WR regions are more preferentially found in interacting/merging galaxies. We estimate the stellar mass function of WR galaxies and the mass-dependent detection rate. The detection rate of WR galaxies is typically ~2%, with weak dependence on stellar mass. This detection rate is about 40 times higher than previous studies with SDSS single fiber data, and by a factor of 2 lower than the CALIFA-based WR catalog. We make comparisons with SDSS and CALIFA studies, and conclude that different detection rates can be explained mainly by three factors: spatial coverage, spectral signal-to-noise ratio, and redshift ranges of the parent sample.

We develop a method to estimate the dust attenuation curve of galaxies from full spectral fitting of their optical spectra. Motivated from previous studies, we separate the small-scale features from the large-scale spectral shape, by performing a moving average method to both the observed spectrum and the simple stellar population model spectra. The intrinsic dust-free model spectrum is then derived by fitting the observed ratio of the small-scale to large-scale (S/L) components with the S/L ratios of the SSP models. The selective dust attenuation curve is then determined by comparing the observed spectrum with the dust-free model spectrum. One important advantage of this method is that the estimated dust attenuation curve is independent of the shape of theoretical dust attenuation curves. We have done a series of tests on a set of mock spectra covering wide ranges of stellar age and metallicity. We show that our method is able to recover the input dust attenuation curve accurately, although the accuracy depends slightly on signal-to-noise ratio of the spectra. We have applied our method to a number of edge-on galaxies with obvious dust lanes from the ongoing MaNGA survey, deriving their dust attenuation curves and $E(B-V)$ maps, as well as dust-free images in $g$, $r$, and $i$ bands. These galaxies show obvious dust lane features in their original images, which largely disappear after we have corrected the effect of dust attenuation. The vertical brightness profiles of these galaxies become axis-symmetric and can well be fitted by a simple model proposed for the disk vertical structure. Comparing the estimated dust attenuation curve with the three commonly-adopted model curves, we find that the Calzetti curve provides the best description of the estimated curves for the inner region of galaxies, while the Milky Way and SMC curves work better for the outer region.

In this paper we develop a new unsupervised machine learning technique comprised of a feature extractor, a convolutional autoencoder (CAE), and a clustering algorithm consisting of a Bayesian Gaussian mixture model (BGM). We apply this technique to visual band space-based simulated imaging data from the Euclid Space Telescope using data from the Strong Gravitational Lenses Finding Challenge. Our technique promisingly captures a variety of lensing features such as Einstein rings with different radii, distorted arc structures, etc, without using predefined labels. After the clustering process, we obtain several classification clusters separated by different visual features which are seen in the images. Our method successfully picks up $\sim$63\ percent of lensing images from all lenses in the training set. With the assumed probability proposed in this study, this technique reaches an accuracy of $77.25\pm 0.48$\% in binary classification using the training set. Additionally, our unsupervised clustering process can be used as the preliminary classification for future surveys of lenses to efficiently select targets and to speed up the labelling process. As the starting point of the astronomical application using this technique, we not only explore the application to gravitationally lensed systems, but also discuss the limitations and potential future uses of this technique.

There are several supervised machine learning methods used for the application of automated morphological classification of galaxies; however, there has not yet been a clear comparison of these different methods using imaging data, or a investigation for maximising their effectiveness. We carry out a comparison between several common machine learning methods for galaxy classification (Convolutional Neural Network (CNN), K-nearest neighbour, Logistic Regression, Support Vector Machine, Random Forest, and Neural Networks) by using Dark Energy Survey (DES) data combined with visual classifications from the Galaxy Zoo 1 project (GZ1). Our goal is to determine the optimal machine learning methods when using imaging data for galaxy classification. We show that CNN is the most successful method of these ten methods in our study. Using a sample of $\sim$2,800 galaxies with visual classification from GZ1, we reach an accuracy of $\sim$0.99 for the morphological classification of Ellipticals and Spirals. The further investigation of the galaxies that have a different ML and visual classification but with high predicted probabilities in our CNN usually reveals an the incorrect classification provided by GZ1. We further find the galaxies having a low probability of being either spirals or ellipticals are visually Lenticulars (S0), demonstrating that supervised learning is able to rediscover that this class of galaxy is distinct from both Es and Spirals. We confirm that $\sim$2.5\% galaxies are misclassified by GZ1 in our study. After correcting these galaxies' labels, we improve our CNN performance to an average accuracy of over 0.99 (accuracy of 0.994 is our best result).

The precise measurement of cosmic-ray antinuclei serves as an important means for identifying the nature of dark matter and other new astrophysical phenomena, and could be used with other cosmic-ray species to understand cosmic-ray production and propagation in the Galaxy. For instance, low-energy antideuterons would provide a "smoking gun" signature of dark matter annihilation or decay, essentially free of astrophysical background. Studies in recent years have emphasized that models for cosmic-ray antideuterons must be considered together with the abundant cosmic antiprotons and any potential observation of antihelium. Therefore, a second dedicated Antideuteron Workshop was organized at UCLA in March 2019, bringing together a community of theorists and experimentalists to review the status of current observations of cosmic-ray antinuclei, the theoretical work towards understanding these signatures, and the potential of upcoming measurements to illuminate ongoing controversies. This review aims to synthesize this recent work and present implications for the upcoming decade of antinuclei observations and searches. This includes discussion of a possible dark matter signature in the AMS-02 antiproton spectrum, the most recent limits from BESS Polar-II on the cosmic antideuteron flux, and reports of candidate antihelium events by AMS-02; recent collider and cosmic-ray measurements relevant for antinuclei production models; the state of cosmic-ray transport models in light of AMS-02 and Voyager data; and the prospects for upcoming experiments, such as GAPS. This provides a roadmap for progress on cosmic antinuclei signatures of dark matter in the coming years.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is acquiring integral-field spectroscopy for the largest sample of galaxies to date. By 2020, the MaNGA Survey --- one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) --- will have observed a statistically representative sample of 10$^4$ galaxies in the local Universe ($z\lesssim0.15$). In addition to a robust data-reduction pipeline (DRP), MaNGA has developed a data-analysis pipeline (DAP) that provides higher-level data products. To accompany the first public release of its code base and data products, we provide an overview of the MaNGA DAP, including its software design, workflow, measurement procedures and algorithms, performance, and output data model. In conjunction with our companion paper Belfiore et al., we also assess the DAP output provided for 4718 observations of 4648 unique galaxies in the recent SDSS Data Release 15 (DR15). These analysis products focus on measurements that are close to the data and require minimal model-based assumptions. Namely, we provide stellar kinematics (velocity and velocity dispersion), emission-line properties (kinematics, fluxes, and equivalent widths), and spectral indices (e.g., D4000 and the Lick indices). We find that the DAP provides robust measurements and errors for the vast majority ($>$99%) of analyzed spectra. We summarize assessments of the precision and accuracy of our measurements as a function of signal-to-noise, and provide specific guidance to users regarding the limitations of the data. The MaNGA DAP software is publicly available and we encourage community involvement in its development.

Forthcoming surveys such as the Large Synoptic Survey Telescope (LSST) and Euclid necessitate automatic and efficient identification methods of strong lensing systems. We present a strong lensing identification approach that utilizes a feature extraction method from computer vision, the Histogram of Oriented Gradients (HOG), to capture edge patterns of arcs. We train a supervised classifier model on the HOG of mock strong galaxy-galaxy lens images similar to observations from the Hubble Space Telescope (HST) and LSST. We assess model performance with the area under the curve (AUC) of a Receiver Operating Characteristic (ROC) curve. Models trained on 10,000 lens and non-lens containing images images exhibit an AUC of 0.975 for an HST-like sample, 0.625 for one exposure of LSST, and 0.809 for 10-year mock LSST observations. Performance appears to continually improve with the training set size. Models trained on fewer images perform better in absence of the lens galaxy light. However, with larger training data sets, information from the lens galaxy actually improves model performance, indicating that HOG captures much of the morphological complexity of the arc finding problem. We test our classifier on data from the Sloan Lens ACS Survey and find that small scale image features reduces the efficiency of our trained model. However, these preliminary tests indicate that some parameterizations of HOG can compensate for differences between observed mock data. One example best-case parameterization results in an AUC of 0.6 in the F814 filter image with other parameterization results equivalent to random performance.

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

We explore the effectiveness of deep learning convolutional neural networks (CNNs) for estimating strong gravitational lens mass model parameters. We have investigated a number of practicalities faced when modelling real image data, such as how network performance depends on the inclusion of lens galaxy light, the addition of colour information and varying signal-to-noise. Our CNN was trained and tested with strong galaxy-galaxy lens images simulated to match the imaging characteristics of the Large Synoptic Survey Telescope (LSST) and Euclid. For images including lens galaxy light, the CNN can recover the lens model parameters with an acceptable accuracy, although a 34 per cent average improvement in accuracy is obtained when lens light is removed. However, the inclusion of colour information can largely compensate for the drop in accuracy resulting from the presence of lens light. While our findings show similar accuracies for single epoch Euclid VIS and LSST r-band datasets, we find a 24 per cent increase in accuracy by adding g- and i-band images to the LSST r-band without lens light and a 20 per cent increase with lens light. The best network performance is obtained when it is trained and tested on images where lens light exactly follows the mass, but when orientation and ellipticity of the light is allowed to differ from those of the mass, the network performs most consistently when trained with a moderate amount of scatter in the difference between the mass and light profiles.

Cosmological cluster-scale strong gravitational lensing probes the mass distribution of the dense cores of massive dark matter halos and the structures along the line of sight from background sources to the observer. It is frequently assumed that the primary lens mass dominates the lensing, with the contribution of secondary masses along the line of sight being neglected. Secondary mass structures may, however, affect both the detectability of strong lensing in a given survey and modify the properties of the lensing that is detected. In this paper, we utilize a large cosmological N-body simulation and a multiple lens plane (and many source planes) ray-tracing technique to quantify the influence of line of sight halos on the detectability of cluster-scale strong lensing in a cluster sample with a mass limit that encompasses current cluster catalogs from the South Pole Telescope. We extract both primary and secondary halos from the "Outer Rim" simulation and consider two strong lensing realizations: one with only the primary halos included, and the other contains all secondary halos down to a mass limit. In both cases, we use the same source information extracted from the Hubble Ultra Deep Field, and create realistic lensed images consistent with moderately deep ground-based imaging. The results demonstrate that down to the mass limit considered the total number of lenses is boosted by about 13-21% when considering the complete multi-halo lightcone. The increment in strong lens counts peaks at lens redshifts of 0.6 approximately with no significant effect at z<0.3. The strongest trends are observed relative to the primary halo mass, with no significant impact in the most massive quintile of the halo sample, but increasingly boosting the observed lens counts toward small primary halo masses, with an enhancement greater than 50% in the least massive quintile of the halo masses considered.

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged)

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 study the spatially resolved star formation history and molecular gas distribution of 58 nearby galaxies, using integral field spectroscopy from the CALIFA survey and CO $J=1\rightarrow 0$ intensity mapping from the CARMA EDGE survey. We use the 4000 \AA\ break (D$_n$4000), the equivalent width of the H$\delta$ absorption line (EW H$\delta_A$), and the equivalent width of the H$\alpha$ emission line (EW H$\alpha$) to measure the recent star formation history (SFH) of these galaxies. We measure radial profiles of the three SFH indicators and molecular gas mass surface density, from which we measure the level of centrally enhanced star formation and the molecular gas concentration. When we separate our galaxies into categories of barred (17 galaxies), unbarred (24 galaxies), and merging/paired (17 galaxies) we find that the galaxies which have centrally-enhanced star formation (19/58) are either barred (13/19) or in mergers/pairs (6/19) with relatively high molecular gas concentrations. A comparison between our barred galaxies and a snapshot of a hydrodynamic $N$-body simulation of a barred galaxy shows that the current theory of bar formation and evolution can qualitatively reproduce the main features of the observed galaxies in our sample, including both the sharp decrease of stellar age in the galactic center and the gradual decrease of age with increasing distance from center. These findings provide substantial evidence for a picture in which cold gas is transported inward by a bar or tidal interaction, which leads to the growth and rejuvenation of star formation in the central region.

As the second paper of a series on studying galaxy-galaxy lensing signals using the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our measurement and modelling of the lensing signals around groups of galaxies. We divide the groups into four halo mass bins, and measure the signals around four different halo-center tracers: brightest central galaxy (BCG), luminosity-weighted center, number-weighted center and X-ray peak position. For X-ray and SDSS DR7 cross identified groups, we further split the groups into low and high X-ray emission subsamples, both of which are assigned with two halo-center tracers, BCGs and X-ray peak positions. The galaxy-galaxy lensing signals show that BCGs, among the four candidates, are the best halo-center tracers. We model the lensing signals using a combination of four contributions: off-centered NFW host halo profile, sub-halo contribution, stellar contribution, and projected 2-halo term. We sample the posterior of 5 parameters i.e., halo mass, concentration, off-centering distance, sub halo mass, and fraction of subhalos via a MCMC package using the galaxy-galaxy lensing signals. After taking into account the sampling effects (e.g. Eddington bias), we found the best fit halo masses obtained from lensing signals are quite consistent with those obtained in the group catalog based on an abundance matching method, except in the lowest mass bin. Subject headings: (cosmology:) gravitational lensing, galaxies: clusters: general

Nuclear clustering describes the appearance of structures resembling smaller nuclei such as alpha particles (4He nuclei) within the interior of a larger nucleus. While clustering is important for several well-known examples, much remains to be discovered about the general nature of clustering in nuclei. In this letter we present lattice Monte Carlo calculations based on chiral effective field theory for the ground states of helium, beryllium, carbon, and oxygen isotopes. By computing model-independent measures that probe three- and four-nucleon correlations at short distances, we determine the shape of the alpha clusters and the entanglement of nucleons comprising each alpha cluster with the outside medium. We also introduce a new computational approach called the pinhole algorithm, which solves a long-standing deficiency of auxiliary-field Monte Carlo simulations in computing density correlations relative to the center of mass. We use the pinhole algorithm to determine the proton and neutron density distributions and the geometry of cluster correlations in 12C, 14C, and 16C. The structural similarities among the carbon isotopes suggest that 14C and 16C have excitations analogous to the well-known Hoyle state resonance in 12C.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

A few dark matter substructures have recently been detected in strong gravitational lenses though their perturbations of highly magnified images. We derive a characteristic scale for lensing perturbations and show that this is significantly larger than the perturber's Einstein radius. We show that the perturber's projected mass enclosed within this radius, scaled by the log-slope of the host galaxy's density profile, can be robustly inferred even if the inferred density profile and tidal radius of the perturber are biased. We demonstrate the validity of our analytic derivation by using several gravitational lens simulations where the tidal radii and the inner log-slopes of the density profile of the perturbing subhalo are allowed to vary. By modeling these simulated data we find that our mass estimator, which we call the effective subhalo lensing mass, is accurate to within about 10\% or smaller in each case, whereas the inferred total subhalo mass can potentially be biased by nearly an order of magnitude. We therefore recommend that the effective subhalo lensing mass be reported in future lensing reconstructions, as this will allow for a more accurate comparison with the results of dark matter simulations.

Galaxy-scale strong gravitational lensing is not only a valuable probe of the dark matter distribution of massive galaxies, but can also provide valuable cosmological constraints, either by studying the population of strong lenses or by measuring time delays in lensed quasars. Due to the rarity of galaxy-scale strongly lensed systems, fast and reliable automated lens finding methods will be essential in the era of large surveys such as LSST, Euclid, and WFIRST. To tackle this challenge, we introduce CMU DeepLens, a new fully automated galaxy-galaxy lens finding method based on Deep Learning. This supervised machine learning approach does not require any tuning after the training step which only requires realistic image simulations of strongly lensed systems. We train and validate our model on a set of 20,000 LSST-like mock observations including a range of lensed systems of various sizes and signal-to-noise ratios (S/N). We find on our simulated data set that for a rejection rate of non-lenses of 99%, a completeness of 90% can be achieved for lenses with Einstein radii larger than 1.4" and S/N larger than 20 on individual $g$-band LSST exposures. Finally, we emphasize the importance of realistically complex simulations for training such machine learning methods by demonstrating that the performance of models of significantly different complexities cannot be distinguished on simpler simulations. We make our code publicly available at https://github.com/McWilliamsCenter/CMUDeepLens .

Sometimes the early star formation can be found in cold and dense molecular clouds, such as infrared dark cloud (IRDC). Considering star formation often occurs in clustered condition, HII regions may be triggering a new generation of star formation, so we can search for initial stage of massive star formation around HII regions. Based on that above, this work is to introduce one method of how to search for initial stage of massive star formation around HII regions. Towards one sample of the HII region G18.2-0.3, multiwavelength observations are carried out to investigate its physical condition. In contrast and analysis, we find three potential initial stages of massive star formation, suggesting that it is feasible to search for initial stage of massive star formation around HII regions.

Gravitational lensing has become one of the most powerful tools available for investigating the 'dark side' of the universe. Cosmological strong gravitational lensing, in particular, probes the properties of the dense cores of dark matter halos over decades in mass and offers the opportunity to study the distant universe at flux levels and spatial resolutions otherwise unavailable. Studies of strongly-lensed variable sources offer yet further scientific opportunities. One of the challenges in realizing the potential of strong lensing is to understand the statistical context of both the individual systems that receive extensive follow-up study, as well as that of the larger samples of strong lenses that are now emerging from survey efforts. Motivated by these challenges, we have developed an image-simulation pipeline, PICS (Pipeline for Images of Cosmological Strong lensing) to generate realistic strong gravitational lensing signals from group and cluster scale lenses. PICS uses a low-noise and unbiased density estimator based on (resampled) Delaunay Tessellations to calculate the density field; lensed images are produced by ray-tracing images of actual galaxies from deep Hubble Space Telescope observations. Other galaxies, similarly sampled, are added to fill in the light cone. The pipeline further adds cluster-member galaxies and foreground stars into the lensed images. The entire image ensemble is then observed using a realistic point spread function which includes appropriate detector artifacts for bright stars. Noise is further added, including such non-Gaussian elements as noise window-paning from mosaiced observations, residual bad pixels, and cosmic rays. The aim is to produced simulated images that appear identical - to the eye (expert or otherwise) - to real observations in various imaging surveys.

Dark energy is investigated from the perspective of quantum cosmology. It is found that, together with an appropriate normal ordering factor $q$, only when there is dark energy then can the cosmological wave function be normalized. This interesting observation may require further attention.

As the first paper in a series on the study of the galaxy-galaxy lensing from Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our image processing pipeline that corrects the systematics primarily introduced by the Point Spread Function (PSF). Using this pipeline, we processed SDSS DR7 imaging data in $r$ band and generated a background galaxy catalog containing the shape information of each galaxy. Based on our own shape measurements of the galaxy images from SDSS DR7, we extract the galaxy-galaxy (GG) lensing signals around foreground spectroscopic galaxies binned in different luminosity and stellar mass. The overall signals are in good agreement with those obtained by \citet{Mandelbaum2005, Mandelbaum2006} from the SDSS DR4. The results in this paper with higher signal to noise ratio is due to the larger survey area than SDSS DR4, confirm that more luminous/massive galaxies bear stronger GG lensing signal. We also divide the foreground galaxies into red/blue and star forming/quenched subsamples and measured their GG lensing signals, respectively. We find that, at a specific stellar mass/luminosity, the red/quenched galaxies have relatively stronger GG lensing signals than their counterparts especially at large radii. These GG lensing signals can be used to probe the galaxy-halo mass relations and their environmental dependences in the halo occupation or conditional luminosity function framework.

Characterization of the morphology of strongly lensed galaxies is challenging because images of such galaxies are typically highly distorted. Lens modeling and source plane reconstruction is one approach that can provide reasonably undistorted images from which morphological measurements can be made, although at the expense of a highly spatially variable telescope PSF when mapped back to the source plane. Unfortunately, modeling the lensing mass is a time and resource intensive process, and in many cases there are too few constraints to precisely model the lensing mass. If, however, useful morphological measurements could be made in the image plane rather than the source plane, it would bypass this issue and obviate the need for a source reconstruction process for some applications. We examine the use of the Gini coefficient as one such measurement. Because it depends on the cumulative distribution of the light of a galaxy, but not the relative spatial positions, the fact that surface brightness is conserved by lensing means that the Gini coefficient may be well-preserved by strong gravitational lensing. Through simulations, we test the extent to which the Gini coefficient is conserved, including by effects due to PSF convolution and pixelization, to determine whether it is invariant enough under lensing to be used as a measurement of galaxy morphology that can be made in the image plane.

We explore the cosmological consequences of interacting dark energy (IDE) models using the SNLS3 supernova samples. In particular, we focus on the impacts of different SNLS3 light-curve fitters (LCF) (corresponding to "SALT2", "SiFTO", and "Combined" sample). Firstly, making use of the three SNLS3 data sets, as well as the Planck distance priors data and the galaxy clustering data, we constrain the parameter spaces of three IDE models. Then, we study the cosmic evolutions of Hubble parameter $H(z)$, deceleration diagram $q(z)$, statefinder hierarchy $S^{(1)}_3(z)$ and $S^{(1)}_4(z)$, and check whether or not these dark energy diagnosis can distinguish the differences among the results of different SNLS3 LCF. At last, we perform high redshift cosmic age test using three old high redshift objects (OHRO), and explore the fate of the Universe. We find that, the impacts of different SNLS3 LCF are rather small, and can not be distinguished by using $H(z)$, $q(z)$, $S^{(1)}_3(z)$, $S^{(1)}_4(z)$, and the age data of OHRO. In addition, we infer, from the current observations, how far we are from a cosmic doomsday in the worst case, and find that the "Combined" sample always gives the largest 2$\sigma$ lower limit of the time interval between "big rip" and today, while the results given by the "SALT2" and the "SiFTO" sample are close to each other. These conclusions are insensitive to a specific form of dark sector interaction. Our method can be used to distinguish the differences among various cosmological observations.

To investigate the environment of HII region Sh2-163 and search for evidence of triggered star formation in this region, we performed a multi-wavelength study of this HII region. Most of our data were taken from large-scale surveys: 2MASS, CGPS, MSX and SCUBA. We also made CO molecular line observations, using the 13.7-m telescope. The ionized region of Sh2-163 is detected by both the optical and radio continuum observations. Sh2-163 is partially bordered by an arc-like photodissociation region (PDR), which is coincident with the strongest optical and radio emissions, indicating interactions between the HII region and the surrounding interstellar medium. Two molecular clouds were discovered on the border of the PDR. The morphology of these two clouds suggests they are compressed by the expansion of Sh2-163. In cloud A, we found two molecular clumps. And it seems star formation in clump A2 is much more active than in clump A1. In cloud B, we found new outflow activities and massive star(s) are forming inside. Using 2MASS photometry, we tried to search for embedded young stellar object (YSO) candidates in this region. The very good agreement between CO emission, infrared shell and YSOs suggest that it is probably a star formation region triggered by the expansion of Sh2-163.We also found the most likely massive protostar related to IRAS 23314+6033.

We have studied the physical and chemical properties of 18 southern Red Midcourse Space Experiment Sources (RMSs), using archival data taken from the Atacama Pathfinder Experiment (APEX) Telescope Large Area Survey of the Galaxy, the Australia Telescope Compact Array, and the Millimeter Astronomy Legacy Team Survey at 90 GHz. Most of our sources have simple cometary/unresolved radio emissions at 4.8 and/or 8.6GHz. The large number of Lyman continuum fluxes (NL) indicates they are probably massive O- or early B-type star formation regions. Archival IRAS infrared data are used to estimate the dust temperature, which is about 30 K of our sources. Then, the H2 column densities and the volume-averaged H2 number densities are estimated using the 0.87 mm dust emissions. Large-scale infall and ionized accretions may be occurring in G345.4881+00.3148. We also attempt to characterize the chemical properties of these RMSs through molecular line (N2H+ (1-0) and HCO+ (1-0)) observations. Most of the detected N2H+ and HCO+ emissions match well with the dust emission, implying a close link to their chemical evolution in the RMSs. We found that the abundance of N2H+ is one order of magnitude lower than that in other surveys of infrared dark clouds, and a positive correlation between the abundances of N2H+ and HCO+. The fractional abundance of N2H+ with respect to H2 seems to decrease as a function of NL. These observed trends could be interpreted as an indication of enhanced destruction of N2H+, either by CO or through dissociative recombination with electrons produced by central UV photons.

In this work, by applying the redshift tomography method to Joint Light-curve Analysis (JLA) supernova sample, we explore the possible redshift-dependence of stretch-luminosity parameter $\alpha$ and color-luminosity parameter $\beta$. The basic idea is to divide the JLA sample into different redshift bins, assuming that $\alpha$ and $\beta$ are piecewise constants. Then, by constraining the $\Lambda$CDM model, we check the consistency of cosmology-fit results given by the SN sample of each redshift bin. We also adopt the same technique to explore the possible evolution of $\beta$ in various subsamples of JLA. Using the full JLA data, we find that $\alpha$ is always consistent with a constant. In contrast, at high redshift $\beta$ has a significant trend of decreasing, at $\sim 3.5\sigma$ confidence level (CL). Moreover, we find that low-$z$ subsample favors a constant $\beta$; in contrast, SDSS and SNLS subsamples favor a decreasing $\beta$ at 2$\sigma$ and $3.3\sigma$ CL, respectively. Besides, by using a binned parameterization of $\beta$, we study the impacts of $\beta$'s evolution on parameter estimation. We find that compared with a constant $\beta$, a varying $\beta$ yields a larger best-fit value of fractional matter density $\Omega_{m0}$, which slightly deviates from the best-fit result given by other cosmological observations. However, for both the varying $\beta$ and the constant $\beta$ cases, the $1\sigma$ regions of $\Omega_{m0}$ are still consistent with the result given by other observations.

Shafieloo ea al. firstly proposed the possibility that the current cosmic acceleration (CA) is slowing down. However, this is rather counterintuitive because a slowing down CA cannot be accommodated in most mainstream cosmological models. In this work, by exploring the evolutionary trajectories of dark energy equation of state $w(z)$ and deceleration parameter $q(z)$, we present a comprehensive investigation on the slowing down of CA from both the theoretical and the observational sides. For the theoretical side, we study the impact of different $w(z)$ by using six parametrization models, and then discuss the effects of spatial curvature. For the observational side, we investigate the effects of different type Ia supernovae (SNe Ia), different baryon acoustic oscillation (BAO), and different cosmic microwave background (CMB) data, respectively. We find that (1) The evolution of CA are insensitive to the specific form of $w(z)$; in contrast, a non-flat Universe more favors a slowing down CA than a flat Universe. (2) SNLS3 SNe Ia datasets favor a slowing down CA at 1$\sigma$ confidence level, while JLA SNe Ia samples prefer an eternal CA; in contrast, the effects of different BAO data are negligible. (3) Compared with CMB distance prior data, full CMB data more favor a slowing down CA. (4) Due to the low significance, the slowing down of CA is still a theoretical possibility that cannot be confirmed by the current observations.

The sample of cosmological strong lensing systems has been steadily growing in recent years and with the advent of the next generation of space-based survey telescopes, the sample will reach into the thousands. The accuracy of strong lens models relies on robust identification of multiple image families of lensed galaxies. For the most massive lenses, often more than one background galaxy is magnified and multiply-imaged, and even in the cases of only a single lensed source, identification of counter images is not always robust. Recently, we have shown that the Gini coefficient in space-telescope-quality imaging is a measurement of galaxy morphology that is relatively well-preserved by strong gravitational lensing. Here, we investigate its usefulness as a diagnostic for the purposes of image family identification and show that it can remove some of the degeneracies encountered when using color as the sole diagnostic, and can do so without the need for additional observations since whenever a color is available, two Gini coefficients are as well.

Inspired by the multiverse scenario, we study a heterotic dark energy model in which there are two parts, the first being the cosmological constant and the second being the holographic dark energy, thus this model is named the $\Lambda$HDE model. By studying the $\Lambda$HDE model theoretically, we find that the parameters $d$ and $\Omega_{hde}$ are divided into a few domains in which the fate of the universe is quite different. We investigate dynamical behaviors of this model, and especially the future evolution of the universe. We perform fitting analysis on the cosmological parameters in the $\Lambda$HDE model by using the recent observational data. We find the model yields $\chi^2_{\rm min}=426.27$ when constrained by $\rm Planck+SNLS3+BAO+HST$, comparable to the results of the HDE model (428.20) and the concordant $\Lambda$CDM model (431.35). At 68.3\% CL, we obtain $-0.07<\Omega_{\Lambda0}<0.68$ and correspondingly $0.04<\Omega_{hde0}<0.79$, implying at present there is considerable degeneracy between the holographic dark energy and cosmological constant components in the $\Lambda$HDE model.

A variation of Affleck-Dine mechanism was proposed to generate the observed baryon asymmetry in [1], in which the inflaton was assumed to be a complex scalar field with a weakly broken $U(1)$ symmetry, and the baryon asymmetry generation was easily unified with the stage of inflation and reheating. We adapt this mechanism to natural inflation scenarios and compare the results with those in chaotic inflation models. We compute the net particle number obtained at the end of inflation and transform it into net baryon number after reheatings. We observed that in natural inflation models, the desired baryon-to-photon ratio can be achieved equally well as in chaotic models.

An inflationary model in the framework of noncommutative space-time may generate a nontrivial running of the scalar spectral index, but usually induces a large tensor-to-scalar ratio simultaneously. With the latest observational data from the Planck mission, we reexamine the inflationary scenarios in a noncommutative space-time. We find that either the running of the spectral index is tiny compared with the recent observational result, or the tensor-to-scalar ratio is too large to allow a sufficient number of $e$-folds. As examples, we show that the chaotic and power-law inflation models with the noncommutative effects are not favored by the current Planck data.

In the 12CO (J=1-0) survey for the 1331 cold IRAS sources 214 sources show profiles with multiple-peak profiles and are selected as cloud-cloud collision candidates. In January 2005, 201 sources are detected with 12CO(1-0), 13CO(1-0), and C18O(1-0) emission by the 13.7m telescope at Purple Mount Observatory. This is the first CO and its isotope lines directed toward possible cloud-cloud collision regions. According to the statistics of the 201 sources in Galactic distribution, the 201 sources show a similar distribution to the parent sample (1331 cold IRAS sources). These sources are located over a wide range of the Galactocentric distances, and are partly associated with the star formation region. Based on preliminary criteria which describe the spectrum properties of the possible cloud-cloud collision region, the 201 sources are classified into four types by the fit of the spectral profiles between the optically thick and thin lines toward each source. The survey is focused on the possible cloud-cloud collision regions, and gives some evidences to help us with selecting the target region. Then we will carry on the mapping and multi-wavelength study for the selected region in future.

Multi-wavelength astronomical studies require cross-identification of detections of the same celestial objects in multiple catalogs based on spherical coordinates and other properties. Because of the large data volumes and spherical geometry, the symmetric N-way association of astronomical detections is a computationally intensive problem, even when sophisticated indexing schemes are used to exclude obviously false candidates. Legacy astronomical catalogs already contain detections of more than a hundred million objects while the ongoing and future surveys will produce catalogs of billions of objects with multiple detections of each at different times. The varying statistical error of position measurements, moving and extended objects, and other physical properties make it necessary to perform the cross-identification using a mathematically correct, proper Bayesian probabilistic algorithm, capable of including various priors. One time, pair-wise cross-identification of these large catalogs is not sufficient for many astronomical scenarios. Consequently, a novel system is necessary that can cross-identify multiple catalogs on-demand, efficiently and reliably. In this paper, we present our solution based on a cluster of commodity servers and ordinary relational databases. The cross-identification problems are formulated in a language based on SQL, but extended with special clauses. These special queries are partitioned spatially by coordinate ranges and compiled into a complex workflow of ordinary SQL queries. Workflows are then executed in a parallel framework using a cluster of servers hosting identical mirrors of the same data sets.

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.

In the context of averaging an inhomogeneous cosmological model, we propose a natural measure identical to the Kullback-Leibler relative information entropy, which expresses the distinguishability of the local inhomogeneous density field from its spatial average on arbitrary compact domains. This measure is expected to be an increasing function in time and thus to play a significant role in studying gravitational entropy. To verify this conjecture, we explore the time evolution of the measure using the linear perturbation theory of a spatially flat FLRW model and a spherically symmetric nonlinear solution. We discuss the generality and conditions for the time-increasing nature of the measure, and also the connection to the backreaction effect caused by inhomogeneities.

Cosmological numerical simulations of galaxy formation have led to the cuspy density profile of a pure cold dark matter halo toward the center, which is in sharp contradiction with the observations of the rotation curves of cold dark matter-dominated dwarf and low surface brightness disk galaxies, with the latter tending to favor mass profiles with a flat central core. Many efforts have been devoted to resolve this cusp-core problem in recent years, among them, baryon-cold dark matter interactions are considered to be the main physical mechanisms erasing the cold dark matter (CDM) cusp into a flat core in the centers of all CDM halos. Clearly, baryon-cold dark matter interactions are not customized only for CDM-dominated disk galaxies, but for all types, including giant ellipticals. We first fit the most recent high resolution observations of rotation curves with the Burkert profile, then use the constrained core size-halo mass relation to calculate the lensing frequency, and compare the predicted results with strong lensing observations. Unfortunately, it turns out that the core size constrained from rotation curves of disk galaxies cannot be extrapolated to giant ellipticals. We conclude that, in the standard cosmological paradigm, baryon-cold dark matter interactions are not universal mechanisms for galaxy formation, and therefore, they cannot be true solutions to the cusp-core problem.

Due to the non-commutation of spatial averaging and temporal evolution, inhomogeneities and anisotropies (cosmic structures) influence the evolution of the averaged Universe via the cosmological backreaction mechanism. We study the backreaction effect as a function of averaging scale in a perturbative approach up to higher orders. We calculate the hierarchy of the critical scales, at which 10% effects show up from averaging at different orders. The dominant contribution comes from the averaged spatial curvature, observable up to scales of 200 Mpc. The cosmic variance of the local Hubble rate is 10% (5%) for spherical regions of radius 40 (60) Mpc. We compare our result to the one from Newtonian cosmology and Hubble Space Telescope Key Project data.

The present standard model of cosmology states that the known particles carry only a tiny fraction of total mass and energy of the Universe. Rather, unknown dark matter and dark energy are the dominant contributions to the cosmic energy budget. We review the logic that leads to the postulated dark energy and present an alternative point of view, in which the puzzle may be solved by properly taking into account the influence of cosmic structures on global observables. We illustrate the effect of averaging on the measurement of the Hubble constant.

Cosmological backreaction has been suggested as an explanation of dark energy and is heavily disputed since. We combine cosmological perturbation theory with Buchert's non-perturbative framework, calculate the relevant averaged observables up to second order in the comoving synchronous gauge, and discuss their gauge dependence. With the help of an integrability condition the leading second order contributions follow from the first order calculation. We focus on the onset of cosmological backreaction, as a perturbative treatment is necessarily restricted to the era when the effect is still small. We demonstrate that the leading contributions to all averaged physical observables are completely specified on the boundary of the averaged domain. For any finite domain these surface terms are nonzero in general and thus backreaction is for real. We map the backreaction effect on an effectively homogeneous and isotropic fluid. The generic effective equation of state is not only time dependent, but also depends on scale.

We discuss the relations between the mixing angles and the mixing matrices of quarks and leptons. With Raidal's numerical relations, we parametrize the lepton mixing (PMNS) matrix with the parameters of the quark mixing (CKM) matrix, and calculate the products of $V_{\mathrm{CKM}}U_{\mathrm{PMNS}}$ and $U_{\mathrm{PMNS}}V_{\mathrm{CKM}}$. Also, under the conjectures $V_{\mathrm{CKM}}U_{\mathrm{PMNS}}=U_{\mathrm{bimax}}$ or $U_{\mathrm{PMNS}}V_{\mathrm{CKM}}=U_{\mathrm{bimax}}$, we get the PMNS matrix naturally, and test Raidal's relations in these two different versions. The similarities and the differences between the different versions are discussed in detail.

We apply Koide's mass relation of charged leptons to neutrinos and quarks, with both the normal and inverted mass schemes of neutrinos discussed. We introduce the parameters $k_{\nu}$, $k_u$ and $k_d$ to describe the deviations of neutrinos and quarks from Koide's relation, and suggest a quark-lepton complementarity of masses such as $ k_{l}+k_{d} \approx k_{\nu}+k_{u} \approx 2$. The masses of neutrinos are determined from the improved relation, and they are strongly hierarchical (with the different orders of magnitude of $10^{-5} eV$, $10^{-3} eV$, and $10^{-2} eV$).

We present a unified parametrization of quark and lepton mixing matrices. By using some simple relations between the mixing angles of quarks and leptons, i.e., the quark-lepton complementarity, we parametrize the lepton mixing matrix with the Wolfenstein parameters $\lambda$ and $A$ of the quark mixing matrix. It is shown that the Wolfenstein parameter $\lambda$ can measure both the deviation of the quark mixing matrix from the unit matrix, and the deviation of the lepton mixing matrix from the exactly bimaximal mixing pattern.

The neutrino mixing matrix is expanded in powers of a small parameter $\lambda$, which approximately equals to 0.1. The meaning of every order of the expansion is discussed respectively, and the range of $\lambda$ is carefully calculated. We also present some applications of this new parametrization, such as to the expression of the Jarlskog parameter $J$, in which the simplicities and advantages of this parametrization are shown.