
Quantum oscillations are usually the manifestation of the underlying physical
nature in condensed matter systems. Here we report a new type of logperiodic
quantum oscillations in ultraquantum threedimensional topological materials.
Beyond the quantum limit (QL), the logperiodic oscillations involving up to
five oscillating cycles (5 peaks and 5 dips) are observed on the
magnetoresistance (MR) of high quality singlecrystal ZrTe5, virtually showing
the clearest feature of discrete scale invariance (DSI). Further theoretical
analyses show that the twobody quasibound states can be responsible for the
DSI feature. Our work provides a new perspective on the ground state of
topological materials beyond the QL.

We study the Combinatorial Pure Exploration problem with Continuous and
Separable reward functions (CPECS) in the stochastic multiarmed bandit
setting. In a CPECS instance, we are given several stochastic arms with
unknown distributions, as well as a collection of possible decisions. Each
decision has a reward according to the distributions of arms. The goal is to
identify the decision with the maximum reward, using as few arm samples as
possible. The problem generalizes the combinatorial pure exploration problem
with linear rewards, which has attracted significant attention in recent years.
In this paper, we propose an adaptive learning algorithm for the CPECS
problem, and analyze its sample complexity. In particular, we introduce a new
hardness measure called the consistent optimality hardness, and give both the
upper and lower bounds of sample complexity. Moreover, we give examples to
demonstrate that our solution has the capacity to deal with nonlinear reward
functions.

The central engine of gammaray bursts (GRBs) is poorly constrained. There
exist two main candidates: a fastrotating black hole and a rapidly spinning
magnetar. Furthermore, the Xray plateaus are widely accepted by the energy
injection into the external shock. In this paper, we systematically analyze the
\emph{Swift}/XRT light curves of 101 GRBs having plateau phases and known
redshifts (before May 2017). Since a maximum energy budget
($\sim2\times10^{52}$ erg) exists for magnetars but not for a black hole, this
provides a good clue to identify the type of GRB central engine. We calculate
the isotropic kinetic energy $E_{\rm K,iso}$ and the isotropic Xray energy
release $E_{\rm X,iso}$ for individual GRB. We identify three categories based
on how likely a black hole harbor at central engine: 'Gold' (9 out of 101, both
$E_{\rm X,iso}$ and $E_{\rm K,iso}$ exceed the energy budget), 'Silver' (69 out
of 101, $E_{\rm X,iso}$ less than the limit but $E_{\rm K,iso}$ is greater than
the limit), and 'Bronze' (23 out of 101, the energies are not above the limit).
We then derive and test the black hole parameters with the Blandford Znajek
mechanism, and find that the observations of the black hole candidate
('Gold'+'Silver') samples are consistent with the expectations of the black
hole model. Furthermore, we also test the magnetar candidate ('Bronze') sample
with the magnetar model, and find that the magnetar surface magnetic led
($B_{p}$) and initial spin period ($P_{0}$) are consistent with the
expectations of the magnetar model. Our analysis indicates that, if the
magnetar wind is isotropic, a magnetar central engine is possible for 20\% of
the analyzed GRBs. For most GRBs a black hole is most likely operating.

The spectrum of gammaray burst (GRB) afterglows can be studied with color
indices. Here we present a large comprehensive catalogue of 70 GRBs with
multiwavelength optical transient data on which we perform a systematic study
to find the temporal evolution of color indices. We categorize them into two
samples based on how well the color indices are evaluated. The Golden sample
includes 25 bursts mostly observed by GROND, and the Silver sample includes 45
bursts observed by other telescopes. For the Golden sample, we find that 95\%
of the color indices do not vary over time. However, the color indices do vary
during short periods in most bursts. The observed variations are consistent
with effects of (i) the cooling frequency crossing the studied energy bands in
a wind medium (43\%) and in a constant density medium (30\%), (ii) early dust
extinction (12\%), (iii) transition from reverse shock to forward shock
emission (5\%), or (iv) an emergent supernova emission (10\%). We also study
the evolutionary properties of the mean color indices for different emission
episodes. We find that 86\% of the color indices in the 70 bursts show
constancy between consecutive ones. The color index variations occur mainly
during the late GRBSN bump, the flare and early reversedshock emission
components. We further perform a statistical analysis of various observational
properties and model parameters (spectral index $\beta_{o}^{\rm CI}$, electron
spectral indices $p^{\rm CI}$, etc.) using color indices. Overall, we conclude
that $\sim$ 90\% of colors are constant in time and can be accounted for by the
simplest external forward shock model, while the varying color indices call for
more detailed modeling.

