
The problem of detecting a few anomalous processes among a large number of
data streams is considered. At each time, aggregated observations can be taken
from a chosen subset of the processes, where the chosen subset conforms to a
given tree structure. The random observations are drawn from a general
distribution that may depend on the size of the chosen subset and the number of
anomalous processes in the subset. We propose a sequential search strategy by
devising an informationdirected random walk (IRW) on the treestructured
observation hierarchy. Subject to a reliability constraint, the proposed policy
is shown to be asymptotically optimal with respect to the detection accuracy.
Furthermore, it achieves the optimal logarithmicorder sample complexity with
respect to the size of the search space provided that the KullbackLeibler
divergence between aggregated observations in the presence and the absence of
anomalous processes are bounded away from zero at all levels of the tree
structure as the size of the search space approaches infinity. Sufficient
conditions on the decaying rate of the aggregated observations to pure noise
under which a sublinear scaling in the size of the search space is preserved
are also identified for the Bernoulli case.

We propose a lattice scale twoband generalized Hubbard model as a caricature
of the electronic structure of twisted bilayer graphene. Various possible
broken symmetry phases can arise, including a nematic phase (which is a form of
orbital ferromagnet) and an orbitaltriplet spinsinglet superconducting phase.
Concerning the mechanism of superconductivity  we propose an analogy with
superconductivity in alkalidoped C$_{60}$ in which a violation of Hund's first
rule plays a central role.

The authors of the comment[Phys. Rev. A 97, 046101 (2018)] raise that the
inconsistency in calculating some common quantumspeedlimit (QSL) bounds,
which is presented in our paper [Phys. Rev. A 95, 052118 (2017)], does not
exist in their paper [Phys. Rev. A 94, 052125 (2016)]. Therefore they insist
that their criticism to some QSL bounds is still valid. We demonstrate all the
QSL bounds mentioned in the comment are similar in essence. We also show the
inconsistency which is presented in their original paper can not be used to
support their criticism. Furthermore, we exhibit although the inconsistency
presented in our paper is not a same one to which is presented by the authors
of the comment, it exists and be unavoidable in their example. Accordingly, we
believe the result in our paper is usable in numerical calculations.

We study the $S$wave $D\Xi$ bound state in the BetheSalpeter formalism in
the ladder and instantaneous approximations. With the kernel generated by the
hadronic effective Lagrangian, two opencharmed bound states, which quantum
numbers are $I=0$, $J^P=(\frac{1}{2})^$ and $I=1$, $J^P=(\frac{1}{2})^$,
respectively, are predicted as new candidates of hadronic pentaquark molecules
in our formalism. If existing, they could contribute to the broad 3188 eV
structure near the five new narrow $\Omega_c$ states observed recently by the
LHCb Collaboration.

Using the Lawson's existence theorem of minimal surfaces and the symmetries
of the Hopf fibration, we will construct symmetric embedded closed minimal
surfaces in the three dimensional sphere. These surfaces contain the Clifford
torus, the Lawson's minimal surfaces, and seven new minimal surfaces with
genera 9, 25, 49, 121, 121, 361 and 841. We will also discuss the relation
between such surfaces and the maximal extendable group actions on subsurfaces
of the three dimensional sphere.

We study the electromagnetic form factors (EMFFs) of $\Lambda_c$ and the
quark and diquark current contributes to the EMFFs of $\Lambda_c$ in the
spacelike (SL) region in the BetheSalpeter equation approach. In this
picture, the heavy baryon $\Lambda_c $ is regarded as composed of a heavy quark
and a scalar diquark. We find that for different values of parameters the quark
and diquark current contribute to the EMFFs of $\Lambda_c$ is very different,
but the total contribute to the EMFFs of $\Lambda_c$ is similarly. The EMFFs of
$\Lambda_c$ are similar to those of other baryons (proton, $\Xi^$, $\Sigma^+$)
with a peak at $\omega =1$ ($\omega=v^\prime \cdot v $ is the velocity transfer
between the initial state (with velocity $v$) and the final state (with
velocity $v^\prime$) of $\Lambda_c$).

