• We present a new millimeter CO-line observation towards supernova remnant (SNR) CTB 87, which was regarded purely as a pulsar wind nebula (PWN), and an optical investigation of a coincident surrounding superbubble. The CO observation shows that the SNR delineated by the radio emission is projectively covered by a molecular cloud (MC) complex at V$_{\rm {LSR}}$ = $-60$-$-54$ km s$^{-1}$. Both the symmetric axis of the radio emission and the trailing X-ray PWN appear projectively to be along a gap between two molecular gas patches at $-58$-$-57$ km s$^{-1}$. Asymmetric broad profiles of $^{12}$CO lines peaked at $-58$ km s$^{-1}$ are found at the eastern and southwestern edges of the radio emission. This represents a kinematic signature consistent with an SNR-MC interaction. We also find that a superbubble, $\sim 37'$ in radius, appears to surround the SNR from HI 21cm (V$_{\rm {LSR}} \sim -61$-$-68$ km s$^{-1}$), WISE mid-IR, and optical extinction data. We build a multi-band photometric stellar sample of stars within the superbubble region and find 82 OB star candidates. The likely peak distance in the stars' distribution seems consistent with the distance previously suggested for CTB 87. We suggest the arc-like radio emission is mainly the relic of the part of blastwave that propagates into the MC complex and is now in a radiative stage while the other part of blastwave has been expanding into the low-density region in the superbubble. This scenario naturally explains the lack of the X-ray emission related to the ejecta and blastwave. The SNR-MC interaction also favours a hadronic contribution to the {\gamma}-ray emission from the CTB 87 region.
  • Mendelian randomization (MR) is an instrumental variable method of estimating the causal effect of risk exposures in epidemiology, where genetic variants are used as instruments. With the increasing availability of large-scale genome-wide association studies, it is now possible to greatly improve the power of MR by using genetic variants that are only weakly relevant. We consider how to increase the efficiency of Mendelian randomization by a genome-wide design where more than a thousand genetic instruments are used. An empirical partially Bayes estimator is proposed, where weaker instruments are shrunken more heavily and thus brings less variation to the MR estimate. This is generally more efficient than the profile-likelihood-based estimator which gives no shrinkage to weak instruments. We apply our method to estimate the causal effect of blood lipids on cardiovascular diseases. We find high-density lipoprotein cholesterol (HDL-c) has a significantly protective effect on heart diseases, while previous MR studies reported null findings.
  • If we can automatically derive compiler optimizations, we might be able to sidestep some of the substantial engineering challenges involved in creating and maintaining a high-quality compiler. We developed Souper, a synthesizing superoptimizer, to see how far these ideas might be pushed in the context of LLVM. Along the way, we discovered that Souper's intermediate representation was sufficiently similar to the one in Microsoft Visual C++ that we applied Souper to that compiler as well. Shipping, or about-to-ship, versions of both compilers contain optimizations suggested by Souper but implemented by hand. Alternately, when Souper is used as a fully automated optimization pass it compiles a Clang compiler binary that is about 3 MB (4.4%) smaller than the one compiled by LLVM.
  • We study the recurrence coefficients of the monic polynomials $P_n(z)$ orthogonal with respect to the deformed (also called semi-classical) Freud weight \begin{equation*} w_{\alpha}(x;s,N)=|x|^{\alpha}{\rm e}^{-N\left[x^{2}+s\left(x^{4}-x^{2}\right)\right]}, ~~x\in\mathbb{R}, \end{equation*} with parameters $\alpha>-1,~N>0,~s\in[0,1]$. We show that the recurrence coefficients $\beta_{n}(s)$ satisfy the first discrete Painlev\'{e} equation (denoted by d${\rm P_{I}}$), a differential-difference equation and a second order nolinear ordinary differential equation (ODE) in $s$. Here $n$ is the order of the Hankel matrix generated by $w_{\alpha}(x;s,N)$. We describe the asymptotic behavior of the recurrence coefficients in three situations, (i) $s\rightarrow0$, $n,N$ finite, (ii) $n\rightarrow\infty$, $N$ finite, (iii) $n, N\rightarrow\infty$, such that the radio $r:=\frac{n}{N}$ is bounded away from $0$ and closed to $1$. We also investigate the existence and uniqueness for the positive solutions of the d${\rm P_{I}}$. Further more, we derive, using the ladder approach, a second order linear ODE satisfied by the polynomials $P_n(z)$. It is found as $n\rightarrow\infty$, the linear ODE turns to be a biconfluent Heun equation. This paper concludes with the study of the Hankel determinant, $D_{n}(s)$, associated with $w_{\alpha}(x;s,N)$ when $n$ tends to infinity.
  • We study the asymptotic behavior of the smallest eigenvalue, $\lambda_{N}$, of the Hankel (or moments) matrix denoted by $\mathcal{H}_{N}=\left(\mu_{m+n}\right)_{0\leq m,n\leq N}$, with respect to the weight $w(x)=x^{\alpha}{\rm e}^{-x^{\beta}},~x\in[0,\infty),~\alpha>-1,~\beta>\frac{1}{2}$. Based on the research by Szeg\"{o}, Chen, etc., we obtain an asymptotic expression of the orthonormal polynomials $\mathcal{P}_{N}(z)$ as $N\rightarrow\infty$, associated with $w(x)$. Using this, we obtain the specific asymptotic formulas of $\lambda_{N}$ in this paper. Applying the parallel algorithm discovered by Emmart, Chen and Weems, we get a variety of numerical results of $\lambda_{N}$ corresponding to our theoretical calculations.
  • We investigate the large $N$ behavior of the smallest eigenvalue, $\lambda_{N}$, of an $\left(N+1\right)\times \left(N+1\right)$ Hankel (or moments) matrix $\mathcal{H}_{N}$, generated by the weight $w(x)=x^{\alpha}(1-x)^{\beta},~x\in[0,1],~ \alpha>-1,~\beta>-1$. By applying the arguments of Szeg\"{o}, Widom and Wilf, we establish the asymptotic formula for the orthonormal polynomials $P_{n}(z),z\in\mathbb{C}\setminus[0,1]$, associated with $w(x)$, which are required in the determination of $\lambda_{N}$. Based on this formula, we produce the expressions for $\lambda_{N}$, for large $N$. Using the parallel algorithm presented by Emmart, Chen and Weems, we show that the theoretical results are in close proximity to the numerical results for sufficiently large $N$.
  • In this paper, we study the gap probability problem of the (symmetric) Jacobi unitary ensemble of Hermitian random matrices, namely the probability that the interval $(-a,a)\:(0<a<1)$ is free of eigenvalues. Using the ladder operator technique for orthogonal polynomials and the associated supplementary conditions, we derive three quantities instrumental in the gap probability, denoted by $H_{n}(a)$, $R_{n}(a)$ and $r_{n}(a)$. We find that each one satisfies a second order differential equation. We show that after a double scaling, the large second order differential equation in the variable $a$ with $n$ as parameter satisfied by $H_{n}(a)$, can be reduced to the Jimbo-Miwa-Okamoto $\sigma$ form of the Painlev\'{e} V equation.
  • This paper studies the generating function of a discontinuous linear statistic in the Gaussian unitary ensemble. It is an extension of Chen and Pruessner \cite{Chen2005}, in which they studied the discontinuous Gaussian weight with a single jump. By using the ladder operator approach, we obtain a series of difference and differential equations to describe the generating function for the single jump case. These equations include the Chazy II equation, continuous and discrete Painlev\'{e} IV. Moreover, we consider the large $n$ behavior of the corresponding orthogonal polynomials and prove that they satisfy the biconfluent Heun equation. We also consider the jump at the edge under a double scaling, from which a Painlev\'{e} XXXIV appeared. Furthermore, we study the linear statistics with two jumps, and show that a quantity related to the generating function satisfies a two variables' generalization of the Jimbo-Miwa-Okamoto $\sigma$ form of the Painlev\'{e} IV.
  • Multiphoton entanglement plays a critical role in quantum information processing, and greatly improves our fundamental understanding of the quantum world. Despite tremendous efforts in either bulk media or fiber-based devices, nonlinear interactions in integrated circuits show great promise as an excellent platform for photon pair generation with its high brightness, stability and scalability \cite{Caspani2017}. Here, we demonstrate the generation of bi- and multiphoton polarization entangled qubits in a single silicon nanowire waveguide, and these qubits directly compatible with the dense wavelength division multiplexing in telecommunication system. Multiphoton interference and quantum state tomography were used to characterize the quality of the entangled states. Four-photon entanglement states among two frequency channels were ascertained with a fidelity of $0.78\pm0.02$. Our work realizes the integrated multiphoton source in a relatively simple pattern and paves a way for the revolution of multiphoton quantum science.
  • An on-chip photonic quantum source, especially an on-chip entangled photon source, is an essential resource for quantum information applications. Here, we report an on-chip transverse-mode entangled photon source, which is realized using spontaneous four-wave mixing processes in a multimode silicon waveguide. Transverse-mode entangled photon pairs are generated and experimentally verified with a bandwidth of $\sim 2\ THz$; a maximally transverse-mode entangled Bell state can also be produced with a fidelity of $0.92\pm0.01$. The demonstrated on-chip entangled photon source provides one of the most important key elements for developing quantum photonics using transverse-mode freedom, which can be used to encode quantum information within a high-dimensional Hilbert space. And the transverse-mode entanglement can be converted coherently to path and polarization entanglement on-chip. This paves the way to realize highly complex quantum photonic circuits with multiple degrees of freedom and plays an important role in the high-dimensional quantum information processing.
  • In this paper, we study the probability density function, $\mathbb{P}(c,\alpha,\beta, n)\,dc$, of the center of mass of the finite $n$ Jacobi unitary ensembles with parameters $\alpha\,>-1$ and $\beta >-1$; that is the probability that ${\rm tr}M_n\in(c, c+dc),$ where $M_n$ are $n\times n$ matrices drawn from the unitary Jacobi ensembles. We first compute the exponential moment generating function of the linear statistics $\sum_{j=1}^{n}\,f(x_j):=\sum_{j=1}^{n}x_j,$ denoted by $\mathcal{M}_f(\lambda,\alpha,\beta,n)$.
  • Precise studies on the Galactic bulge, globular cluster, Galactic halo and Galactic thick disk require stellar models with alpha enhancement and various values of helium content. These models are also important for extra-Galactic population synthesis studies. For this purpose, we complement the existing PARSEC models, which are based on the solar partition of heavy elements, with alpha-enhanced partitions. We collect detailed measurements on the metal mixture and helium abundance for the two populations of 47 Tuc (NGC 104) from the literature, and calculate stellar tracks and isochrones with these alpha-enhanced compositions. By fitting the precise color-magnitude diagram with HST ACS/WFC data, from low main sequence till horizontal branch, we calibrate some free parameters that are important for the evolution of low mass stars like the mixing at the bottom of the convective envelope. This new calibration significantly improves the prediction of the RGB bump brightness. Comparison with the observed RGB and HB luminosity functions also shows that the evolutionary lifetimes are correctly predicted. As a further result of this calibration process, we derive the age, distance modulus, reddening, and the red giant branch mass loss for 47 Tuc. We apply the new calibration and alpha-enhanced mixtures of the two 47 Tuc populations ( [alpha/Fe] ~0.4 and 0.2) to other metallicities. The new models reproduce the RGB bump observations much better than previous models. This new PARSEC database, with the newly updated alpha-enhanced stellar evolutionary tracks and isochrones, will also be part of the new stellar products for Gaia.
  • The cross-site linking function is widely adopted by online social networks (OSNs). This function allows a user to link her account on one OSN to her accounts on other OSNs. Thus, users are able to sign in with the linked accounts, share contents among these accounts and import friends from them. It leads to the service integration of different OSNs. This integration not only provides convenience for users to manage accounts of different OSNs, but also introduces usefulness to OSNs that adopt the cross-site linking function. In this paper, we investigate this usefulness based on users' data collected from a popular OSN called Medium. We conduct a thorough analysis on its social graph, and find that the service integration brought by the cross-site linking function is able to change Medium's social graph structure and attract a large number of new users. However, almost none of the new users would become high PageRank users (PageRank is used to measure a user's influence in an OSN). To solve this problem, we build a machine-learning-based model to predict high PageRank users in Medium based on their Twitter data only. This model achieves a high F1-score of 0.942 and a high area under the curve (AUC) of 0.986. Based on it, we design a system to assist new OSNs to identify and attract high PageRank users from other well-established OSNs through the cross-site linking function.
  • In this paper, we address a class of problems in unitary ensembles. Specifically, we study the probability that a gap symmetric about 0, i.e. $(-a,a)$ is found in the Gaussian unitary ensembles (GUE) and the Jacobi unitary ensembles (JUE) (where in the JUE, we take the parameters $\alpha=\beta$). By exploiting the even parity of the weight, a doubling of the interval to $(a^2,\infty)$ for the GUE, and $(a^2,1)$, for the (symmetric) JUE, shows that the gap probabilities maybe determined as the product of the smallest eigenvalue distributions of the LUE with parameter $\alpha=-1/2,$ and $\alpha=1/2$ and the (shifted) JUE with weights $x^{1/2}(1-x)^{\beta}$ and $x^{-1/2}(1-x)^{\beta}$ The $\sigma$ function, namely, the derivative of the log of the smallest eigenvalue distributions of the finite-$n$ LUE or the JUE, satisfies the Jimbo-Miwa-Okamoto $\sigma$ form of $P_{V}$ and $P_{VI}$, although in the shift Jacobi case, with the weight $x^{\alpha}(1-x)^{\beta},$ the $\beta$ parameter does not show up in the equation. We also obtain the asymptotic expansions for the smallest eigenvalue distributions of the Laguerre unitary and Jacobi unitary ensembles after appropriate double scalings, and obtained the constants in the asymptotic expansion of the gap probablities, expressed in term of the Barnes $G-$ function valuated at special point.
  • HESS J1912+101 is a shell-like TeV source that has no clear counterpart in multiwavelength. Using CO and H i data, we reveal that VLSR~+60 km/s molecular clouds (MCs), together with shocked molecular gas and high-velocity neutral atomic shells, are concentrated toward HESS J1912+101. The prominent wing profiles up to VLSR~+80 km/s seen in 12CO (J=1-0 and J=3-2) data, as well as the high-velocity expanding H i shells up to VLSR~ +100 km/s, exhibit striking redshifted-broadening relative to the quiescent gas. These features provide compelling evidences for large-scale perturbation in the region. We argue that the shocked MCs and the high-velocity Hi shells may originate from an old supernova remnant (SNR). The distance to the SNR is estimated to be ~4.1 kpc based on the Hi self-absorption method, which leads to a physical radius of 29.0 pc for the ~(0.7-2.0)x10e5 years old remnant with an expansion velocity of >40 km/s. The +60 km/s MCs and the disturbed gas are indeed found to coincide with the bright TeV emission, supporting the physical association between them. Naturally, the shell-like TeV emission comes from the decay of neutral pions produced by interactions between the accelerated hadrons from the SNR and the surrounding high-density molecular gas.
  • Given the potential X-ray radiation risk to the patient, low-dose CT has attracted a considerable interest in the medical imaging field. The current main stream low-dose CT methods include vendor-specific sinogram domain filtration and iterative reconstruction, but they need to access original raw data whose formats are not transparent to most users. Due to the difficulty of modeling the statistical characteristics in the image domain, the existing methods for directly processing reconstructed images cannot eliminate image noise very well while keeping structural details. Inspired by the idea of deep learning, here we combine the autoencoder, the deconvolution network, and shortcut connections into the residual encoder-decoder convolutional neural network (RED-CNN) for low-dose CT imaging. After patch-based training, the proposed RED-CNN achieves a competitive performance relative to the-state-of-art methods in both simulated and clinical cases. Especially, our method has been favorably evaluated in terms of noise suppression, structural preservation and lesion detection.
  • The present work carries out molecular dynamics simulations to compute the thermal conductivity of the borophene nanoribbon and the borophene nanotube using the Muller-Plathe approach. We investigate the thermal conductivity of the armchair and zigzag borophenes, and show the strong anisotropic thermal conductivity property of borophene. We compare the results of the borophene nanoribbon and the borophene nanotube, and find the thermal conductivity of the borophene is structure dependent.
  • At present, multi-electrode array (MEA) approach and optical recording allow us to acquire plant electrical activity with higher spatio-temporal resolution. To understand the dynamic information flow of the electrical signaling system and estimate the effective connectivity, we proposed a solution to combine the two casualty analysis approaches, i.e. Granger causality and transfer entropy, which they complement each other to measure dynamics effective connectivity of the complex system. Our findings in three qualitatively different levels of plant bioelectrical activities revealed direction of information flow and dynamic complex causal connectives by using the two causal analysis approaches, especially indicated that the direction of information flow is not only along the longitudinal section but also spreading in transection.
  • In this paper, we propose a new simple and learning-free deep learning network named MomentsNet, whose convolution layer, nonlinear processing layer and pooling layer are constructed by Moments kernels, binary hashing and block-wise histogram, respectively. Twelve typical moments (including geometrical moment, Zernike moment, Tchebichef moment, etc.) are used to construct the MomentsNet whose recognition performance for binary image is studied. The results reveal that MomentsNet has better recognition performance than its corresponding moments in almost all cases and ZernikeNet achieves the best recognition performance among MomentsNet constructed by twelve moments. ZernikeNet also shows better recognition performance on binary image database than that of PCANet, which is a learning-based deep learning network.
  • As an alternative option to kinetic electrons, the gyrokinetic total-f particle-in-cell (PIC) code XGC1 has been extended to the MHD/fluid type electromagnetic regime by combining gyrokinetic PIC ions with massless drift-fluid electrons analogous to Chen and Parker, Physics of Plasmas 8, 441 (2001). Two representative long wavelength modes, shear Alfv\'en waves and resistive tearing modes, are verified in cylindrical and toroidal magnetic field geometries.
  • We present CO observations toward three large supernova remnants (SNRs) in the third Galactic quadrant using the Purple Mountain Observatory Delingha 13.7m radio telescope. The observations are part of the high-resolution CO survey of the Galactic plane between Galactic longitudes l=-10deg to 250deg and latitudes b=-5deg to 5d. CO emission was detected toward the three SNRs: G205.5+0.5 (Monoceros Nebula), G206.9+2.3 (PKS 0646+06), and G213.0-0.6. Both of SNRs G205.5+0.5 and G213.0-0.6 exhibit the morphological agreement (or spatial correspondences) between the remnant and the surrounding molecular clouds (MCs), as well as kinematic signatures of shock perturbation in the molecular gas. We confirm that the two SNRs are physically associated with their ambient MCs and the shock of SNRs is interacting with the dense, clumpy molecular gas. SNR G206.9+2.3, which is close to the northeastern edge of the Monoceros Nebula, displays the spatial coincidence with molecular partial shell structures at VLSR~15km/s. While no significant line broadening has been detected within or near the remnant, the strong morphological correspondence between the SNR and the molecular cavity implies that SNR G206.9+2.3 is probably associated with these CO gas and is evolving in the low-density environment. The physical features of individual SNRs, together with the relationship between SNRs and their nearby objects, are also discussed.
  • Hyperspectral image (HSI) denoising has been attracting much research attention in remote sensing area due to its importance in improving the HSI qualities. The existing HSI denoising methods mainly focus on specific spectral and spatial prior knowledge in HSIs, and share a common underlying assumption that the embedded noise in HSI is independent and identically distributed (i.i.d.). In real scenarios, however, the noise existed in a natural HSI is always with much more complicated non-i.i.d. statistical structures and the under-estimation to this noise complexity often tends to evidently degenerate the robustness of current methods. To alleviate this issue, this paper attempts the first effort to model the HSI noise using a non-i.i.d. mixture of Gaussians (NMoG) noise assumption, which is finely in accordance with the noise characteristics possessed by a natural HSI and thus is capable of adapting various noise shapes encountered in real applications. Then we integrate such noise modeling strategy into the low-rank matrix factorization (LRMF) model and propose a NMoG-LRMF model in the Bayesian framework. A variational Bayes algorithm is designed to infer the posterior of the proposed model. All involved parameters can be recursively updated in closed-form. Compared with the current techniques, the proposed method performs more robust beyond the state-of-the-arts, as substantiated by our experiments implemented on synthetic and real noisy HSIs.
  • We introduce a new generation of PARSEC-COLIBRI stellar isochrones that include a detailed treatment of the thermally-pulsing asymptotic giant branch (TP-AGB) phase, and covering a wide range of initial metallicities (0.0001<Zi<0.06). Compared to previous releases, the main novelties and improvements are: use of new TP-AGB tracks and related atmosphere models and spectra for M and C-type stars; inclusion of the surface H+He+CNO abundances in the isochrone tables, accounting for the effects of diffusion, dredge-up episodes and hot-bottom burning; inclusion of complete thermal pulse cycles, with a complete description of the in-cycle changes in the stellar parameters; new pulsation models to describe the long-period variability in the fundamental and first overtone modes; new dust models that follow the growth of the grains during the AGB evolution, in combination with radiative transfer calculations for the reprocessing of the photospheric emission. Overall, these improvements are expected to lead to a more consistent and detailed description of properties of TP-AGB stars expected in resolved stellar populations, especially in regard to their mean photometric properties from optical to mid-infrared wavelengths. We illustrate the expected numbers of TP-AGB stars of different types in stellar populations covering a wide range of ages and initial metallicities, providing further details on the C-star island that appears at intermediate values of age and metallicity, and about the AGB-boosting effect that occurs at ages close to 1.6 Gyr for populations of all metallicities. The isochrones are available through a new dedicated web server.
  • Strong light-matter interaction and high-efficiency optical collection of fluorescence from quantum emitters are crucial topics in quantum and nanophotonic fields. High-quality cavities, dispersive photonic crystal waveguides and even plasmonic structures have been used to enhance the interaction with quantum emitters, thus realize efficient collection of the fluorescence. In this work, a new method is proposed to collect the fluorescence of quantum dots (QDs) with a fiber-integrated multimode silver nanowire (AgNW) waveguide. Fluorescence lifetime measurement is performed to investigate the coupling between QDs and different plasmonic modes. Compared with far-field collection method, the AgNW-fiber probe can realize near-unity collection efficiency theoretically. This fiber-integrated plasmonic probe may be useful in the area of nanophotonics and also promising for quantum information devices.
  • Statistical distances, divergences, and similar quantities have a large history and play a fundamental role in statistics, machine learning and associated scientific disciplines. However, within the statistical literature, this extensive role has too often been played out behind the scenes, with other aspects of the statistical problems being viewed as more central, more interesting, or more important. The behind the scenes role of statistical distances shows up in estimation, where we often use estimators based on minimizing a distance, explicitly or implicitly, but rarely studying how the properties of a distance determine the properties of the estimators. Distances are also prominent in goodness-of-fit, but the usual question we ask is "how powerful is this method against a set of interesting alternatives" not "what aspect of the distance between the hypothetical model and the alternative are we measuring?" Our focus is on describing the statistical properties of some of the distance measures we have found to be most important and most visible. We illustrate the robust nature of Neyman's chi-squared and the non-robust nature of Pearson's chi-squared statistics and discuss the concept of discretization robustness.