Nematic order often breaks the tetragonal symmetry of ironbased
superconductors. It arises from regular structural transition or electronic
instability in the normal phase. Here, we report the observation of a nematic
superconducting state, by measuring the angular dependence of the inplane and
outofplane magnetoresistivity of Ba0.5K0.5Fe2As2 single crystals. We find
large twofold oscillations in the vicinity of the superconducting transition,
when the direction of applied magnetic field is rotated within the basal plane.
To avoid the influences from sample geometry or current flow direction, the
sample was designed as Corbinoshape for inplane and mesashape for
outofplane measurements. Theoretical analysis shows that the nematic
superconductivity arises from the weak mixture of the quasidegenerate swave
and dwave components of the superconducting condensate, most probably induced
by a weak anisotropy of stresses inherent to single crystals.

Photon channel perspective on high harmonic generation (HHG) is proposed by
quantizing both the driving laser and high harmonics. It is shown that the HHG
yield can be expressed as a sum of the contribution of all the photon channels.
From this perspective, the contribution of a specific photon channel follows a
simple scaling law and the competition between the channels is well
interpreted. Our prediction is shown to be in good agreement with the
simulations by solving the timedependent Schrodinger equation. It also can
well explains the experimental results of the HHG in the noncollinear twocolor
field and bicicular laser field.

We report a new scenario of timeofflight (TOF) technique in which fast
neutrons and delayed gammaray signals were both recorded in a millisecond time
window in harsh environments induced by highintensity lasers. The delayed
gamma signals, arriving far later than the original fast neutron and often
being ignored previously, were identified to be the results of radiative
captures of thermalized neutrons. The linear correlation between gamma photon
number and the fast neutron yield shows that these delayed gamma events can be
employed for neutron diagnosis. This method can reduce the detecting efficiency
dropping problem caused by prompt highflux gamma radiation, and provides a new
way for neutron diagnosing in highintensity lasertarget interaction
experiments.

A tight criterion under which the abstract version Lov\'asz Local Lemma
(abstractLLL) holds was given by Shearer decades ago. However, little is known
about that of the variable version LLL (variableLLL) where events are
generated by independent random variables, though this model of events is
applicable to almost all applications of LLL. We introduce a necessary and
sufficient criterion for variableLLL, in terms of the probabilities of the
events and the eventvariable graph specifying the dependency among the events.
Based on this new criterion, we obtain boundaries for two families of
eventvariable graphs, namely, cyclic and treelike bigraphs. These are the
first two nontrivial cases where the variableLLL boundary is fully
determined. As a byproduct, we also provide a universal constructive method to
find a set of events whose union has the maximum probability, given the
probability vector and the eventvariable graph. Though it is #Phard in
general to determine variableLLL boundaries, we can to some extent decide
whether a gap exists between a variableLLL boundary and the corresponding
abstractLLL boundary. In particular, we show that the gap existence can be
decided without solving Shearer's conditions or checking our variableLLL
criterion. Equipped with this powerful theorem, we show that there is no gap if
the base graph of the eventvariable graph is a tree, while gap appears if the
base graph has an induced cycle of length at least 4. The problem is almost
completely solved except when the base graph has only 3cliques, in which case
we also get partial solutions. A set of reduction rules are established that
facilitate to infer gap existence of an eventvariable graph from known ones.
As an application, various eventvariable graphs, in particular combinatorial
ones, are shown to be gapful/gapless.