We consider the highresolution imaging problem of 3D point source image
recovery from 2D data using a method based on point spread function (PSF)
engineering. The method involves a new technique, recently proposed by
S.~Prasad, based on the use of a rotating PSF with a single lobe to obtain
depth from defocus. The amount of rotation of the PSF encodes the depth
position of the point source. Applications include highresolution single
molecule localization microscopy as well as the problem addressed in this paper
on localization of space debris using a spacebased telescope. The localization
problem is discretized on a cubical lattice where the coordinates of nonzero
entries represent the 3D locations and the values of these entries the fluxes
of the point sources. Finding the locations and fluxes of the point sources is
a largescale sparse 3D inverse problem. A new nonconvex regularization method
with a datafitting term based on KullbackLeibler (KL) divergence is proposed
for 3D localization for the Poisson noise model. In addition, we propose a new
scheme of estimation of the source fluxes from the KL datafitting term.
Numerical experiments illustrate the efficiency and stability of the algorithms
that are trained on a random subset of image data before being applied to other
images. Our 3D localization algorithms can be readily applied to other kinds of
depthencoding PSFs as well.

If a finite group of orientationpreserving diffeomorphisms of the
3dimensional torus leaves invariant an oriented, closed, embedded surface of
genus g>1 and preserves the orientation of the surface, then its order is
bounded from above by 12(g1). In the present paper we classify (up to
conjugation) all such group actions and surfaces for which the maximal possible
order 12(g1) is achieved, and note that the unknotted surfaces can be realized
by equivariant minimal surfaces in a 3torus.

The naturalness of warps is gaining extensive attentions in image stitching.
Recent warps such as SPHP and AANAP, use global similarity warps to mitigate
projective distortion (which enlarges regions), however, they necessarily bring
in perspective distortion (which generates inconsistencies). In this paper, we
propose a novel quasihomography warp, which effectively balances the
perspective distortion against the projective distortion in the nonoverlapping
region to create a more naturallooking panorama. Our approach formulates the
warp as the solution of a bivariate system, where perspective distortion and
projective distortion are characterized as slope preservation and scale
linearization respectively. Because our proposed warp only relies on a global
homography, thus it is totally parameterfree. A comprehensive experiment shows
that a quasihomography warp outperforms some stateoftheart warps in urban
scenes, including homography, AutoStitch and SPHP. A user study demonstrates
that it wins most users' favor, comparing to homography and SPHP.

An ultrafast singlepixel optical 2D imaging system using a single multimode
fiber (MF) is proposed. The MF acted as the alloptical random pattern
generator. Light with different wavelengths pass through a single MF will
generator alloptical random speckle patterns, which have a low correlation of
0.074 with 0.1nm wavelength step from 1518.0nm to 1567.9nm. The alloptical
random speckle patterns are perfect for compressive sensing (CS) imaging with
the advantage of low cost in comparison with the conventional expensive
pseudorandom binary sequence (PRBS). Besides, with the employment of photonic
time stretch (PTS), light of different wavelengths will go through a single
capsuled MF in time serial within a short pulse time, which makes ultrafast
singlepixel alloptical CS imaging possible. In our work, the alloptical
random speckle patterns are analyzed and used to perform CS imaging in our
proposed system and the results shows a singlepixel photodetector can be
employed in CS imaging system and a 27 by 27 pixels image is reconstructed
within 500 measurements. In our proposed imaging system, the fast Fourier
transform (FFT) spatial resolution, which is a combination of multiple
Gaussians, is analyzed. Considering 4 optical speckle patterns, the FFT spatial
resolution is 50 by 50 pixels. This resolution limit has been obtained by
removing the central low frequency components and observing the significant
spectral power along all the radial directions.

