
The statistical problem of using an initial sample to estimate the number of
species in a larger sample has found important applications in fields far
removed from ecology. Here we address the general problem of estimating the
number of species that will be represented by at least a number r of
observations in a future sample. The number r indicates species with sufficient
observations, which are commonly used as a necessary condition for any robust
statistical inference. We derive a procedure to construct consistent estimators
that apply universally for a given population: once constructed, they can be
evaluated as a simple function of r. Our approach is based on a relation
between the number of species represented at least r times and the higher
derivatives of the expected number of species discovered per unit of time.
Combining this relation with a rational function approximation, we propose
nonparametric estimators that are accurate for both large values of r and
longrange extrapolations. We further show that our estimators retain
asymptotic behaviors that are essential for applications on largescale
datasets. We evaluate the performance of this approach by both simulation and
real data applications for inferences of the vocabulary of Shakespeare and
Dickens, the topology of a Twitter social network, and molecular diversity in
DNA sequencing data.

In the minimal supersymmetric standard model (MSSM) the binowino
coannihilation provides a feasible way to accommodate the observed cosmological
dark matter (DM) relic density. However, such a scenario usually predicts a
very small DMnucleon scattering cross section that is below the neutrino
floor, and can not be tested by DM direct detection experiments. In this work,
we investigate the discovery potential of this binowino coannihilation region
by searching for the soft dilepton events from the process $pp \to \chi^0_2
(\to \ell^+\ell^ \chi^0_1) \chi^\pm_1+jets$ at the LHC. We find that the mass
of the winolike $\chi^0_2$ can be probed up to about 310 (230) GeV at
$2\sigma$ ($5\sigma$) level for an integrated luminosity ${\cal L}=300$
fb$^{1}$. In the future HLLHC with 3000 fb$^{1}$ luminosity, the
corresponding mass limits can be pushed up to 430 (330) GeV.

We explore the fidelity susceptibility and the quantum coherence along with
the entanglement entropy in the groundstate of onedimensional spin1 XXZ
chains with the rhombic singleion anisotropy. By using the techniques of
density matrix renormalization group, effects of the rhombic singleion
anisotropy on a few information theoretical measures are investigated, such as
the fidelity susceptibility, the quantum coherence and the entanglement
entropy. Their relations with the quantum phase transitions are also analyzed.
The phase transitions from the YN\'{e}el phase to the Large$E_x$ or the
Haldane phase can be well characterized by the fidelity susceptibility. The
secondorder derivative of the groundstate energy indicates all the
transitions are of second order. We also find that the quantum coherence, the
entanglement entropy, the Schmidt gap can be used to detect the critical points
of quantum phase transitions. Conclusions drawn from these quantum information
observables agree well with each other. Finally we provide a groundstate phase
diagram as functions of the exchange anisotropy $\Delta$ and the rhombic
singleion anisotropy $E$.

The system of cold twolevel atoms inside an optical lattice potential has
been used to simulate various models encountered in fundamental and
condensedmatter physics. When the optical lattice is periodically shaken, some
interesting dynamical effects were observed. In this work, by utilizing the
ultranarrow linewidth of the dipoleforbidden transition in $^{87}$Sr atoms
loaded in a shaken optical lattice, we reach a new operation regime where the
shaking rate is much larger than linewidth of the twolevel transition but much
smaller than the frequency interval of the harmonic lattice potential, which
was not realizable before. Under such unique conditions, an interesting quantum
dynamical effect, the highorder sideband effect (HSE), has been experimentally
observed. All the results can be well explained by a developed theoretical
model with analytical solution. This HSE is governed by different mechanism as
the similar phenomenon observed in semiconductors and has unique
characteristics.

Inspired by the topological insulator circuit proposed and experimentally
verified by Jia., et al. \cite{1}, we theoretically realized the topological
Lieb lattice, a line centered square lattice with rich topological properties,
in a radiofrequency circuit. We open the topological nontrivial bandgap
through specific capacitorinductor network, which resembles adding intrinsic
spin orbit coupling term into the tight binding model. Finally, we discuss the
extension of the $\phi=\pi/2$ phase change of hopping between sites to
arbitrary value, and investigate the topological phase transition of the band
structure vary with capacitance, thereby paving the way for designing tunable
lattices using the presented framework.