This paper focuses on quantifying and estimating the predictive accuracy of
prognostic models for timetoevent outcomes with competing events. We consider
the timedependent discrimination and calibration metrics, including the
receiver operating characteristics curve and the Brier score, in the context of
competing risks. To address censoring, we propose a unified nonparametric
estimation framework for both discrimination and calibration measures, by
weighting the censored subjects with the conditional probability of the event
of interest given the observed data. We demonstrate through simulations that
the proposed estimator is unbiased, efficient and robust against model
misspecification in comparison to other methods published in the literature. In
addition, the proposed method can be extended to timedependent predictive
accuracy metrics constructed from a general class of loss functions. We apply
the methodology to a data set from the African American Study of Kidney Disease
and Hypertension to evaluate the predictive accuracy of a prognostic risk score
in predicting endstage renal disease (ESRD), accounting for the competing risk
of preESRD death.

The spacetime structure of the Higgs boson decay are carefully studied with
the HBT correlation method using e$^+$e$^$ collision events produced through
Monte Carlo generator PYTHIA 8.2 at $\sqrt{s}$=250GeV. The Higgs boson jets
(Higgsjets) are identified by Htag tracing. The measurement of the Higgs
boson radius and decay lifetime are derived from HBT correlation of its decay
final state pions inside Higgsjets in the e$^+$e$^$ collisions events with an
upper bound of $R_H \le 1.03\pm 0.05$ fm and $\tau_H \le (1.29\pm0.15)\times
10^{7}$ fs. This result is consistent with CMS data.

In many clinical trials studying neurodegenerative diseases such as
Parkinson's disease (PD), multiple longitudinal outcomes are collected to fully
explore the multidimensional impairment caused by this disease. If the outcomes
deteriorate rapidly, patients may reach a level of functional disability
sufficient to initiate levodopa therapy for ameliorating disease symptoms. An
accurate prediction of the time to functional disability is helpful for
clinicians to monitor patients' disease progression and make informative
medical decisions. In this article, we first propose a joint model that
consists of a semiparametric multilevel latent trait model (MLLTM) for the
multiple longitudinal outcomes, and a survival model for event time. The two
submodels are linked together by an underlying latent variable. We develop a
Bayesian approach for parameter estimation and a dynamic prediction framework
for predicting target patients' future outcome trajectories and risk of a
survival event, based on their multivariate longitudinal measurements. Our
proposed model is evaluated by simulation studies and is applied to the DATATOP
study, a motivating clinical trial assessing the effect of deprenyl among
patients with early PD.

Weyl nodes are topological objects in threedimensional metals. Their
topological property can be revealed by studying the highfield transport
properties of a Weyl semimetal. While the energy of the lowest Landau band
(LLB) of a conventional Fermi pocket always increases with magnetic field due
to the zero point energy, the LLB of Weyl cones remains at zero energy unless a
strong magnetic field couples the Weyl fermions of opposite chirality. In the
Weyl semimetal TaP, we achieve such a magnetic coupling between the
electronlike Fermi pockets arising from the W1 Weyl fermions. As a result,
their LLBs move above chemical potential, leading to a sharp sign reversal in
the Hall resistivity at a specific magnetic field corresponding to the W1 Weyl
node separation. By contrast, despite having almost identical carrier density,
the annihilation is unobserved for the holelike pockets because the W2 Weyl
nodes are much further separated. These key findings, corroborated by other
systematic analyses, reveal the nontrivial topology of Weyl fermions in
highfield measurements.

Semantic parsing has emerged as a significant and powerful paradigm for
natural language interface and question answering systems. Traditional methods
of building a semantic parser rely on highquality lexicons, handcrafted
grammars and linguistic features which are limited by applied domain or
representation. In this paper, we propose a general approach to learn from
denotations based on Seq2Seq model augmented with attention mechanism. We
encode input sequence into vectors and use dynamic programming to infer
candidate logical forms. We utilize the fact that similar utterances should
have similar logical forms to help reduce the searching space. Under our
learning policy, the Seq2Seq model can learn mappings gradually with noises.
Curriculum learning is adopted to make the learning smoother. We test our
method on the arithmetic domain which shows our model can successfully infer
the correct logical forms and learn the word meanings, compositionality and
operation orders simultaneously.