The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey is proposed
to obtain deep optical $ugri$ imaging of the CDFS and ES1 fields using the VLT
Survey Telescope (VST). At present, the observations for the CDFS field have
been completed, and comprise in total about 4.9 deg$^2$ down to
$r_\mathrm{AB}$$\sim$26mag. Weak lensing shear is measured in a companion paper
for $r$band images with seeing $\leq$ 0.9 arcsec. In this work, we perform
image simulations to calibrate possible biases of the measured shear signals.
Statistically, the properties of the simulated point spread function (PSF) and
galaxies show good agreement with those of observations. The multiplicative
bias is calibrated to reach an accuracy of 3.0%. We study the bias
sensitivities to undetected faint galaxies and to the existence of neighboring
galaxies. We find that undetected galaxies contribute to 0.3% multiplicative
bias. Further analysis shows that galaxies with lower signaltonoise ratio
(SNR) suffer from more significant impact due to the skewed background noise
caused by undetected galaxies. Meanwhile, we find that although most of the
neighbors have been rejected in the shape measurement procedure, about one
third of them still remain in the final shear sample. They show larger
ellipticity dispersion and contribute to 0.2% multiplicative bias. Simply
rejecting these galaxies can reduce the effective number density of galaxies
for deep surveys, such as VOICE. Therefore efficient methods should be
developed for future weak lensing surveys in order to decompose their
overlapped surface brightness distributions.

In the heavy quark limit, we study the $X(5568)$ state as a $B\bar{K}$
molecule in the BetheSalpeter equation approach. With the kernel containing
oneparticleexchange diagrams, we solve the BetheSalpeter equation
numerically in the covariant instantaneous approximation and find that the
$X(5568)$ can exist as a $B\bar{K}$ molecular state with quantum numbers $I
(J^P) = 1(0^+)$. In this picture we calculate the strong decay width of
$X(5568)\rightarrow B_s^0\pi^+$ and find it to be in the range 19.83  22.45
MeV, which is consistent with the experimental data from the D0 Collaboration.

We report combined magnetic susceptibility, dielectric constant, nuclear
quadruple resonance (NQR) and zerofield nuclear magnetic resonance (NMR)
measurements on single crystals of multiferroics CuBr$_2$. High quality of the
sample is demonstrated by the sharp magnetic and magneticdriven ferroelectric
transition at $T_N=T_C\approx$ 74~K. The zerofield $^{79}$Br and $^{81}$Br NMR
are resolved below $T_N$. The spinlattice relaxation rates reveal charge
fluctuations when cooled below 60~K. Evidences of an increase of NMR linewidth,
a reduction of dielectric constant, and an increase of magnetic susceptibility
are also seen at low temperatures. These data suggest an emergent instability
which competes with the spiral magnetic ordering and the ferroelectricity.
Candidate mechanisms are discussed based on the quasionedimensional (1D)
nature of the magnetic system.

Semantic layouts based Image synthesizing, which has benefited from the
success of Generative Adversarial Network (GAN), has drawn much attention in
these days. How to enhance the synthesis image equality while keeping the
stochasticity of the GAN is still a challenge. We propose a novel denoising
framework to handle this problem. The overlapped objects generation is another
challenging task when synthesizing images from a semantic layout to a realistic
RGB photo. To overcome this deficiency, we include a onehot semantic label map
to force the generator paying more attention on the overlapped objects
generation. Furthermore, we improve the loss function of the discriminator by
considering perturb loss and cascade layer loss to guide the generation
process. We applied our methods on the Cityscapes, Facades and NYU datasets and
demonstrate the image generation ability of our model.

Using the nonequilibrium Green function method, we study the Andreev
reflection in a Yshaped graphenesuperconductor device by tightbinding model.
Considering both the zigzag and armchair terminals, we confirm that the zigzag
terminals are the better choice for detecting the Andreev reflection without no
external field. Due to scattering from the boundaries of the finitesize centre
region, the difference between Andreev retroreflection and specular reflection
is hard to be distinguished. Although adjusting the size of the device makes
the difference visible, to distinguish them quantitatively is still impossible
through the transport conductance. The problem is circumvented when applying a
perpendicular magnetic field on the centre region, which makes the incident
electrons and the reflected holes propagate along the edge or the interface. In
this case, the retroreflected and specular reflected holes from the different
bands have opposite effective masses, therefore the moving direction of one is
opposite to the other. Which external terminal the reflected holes flow into
depends entirely on the kind of the Andreev reflection. Therefore, the specular
Andreev reflection can be clearly distinguished from the retroreflected one in
the presence of strong magnetic field, even for the device with finite size.

