
This paper is concerned with the hard thresholding operator which sets all
but the $k$ largest absolute elements of a vector to zero. We establish a {\em
tight} bound to quantitatively characterize the deviation of the thresholded
solution from a given signal. Our theoretical result is universal in the sense
that it holds for all choices of parameters, and the underlying analysis
depends only on fundamental arguments in mathematical optimization. We discuss
the implications for two domains:
Compressed Sensing. On account of the crucial estimate, we bridge the
connection between the restricted isometry property (RIP) and the sparsity
parameter for a vast volume of hard thresholding based algorithms, which
renders an improvement on the RIP condition especially when the true sparsity
is unknown. This suggests that in essence, many more kinds of sensing matrices
or fewer measurements are admissible for the data acquisition procedure.
Machine Learning. In terms of largescale machine learning, a significant yet
challenging problem is learning accurate sparse models in an efficient manner.
In stark contrast to prior work that attempted the $\ell_1$relaxation for
promoting sparsity, we present a novel stochastic algorithm which performs hard
thresholding in each iteration, hence ensuring such parsimonious solutions.
Equipped with the developed bound, we prove the {\em global linear convergence}
for a number of prevalent statistical models under mild assumptions, even
though the problem turns out to be nonconvex.

In the context of HumanRobot Interaction (HRI), face ReIdentification (face
ReID) aims to verify if certain detected faces have already been observed by
robots. The ability of distinguishing between different users is crucial in
social robots as it will enable the robot to tailor the interaction strategy
toward the users' individual preferences. So far face recognition research has
achieved great success, however little attention has been paid to the realistic
applications of Face ReID in social robots. In this paper, we present an
effective and unsupervised face ReID system which simultaneously reidentifies
multiple faces for HRI. This ReID system employs Deep Convolutional Neural
Networks to extract features, and an online clustering algorithm to determine
the face's ID. Its performance is evaluated on two datasets: the TERESA video
dataset collected by the TERESA robot, and the YouTube Face Dataset (YTF
Dataset). We demonstrate that the optimised combination of techniques achieves
an overall 93.55% accuracy on TERESA dataset and an overall 90.41% accuracy on
YTF dataset. We have implemented the proposed method into a software module in
the HCI^2 Framework for it to be further integrated into the TERESA robot, and
has achieved realtime performance at 10~26 Frames per second.

The number of electrons on isolated small metallic islands and semiconductor
quantum dots is quantized. When tunnelling is enabled via opaque barriers this
number can change by an integer.Every extra electron adds an elementary charge,
e, at an energy cost, a charging energy, which at low temperatures regulates
the electron flow as onebyone, single electron tunnelling.In superconductors
the addition is in units of 2e charges, reflecting that the Cooper pair
condensate must have an even parity [3]. This evenparity ground state is
foundational for all superconducting qubit devices. Here, we study a hybrid
superconducting (aluminium)semiconducting (InSb) island and find that a
magnetic field can induce an even to odd parity transition in the
superconducting ground state. This parity transition can occur when a single,
spinresolved subgap state (i.e. an Andreev bound state, ABS) has crossed zero
energy.In addition, we also find that the magnetic field can cause a change
from 2e to 1e charge quantization while the aluminium remains superconducting.
This observation is compatible with ABS at zero energy or the presence of
Majorana zero modes (MZMs).

We consider a family of random locations, called intrinsic location
functionals, of periodic stationary processes. This family includes but is not
limited to the location of the path supremum and first/last hitting times. We
first show that the set of all possible distributions of intrinsic location
functionals for periodic stationary processes is the convex hull generated by a
specific group of distributions. We then focus on two special subclasses of
these random locations. For the first subclass, the density has a uniform lower
bound; for the second subclass, the possible distributions are closely related
to the concept of joint mixability.