Recently, two families of HSSbased iteration methods are constructed for
solving the system of absolute value equations (AVEs), which is a class of
nondifferentiable NPhard problems. In this study, we establish the
PicardCSCS iteration method and the nonlinear CSCSlike iteration method for
AVEs involving the Toeplitz matrix. Then, we analyze the convergence of the
PicardCSCS iteration method for solving AVEs. By using the theory about
nonsmooth analysis, we particularly prove the convergence of the nonlinear
CSCSlike iterationsolver for AVEs. The advantage of these methods is that they
do not require the storage of coefficient matrices at all, and the subsystem
of linear equations can be solved efficiently via the fast Fourier transforms
(FFTs). Therefore, computational cost and storage can be saved in practical
implementations. Numerical examples including numerical solutions of nonlinear
fractional diffusion equations are reported to show the effectiveness of the
proposed methods in comparison with some existing methods.

We report measurements of the upper critical field $H_{\mathrm{c2}}$ as
functions of temperature $T$, polar angle $\theta$ (of the field direction with
respect to the crystallographic $c$ axis), and azimuthal angle $\phi$ (of the
field direction relative to the $a$ axis within the $ab$ plane) for the
Crbased superconductor K$_{2}$Cr$_{3}$As$_{3}$ with a quasionedimensional
and noncentrosymmetric crystal structure. We confirm that the anisotropy in
$H_{\mathrm{c2}}(T)$ becomes inverse with decreasing temperature. At low
temperatures, $H_{\mathrm{c2}}(\theta)$ data are featured by two maxima at
$\theta$ = 0 ($\mathbf{H}\parallel c$) and $\pi/2$ ($\mathbf{H}\bot c$), which
can be quantitatively understood only if uniaxial effectivemass anisotropy and
absence of Pauli paramagnetic effect for $\mathbf{H}\bot c$ are taken
simultaneously into consideration. The inplane $H_{\mathrm{c2}}(\phi)$ profile
shows a unique threefold modulation especially at low temperatures. Overall,
the characteristic of the $H_{\mathrm{c2}}(\theta, \phi, T)$ data mostly
resemble those of the heavyfermion superconductor UPt$_3$, and we argue in
favor of a dominant spintriplet superconductivity with odd parity in
K$_{2}$Cr$_{3}$As$_{3}$.

Public opinion polling is usually done by random sampling from the entire
population, treating individual opinions as independent. In the real world,
individuals' opinions are often correlated, e.g., among friends in a social
network. In this paper, we explore the idea of partitioned sampling, which
partitions individuals with high opinion similarities into groups and then
samples every group separately to obtain an accurate estimate of the population
opinion. We rigorously formulate the above idea as an optimization problem. We
then show that the simple partitions which contain only one sample in each
group are always better, and reduce finding the optimal simple partition to a
wellstudied MinrPartition problem. We adapt an approximation algorithm and a
heuristic algorithm to solve the optimization problem. Moreover, to obtain
opinion similarity efficiently, we adapt a wellknown opinion evolution model
to characterize social interactions, and provide an exact computation of
opinion similarities based on the model. We use both synthetic and realworld
datasets to demonstrate that the partitioned sampling method results in
significant improvement in sampling quality and it is robust when some opinion
similarities are inaccurate or even missing.

We propose Hybridizable Discontinuous Galerkin (HDG) methods for solving the
frequencydomain Maxwell's equations coupled to the Nonlocal Hydrodynamic Drude
(NHD) and Generalized Nonlocal Optical Response (GNOR) models, which are
employed to describe the optical properties of nanoplasmonic scatterers and
waveguides. Brief derivations for both the NHD model and the GNOR model are
presented. The formulations of the HDG method are given, in which we introduce
two hybrid variables living only on the skeleton of the mesh. The local field
solutions are expressed in terms of the hybrid variables in each element. Two
conservativity conditions are globally enforced to make the problem solvable
and to guarantee the continuity of the tangential component of the electric
field and the normal component of the current density. Numerical results show
that the proposed HDG methods converge at optimal rate. We benchmark our
implementation and demonstrate that the HDG method has the potential to solve
complex nanophotonic problems.