The idea of statistical transmutation plays a crucial role in descriptions of
the fractional quantum Hall effect. However, a recently conjectured duality
between a critical boson and a massless 2component Dirac fermion extends this
notion to gapless systems. This duality sheds light on highly nontrivial
problems such as the halffilled Landau level, the superconductorinsulator
transition, and surface states of strongly coupled topological insulators.
Although this bosonfermion duality has undergone many consistency checks, it
has remained unproven. We describe the duality in a nonperturbative fashion
using an exact UV mapping of partition functions on a 3D Euclidean lattice.

We study the threebody radiative decays $B\to \phi(\rho) K\gamma$ induced by
a flavorchanging neutral current in the perturbative QCD approach.
Pseudoscalarvector ($PV$) distribution amplitudes (DAs) are introduced for the
finalstate $\phi K$ ($\rho K$) pair to capture important infrared dynamics in
the region with a small $PV$pair invariant mass. The dependence of these $PV$
DAs on the parton momentum fraction is parametrized in terms of the Gegenbauer
polynomials, and the dependence on the meson momentum fraction is derived
through their normalizations to timelike $PV$ form factors. In addition to the
dominant electromagnetic penguin, the subleading chromomagnetic penguin,
quarkloop and annihilation diagrams are also calculated. After determining the
$PV$ DAs from relevant branchingratio data, the direct $CP$ asymmetries and
decay spectra in the $PV$pair invariant mass are predicted for each $B\to
\phi(\rho) K\gamma$ mode.

We discuss the possibility that the X(3872) can be a $D\bar{D}^*$ molecular
bound state in the BetheSalpeter equation approach in the ladder and
instantaneous approximations. We show that the $D\bar{D}^*$ bound state with
quantum numbers $J^{PC}=1^{++}$ exists. We also calculate the decay width of
$X(3872) \rightarrow \gamma J/\psi$ channel and compare our result with those
from previous calculations.

We propose a memorymodelaware static program analysis method for accurately
analyzing the behavior of concurrent software running on processors with weak
consistency models such as x86TSO, SPARCPSO, and SPARCRMO. At the center of
our method is a unified framework for deciding the feasibility of interthread
interferences to avoid propagating spurious data flows during static analysis
and thus boost the performance of the static analyzer. We formulate the
checking of interference feasibility as a set of Datalog rules which are both
efficiently solvable and general enough to capture a range of hardwarelevel
memory models. Compared to existing techniques, our method can significantly
reduce the number of bogus alarms as well as unsound proofs. We implemented the
method and evaluated it on a large set of multithreaded C programs. Our
experiments showthe method significantly outperforms stateoftheart
techniques in terms of accuracy with only moderate runtime overhead.

Interrupts have been widely used in safetycritical computer systems to
handle outside stimuli and interact with the hardware, but reasoning about
interruptdriven software remains a difficult task. Although a number of static
verification techniques have been proposed for interruptdriven software, they
often rely on constructing a monolithic verification model. Furthermore, they
do not precisely capture the complete execution semantics of interrupts such as
nested invocations of interrupt handlers. To overcome these limitations, we
propose an abstract interpretation framework for static verification of
interruptdriven software that first analyzes each interrupt handler in
isolation as if it were a sequential program, and then propagates the result to
other interrupt handlers. This iterative process continues until results from
all interrupt handlers reach a fixed point. Since our method never constructs
the global model, it avoids the upfront blowup in model construction that
hampers existing, nonmodular, verification techniques. We have evaluated our
method on 35 interruptdriven applications with a total of 22,541 lines of
code. Our results show the method is able to quickly and more accurately
analyze the behavior of interrupts.