Silent speech interfaces have been recently proposed as a way to enable
communication when the acoustic signal is not available. This introduces the
need to build visual speech recognition systems for silent and whispered
speech. However, almost all the recently proposed systems have been trained on
vocalised data only. This is in contrast with evidence in the literature which
suggests that lip movements change depending on the speech mode. In this work,
we introduce a new audiovisual database which is publicly available and
contains normal, whispered and silent speech. To the best of our knowledge,
this is the first study which investigates the differences between the three
speech modes using the visual modality only. We show that an absolute decrease
in classification rate of up to 3.7% is observed when training and testing on
normal and whispered, respectively, and vice versa. An even higher decrease of
up to 8.5% is reported when the models are tested on silent speech. This
reveals that there are indeed visual differences between the 3 speech modes and
the common assumption that vocalized training data can be used directly to
train a silent speech recognition system may not be true.

This paper presents the Graph Analytics Repository for Designing
Nextgeneration Accelerators (GARDENIA), a benchmark suite for studying
irregular algorithms on massively parallel accelerators. Existing generic
benchmarks for accelerators have mainly focused on high performance computing
(HPC) applications with limited control and data irregularity, while available
graph analytics benchmarks do not apply stateoftheart algorithms and/or
optimization techniques. GARDENIA includes emerging irregular applications in
bigdata and machine learning domains which mimic massively multithreaded
commercial programs running on modern largescale datacenters. Our
characterization shows that GARDENIA exhibits irregular microarchitectural
behavior which is quite different from structured workloads and
straightforwardimplemented graph benchmarks.

We study asymptotically and numerically the fundamental gap  the difference
between the first two smallest (and distinct) eigenvalues  of the fractional
Schr\"{o}dinger operator (FSO) and formulate a gap conjecture on the
fundamental gap of the FSO. We begin with an introduction of the FSO on bounded
domains with homogeneous Dirichlet boundary conditions, while the fractional
Laplacian operator defined either via the local fractional Laplacian (i.e. via
the eigenfunctions decomposition of the Laplacian operator) or via the
classical fractional Laplacian (i.e. zero extension of the eigenfunctions
outside the bounded domains and then via the Fourier transform). For the FSO on
bounded domains with either the local fractional Laplacian or the classical
fractional Laplacian, we obtain the fundamental gap of the FSO analytically on
simple geometry without potential and numerically on complicated geometries
and/or with different convex potentials. Based on the asymptotic and extensive
numerical results, a gap conjecture on the fundamental gap of the FSO is
formulated. Surprisingly, for two and higher dimensions, the lower bound of the
fundamental gap depends not only on the diameter of the domain, but also the
diameter of the largest inscribed ball of the domain, which is completely
different from the case of the Schr\"{o}dinger operator. Extensions of these
results for the FSO in the whole space and on bounded domains with periodic
boundary conditions are presented.

We propose a new numerical technique to deal with nonlinear terms in gradient
flows. By introducing a scalar auxiliary variable (SAV), we construct efficient
and robust energy stable schemes for a large class of gradient flows. The SAV
approach is not restricted to specific forms of the nonlinear part of the free
energy, and only requires to solve {\it decoupled} linear equations with {\it
constant coefficients}. We use this technique to deal with several challenging
applications which can not be easily handled by existing approaches, and
present convincing numerical results to show that our schemes are not only much
more efficient and easy to implement, but can also better capture the physical
properties in these models. Based on this SAV approach, we can construct
unconditionally secondorder energy stable schemes; and we can easily construct
even third or fourth order BDF schemes, although not unconditionally stable,
which are very robust in practice. In particular, when coupled with an adaptive
time stepping strategy, the SAV approach can be extremely efficient and
accurate.