Nuclear fusion reactions are the most important processes in nature to power
stars and produce new elements, and lie at the center of the understanding of
nucleosynthesis in the universe. It is critically important to study the
reactions in full plasma environments that are close to true astrophysical
conditions. By using laserdriven counterstreaming collisionless plasmas, we
studied the fusion D$+$D$\rightarrow n +^3$He in a Gamowlike window around 27
keV. The results show that astrophysical nuclear reaction yield can be
modulated significantly by the selfgenerated electromagnetic fields and the
collective motion of the plasma. This plasmaversion minicollider may provide
a novel tool for studies of astrophysicsinterested nuclear reactions in plasma
with tunable energies in earthbased laboratories.

Emergency of superconductivity at the instabilities of antiferromagnetism has
been widely recognized in unconventional superconductors. In copperoxide
superconductors, spin fluctuations play a predominant role in electron pairing
with electron dopants yet composite orders veil the nature of superconductivity
for holedoped family. However, in electrondoped copper oxide superconductors
(cuprates) the AFM critical end point is still in controversy for different
probes, demonstrating high sensitivity to oxygen content. Here, by carefully
tuning the oxygen content, a systematic study of the Hall signal and
magnetoresistivity up to 58 Tesla on LCCO thin films identifies two
characteristic temperatures. The former is quite robust, whereas the latter
becomes flexible with increasing magnetic field, thereby linking respectively
to two and threedimensional AFM, evident from the multidimensional phase
diagram as a function of oxygen and Ce dopants. A rigorous theoretical analysis
of the presented data suggest the existence of conductive nanofilamentary
structures that effectively corroborate all previously reported field studies.
The new findings provide a uniquely consistent alternative picture in
understanding the interactions between AFM and superconductivity in
electrondoped cuprates and offer a consolidating interpretation to the
pioneering scaling law in cuprates recently established by Bozovic et al.
(Nature, 2016)

We investigate the polarization properties of high harmonics generated with
the bichromatic counterrotating circularly polarized (BCCP) laser fields by
numerically solving timedependent Schr\"odinger equation (TDSE). It is found
that, the helicity of the elliptically polarized harmonic emission is reversed
at particular harmonic orders. Based on the timefrequency analysis and the
classical threestep model, the correspondence between the positions of
helicity reversions and the classical trajectories of continuum electrons is
established. It is shown that, the electrons ionized at one lobe of laser field
can be divided into different groups based on the different lobes they
recombine at, and the harmonics generated by adjacent groups have opposite
helicities. Our study performs a detailed analysis of high harmonics in terms
of electron trajectories and depicts a clear and intuitive physical picture of
the HHG process in BCCP laser field.

The selection rules of high harmonic generation (HHG) are investigated using
threedimensional timedependent density functional theory (TDDFT). From the
harmonic spectra obtained with various real molecules and different forms of
laser fields, several factors that contribute to selection rules are revealed.
Extending the targets to stereoscopic molecules, it is shown that the allowed
harmonics are dependent on the symmetries of the projections of molecules. For
laser fields, the symmetries contributing to the selection rules are discussed
according to Lissajous figures and their dynamical directivities. All the
phenomena are explained by the symmetry of the full timedependent Hamiltonian
under a combined transformation. We present a systematic study on the selection
rules and propose an intuitive method for the judgment of allowed harmonic
orders, which can be extended to more complex molecules and various forms of
laser pulses.