Genome and metagenome comparisons based on large amounts of nextgeneration
sequencing (NGS) data pose significant challenges for alignmentbased
approaches due to the huge data size and the relatively short length of the
reads. Alignmentfree approaches based on the counts of word patterns in NGS
data do not depend on the complete genome and are generally computationally
efficient. Thus, they contribute significantly to genome and metagenome
comparison. Recently, novel statistical approaches have been developed for the
comparison of both long and shotgun sequences. These approaches have been
applied to many problems including the comparison of gene regulatory regions,
genome sequences, metagenomes, binning contigs in metagenomic data,
identification of virushost interactions, and detection of horizontal gene
transfers. We provide an updated review of these applications and other related
developments of wordcount based approaches for alignmentfree sequence
analysis.

By exploring the extraordinary property of exceptional points (EPs) in
nonHermitian systems, we here demonstrate that losses can play constructive
roles in controlling bistable states. We experimentally realize the EP in a
nonHermitian metamolecule of coupled resonators in microwave regime. By
increasing the loss, we first observe the bistable state suppression at the
weakdissipative regime, but then the bistable state recovery in
strongdissipative regime. Both the experimental and theoretical analysis
demonstrate that the revival of bistable states results from the revival of the
field intensity after the system encounters EPs, in spite of the increasing
loss. Our results provide an alternative way to controlling and manifesting
bistable systems, so as to achieve flexible photonic devices, not limited to
the microwave regime.

Zak phase, which refers to the Berry's phase picked up by a particle moving
across the Brillouin zone, characterizes the topological properties of Bloch
bands in onedimensional periodic system. Here the Zak phase in dimerized
onedimensional locally resonant metamaterials is investigated. It is found
that there are some singular points in the bulk band across which the Bloch
states contribute {\pi} to the Zak phase, whereas while in the rest of the band
the contribution is nearly zero. These singular points associated with zero
reflection are caused by two different mechanisms: the dimerizationindependent
antiresonating of each branch, and the dimerizationdependent destructive
interference in multiple backscattering. The structure undergoes a topological
transition point in the band structure where the band inverts and the Zak
phase, which is determined by the numbers of singular points in the bulk band,
changes following a shift in dimerization parameter. Finally, the interface
state between two dimerized metamaterial structures with different topological
property in the first band gap is demonstrated experimentally. The
quasionedimensional configuration of the system allows one to explore
topologyinspired new methods and applications in the subwavelength scale.

Investigating wave propagation in fluid enables a variety of important
applications in underwater communications, object detections and unmanned robot
control. Conventionally, momentum and spin reveal fundamental physical
properties about propagating waves. Yet, vast previous research focused on the
orbital angular momentum of acoustics without thinking about the existence
possibility of spin due to the longitudinal wave nature. Here, we show that
underwater acoustic wave processes the nontrivial spin angular momentum
intrinsically, which is associated with its special spinorbital coupling
relation for longitudinal waves. Furthermore, we demonstrate that this
intrinsic spin, although unobservable in plane wave form, can be detected by
four approaches: wave interference, Gaussian exponential decay form, boundary
evanescent wave, and waveguide mode. We further show that the strong
spinorbital coupling can be exploited to achieve unidirectional excitation and
backscattering immune transport. We hope the present results can improve the
geometric and topological understanding about underwater acoustic wave and pave
the way on the spinrelated underwater applications.

Tellurium is a highperformance thermoelectric material due to its superior
electronic transport and low lattice thermal conductivity ($\kappa_L$). Here,
we report the ultralow $\kappa_L$ in the monolayer tellurium, i.e., tellurene,
which has been successfully synthesized in recent experiments. We find
tellurene has a compellingly low room temperature $\kappa_L$ of 2.16 and 4.08 W
m$^{1}$ K$^{1}$ along the armchair and zigzag directions, respectively, which
is lower than any reported values for other 2D materials. We attribute this
unusually low $\kappa_L$ to the soft acoustic modes, extremely lowenergy
optical modes and the strong scattering among opticalacoustic phonons, which
place tellurene as a potential novel thermoelectric material. Finally, we
disclose that $\kappa_L$ is proportional to the largest acoustic phonon
frequency ($\omega_{D}^{a}$) and the lowest optical phonon frequency at
$\Gamma$ point ($\omega_{\Gamma}^{o}$) in 2D materials, which reflect both
harmonic and anharmonic thermal properties respectively.