In this paper, we characterize compatibility of distributions and probability
measures on a measurable space. For a set of indices $\mathcal J$, we say that
the tuples of probability measures $(Q_i)_{i\in \mathcal J} $ and distributions
$(F_i)_{i\in \mathcal J} $ are {compatible} if there exists a random variable
having distribution $F_i$ under $Q_i$ for each $i\in \mathcal J$. We first
establish an equivalent condition using conditional expectations for general
(possibly uncountable) $\mathcal J$. For a finite $n$, it turns out that
compatibility of $(Q_1,\dots,Q_n)$ and $(F_1,\dots,F_n)$ depends on the
heterogeneity among $Q_1,\dots,Q_n$ compared with that among $F_1,\dots,F_n$.
We show that, under an assumption that the measurable space is rich enough,
$(Q_1,\dots,Q_n)$ and $(F_1,\dots,F_n)$ are compatible if and only if
$(Q_1,\dots,Q_n)$ dominates $(F_1,\dots,F_n)$ in a notion of heterogeneity
order, defined via multivariate convex order between the RadonNikodym
derivatives of $(Q_1,\dots,Q_n)$ and $(F_1,\dots,F_n)$ with respect to some
reference measures.
We then proceed to generalize our results to stochastic processes, and
conclude the paper with an application to portfolio selection problems under
multiple constraints.

Object proposals are an ensemble of bounding boxes with high potential to
contain objects. In order to determine a small set of proposals with a high
recall, a common scheme is extracting multiple features followed by a ranking
algorithm which however, incurs two major challenges: {\bf 1)} The ranking
model often imposes pairwise constraints between each proposal, rendering the
problem away from an efficient training/testing phase; {\bf 2)} Linear kernels
are utilized due to the computational and memory bottleneck of training a
kernelized model.
In this paper, we remedy these two issues by suggesting a {\em kernelized
partial ranking model}. In particular, we demonstrate that {\bf i)} our partial
ranking model reduces the number of constraints from $O(n^2)$ to $O(nk)$ where
$n$ is the number of all potential proposals for an image but we are only
interested in top$k$ of them that has the largest overlap with the ground
truth; {\bf ii)} we permit nonlinear kernels in our model which is often
superior to the linear classifier in terms of accuracy. For the sake of
mitigating the computational and memory issues, we introduce a consistent
weighted sampling~(CWS) paradigm that approximates the nonlinear kernel as
well as facilitates an efficient learning. In fact, as we will show, training a
linear CWS model amounts to learning a kernelized model. Extensive experiments
demonstrate that equipped with the nonlinear kernel and the partial ranking
algorithm, recall at top$k$ proposals can be substantially improved.

Semiconductor nanowires provide an ideal platform for various lowdimensional
quantum devices. In particular, topological phases of matter hosting
nonAbelian quasiparticles can emerge when a semiconductor nanowire with
strong spinorbit coupling is brought in contact with a superconductor. To
fully exploit the potential of nonAbelian anyons for topological quantum
computing, they need to be exchanged in a wellcontrolled braiding operation.
Essential hardware for braiding is a network of singlecrystalline nanowires
coupled to superconducting islands. Here, we demonstrate a technique for
generic bottomup synthesis of complex quantum devices with a special focus on
nanowire networks having a predefined number of superconducting islands.
Structural analysis confirms the high crystalline quality of the nanowire
junctions, as well as an epitaxial superconductorsemiconductor interface.
Quantum transport measurements of nanowire "hashtags" reveal AharonovBohm and
weakantilocalization effects, indicating a phase coherent system with strong
spinorbit coupling. In addition, a proximityinduced hard superconducting gap
is demonstrated in these hybrid superconductorsemiconductor nanowires,
highlighting the successful materials development necessary for a first
braiding experiment. Our approach opens new avenues for the realization of
epitaxial 3dimensional quantum device architectures.

Tempered fractional derivatives originated from the tempered fractional
diffusion equations (TFDEs) modeled on the whole space R (see [23]). For
numerically solving TFDEs, two kinds of generalized Laguerre functions were
defined and some important properties were proposed to establish the
approximate theory. The related prototype tempered fractional differ ential
problems was proposed and solved as the guidance. TFDEs are numerically solved
by two domains Laguerre spectral method and the numerical experiments show some
properties of the TFDEs and verify the efficiency of the spectral scheme.

