• 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 information-directed random walk (IRW) on the tree-structured 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 logarithmic-order sample complexity with respect to the size of the search space provided that the Kullback-Leibler 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 two-band 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 orbital-triplet spin-singlet superconducting phase. Concerning the mechanism of superconductivity -- we propose an analogy with superconductivity in alkali-doped 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 quantum-speed-limit (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 Bethe-Salpeter formalism in the ladder and instantaneous approximations. With the kernel generated by the hadronic effective Lagrangian, two open-charmed 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 space-like (SL) region in the Bethe-Salpeter 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 high-resolution 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 high-resolution single molecule localization microscopy as well as the problem addressed in this paper on localization of space debris using a space-based 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 large-scale sparse 3D inverse problem. A new nonconvex regularization method with a data-fitting term based on Kullback-Leibler (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 data-fitting 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 depth-encoding PSFs as well.
  • If a finite group of orientation-preserving diffeomorphisms of the 3-dimensional 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(g-1). In the present paper we classify (up to conjugation) all such group actions and surfaces for which the maximal possible order 12(g-1) is achieved, and note that the unknotted surfaces can be realized by equivariant minimal surfaces in a 3-torus.
  • 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 quasi-homography warp, which effectively balances the perspective distortion against the projective distortion in the non-overlapping region to create a more natural-looking 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 parameter-free. A comprehensive experiment shows that a quasi-homography warp outperforms some state-of-the-art 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 single-pixel optical 2D imaging system using a single multimode fiber (MF) is proposed. The MF acted as the all-optical random pattern generator. Light with different wavelengths pass through a single MF will generator all-optical random speckle patterns, which have a low correlation of 0.074 with 0.1nm wavelength step from 1518.0nm to 1567.9nm. The all-optical 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 single-pixel all-optical CS imaging possible. In our work, the all-optical random speckle patterns are analyzed and used to perform CS imaging in our proposed system and the results shows a single-pixel photo-detector 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 signal-to-noise 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 Bethe-Salpeter equation approach. With the kernel containing one-particle-exchange diagrams, we solve the Bethe-Salpeter 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 zero-field nuclear magnetic resonance (NMR) measurements on single crystals of multiferroics CuBr$_2$. High quality of the sample is demonstrated by the sharp magnetic and magnetic-driven ferroelectric transition at $T_N=T_C\approx$ 74~K. The zero-field $^{79}$Br and $^{81}$Br NMR are resolved below $T_N$. The spin-lattice 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 quasi-one-dimensional (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 one-hot 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 non-equilibrium Green function method, we study the Andreev reflection in a Y-shaped graphene-superconductor device by tight-binding 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 finite-size 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 2-component Dirac fermion extends this notion to gapless systems. This duality sheds light on highly non-trivial problems such as the half-filled Landau level, the superconductor-insulator transition, and surface states of strongly coupled topological insulators. Although this boson-fermion duality has undergone many consistency checks, it has remained unproven. We describe the duality in a non-perturbative fashion using an exact UV mapping of partition functions on a 3D Euclidean lattice.
  • We study the three-body radiative decays $B\to \phi(\rho) K\gamma$ induced by a flavor-changing neutral current in the perturbative QCD approach. Pseudoscalar-vector ($PV$) distribution amplitudes (DAs) are introduced for the final-state $\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 time-like $PV$ form factors. In addition to the dominant electromagnetic penguin, the subleading chromomagnetic penguin, quark-loop and annihilation diagrams are also calculated. After determining the $PV$ DAs from relevant branching-ratio 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 Bethe-Salpeter 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 memory-model-aware static program analysis method for accurately analyzing the behavior of concurrent software running on processors with weak consistency models such as x86-TSO, SPARC-PSO, and SPARC-RMO. At the center of our method is a unified framework for deciding the feasibility of inter-thread 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 hardware-level 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 state-of-the-art techniques in terms of accuracy with only moderate run-time overhead.
  • Interrupts have been widely used in safety-critical computer systems to handle outside stimuli and interact with the hardware, but reasoning about interrupt-driven software remains a difficult task. Although a number of static verification techniques have been proposed for interrupt-driven 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 interrupt-driven 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 up-front blowup in model construction that hampers existing, non-modular, verification techniques. We have evaluated our method on 35 interrupt-driven 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 constraint-based flow-sensitive static analysis for concurrent programs by iteratively composing thread-modular 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 flow-sensitive 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 off-the-shelf Datalog engine. Our new method has the advantage of being more accurate than existing, flow-insensitive, 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 thread-modular 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 fixed-point iterative method is proposed to retrieve the absorption coefficient for its cheap computational cost. And we prove the convergence. To retrieve two coefficients simultaneously, Barzilai-Borwein (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 fixed-point iterative method and the BB method are the comparative methods for QPAT in two cases.
  • We study the $s$-wave kaon-nucleon bound state with the strangeness $S=-1$ in the Bethe-Salpeter formalism in the ladder and instantaneous approximations. We solve the Bethe-Salpeter equation of the bound state and obtain the Bethe-Salpeter 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 Bethe-Salpeter 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 non-zero 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.
  • Large-scale 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 block-circulant 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 6-102X energy efficiency improvements compared with the best state-of-the-art results.