Silicon dioxide or silica, normally existing in various bulk crystalline and
amorphous forms, is recently found to possess a twodimensional structure. In
this work, we use ab initio calculation and evolutionary algorithm to unveil
three new 2D silica structures whose themal, dynamical and mechanical
stabilities are compared with many typical bulk silica. In particular, we find
that all these three 2D silica have large inplane negative Poisson's ratios
with the largest one being double of pentagraphene and three times of
borophenes. The negative Poisson's ratio originates from the interplay of
lattice symmetry and SiO tetrahedron symmetry. Slab silica is also an
insulating 2D material, with the highest electronic band gap (> 7 eV) among
reported 2D structures. These exotic 2D silica with inplane negative Poisson's
ratios and widest band gaps are expected to have great potential applications
in nanomechanics and nanoelectronics.

Unveiling intrinsic spins of propagating waves usually offers people a
fundamental understanding of the geometrical and topological properties of
waves from classical to quantum aspects. A great variety of research has shown
that transverse waves can possess nontrivial quantum spins and topology
without help of strong wavematter interaction. However, until now we still
lack essential physical insights about the spin and topological nature of
longitudinal waves. Here, demonstrated by elastic waves we uncover unique
quantum spins for longitudinal waves and the mixed longitudinaltransverse
waves that play essential roles in topological spinmomentum locking. Based on
this quantum spin perspective, several abnormal phenomena beyond pure
transverse waves are attributed to the hybrid spin induced by mixed
longitudinaltransverse waves. The intrinsic hybrid spin reveals the complex
spin essence in elastic waves and advances our understanding about their
fundamental topological properties. We also show these spindependent phenomena
can be exploited to control the wave propagation, such as nonsymmetric elastic
wave excitation by spin pairs, unidirectional Rayleigh wave and spinselected
elastic wave routing. These findings are generally applicable for arbitrary
waves with longitudinal and transverse components.

We carry out a direct comparison of experimental and numerical realizations
of the exact same granular system as it undergoes shear jamming. We adjust the
numerical methods used to optimally represent the experimental settings and
outcomes up to microscopic contact force dynamics. Measures presented here
range form microscopic, through mesoscopic to systemwide characteristics of
the system. Topological properties of the mesoscopic force networks provide a
key link between micro and macro scales. We report two main findings: the
number of particles in the packing that have at least two contacts is a good
predictor for the mechanical state of the system, regardless of strain history
and packing density. All measures explored in both experiments and numerics,
including stress tensor derived measures and contact numbers depend in a
universal manner on the fraction of nonrattler particles, $f_{NR}$. The force
network topology also tends to show this universality, yet the shape of the
master curve depends much more on the details of the numerical simulations. In
particular we show that adding force noise to the numerical data set can
significantly alter the topological features in the data. We conclude that both
$f_{NR}$ and topological metrics are useful measures to consider when
quantifying the state of a granular system.

Directional design of functional materials with multiobjective constraints
is a big challenge, whose performance and stability are determined by different
physics factors entangled with each other complicatedly. In this work, we apply
the multiobjective optimization based on the Pareto Efficiency and
ParticleSwarm Optimization methods to design new functional materials
directionally. As a demonstration, we achieve the thermoelectric design of 2D
SnSe materials through the methods. We identify several novel metastable 2D
SnSe structures with simultaneously lower free energy and better thermoelectric
performance over the experimentallyreported monolayer structures. We hope our
results about the multiobjective Pareto Optimization method can make a step
towards the integrative design of multiobjective and multifunctional
materials in the future.