How to develop efficient numerical schemes while preserving the energy
stability at the discrete level is a challenging issue for the three component
CahnHilliard phasefield model. In this paper, we develop first and second
order temporal approximation schemes based on the "Invariant Energy
Quadratization" approach, where all nonlinear terms are treated
semiexplicitly. Consequently, the resulting numerical schemes lead to a
wellposed linear system with the symmetric positive definite operator to be
solved at each time step. We rigorously prove that the proposed schemes are
unconditionally energy stable. Various 2D and 3D numerical simulations are
presented to demonstrate the stability and the accuracy of the schemes.

In this paper, we consider numerical approximations of a hydrodynamically
coupled phase field diblock copolymer model, in which the free energy contains
a kinetic potential, a gradient entropy, a GinzburgLandau double well
potential, and a long range nonlocal type potential. We develop a set of second
order time marching schemes for this system using the "Invariant Energy
Quadratization" approach for the double well potential, the projection method
for the NavierStokes equation, and a subtle implicitexplicit treatment for
the stress and convective term. The resulting schemes are linear and lead to
symmetric positive definite systems at each time step, thus they can be
efficiently solved. We further prove that these schemes are unconditionally
energy stable. Various numerical experiments are performed to validate the
accuracy and energy stability of the proposed schemes.

Because they may provide ultrathin, highflux, and energyefficient membranes
for precise ionic and molecular sieving in aqueous solution, GO membranes
(partially oxidized, stacked sheets of graphene) have shown great potential in
water desalination and purification, gas and ion separation, biosensors, proton
conductors, lithiumbased batteries and supercapacitors. Unlike carbon
nanotube (CNT) membranes, in which the nanotube pores have fixed sizes, the
pores of GO membranes  the interlayer spacing between GO sheets  are of
variable size. This presents a challenge for using GO membranes for filtration.
Despite the great efforts to tune and fix the interlayer spacing, it remains
difficult both to reduce the interlayer spacing sufficiently to exclude small
ions while keeping this separation constant against the tendency of GO
membranes to swell when immersed in aqueous solution, which greatly affects the
applications of GO membranes. Here, we demonstrate experimentally that highly
efficient and selective ion rejection by GO membranes can be readily achieved
by controlling the interlayer spacing of GO membranes using cations (K+, Na+,
Ca2+, Li+ and Mg2+) themselves. The interspacing can be controlled with
precision as small as 1 A, and GO membranes controlled by one kind of cation
can exclude other cations with a larger hydrated volume, which can only be
accommodated with a larger interlayer spacing. Firstprinciples calculations
reveal that the strong noncovalent cationpi interactions between hydrated
cations in solution and aromatic ring structures in GO are the cause of this
unexpected behavior. These findings open up new avenues for using GO membranes
for water desalination and purification, lithiumbased batteries and
supercapacitors, molecular sieves for separating ions or molecules, and many
other applications.

The quantized version of the anomalous Hall effect has been predicted to
occur in magnetic topological insulators, but the experimental realization has
been challenging. Here, we report the observation of the quantum anomalous Hall
(QAH) effect in thin films of Crdoped (Bi,Sb)2Te3, a magnetic topological
insulator. At zero magnetic field, the gatetuned anomalous Hall resistance
reaches the predicted quantized value of h/e^2,accompanied by a considerable
drop of the longitudinal resistance. Under a strong magnetic field, the
longitudinal resistance vanishes whereas the Hall resistance remains at the
quantized value. The realization of the QAH effect may lead to the development
of lowpowerconsumption electronics.

