
Graph based semisupervised learning (GSSL) has intuitive representation and
can be improved by exploiting the matrix calculation. However, it has to
perform iterative optimization to achieve a preset objective, which usually
leads to low efficiency. Another inconvenience lying in GSSL is that when new
data come, the graph construction and the optimization have to be conducted all
over again. We propose a sound assumption, arguing that: the neighboring data
points are not in peertopeer relation, but in a partialordered relation
induced by the local density and distance between the data; and the label of a
center can be regarded as the contribution of its followers. Starting from the
assumption, we develop a highly efficient noniterative label propagation
algorithm based on a novel data structure named as optimal leading forest
(LaPOLeaF). The major weaknesses of the traditional GSSL are addressed by this
study. We further scale LaPOLeaF to accommodate big data by utilizing block
distance matrix technique, parallel computing, and LocalitySensitive Hashing
(LSH). Experiments on large datasets have shown the promising results of the
proposed methods.

We investigate the leading twist lightcone distribution amplitudes (LCDAs)
of vector meson in the framework of large momentum effective theory. We derive
the matching equation for the LCDAs and quasi distribution amplitudes. The
matching coefficients are determined to one loop accuracy, both in the
ultraviolet cutoff and dimensional regularization schemes. This calculation
provides the possibility of studying the full $x$ behavior of LCDAs and
extracting LCDAs of vector mesons from lattice simulations.

After the experimental establishment of doubly heavy baryons, baryons with
three quarks are the last missing pieces of the lowestlying baryon multiplets
in quark model. In this work we study semileptonic and nonleptonic weak decays
of triply heavy baryons, $\Omega_{ccc}^{++}, \Omega_{ccb}^{+},
\Omega_{cbb}^{0}, \Omega_{bbb}^{}$. Decay amplitudes for various channels are
parametrized in terms of a few SU(3) irreducible amplitudes. We point out that
branching fractions for Cabibbo allowed processes, $\Omega_{ccc}\to
(\Xi_{cc}^{++} \overline K^0, \Xi_{cc}^{++}K^\pi^+, \Omega_{cc}^{+}\pi^+,
\Xi_{c}^+ D^+, \Xi_{c}^{\prime} D^+, \Lambda_c D^+\overline K^0, \Xi_{c}^+ D^0
\pi^+, \Xi_{c}^0 D^+\pi^+)$ may reach a few percents. We suggest our
experimental colleagues to perform a search at hadron colliders and the
electron and positron collisions in future, which will presumably lead to
discoveries of triply heavy baryons and complete the baryon multiplets. Using
the expanded amplitudes, we derive a number of relations for the partial widths
which can be examined in future.

We consider an $\ell_2$regularized nonconvex optimization problem for
recovering signals from their noisy phaseless observations. We design and study
the performance of a message passing algorithm that aims to solve this
optimization problem. We consider the asymptotic setting $m,n \rightarrow
\infty$, $m/n \rightarrow \delta$ and obtain sharp performance bounds, where
$m$ is the number of measurements and $n$ is the signal dimension. We show that
for complex signals the algorithm can perform accurate recovery with only $m=
\left(\frac{64}{\pi^2}4\right)n \approx 2.5n$ measurements. Also, we provide
sharp analysis on the sensitivity of the algorithm to noise. We highlight the
following facts about our message passing algorithm: (i) Adding $\ell_2$
regularization to the nonconvex loss function can be beneficial. (ii) Spectral
initialization has marginal impact on the performance of the algorithm. The
sharp analyses in this paper, not only enable us to compare the performance of
our method with other phase recovery schemes, but also shed light on designing
better iterative algorithms for other nonconvex optimization problems.

The newlydiscovered $\Xi_{cc}^{++}$ decays into the $ \Lambda_{c}^+
K^\pi^+\pi^+$, but the experimental data has indicated that this decay is not
saturated by any twobody intermediate state. In this work, we analyze the
multibody weak decays of doubly heavy baryons $\Xi_{cc}$, $\Omega_{cc}$,
$\Xi_{bc}$, $\Omega_{bc}$, $\Xi_{bb}$ and $\Omega_{bb}$, in particular the
threebody nonleptonic decays and fourbody semileptonic decays. We classify
various decay modes according to the quarklevel transitions and present an
estimate of the typical branching fractions for a few golden decay channels.
Decay amplitudes are then parametrized in terms of a few SU(3) irreducible
amplitudes. With these amplitudes, we find a number of relations for decay
widths, which can be examined in future.