Limited angle problem is a challenging issue in xray computed tomography
(CT) field. Iterative reconstruction methods that utilize the additional prior
can suppress artifacts and improve image quality, but unfortunately require
increased computation time. An interesting way is to restrain the artifacts in
the images reconstructed from the practical filtered back projection (FBP)
method. Frikel and Quinto have proved that the streak artifacts in FBP results
could be characterized. It indicates that the artifacts created by FBP method
have specific and similar characteristics in a stationary limitedangle
scanning configuration. Based on this understanding, this work aims at
developing a method to extract and suppress specific artifacts of FBP
reconstructions for limitedangle tomography. A datadriven learningbased
method is proposed based on a deep convolutional neural network. An endtoend
mapping between the FBP and artifactfree images is learned and the implicit
features involving artifacts will be extracted and suppressed via nonlinear
mapping. The qualitative and quantitative evaluations of experimental results
indicate that the proposed method show a stable and prospective performance on
artifacts reduction and detail recovery for limited angle tomography. The
presented strategy provides a simple and efficient approach for improving image
quality of the reconstruction results from limited projection data.

In this paper, we study information cascades on graphs. In this setting, each
node in the graph represents a person. One after another, each person has to
take a decision based on a private signal as well as the decisions made by
earlier neighboring nodes. Such information cascades commonly occur in practice
and have been studied in complete graphs where everyone can overhear the
decisions of every other player. It is known that information cascades can be
fragile and based on very little information, and that they have a high
likelihood of being wrong.
Generalizing the problem to arbitrary graphs reveals interesting insights. In
particular, we show that in a random graph $G(n,q)$, for the right value of
$q$, the number of nodes making a wrong decision is logarithmic in $n$. That
is, in the limit for large $n$, the fraction of players that make a wrong
decision tends to zero. This is intriguing because it contrasts to the two
natural corner cases: empty graph (everyone decides independently based on his
private signal) and complete graph (all decisions are heard by all nodes). In
both of these cases a constant fraction of nodes make a wrong decision in
expectation. Thus, our result shows that while both too little and too much
information sharing causes nodes to take wrong decisions, for exactly the right
amount of information sharing, asymptotically everyone can be right. We further
show that this result in random graphs is asymptotically optimal for any
topology, even if nodes follow a globally optimal algorithmic strategy. Based
on the analysis of random graphs, we explore how topology impacts global
performance and construct an optimal deterministic topology among layer graphs.

We apply the concept of artificial and controlled interference in a
triangular relay network with an untrusted relay aiming at enhancing the
wireless communication secrecy between the source and the destination node. In
order to shield the square quadrature amplitude modulated (QAM) signals
transmitted from the source node to the relay, the destination node designs and
transmits artificial noise (AN) symbols to jam the relay reception. The
objective of our considered AN design is to degrade the error probability
performance at the untrusted relay, for different types of channel state
information (CSI) at the destination. By considering perfect knowledge of the
instantaneous CSI of the sourcetorelay and relaytodestination links, we
first present an analytical expression for the symbol error rate (SER)
performance at the relay. Based on the assumption of an average power
constraint at the destination node, we then derive the optimal phase and power
distribution of the AN that maximizes the SER at the relay. Furthermore, we
obtain the optimal AN design for the case where only statistical CSI is
available at the destination node. For both cases, our study reveals that the
Gaussian distribution is generally not optimal to generate AN symbols. In
particular, the conducted numerical results show that the Gaussian distributed
AN is far from optimal. The presented AN design takes into account practical
parameters for the communication links, such as QAM signaling and maximum
likelihood decoding.

HfTe5 is predicted to be a promising platform for studying topological
phases. Here through an electrical transport study, we present the first
observation of chiral anomaly and ultrahigh mobility in HfTe5 crystals.
Negative magnetoresistivity in HfTe5 is observed when the external magnetic
and electrical fields are parallel (B//E) and quickly disappears once B
deviates from the direction of E. Quantitative fitting further confirms the
chiral anomaly as the underlying physics. Moreover, by analyzing the
conductivity tensors of longitudinal and Hall traces, ultrahigh mobility and
ultralow carrier density are revealed in HfTe5, which paves the way for
potential electronic applications.