Web browsing is an activity that billions of mobile users perform on a daily
basis. Battery life is a primary concern to many mobile users who often find
their phone has died at most inconvenient times. The heterogeneous multicore
architecture is a solution for energyefficient processing. However, the
current mobile web browsers rely on the operating system to exploit the
underlying hardware, which has no knowledge of individual web contents and
often leads to poor energy efficiency. This paper describes an automatic
approach to render mobile web workloads for performance and energy efficiency.
It achieves this by developing a machine learning based approach to predict
which processor to use to run the web rendering engine and at what frequencies
the processors should operate. Our predictor learns offline from a set of
training web workloads. The built predictor is then integrated into the browser
to predict the optimal processor configuration at runtime, taking into account
the web workload characteristics and the optimisation goal: whether it is load
time, energy consumption or a tradeoff between them. We evaluate our approach
on a representative ARM big.LITTLE mobile architecture using the hottest 500
webpages. Our approach achieves 80% of the performance delivered by an ideal
predictor. We obtain, on average, 45%, 63.5% and 81% improvement respectively
for load time, energy consumption and the energy delay product, when compared
to the Linux heterogeneous multiprocessing scheduler.

We investigate the rectification and negative differential effects of
fullcounting statistics of conjugate thermal spin transport in a spin Seebeck
engine. The engine is made by an electron/magnon interface diode driven by a
temperature bias, of which the scaled cumulant generating function is
formulated in the framework of nonequilibrium Green's function. The strongly
fluctuating interfacial electron density of states induced by quantum dot is
responsible for these intriguing properties. This work is relevant for
designing efficient spin caloritronic engine and can broaden our view in
nonequilibrium thermodynamics and the nonlinear phenomenon in quantum transport
system.

Transient dynamics of heat conduction in isotropic fractal metamaterials is
investigated. By using the Laplacian operator in noninteger dimension, we
analytically and numerically study the effect of fractal dimensionality on the
evolution of the temperature profile, heat flux and excess energy under certain
initial and boundary conditions. Particularly, with randomly distributed
absorbing heat sinks in the fractal metamaterials, we obtain an anomalous
nonexponential decay behavior of the heat pulse diffusion. and an optimal
dimension for efficient heat absorption as a function of sink concentrations.
Our results may have potential applications in controlling transient heat
conduction in fractal media, which will be ubiquitous as porous, composite,
networked, hierarchical metamaterials.

Investigation of wellmotivated parameter space in the theories of Beyond the
Standard Model (BSM) plays an important role in new physics discoveries.
However, a largescale exploration of models with multiparameter or equivalent
solutions with a finite separation, such as supersymmetric models, is typically
a timeconsuming and challenging task. In this paper, we propose a
selfexploration method, named Machine Learning Scan (MLS), to achieve an
efficient test of models. As a proofofconcept, we apply MLS to investigate
the subspace of MSSM and CMSSM and find that such a method can reduce the
computational cost and may be helpful for accelerating the exploration of
supersymmetry.

Deep learning, e.g., convolutional neural networks (CNNs), has achieved great
success in image processing and computer vision especially in high level vision
applications such as recognition and understanding. However, it is rarely used
to solve lowlevel vision problems such as image compression studied in this
paper. Here, we move forward a step and propose a novel compression framework
based on CNNs. To achieve highquality image compression at low bit rates, two
CNNs are seamlessly integrated into an endtoend compression framework. The
first CNN, named compact convolutional neural network (ComCNN), learns an
optimal compact representation from an input image, which preserves the
structural information and is then encoded using an image codec (e.g., JPEG,
JPEG2000 or BPG). The second CNN, named reconstruction convolutional neural
network (RecCNN), is used to reconstruct the decoded image with highquality in
the decoding end. To make two CNNs effectively collaborate, we develop a
unified endtoend learning algorithm to simultaneously learn ComCNN and
RecCNN, which facilitates the accurate reconstruction of the decoded image
using RecCNN. Such a design also makes the proposed compression framework
compatible with existing image coding standards. Experimental results validate
that the proposed compression framework greatly outperforms several compression
frameworks that use existing image coding standards with stateoftheart
deblocking or denoising postprocessing methods.

The Lieb Lattice exhibits intriguing properties that are of general interest
in both the fundamental physics and practical applications. Here, we
investigate the topological LandauZener Bloch oscillation in a photonic
Floquet Lieb lattice, where the dimerized helical waveguides is constructed to
realize the synthetic spinorbital interaction through the Floquet mechanism,
rendering us to study the impacts of topological transition from trivial gaps
to nontrivial ones. The compact localized states of flat bands supported by
the local symmetry of Lieb lattice will be associated with other bands by
topological invariants, Chern number, and involved into LandauZener transition
during Bloch oscillation. Importantly, the nontrivial geometrical phases after
topological transitions will be taken into account for constructive and
destructive interferences of wave functions. The numerical calculations of
continuum photonic medium demonstrate reasonable agreements with theoretical
tightbinding model. Our results provide an ongoing effort to realize designed
quantum materials with tailored properties.