Long ShortTerm Memory (LSTM) is the primary recurrent neural networks
architecture for acoustic modeling in automatic speech recognition systems.
Residual learning is an efficient method to help neural networks converge
easier and faster. In this paper, we propose several types of residual LSTM
methods for our acoustic modeling. Our experiments indicate that, compared with
classic LSTM, our architecture shows more than 8% relative reduction in Phone
Error Rate (PER) on TIMIT tasks. At the same time, our residual fast LSTM
approach shows 4% relative reduction in PER on the same task. Besides, we find
that all this architecture could have good results on THCHS30, Librispeech and
Switchboard corpora.

Motivated by the recent LHCb observation of doublycharmed baryon
$\Xi_{cc}^{++}$ in the $\Lambda_c^+ K^\pi^+\pi^+$ final state, we analyze the
weak decays of doubly heavy baryons $\Xi_{cc}$, $\Omega_{cc}$, $\Xi_{bc}$,
$\Omega_{bc}$, $\Xi_{bb}$ and $\Omega_{bb}$ under the flavor SU(3) symmetry.
Decay amplitudes for various semileptonic and nonleptonic decays are
parametrized in terms of a few SU(3) irreducible amplitudes. We find a number
of relations or sum rules between decay widths and CP asymmetries, which can be
examined in future measurements at experimental facilities like LHC, Belle II
and CEPC. Moreover once a few decay branching fractions are measured in future,
some of these relations may provide hints for exploration of new decay modes.

Fluctuations of conserved quantities, such as baryon, electric charge and
strangeness number, are sensitive observables in heavyion collisions to search
for the QCD phase transition and critical point. In this paper, we performed a
systematical analysis on the various cumulants and cumulant ratios of
eventbyevent netstrangeness distributions in Au+Au collisions at
$\sqrt{s_{NN}}$=7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV from UrQMD model. We
performed a systematical study on the contributions from various strange
baryons and mesons to the netstrangeness fluctuations. The results demonstrate
that the cumulants and cumulant ratios of netstrangeness distributions
extracted from different strange particles show very different centrality and
energy dependence behavior. By comparing with the netkaon fluctuations, we
found that the strange baryons play an important role in the fluctuations of
netstrangeness. This study can provide useful baselines to study the QCD phase
transition and search for the QCD critical point by using the fluctuations of
netstrangeness in heavyion collisions experiment. It can help us to
understand noncritical physics contributions to the fluctuations of
netstrangeness.

In collinear factorization, lightcone distribution amplitudes (LCDAs) are
key ingredients to calculate the production rate of a hadron in high energy
exclusive processes. For a doublyheavy meson system (such as $B_c, J/\psi,
\Upsilon$ etc), the LCDAs contain perturbative scales that can be integrated
out and then are refactorized into products of perturbatively calculable
distribution parts and nonrelativistic QCD matrix elements. In this
refactorization scheme, the LCDAs are known at nexttoleading order in the
strong coupling constant $\alpha_s$ and at leading order in the velocity
expansion. In this work, we calculate the ${\cal O}( { v}^2)$ corrections to
twist2 LCDAs of Swave $B_c$ mesons. These results are applicable to heavy
quarkonia like $\eta_{c,b}$, $J/\psi$ and $\Upsilon$ by setting $m_b=m_c$. We
apply these relativistically corrected LCDAs to study their inverse moments and
a few Gegenbauer moments which are important for phenomenological study. We
point out that the relativistic corrections are sizable, and comparable with
the nexttoleading order radiative corrections. These results for LCDAs are
useful in future theoretical analyses of the productions of heavy quarkonia and
$B_c$ mesons.