LowRank Representation~(LRR) has been a significant method for segmenting
data that are generated from a union of subspaces. It is also known that
solving LRR is challenging in terms of time complexity and memory footprint, in
that the size of the nuclear norm regularized matrix is $n$by$n$ (where $n$
is the number of samples). In this paper, we thereby develop a novel online
implementation of LRR that reduces the memory cost from $O(n^2)$ to $O(pd)$,
with $p$ being the ambient dimension and $d$ being some estimated rank~($d < p
\ll n$). We also establish the theoretical guarantee that the sequence of
solutions produced by our algorithm converges to a stationary point of the
expected loss function asymptotically. Extensive experiments on synthetic and
realistic datasets further substantiate that our algorithm is fast, robust and
memory efficient.

Maxnorm regularizer has been extensively studied in the last decade as it
promotes an effective lowrank estimation for the underlying data. However,
such maxnorm regularized problems are typically formulated and solved in a
batch manner, which prevents it from processing big data due to possible memory
budget. In this paper, hence, we propose an online algorithm that is scalable
to largescale setting. Particularly, we consider the matrix decomposition
problem as an example, although a simple variant of the algorithm and analysis
can be adapted to other important problems such as matrix completion. The
crucial technique in our implementation is to reformulating the maxnorm to an
equivalent matrix factorization form, where the factors consist of a (possibly
overcomplete) basis component and a coefficients one. In this way, we may
maintain the basis component in the memory and optimize over it and the
coefficients for each sample alternatively. Since the memory footprint of the
basis component is independent of the sample size, our algorithm is appealing
when manipulating a large collection of samples. We prove that the sequence of
the solutions (i.e., the basis component) produced by our algorithm converges
to a stationary point of the expected loss function asymptotically. Numerical
study demonstrates encouraging results for the efficacy and robustness of our
algorithm compared to the widely used nuclear norm solvers.

For topological insulators and topological crystalline insulators (TCIs),
their exotic surface states are promising for fundamental condensed matter
physics research as well as future electronics such as lowdissipation
electronics and spintronics. However, the high bulk carrier density that often
dominates the transport property is the major materials challenge, critically
hindering our ability to study and manipulate the surface states. In this
manuscript, we demonstrate an alloying strategy, SnSexTe1x, to effectively
reduce the bulk carrier density. As long as SnSexTe1x remains in the cubic
crystal structure, it is predicted to be a TCI. We show systematic decrease of
the bulk carrier density with the increasing Se concentration, demonstrating
that the alloying principle works. In addition, we map out the phase diagram of
the cubic to the orthorhombic structural transition as a function of the Se
concentration. This was made possible by studying alloy nanoplates which remain
singlecrystalline and is either in the cubic or the orthorhombic phase, in
contrast to bulk alloys that would exhibit polycrystalline grains. Lastly, we
investigate systematically the ferroelectric transition associated with the
structural transition from the cubic to the rhombohedral phase for SnSexTe1x.
This is the first ferroelectric transition study of the alloy system
SnSexTe1x.

Recently, much attention has been paid to search for Majorana fermions in
solidstate systems. Among various proposals there is one based on
radiofrequency superconducting quantum interference devices (rfSQUIDs), in
which the appearance of 4$\pi$period energyphase relations is regarded as
smokinggun evidence of Majorana fermion states. Here we report the observation
of truncated 4$\pi$period (i.e., 2$\pi$period but fully skewed) oscillatory
patterns of contact resistance on rfSQUIDs constructed on the surface of
threedimensional topological insulator Bi$_2$Te$_3$. The results reveal the
existence of 1/2 fractional modes of Cooper pairs and the occurrence of parity
switchings, both of which are necessary signatures accompanied with the
formation of Majorana fermion states.

Social trust prediction addresses the significant problem of exploring
interactions among users in social networks. Naturally, this problem can be
formulated in the matrix completion framework, with each entry indicating the
trustness or distrustness. However, there are two challenges for the social
trust problem: 1) the observed data are with sign (1bit) measurements; 2) they
are typically sampled nonuniformly. Most of the previous matrix completion
methods do not well handle the two issues. Motivated by the recent progress of
maxnorm, we propose to solve the problem with a 1bit maxnorm constrained
formulation. Since maxnorm is not easy to optimize, we utilize a reformulation
of maxnorm which facilitates an efficient projected gradient decent algorithm.
We demonstrate the superiority of our formulation on two benchmark datasets.

