
We present a new millimeter COline 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 Xray
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 SNRMC
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 midIR, and optical extinction data. We build a multiband 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
arclike 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 lowdensity region in the
superbubble. This scenario naturally explains the lack of the Xray emission
related to the ejecta and blastwave. The SNRMC 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 largescale
genomewide 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
genomewide 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 profilelikelihoodbased 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 highdensity
lipoprotein cholesterol (HDLc) 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 highquality 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 abouttoship, 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 semiclassical) 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 differentialdifference 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}(1x)^{\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
JimboMiwaOkamoto $\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 JimboMiwaOkamoto $\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 fiberbased 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. Fourphoton 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 onchip photonic quantum source, especially an onchip entangled photon
source, is an essential resource for quantum information applications. Here, we
report an onchip transversemode entangled photon source, which is realized
using spontaneous fourwave mixing processes in a multimode silicon waveguide.
Transversemode entangled photon pairs are generated and experimentally
verified with a bandwidth of $\sim 2\ THz$; a maximally transversemode
entangled Bell state can also be produced with a fidelity of $0.92\pm0.01$. The
demonstrated onchip entangled photon source provides one of the most important
key elements for developing quantum photonics using transversemode freedom,
which can be used to encode quantum information within a highdimensional
Hilbert space. And the transversemode entanglement can be converted coherently
to path and polarization entanglement onchip. This paves the way to realize
highly complex quantum photonic circuits with multiple degrees of freedom and
plays an important role in the highdimensional 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 extraGalactic
population synthesis studies. For this purpose, we complement the existing
PARSEC models, which are based on the solar partition of heavy elements, with
alphaenhanced 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
alphaenhanced compositions. By fitting the precise colormagnitude 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
alphaenhanced 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 alphaenhanced stellar evolutionary tracks and isochrones, will also be
part of the new stellar products for Gaia.

The crosssite 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 crosssite 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 crosssite 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 machinelearningbased model to predict high
PageRank users in Medium based on their Twitter data only. This model achieves
a high F1score 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 wellestablished OSNs through the crosssite 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}(1x)^{\beta}$ and
$x^{1/2}(1x)^{\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 JimboMiwaOkamoto $\sigma$ form of $P_{V}$ and $P_{VI}$,
although in the shift Jacobi case, with the weight $x^{\alpha}(1x)^{\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 shelllike 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 highvelocity neutral
atomic shells, are concentrated toward HESS J1912+101. The prominent wing
profiles up to VLSR~+80 km/s seen in 12CO (J=10 and J=32) data, as well as
the highvelocity expanding H i shells up to VLSR~ +100 km/s, exhibit striking
redshiftedbroadening relative to the quiescent gas. These features provide
compelling evidences for largescale perturbation in the region. We argue that
the shocked MCs and the highvelocity 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 selfabsorption method, which leads to a physical radius of
29.0 pc for the ~(0.72.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 shelllike TeV emission comes from the decay of neutral pions
produced by interactions between the accelerated hadrons from the SNR and the
surrounding highdensity molecular gas.

Given the potential Xray radiation risk to the patient, lowdose CT has
attracted a considerable interest in the medical imaging field. The current
main stream lowdose CT methods include vendorspecific 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 encoderdecoder convolutional neural
network (REDCNN) for lowdose CT imaging. After patchbased training, the
proposed REDCNN achieves a competitive performance relative to
thestateofart 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 MullerPlathe 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, multielectrode array (MEA) approach and optical recording allow
us to acquire plant electrical activity with higher spatiotemporal 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 learningfree deep learning
network named MomentsNet, whose convolution layer, nonlinear processing layer
and pooling layer are constructed by Moments kernels, binary hashing and
blockwise 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 learningbased deep learning network.

As an alternative option to kinetic electrons, the gyrokinetic totalf
particleincell (PIC) code XGC1 has been extended to the MHD/fluid type
electromagnetic regime by combining gyrokinetic PIC ions with massless
driftfluid 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 highresolution 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.00.6. Both
of SNRs G205.5+0.5 and G213.00.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 lowdensity
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 noni.i.d. statistical structures and the
underestimation 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 noni.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 lowrank matrix factorization (LRMF) model and
propose a NMoGLRMF 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 closedform. Compared with
the current techniques, the proposed method performs more robust beyond the
stateofthearts, as substantiated by our experiments implemented on synthetic
and real noisy HSIs.

We introduce a new generation of PARSECCOLIBRI stellar isochrones that
include a detailed treatment of the thermallypulsing asymptotic giant branch
(TPAGB) 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 TPAGB tracks and related atmosphere models and
spectra for M and Ctype stars; inclusion of the surface H+He+CNO abundances in
the isochrone tables, accounting for the effects of diffusion, dredgeup
episodes and hotbottom burning; inclusion of complete thermal pulse cycles,
with a complete description of the incycle changes in the stellar parameters;
new pulsation models to describe the longperiod 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 TPAGB stars expected in resolved stellar populations, especially
in regard to their mean photometric properties from optical to midinfrared
wavelengths. We illustrate the expected numbers of TPAGB stars of different
types in stellar populations covering a wide range of ages and initial
metallicities, providing further details on the Cstar island that appears at
intermediate values of age and metallicity, and about the AGBboosting 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 lightmatter interaction and highefficiency optical collection of
fluorescence from quantum emitters are crucial topics in quantum and
nanophotonic fields. Highquality 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 fiberintegrated multimode silver
nanowire (AgNW) waveguide. Fluorescence lifetime measurement is performed to
investigate the coupling between QDs and different plasmonic modes. Compared
with farfield collection method, the AgNWfiber probe can realize nearunity
collection efficiency theoretically. This fiberintegrated 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 goodnessoffit, 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 chisquared and the nonrobust nature of Pearson's
chisquared statistics and discuss the concept of discretization robustness.