Building on [1], we examine a holographic model in which a U(1) symmetry and
translational invariance are broken spontaneously at the same time. The
symmetry breaking is realized through the St\"{u}ckelberg mechanism, and leads
to a scalar condensate and a charge density which are spatially modulated and
exhibit unidirectional stripe order. Depending on the choice of parameters, the
oscillations of the scalar condensate can average out to zero, with a frequency
which is half of that of the charge density. In this case the system realizes
some of the key features of pair density wave order. The model also admits a
phase with coexisting superconducting and charge density wave orders, in which
the scalar condensate has a uniform component. In our construction the various
orders are intertwined with each other and have a common origin. The fully
backreacted geometry is computed numerically, including for the case in which
the theory contains axions. The latter can be added to explicitly break
translational symmetry and mimic latticetype effects.

Solids are distinguished from fluids by their ability to resist shear. In
traditional solids, the resistance to shear is associated with the emergence of
broken translational symmetry as exhibited by a nonuniform density pattern,
which results from either minimizing the energy cost or maximizing the entropy
or both. In this work, we focus on a class of systems, where this paradigm is
challenged. We show that sheardriven jamming in dry granular materials is a
collective process controlled solely by the constraints of mechanical
equilibrium. We argue that these constraints lead to a broken translational
symmetry in a dual space that encodes the statistics of contact forces and the
topology of the contact network. The shearjamming transition is marked by the
appearance of this broken symmetry. We extend our earlier work, by comparing
and contrasting real space measures of rheology with those obtained from the
dual space. We investigate the structure and behavior of the dual space as the
system evolves through the rigidity transition in two different shear
protocols. We analyze the robustness of the shearjamming scenario with respect
to protocol and packing fraction, and demonstrate that it is possible to define
a protocolindependent order parameter in this dual space, which signals the
onset of rigidity.

A key issue in the control of distributed discrete systems modeled as Markov
decisions processes, is that often the state of the system is not directly
observable at any single location in the system. The participants in the
control scheme must share information with one another regarding the state of
the system in order to collectively make informed control decisions, but this
information sharing can be costly. Harnessing recent results from information
theory regarding distributed function computation, in this paper we derive, for
several information sharing model structures, the minimum amount of control
information that must be exchanged to enable local participants to derive the
same control decisions as an imaginary omniscient controller having full
knowledge of the global state. Incorporating consideration for this amount of
information that must be exchanged into the reward enables one to trade the
competing objectives of minimizing this control information exchange and
maximizing the performance of the controller. An alternating optimization
framework is then provided to help find the efficient controllers and messaging
schemes. A series of running examples from wireless resource allocation
illustrate the ideas and design tradeoffs.

We report the first result on Ge76 neutrinoless double beta decay from
CDEX1 experiment at China Jinping Underground Laboratory. A mass of 994 g
ptype pointcontact high purity germanium detector has been installed to
search the neutrinoless double beta decay events, as well as to directly detect
dark matter particles. An exposure of 304 kg*day has been analyzed. The
wideband spectrum from 500 keV to 3 MeV was obtained and the average event rate
at the 2.039 MeV energy range is about 0.012 count per keV per kg per day. The
halflife of Ge76 neutrinoless double beta decay has been derived based on
this result as: T 1/2 > 6.4*10^22 yr (90% C.L.). An upper limit on the
effective Majorananeutrino mass of 5.0 eV has been achieved. The possible
methods to further decrease the background level have been discussed and will
be pursued in the next stage of CDEX experiment.

We have developed a method to improve the doping computation efficiency, this
method is based on first principles calculations and cluster expansion. First
principles codes produce highly accurate total energies and optimized
geometries for any given structure. Cluster expansion method constructs a
cluster expansion using partial first principles results and computes the
energies for other structures derived from a parent lattice. Using this method,
energies for multiple doping structures can be predicted quickly without series
of first principles calculations. This method has been packaged into a tool
named as CE Screen and integrated into MatCloud (A highthroughput first
principles calculation platform). This makes the tool simple and easy for all
the users.