The main aim of this paper is to describe the dynamic transitions in flows
described by the twodimensional, barotropic vorticity equation in a periodic
zonal channel. In \cite{CGSW03}, the existence of a Hopf bifurcation in this
model as the Reynolds number crosses a critical value was proven. In this
paper, we extend the results in \cite{CGSW03} by addressing the stability
problem of the bifurcated periodic solutions. Our main result is the explicit
expression of a nondimensional number $\gamma$ which controls the transition
behavior. We prove that depending on $\gamma$, the modeled flow exhibits either
a continuous (Type I) or catastrophic (Type II) transition. Numerical
evaluation of $\gamma$ for a physically realistic region of parameter space
suggest that a catastrophic transition is preferred in this flow.

Indium (In) doping in topological crystalline insulator SnTe induces
superconductivity, making Indoped SnTe a candidate for a topological
superconductor. SnTe nanostructures offer welldefined nanoscale morphology and
high surfacetovolume ratios to enhance surface effects. Here, we study
Indoped SnTe nanoplates, InxSn1xTe, with x ranging from 0 to 0.1 and show
they superconduct. More importantly, we show that In doping reduces the bulk
mobility of InxSn1xTe such that the surface states are revealed in
magnetotransport despite the high bulk carrier density. This is manifested by
twodimensional linear magnetoresistance in high magnetic fields, which is
independent of temperature up to 10 K. Aging experiments show that the linear
magnetoresistance is sensitive to ambient conditions, further confirming its
surface origin. We also show that the weak antilocalization observed in
InxSn1xTe nanoplates is a bulk effect. Thus, we show that nanostructures and
reducing the bulk mobility are effective strategies to reveal the surface
states and test for topological superconductors.

As a fundamental technique that concerns several vision tasks such as image
parsing, action recognition and clothing retrieval, human pose estimation (HPE)
has been extensively investigated in recent years. To achieve accurate and
reliable estimation of the human pose, it is wellrecognized that the clothing
attributes are useful and should be utilized properly. Most previous
approaches, however, require to manually annotate the clothing attributes and
are therefore very costly. In this paper, we shall propose and explore a
\emph{latent} clothing attribute approach for HPE. Unlike previous approaches,
our approach models the clothing attributes as latent variables and thus
requires no explicit labeling for the clothing attributes. The inference of the
latent variables are accomplished by utilizing the framework of latent
structured support vector machines (LSSVM). We employ the strategy of
\emph{alternating direction} to train the LSSVM model: In each iteration, one
kind of variables (e.g., human pose or clothing attribute) are fixed and the
others are optimized. Our extensive experiments on two realworld benchmarks
show the stateoftheart performance of our proposed approach.

Topological crystalline insulators are topological insulators whose surface
states are protected by the crystalline symmetry, instead of the time reversal
symmetry. Similar to the first generation of threedimensional topological
insulators such as Bi2Se3 and Bi2Te3, topological crystalline insulators also
possess surface states with exotic electronic properties such as spinmomentum
locking and Dirac dispersion. Experimentally verified topological crystalline
insulators to date are SnTe, Pb1xSnxSe, and Pb1xSnxTe. Because topological
protection comes from the crystal symmetry, magnetic impurities or inplane
magnetic fields are not expected to open a gap in the surface states in
topological crystalline insulators. Additionally, because they are cubic
structure instead of layered structure, branched structures or strong coupling
with other materials for large proximity effects are possible, which are
difficult with layered Bi2Se3 and Bi2Te3. Thus, additional fundamental
phenomena inaccessible in threedimensional topological insulators can be
pursued. In this review, topological crystalline insulator SnTe nanostructures
will be discussed. For comparison, experimental results based on SnTe thin
films will be covered. Surface state properties of topological crystalline
insulators will be discussed briefly.