We propose a constraintbased flowsensitive static analysis for concurrent
programs by iteratively composing threadmodular abstract interpreters via the
use of a system of lightweight constraints. Our method is compositional in that
it first applies sequential abstract interpreters to individual threads and
then composes their results. It is flowsensitive in that the causality
ordering of interferences (flow of data from global writes to reads) is modeled
by a system of constraints. These interference constraints are lightweight
since they only refer to the execution order of program statements as opposed
to their numerical properties: they can be decided efficiently using an
offtheshelf Datalog engine. Our new method has the advantage of being more
accurate than existing, flowinsensitive, static analyzers while remaining
scalable and providing the expected soundness and termination guarantees even
for programs with unbounded data. We implemented our method and evaluated it on
a large number of benchmarks, demonstrating its effectiveness at increasing the
accuracy of threadmodular abstract interpretation.

In this paper, we describe the numerical reconstruction method for
quantitative photoacoustic tomography (QPAT) based on the radiative transfer
equation (RTE), which models light propagation more accurately than diffusion
approximation (DA). We investigate the reconstruction of absorption coefficient
and/or scattering coefficient of biological tissues. Given the scattering
coefficient, an improved fixedpoint iterative method is proposed to retrieve
the absorption coefficient for its cheap computational cost. And we prove the
convergence. To retrieve two coefficients simultaneously, BarzilaiBorwein (BB)
method is applied. Since the reconstruction of optical coefficients involves
the solution of original and adjoint RTEs in the framework of optimization, an
efficient solver with high accuracy is improved from~\cite{Gao}. Simulation
experiments illustrate that the improved fixedpoint iterative method and the
BB method are the comparative methods for QPAT in two cases.

We study the $s$wave kaonnucleon bound state with the strangeness $S=1$ in
the BetheSalpeter formalism in the ladder and instantaneous approximations. We
solve the BetheSalpeter equation of the bound state and obtain the
BetheSalpeter amplitude. It is shown that the $K^p$ bound state exists in
this formalism. We also study the decay width of the bound state based on the
BetheSalpeter techniques. The mass of this bound state is 1422 MeV and its
decay width is obviously smaller than that of $\Lambda(1405)$. These results
indicate that there may be some other structures in the observed resonance.

We study the two dimensional water waves problem with surface tension in the
case when there is a nonzero contact angle between the free surface and the
bottom. In the presence of surface tension, dissipations take place at the
contact point. Moreover, when the contact angle is less than $\pi/6$, no
singularity appears in our settings. Using elliptic estimates in corner domains
and a geometric approach, we prove an a priori estimate for the water waves
problem.

Largescale deep neural networks (DNNs) are both compute and memory
intensive. As the size of DNNs continues to grow, it is critical to improve the
energy efficiency and performance while maintaining accuracy. For DNNs, the
model size is an important factor affecting performance, scalability and energy
efficiency. Weight pruning achieves good compression ratios but suffers from
three drawbacks: 1) the irregular network structure after pruning; 2) the
increased training complexity; and 3) the lack of rigorous guarantee of
compression ratio and inference accuracy. To overcome these limitations, this
paper proposes CirCNN, a principled approach to represent weights and process
neural networks using blockcirculant matrices. CirCNN utilizes the Fast
Fourier Transform (FFT)based fast multiplication, simultaneously reducing the
computational complexity (both in inference and training) from O(n2) to
O(nlogn) and the storage complexity from O(n2) to O(n), with negligible
accuracy loss. Compared to other approaches, CirCNN is distinct due to its
mathematical rigor: it can converge to the same effectiveness as DNNs without
compression. The CirCNN architecture, a universal DNN inference engine that can
be implemented on various hardware/software platforms with configurable network
architecture. To demonstrate the performance and energy efficiency, we test
CirCNN in FPGA, ASIC and embedded processors. Our results show that CirCNN
architecture achieves very high energy efficiency and performance with a small
hardware footprint. Based on the FPGA implementation and ASIC synthesis
results, CirCNN achieves 6102X energy efficiency improvements compared with
the best stateoftheart results.