By extending local U(1) gauge symmetry to discontinuous case, we find that
under one special discontinuous U(1) gauge transformation the symmetric and
antisymmetric wave functions can transform into each other in one dimensional
quantum mechanics. The free spinless fermionic system and bosonic system with
$\delta$type vector gauge potential are proved to be equivalent. The relation
also holds in higher spacetime dimensions.

Fluctuations of conserved quantities such as baryon number (B), electric
charge number (Q), and strangeness number (S), are sensitive to the correlation
length and can be used to probe nongaussian fluctuations near the critical
point. Experimentally, higher moments of the multiplicity distributions have
been used to search for the QCD critical point in heavyion collisions. In this
paper, we report the efficiencycorrected cumulants and their ratios of mid
rapidity (y < 0.5) netkaon multiplicity distributions in Au+Au collisions at
7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV collected in 2010, 2011, and
2014 with STAR at RHIC. The centrality and energy dependence of the cumulants
and their ratios, are presented. Furthermore, the comparisons with baseline
calculations (Poisson) and noncriticalpoint models (UrQMD) are also
discussed.

One of the main goals of the RHIC Beam Energy Scan (BES) program is to study
the QCD phase structure, which includes the search for the QCD critical point,
over a wide range of chemical potential. Theoretical calculations predict that
fluctuations of conserved quantities, such as baryon number (B), charge (Q),
and strangeness (S), are sensitive to the correlation length of the dynamical
system. Experimentally, higher moments of multiplicity distributions have been
utilized to search for the QCD critical point in heavyion collisions. In this
paper, we report recent efficiencycorrected cumulants and cumulants ratios of
the net proton, netkaon, and netcharge multiplicity distributions in Au+Au
collisions at 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV collected in the
years 2010, 2011, and 2014 with STAR at RHIC. The centrality and energy
dependence of the cumulants up to the fourth order, as well as their ratios,
are presented. Furthermore, the comparisons with baseline calculations
(Poisson) and noncriticalpoint models (UrQMD) will also be discussed.

Expectation Maximization (EM) is among the most popular algorithms for
estimating parameters of statistical models. However, EM, which is an iterative
algorithm based on the maximum likelihood principle, is generally only
guaranteed to find stationary points of the likelihood objective, and these
points may be far from any maximizer. This article addresses this disconnect
between the statistical principles behind EM and its algorithmic properties.
Specifically, it provides a global analysis of EM for specific models in which
the observations comprise an i.i.d. sample from a mixture of two Gaussians.
This is achieved by (i) studying the sequence of parameters from idealized
execution of EM in the infinite sample limit, and fully characterizing the
limit points of the sequence in terms of the initial parameters; and then (ii)
based on this convergence analysis, establishing statistical consistency (or
lack thereof) for the actual sequence of parameters produced by EM.

The lightcone distribution amplitudes (LCDAs) serve as important
nonperturbative inputs for the study of hard exclusive processes. In this
paper, we calculate ten LCDAs at twist2 for the Swave and Pwave $B_c$ mesons
up to the nexttoleading order (NLO) of the strong coupling $\alpha_s$ and
leading order of the velocity expansion. Each one of these ten LCDAs is
expressed as a product of a perturbatively calculable distribution and a
universal NRQCD matrixelement. By use of the spin symmetry, only two NRQCD
matrixelements will be involved. The reduction of the number of
nonperturbative inputs will improve the predictive power of collinear
factorization.

Fluctuations of conserved quantities are sensitive observables to probe the
signature of QCD phase transition and critical point in heavyion collisions.
With the UrQMD model, we have studied the centrality and energy dependence of
various order cumulants and cumulant ratios (up to fourth order) of
netproton,netcharge and netkaon multiplicity distributions in Au+Au
collisions at $\sqrt{s_{NN}}$= 7.7, 11.5, 19.6, 27, 39, 62.4, 200 GeV. The
model results show that the production mechanism of the particles and
antiparticles have significant impacts on the cumulants of netparticles
multiplicity distributions and show strong energy dependence. We also made
comparisons between model calculations and experimental data measured in the
first phase of the beam energy scan (BES) program by the STAR experiment at
RHIC. The comparisons indicate that the baryon conservation effect strongly
suppress the cumulants of netproton distributions at low energies and the
nonmonotonic energy dependence for the netproton {\KV} at the most central
Au+Au collisions measured by the STAR experiment can not be described by the
UrQMD model. Since there has no physics of QCD phase transition and QCD
critical point implemented in the UrQMD, the model results provide us baselines
and qualitative estimates about the noncritical background contributions to
the fluctuations observables in heavyion collisions.

This paper reveals the tree structure as an intermediate result of clustering
by fast search and find of density peaks (DPCLUS), and explores the power of
using this tree to perform hierarchical clustering. The array used to hold the
index of the nearest higherdensitied object for each object can be transformed
into a Leading Tree (LT), in which each parent node P leads its child nodes to
join the same cluster as P itself, and the child nodes are sorted by their
gamma values in descendant order to accelerate the disconnecting of root in
each subtree. There are two major advantages with the LT: One is dramatically
reducing the running time of assigning noncenter data points to their cluster
ID, because the assigning process is turned into just disconnecting the links
from each center to its parent. The other is that the tree model for
representing clusters is more informative. Because we can check which objects
are more likely to be selected as centers in finer grained clustering, or which
objects reach to its center via less jumps. Experiment results and analysis
show the effectiveness and efficiency of the assigning process with an LT.

Neither of the two prevalent theories, namely thermodynamic stability and
kinetic stability, provides a comprehensive understanding of protein folding.
The thermodynamic theory is misleading because it assumes that free energy is
the exclusive dominant mechanism of protein folding, and attributes the
structural transition from one characteristic state to another to energy
barriers. Conversely, the concept of kinetic stability overemphasizes dominant
mechanisms that are related to kinetic factors. This article explores the
stability condition of protein structures from the viewpoint of mesoscience,
paying attention to the compromise in the competition between minimum free
energy and other dominant mechanisms. Based on our study of complex systems, we
propose that protein folding is a mesoscale, dissipative, nonlinear and
nonequilibrium process that is dominated by the compromise between free energy
and other dominant mechanisms such as environmental factors. Consequently, a
protein shows dynamic structures, featuring characteristic states that appear
alternately and dynamically, only one of which is the state with minimum free
energy. To provide evidence for this concept, we analyzed the time series of
energetic and structural changes of three simulations of protein
folding/unfolding. Our results indicate that thorough consideration of the
multiple dynamic characteristic structures generated by multiple mechanisms may
be the key to understanding protein folding.

Moments (Variance ($\sigma^2$), Skewness($S$), Kurtosis($\kappa$)) of
multiplicity distributions of conserved quantities, such as
netbaryon,netcharge and netstrangeness, are predicted to be sensitive to the
correlation length of the system and connected to the thermodynamic
susceptibilities computed in Lattice QCD and Hadron Resonance Gas (HRG) model.
In this paper, we present several measurement artifacts that could lead to
volume fluctuation and autocorrelation effects in the moment analysis of
netproton multiplicity distributions in heavyion collisions using the UrQMD
model. We discuss methods to overcome these artifacts so that the extracted
moments could be used to obtain physical conclusions. In addition we present
methods to properly estimate the statistical errors in moment analysis.

Molecular dynamics (MD) simulation is a powerful computational tool to study
the behavior of macromolecular systems. But many simulations of this field are
limited in spatial or temporal scale by the available computational resource.
In recent years, graphics processing unit (GPU) provides unprecedented
computational power for scientific applications. Many MD algorithms suit with
the multithread nature of GPU. In this paper, MD algorithms for macromolecular
systems that run entirely on GPU are presented. Compared to the MD simulation
with free software GROMACS on a single CPU core, our codes achieve about 10
times speedup on a single GPU. For validation, we have performed MD
simulations of polymer crystallization on GPU, and the results observed
perfectly agree with computations on CPU. Therefore, our single GPU codes have
already provided an inexpensive alternative for macromolecular simulations on
traditional CPU clusters and they can also be used as a basis to develop
parallel GPU programs to further speedup the computations.