title: BindsNET: A machine learning-oriented spiking neural networks library in Python

abstract: The development of spiking neural network simulation software is a critical component enabling the modeling of neural systems and the development of biologically inspired algorithms. Existing software frameworks support a wide range of neural functionality, software abstraction levels, and hardware devices, yet are typically not suitable for rapid prototyping or application to problems in the domain of machine learning. In this paper, we describe a new Python package for the simulation of spiking neural networks, specifically geared towards machine learning and reinforcement learning. Our software, called BindsNET, enables rapid building and simulation of spiking networks and features user-friendly, concise syntax. BindsNET is built on top of the PyTorch deep neural networks library, enabling fast CPU and GPU computation for large spiking networks. The BindsNET framework can be adjusted to meet the needs of other existing computing and hardware environments, e.g., TensorFlow. We also provide an interface into the OpenAI gym library, allowing for training and evaluation of spiking networks on reinforcement learning problems. We argue that this package facilitates the use of spiking networks for large-scale machine learning experimentation, and show some simple examples of how we envision BindsNET can be used in practice. BindsNET code is available at https://github.com/Hananel-Hazan/bindsnet


title: Period integrals in nonpositively curved manifolds

abstract: We provide an improvement of a half power of log to standard bounds on integrals of Laplace eigenfunctions over submanifolds of codimension 2 or more, where the ambient space is a compact Riemannian manifold with negative sectional curvature. We provide the same improvement for hypersurfaces whose second fundamental form differs sufficiently from that of spheres of infinite radius. This result extends previous ones obtained in the 2-dimensional setting by Chen and Sogge; Sogge, Xi, and Zhang; and the author.


title: A Comparative Analysis of Contact Models in Trajectory Optimization for Manipulation

abstract: In this paper, we analyze the effects of contact models on contact-implicit trajectory optimization for manipulation. We consider three different approaches: (1) a contact model that is based on complementarity constraints, (2) a smooth contact model, and our proposed method (3) a variable smooth contact model. We compare these models in simulation in terms of physical accuracy, quality of motions, and computation time. In each case, the optimization process is initialized by setting all torque variables to zero, namely, without a meaningful initial guess. For simulations, we consider a pushing task with varying complexity for a 7 degrees-of-freedom robot arm. Our results demonstrate that the optimization based on the proposed variable smooth contact model provides a good trade-off between the physical fidelity and quality of motions at the cost of increased computation time.


title: Composite Marginal Likelihood Methods for Random Utility Models

abstract: We propose a novel and flexible rank-breaking-then-composite-marginal-likelihood (RBCML) framework for learning random utility models (RUMs), which include the Plackett-Luce model. We characterize conditions for the objective function of RBCML to be strictly log-concave by proving that strict log-concavity is preserved under convolution and marginalization. We characterize necessary and sufficient conditions for RBCML to satisfy consistency and asymptotic normality. Experiments on synthetic data show that RBCML for Gaussian RUMs achieves better statistical efficiency and computational efficiency than the state-of-the-art algorithm and our RBCML for the Plackett-Luce model provides flexible tradeoffs between running time and statistical efficiency.


title: Analysis of DAWNBench, a Time-to-Accuracy Machine Learning Performance Benchmark

abstract: Researchers have proposed hardware, software, and algorithmic optimizations to improve the computational performance of deep learning. While some of these optimizations perform the same operations faster (e.g., increasing GPU clock speed), many others modify the semantics of the training procedure (e.g., reduced precision), and can impact the final model's accuracy on unseen data. Due to a lack of standard evaluation criteria that considers these trade-offs, it is difficult to directly compare these optimizations. To address this problem, we recently introduced DAWNBench, a benchmark competition focused on end-to-end training time to achieve near-state-of-the-art accuracy on an unseen dataset---a combined metric called time-to-accuracy (TTA). In this work, we analyze the entries from DAWNBench, which received optimized submissions from multiple industrial groups, to investigate the behavior of TTA as a metric as well as trends in the best-performing entries. We show that TTA has a low coefficient of variation and that models optimized for TTA generalize nearly as well as those trained using standard methods. Additionally, even though DAWNBench entries were able to train ImageNet models in under 3 minutes, we find they still underutilize hardware capabilities such as Tensor Cores. Furthermore, we find that distributed entries can spend more than half of their time on communication. We show similar findings with entries to the MLPERF v0.5 benchmark.


title: Geometric distance between positive definite matrices of different dimensions

abstract: We show how the Riemannian distance on $\mathbb{S}^n_{++}$, the cone of $n\times n$ real symmetric or complex Hermitian positive definite matrices, may be used to naturally define a distance between two such matrices of different dimensions. Given that $\mathbb{S}^n_{++}$ also parameterizes $n$-dimensional ellipsoids, and inner products on $\mathbb{R}^n$, $n \times n$ covariance matrices of nondegenerate probability distributions, this gives us a natural way to define a geometric distance between a pair of such objects of different dimensions.


title: Elastic capsules at liquid-liquid interfaces

abstract: We investigate the deformation of elastic microcapsules adsorbed at liquid-liquid interfaces. An initially spherical elastic capsule at a liquid-liquid interface undergoes circumferential stretching due to the liquid-liquid surface tension and becomes lens- or discus-shaped, depending on its bending rigidity. The resulting elastic capsule deformation is qualitatively similar, but distinct from the deformation of a liquid droplet into a liquid lens at a liquid-liquid interface. We discuss the deformed shapes of droplets and capsules adsorbed at liquid-liquid interfaces for a whole range of different surface elasticities: from droplets (only surface tension) deforming into liquid lenses, droplets with a Hookean membrane (finite stretching modulus, zero bending modulus) deforming into elastic lenses, to microcapsules (finite stretching and bending modulus) deforming into rounded elastic lenses. We calculate capsule shapes at liquid-liquid interfaces numerically using shape equations from nonlinear elastic shell theory. Finally, we present theoretical results for the contact angle (or the capsule height) and the maximal capsule curvature at the three phase contact line. These results can be used to infer information about the elastic moduli from optical measurements. During capsule deformation into a lens-like shape, surface energy of the liquid-liquid interface is converted into elastic energy of the capsule shell giving rise to an overall adsorption energy gain by deformation. Soft hollow capsules exhibit a pronounced increase of the adsorption energy as compared to filled soft particles and, thus, are attractive candidates as foam and emulsion stabilizers.


title: Study of Automatic GPU Offloading Technology for Open IoT

abstract: IoT technologies have been progressed. Now Open IoT concept has attracted attentions which achieve various IoT services by integrating horizontal separated devices and services. For Open IoT era, we have proposed the Tacit Computing technology to discover the devices with necessary data for users on demand and use them dynamically. However, existing Tacit Computing does not care about performance and operation cost. Therefore in this paper, we propose an automatic GPU offloading technology as an elementary technology of Tacit Computing which uses Genetic Algorithm to extract appropriate offload loop statements to improve performances. We evaluate a C/C++ matrix manipulation to verify effectiveness of GPU offloading and confirm more than 35 times performances within 1 hour tuning time.


title: A Uniform-in-$P$ Edgeworth Expansion under Weak Cram\'{e}r Conditions

abstract: This paper provides a finite sample bound for the error term in the Edgeworth expansion for a sum of independent, potentially discrete, nonlattice random vectors, using a uniform-in-$P$ version of the weaker Cram\'{e}r condition in Angst and Poly (2017). This finite sample bound can be used to derive an Edgeworth expansion that is uniform over the distributions of the random vectors. Using this result, we derive a uniform-in-$P$ higher order expansion of resampling-based distributions.


title: Pressure tuned superconductivity and normal state behavior in Ba(Fe$_{0.943}$Co$_{0.057}$)$_2$As$_2$ near the antiferromagnetic boundary

abstract: Superconductivity in iron pnictides is unconventional and pairing may be mediated by magnetic fluctuations in the Fe-sublattice. Pressure is a clean method to explore superconductivity in iron based superconductors by tuning the ground state continuously without introducing disorder. Here we present a systematic high pressure transport study in Ba(Fe$_{1-x}$Co$_{x}$)$_2$As$_2$ single crystals with $x=$ 0.057, which is near the antiferromagnetic instability. Resistivity $\rho=\rho_0+AT^{n}$ was studied under applied pressure up to 7.90 GPa. The parameter $n$ approaches a minimum value of $n\approx 1$ at a critical pressure $P_c =$ 3.65 GPa. Near $P_c$, the superconducting transition temperature $T_c$ reaches a maximum value of 25.8 K. In addition, the superconducting diamagnetism at 2 K shows a sudden change around the same critical pressure. These results may be associated with a possible quantum critical point hidden inside the superconducting dome, near optimum $T_c$.


title: On Computing the Multiplicity of Cycles in Bipartite Graphs Using the Degree Distribution and the Spectrum of the Graph

abstract: Counting short cycles in bipartite graphs is a fundamental problem of interest in the analysis and design of low-density parity-check (LDPC) codes. The vast majority of research in this area is focused on algorithmic techniques. Most recently, Blake and Lin proposed a computational technique to count the number of cycles of length $g$ in a bi-regular bipartite graph, where $g$ is the girth of the graph. The information required for the computation is the node degree and the multiplicity of the nodes on both sides of the partition, as well as the eigenvalues of the adjacency matrix of the graph (graph spectrum). In this paper, the result of Blake and Lin is extended to compute the number of cycles of length $g+2, \ldots, 2g-2$, for bi-regular bipartite graphs, as well as the number of $4$-cycles and $6$-cycles in irregular and half-regular bipartite graphs, with $g \geq 4$ and $g \geq 6$, respectively.


title: Dynamical quadrupole structure factor of frustrated ferromagnetic chain

abstract: We investigate the dynamical quadrupole structure factor of a spin-1/2 $J_{1}$-$J_{2}$ Heisenberg chain with competing ferromagnetic $J_{1}$ and antiferromagnetic $J_{2}$ in a magnetic field by exploiting density-matrix renormalization group techniques. In a field-induced spin nematic regime, we observe gapless excitations at $q=\pi$ according to quasi-long-range antiferro-quadrupole correlations. The gapless excitation mode has a quadratic form at the saturation, while it changes into a linear dispersion as the magnetization decreases.


title: First study of small-cell 3D Silicon Pixel Detectors for the High Luminosity LHC

abstract: A study of 3D pixel sensors of cell size 50 {\mu}m x 50 {\mu}m fabricated at IMB-CNM using double-sided n-on-p 3D technology is presented. Sensors were bump-bonded to the ROC4SENS readout chip. For the first time in such a small-pitch hybrid assembly, the sensor response to ionizing radiation in a test beam of 5.6 GeV electrons was studied. Results for non-irradiated sensors are presented, including efficiency, charge sharing, signal-to-noise, and resolution for different incidence angles.


title: The problem of dispersion-free probabilities in Gleason-type theorems for a two-dimensional Hilbert space

abstract: As it is known, Gleason's theorem is not applicable for a two-dimensional Hilbert space since in this situation Gleason's axioms are not strong enough to imply Born's rule thus leaving room for a dispersion-free probability measure i.e., one that has only values 0 and 1. To strengthen Gleason's axioms one must add at least one more assumption. But, as it is argued in the present paper, alternatively one can give up the lattice condition lying in the foundation of Gleason's theorem. Particularly, the lattice structure based on the closed linear subspaces in the Hilbert space could be weakened by the requirement for the meet operation to exist only for the subspaces belonging to commutable projection operators. The paper demonstrates that this weakening can resolve the problem of the dispersion-free probability measure in the case of a qubit.


title: PETSc/TS: A Modern Scalable ODE/DAE Solver Library

abstract: High-quality ordinary differential equation (ODE) solver libraries have a long history, going back to the 1970s. Over the past several years we have implemented, on top of the PETSc linear and nonlinear solver package, a new general-purpose, extensive, extensible library for solving ODEs and differential algebraic equations (DAEs). Package includes support for both forward and adjoint sensitivities that can be easily utilized by the TAO optimization package, which is also part of PETSc. The ODE/DAE integrator library strives to be highly scalable but also to deliver high efficiency for modest-sized problems. The library includes explicit solvers, implicit solvers, and a collection of implicit-explicit solvers, all with a common user interface and runtime selection of solver types, adaptive error control, and monitoring of solution progress. The library also offers enormous flexibility in selection of nonlinear and linear solvers, including the entire suite of PETSc iterative solvers, as well as several parallel direct solvers.


title: Hybrid lattices and thin subgroups of Picard modular groups

abstract: We consider a certain hybridization construction which produces a subgroup of ${\rm PU}(n,1)$ from a pair of lattices in ${\rm PU}(n-1,1)$. Among the Picard modular groups ${\rm PU}(2,1,\mathcal{O}_d)$, we show that the hybrid of pairs of Fuchsian subgroups ${\rm PU}(1,1,\mathcal{O}_d)$ is a lattice when $d=1$ and $d=7$, and a geometrically infinite thin subgroup when $d=3$, that is an infinite-index subgroup with the same Zariski-closure as the full lattice.


title: Effect of finite disk-thickness on swing amplification of non-axisymmetric perturbations in a sheared galactic disk

abstract: A typical galactic disk is observed to have a finite thickness. Here, we present the study of the physical effect of introduction of finite thickness on the generation of small-scale spiral arms by swing amplification in a differentially rotating galactic disk. The galactic disk is modelled first as a one-fluid system, and then as a gravitationally-coupled two-fluid (stars and gas) system where each fluid is taken as isothermal, and corotating with each other. We derived the equations governing the evolution of growth of the non-axisymmetric perturbations in a sheared frame of reference while incorporating the effect of finite thickness of a galactic disk. We found that the finite thickness of a galactic disk has a generic trend of suppressing the growth of the non-axisymmetric perturbations via swing amplification. Moreover, even the observed range of disk-thickness values (~ 300-500 pc) can lead to a complete suppression of swing amplification for Q ~ 1.7, whereas for an infinitesimally-thin disk, the corresponding critical value is Q ~ 2. For a two-fluid (stars and gas) system, the net amplification is shown to be set by the mutual interplay of the effect of interstellar gas in promoting the spiral features and the effect of finite thickness in preventing the spiral arms. The coexistence of these two opposite effects is shown to be capable of giving rise to diverse and complex dynamical behaviour.


title: Measuring the spectral running from cosmic microwave background and primordial black holes

abstract: We constrain the spectral running by combining cosmic microwave background (CMB) data, baryon acoustic oscillation (BAO) data and the constraint from primordial black holes (PBHs). We find that the constraint from PBHs has a significant impact on the running of running of scalar spectral index, and a power-law scalar power spectrum without running is consistent with observational data once the constraint from PBHs is taken into account. In addition, from the constraints on the slow-roll parameters, the derived tensor spectral index in the single-field slow-roll inflation model is quite small, namely $|n_t|\lesssim 9.3\times 10^{-3}$ which will be very difficult to be measured by CMB data only in the future, and the absolute value of derived running of tensor spectral index is not larger than $2.1\times 10^{-4}$ at 95% confidence level.


title: Existence, uniqueness, estimation and continuous dependence of the solutions of a nonlinear integral and an integrodifferential equations of fractional order

abstract: By means of two fractional order integral inequalities we investigate the existence and uniqueness of the solutions of the fractional nonlinear Volterra integral equation and a fractional nonlinear integrodifferential equation in Banach space $C_{\xi}$, using an adequate norm, $||\cdot||_{\xi,\infty}$. We estimate the solutions and investigate their continuous dependence.


title: On the existence and stability for non-instantaneuos impulsive fractional integrodifferential equation

abstract: In this paper, by means of Banach fixed point theorem, we investigate the existence and Ulam--Hyers--Rassias stability of the non-instantaneous impulsive integrodifferential equation by means of $\psi$-Hilfer fractional derivative. In this sense, some examples are presented, in order to consolidate the results obtained.


title: Analysis and Design of 8-Bit CMOS Priority Encoders

abstract: A comprehensive review and fair comparison of previous priority encoder (PE) designs over the past one and a half decades are presented using a 45 nm technology. Further, potential limitations of existed PEs are identified, based on which we propose a robust PE design. The new PE is able to eliminate race condition and charge sharing problem which are suffered by almost all the previous designs. Besides, the proposed PE can also be used in comprising higher order PEs by incorporating a carefully designed look-ahead structure. Simulation results demonstrate that our design can achieve one of the best power and delay performance among previous PEs and are free from potential risks.


title: NDN, CoAP, and MQTT: A Comparative Measurement Study in the IoT

abstract: This paper takes a comprehensive view on the protocol stacks that are under debate for a future Internet of Things (IoT). It addresses the holistic question of which solution is beneficial for common IoT use cases. We deploy NDN and the two popular IP-based application protocols, CoAP and MQTT, in its different variants on a large-scale IoT testbed in single- and multi-hop scenarios. We analyze the use cases of scheduled periodic and unscheduled traffic under varying loads. Our findings indicate that (a) NDN admits the most resource-friendly deployment on nodes, and (b) shows superior robustness and resilience in multi-hop scenarios, while (c) the IP protocols operate at less overhead and higher speed in single-hop deployments. Most strikingly we find that NDN-based protocols are in significantly better flow balance than the UDP-based IP protocols and require less corrective actions.


title: Embedding Logical Queries on Knowledge Graphs

abstract: Learning low-dimensional embeddings of knowledge graphs is a powerful approach used to predict unobserved or missing edges between entities. However, an open challenge in this area is developing techniques that can go beyond simple edge prediction and handle more complex logical queries, which might involve multiple unobserved edges, entities, and variables. For instance, given an incomplete biological knowledge graph, we might want to predict "em what drugs are likely to target proteins involved with both diseases X and Y?" -- a query that requires reasoning about all possible proteins that {\em might} interact with diseases X and Y. Here we introduce a framework to efficiently make predictions about conjunctive logical queries -- a flexible but tractable subset of first-order logic -- on incomplete knowledge graphs. In our approach, we embed graph nodes in a low-dimensional space and represent logical operators as learned geometric operations (e.g., translation, rotation) in this embedding space. By performing logical operations within a low-dimensional embedding space, our approach achieves a time complexity that is linear in the number of query variables, compared to the exponential complexity required by a naive enumeration-based approach. We demonstrate the utility of this framework in two application studies on real-world datasets with millions of relations: predicting logical relationships in a network of drug-gene-disease interactions and in a graph-based representation of social interactions derived from a popular web forum.


title: Two-photon absorption coefficient determination with a differential F-scan technique

abstract: In this paper we present a modification to the recently proposed transmission F-scan technique, the differential F-scan technique. In differential F-scan technique the programmed focal distance in the electronic-tunable lens oscillates, allowing the light detector of the setup to record a signal proportional to the derivative of the signal recorded with an F-scan. As for the differential Z-scan a background-free signal is obtained, but also the optical setup is simplified and the available laser power is double. We also present and validate a new fitting-procedure protocol that increments the accuracy of the technique. Finally, we show that fitting a signal from differential F-scan or the derivate of the signal of transmission F-scan is more accurate than simply fitting the signal from F-scan directly. Results from two-photon absorption at 790 nm of CdS, ZeSe and CdSe are presented.


title: Theory of superconductivity in hole-doped monolayer MoS$_{2}$

abstract: We theoretically investigate the Cooper-pair symmetry to be realized in hole-doped monolayer MoS$_2$ by solving linearized BCS gap equations on the three-orbital attractive Hubbard-like model in the presence of the atomic spin-orbit coupling. In hole-doped monolayer MoS$_2$, both spin-orbit coupling and the multi-orbital effects are more prominent than those of electron-doped system. Near the valence band edge, the Fermi surfaces are composed of three different types of hole pockets, namely, one mainly consisting of the almost spin-degenerate $|d_{z^{2}}\rangle$ orbital near $\Gamma$ point, and the others of the spin-split upper and lower bands near ${\rm K}$ and ${\rm K}'$ points arising from the $|d_{x^{2}-y^{2}}\rangle$ and $|d_{xy}\rangle$ orbitals. The number of relevant Fermi pockets increases with increase of the doping. At very low doping, the upper split bands of $|d_{x^{2}-y^{2}}\rangle$ and $|d_{xy}\rangle$ are concerned, yielding extremely low $T_{\rm c}$ due to small density of states of the split bands. For further doping, the conventional spin-singlet state (SS) appears in the $\Gamma$ pocket, which has a mixture of the spin-triplet (orbital-singlet) (ST-OS) and spin-singlet (orbital-triplet) (SS-OT) states in the K and K$'$ pockets. The ratio of the mixture depends on the relative strength of the interactions, and the sign of the exchange interactions. Moderately strong ferromagnetic exchange interactions even lead to the pairing state with the dominant ST-OS state over the conventional SS one. With these observations, we expect that the fascinating pairing with relatively high $T_{\rm c}$ emerges at high doping that involves all the three Fermi pockets.


title: Coexistence of extended and localized states in one-dimensional systems

abstract: Mobility edge transitions from localized to extended states have been observed in two and three dimensional systems, for which sound theoretical explanations have also been derived. One-dimensional lattice models have failed to predict their emergence, offering no clues on how to actually probe this phenomenon in lower dimensions. This work reports results for a class of tight-binding models with electron-mass position dependence, for which localized-extended wave function transitions can be identified. We show that it is possible to control the density of localized and extended states by tuning the transition-related parameter for a continuous range of energy values. Mathematically exact results for extended or localized states are derived in two extreme conditions of this parameter, as well as an exact energy value for the mobility edge transition in the intermediate regime. Our framework provides a clear point of view on the phenomena and can also be harnessed for setting up experiments to probe to precisely evaluate the associated mobility edges using state-of-the-art technology.


title: The $\varphi$-Dimension of Cyclic Nakayama Algebras

abstract: K. Igusa and G. Todorov introduced the $\varphi$ function which generalizes the notion of projective dimension. We study the behavior of the $\varphi$ function for cyclic Nakayama algebras of infinite global dimension. We prove that the supremum of values of $\varphi$ is always an even number. In particular we show that the $\varphi$-dimension is $2$ if and only if the algebra satisfies certain symmetry conditions. Also we give a sharp upper bound for $\varphi$-dimension in terms of the number of monomial relations which describes the algebra.


title: Asymptotic Refinements of a Misspecification-Robust Bootstrap for Generalized Method of Moments Estimators

abstract: I propose a nonparametric iid bootstrap that achieves asymptotic refinements for t tests and confidence intervals based on GMM estimators even when the model is misspecified. In addition, my bootstrap does not require recentering the moment function, which has been considered as critical for GMM. Regardless of model misspecification, the proposed bootstrap achieves the same sharp magnitude of refinements as the conventional bootstrap methods which establish asymptotic refinements by recentering in the absence of misspecification. The key idea is to link the misspecified bootstrap moment condition to the large sample theory of GMM under misspecification of Hall and Inoue (2003). Two examples are provided: Combining data sets and invalid instrumental variables.


title: Nonlinear imaging through a golden spiral multicore fiber

abstract: We report two-photon lensless imaging through a novel golden spiral multicore fiber. This unique layout optimizes the sidelobe levels, field of view, cross-talk, group delay and mode density to achieve a sidelobe contrast of atleast 10.9 dB. We demonstrate experimentally the ability to generate and scan a focal point with a femtosecond pulse and perform two-photon imaging.


title: Breaking Taylor-Proudman balance by magnetic field in stellar convection zone

abstract: We carry out high-resolution calculations for the stellar convection zone. The main purpose of this study is to investigate the effect of a small-scale dynamo on the differential rotation. The solar differential rotation deviates from the Taylor-Proudman state in which the angular velocity does not change along the rotational axis. To break the Taylor-Proudman state deep in the convection zone, it is thought that a latitudinal entropy gradient is required. In this study, we find that the small-scale dynamo has three roles in the deviation of the stellar differential rotation from the Taylor-Proudman state. 1. The shear of the angular velocity is suppressed. This leads to a situation where the latitudinal entropy gradient efficiently breaks the Taylor-Proudman state. 2. The perturbation of the entropy is increased with suppressing the turbulent velocity between up- and downflows. 3. The convection velocity is reduced. This increases the effect of the rotation on the convection. The second and third factors increase the latitudinal entropy gradient and break the Taylor-Proudman state. We find that the efficient small-scale dynamo has a significant impact on the stellar differential rotation.


title: Calibration for computer experiments with binary responses and application to cell adhesion study

abstract: Calibration refers to the estimation of unknown parameters which are present in computer experiments but not available in physical experiments. An accurate estimation of these parameters is important because it provides a scientific understanding of the underlying system which is not available in physical experiments. Most of the work in the literature is limited to the analysis of continuous responses. Motivated by a study of cell adhesion experiments, we propose a new calibration framework for binary responses. Its application to the T cell adhesion data provides insight into the unknown values of the kinetic parameters which are difficult to determine by physical experiments due to the limitation of the existing experimental techniques.


title: On the relative velocity distribution for general statistics and an application to big-bang nucleosynthesis under Tsallis statistics

abstract: The distribution function of the relative velocity in a two-body reaction of nonrelativistic uncorrelated particles is derived for general cases of given distribution functions of single particle velocities. The distribution function is then used in calculations of thermonuclear reaction rates. As an example, we take the Tsallis non-Maxwellian distribution, and show that the distribution function of the relative velocity is different from the Tsallis distribution. We identify an inconsistency in previous studies of nuclear reaction rates within Tsallis statistics, and derive revised nuclear reaction rates. Utilizing the revised rates, accurate results of big bang nucleosynthesis are obtained for the Tsallis statistics. For this application it is more difficult to reduce the primordial 7Li abundance while keeping other nuclear abundances within the observational constraints. A small deviation from a Maxwell-Boltzmann distribution can increase the D abundance and slightly reduce 7Li abundance. Although it is impossible to realize a 7Li abundance at the level observed in metal-poor stars, a significant decrease is possible while maintaining a consistency with the observed D abundance.


title: EigenNetworks

abstract: Many applications donot have the benefit of the laws of physics to derive succinct descriptive models for observed data. In alternative, interdependencies among $N$ time series $\{ x_{nk}, k>0 \}_{n=1}^{N}$ are nowadays often captured by a graph or network $G$ that in practice may be very large. The network itself may change over time as well (i.e., as $G_k$). Tracking brute force the changes of time varying networks presents major challenges, including the associated computational problems. Further, a large set of networks may not lend itself to useful analysis. This paper approximates the time varying networks $\left\{G_k\right\}$ as weighted linear combinations of eigennetworks. The eigennetworks are fixed building blocks that are estimated by first learning the time series of graphs $G_k$ from the data $\{ x_{nk}, k>0 \}_{n=1}^{N}$, followed by a Principal Network Analysis procedure. The weights of the eigennetwork representation are eigenfeatures and the time varying networks $\left\{G_k\right\}$ describe a trajectory in eigennetwork space. These eigentrajectories should be smooth since the networks $G_k$ vary at a much slower rate than the data $x_{nk}$, except when structural network shifts occur reflecting potentially an abrupt change in the underlying application and sources of the data. Algorithms for learning the time series of graphs $\left\{G_k\right\}$, deriving the eigennetworks, eigenfeatures and eigentrajectories, and detecting changepoints are presented. Experiments on simulated data and with two real time series data (a voting record of the US senate and genetic expression data for the \textit{Drosophila Melanogaster} as it goes through its life cycle) demonstrate the performance of the learning and provide interesting interpretations of the eigennetworks.


title: A verifiable framework of entanglement-free quantum secret sharing with information-theoretical security

abstract: Quantum secret sharing (QSS) schemes without entanglement have huge advantages in scalability and are easier to realize as they only require sequential communications of a single quantum system. However, these schemes often come with drawbacks such as exact ($n, n$) structure, security flaws and absences of effective cheating detections. To address these problems, we propose a verifiable framework by utilizing entanglement-free states to construct ($t, n$)-QSS schemes. Our work is the heuristic step towards information-theoretical security in entanglement-free QSS, and it sheds light on how to establish effective verification mechanism against cheating. As a result, the proposed framework has a significant importance in constructing QSS schemes for versatile applications in quantum networks due to its intrinsic scalability, flexibility and information theoretical security.


title: A Consistent Variance Estimator for 2SLS When Instruments Identify Different LATEs

abstract: Under treatment effect heterogeneity, an instrument identifies the instrument-specific local average treatment effect (LATE). With multiple instruments, two-stage least squares (2SLS) estimand is a weighted average of different LATEs. What is often overlooked in the literature is that the postulated moment condition evaluated at the 2SLS estimand does not hold unless those LATEs are the same. If so, the conventional heteroskedasticity-robust variance estimator would be inconsistent, and 2SLS standard errors based on such estimators would be incorrect. I derive the correct asymptotic distribution, and propose a consistent asymptotic variance estimator by using the result of Hall and Inoue (2003, Journal of Econometrics) on misspecified moment condition models. This can be used to correctly calculate the standard errors regardless of whether there is more than one LATE or not.


title: The Value of Information in Retrospect

abstract: In the course of any statistical analysis, it is necessary to consider issues of data quality and model appropriateness. Value of information methods were initially put forward in the middle of the twentieth century in order to provide a framework for choosing between potential sources of information. However, since their genesis, value of information methods have been largely neglected by statisticians. In this paper we review and extend existing value of information methods and recommend the use of three quantities for identifying influential and outlying data: an influence measure previously suggested by \cite{kempthorne1986}, a related quantity known as the expected value of sample information that is used to gauge how much influence we would expect a portion of the data to have, and the ratio of these two quantities which serves as a comparison between observed influence and expected influence. We study the basic theoretical properties of those quantities and illustrate our proposed approach using two datasets. A data set containing employment rates and other economic factors in U.S. first presented by \cite{longley} is used to provide an example in the case of linear regression. HIV surveillance data collected from prenatal clinics have been the main source of information for monitoring the HIV epidemic in low and middle income countries. A data set providing information about HIV prevalence in Swaziland is used as an example in the case of generalized linear mixed models.


title: The Effect of Multiplicity Cut on Particle Ratios at 7 TeV Proton-Proton Collisions

abstract: The final state hadronic composition of proton-proton collisions at 7 TeV as a function of charged particle multiplicity density is studied. The thermal model is adjusted to match the experimental conditions in order to understand the effect of the experimental cuts on rapidity and multiplicity. This study indicates that the measured hadronic composition under a simultaneous cut on rapidity and multiplicity does not coincide with the fireball hadronic composition. Excluding proton to pion ratio, the adjusted model can explain the measured particle ratios. Therefore, the observed strangeness enhancement as a function of charged particle multiplicity density could be an effect of the experimental cuts.


title: Dynamic Function-on-Scalars Regression

abstract: We develop a modeling framework for dynamic function-on-scalars regression, in which a time series of functional data is regressed on a time series of scalar predictors. The regression coefficient function for each predictor is allowed to be dynamic, which is essential for applications where the association between predictors and a (functional) response is time-varying. For greater modeling flexibility, we design a nonparametric reduced-rank functional data model with an unknown functional basis expansion, which is data-adaptive and, unlike most existing methods, modeled as unknown for appropriate uncertainty quantification. Within a Bayesian framework, we introduce shrinkage priors that simultaneously (i) regularize time-varying regression coefficient functions to be locally static, (ii) effectively remove unimportant predictor variables from the model, and (iii) reduce sensitivity to the dimension of the functional basis. A simulation analysis confirms the importance of these shrinkage priors, with notable improvements over existing alternatives. We develop a novel projection-based Gibbs sampling algorithm, which offers unrivaled computational scalability for fully Bayesian functional regression. We apply the proposed methodology (i) to analyze the time-varying impact of macroeconomic variables on the U.S. yield curve and (ii) to characterize the effects of socioeconomic and demographic predictors on age-specific fertility rates in South and Southeast Asia.


title: Search for Dust Emission from (24) Themis Using the Gemini-North Telescope

abstract: We report the results of a search for a dust trail aligned with the orbit plane of the large main-belt asteroid (24) Themis, which has been reported to have water ice frost on its surface. Observations were obtained with the GMOS instrument on the Gemini-North Observatory in imaging mode, where we used a chip gap to block much of the light from the asteroid, allowing us to take long exposures while avoiding saturation by the object. No dust trail is detected within 2' of Themis to a 3-sigma limiting surface brightness magnitude of 29.7 mag/arcsec^2, as measured along the expected direction of the dust trail. Detailed consideration of dust ejection physics indicates that particles large enough to form a detectable dust trail were unlikely to be ejected as a result of sublimation from an object as large as Themis. We nonetheless demonstrate that our observations would have been capable of detecting faint dust emission as close as 20" from the object, even in a crowded star field. This approach could be used to conduct future searches for sublimation-generated dust emission from Themis or other large asteroids closer to perihelion than was done in this work. It would also be useful for deep imaging of collisionally generated dust emission from large asteroids at times when the visibility of dust features are expected to be maximized, such as during orbit plane crossings, during close approaches to the Earth, or following detected impact events.


title: Global attractor for a nonlinear one-dimensional compressible viscous micropolar fluid

abstract: This paper considers the dynamical behavior of solutions of constitutive systems for 1D compressible viscous and heat-conducting micropolar fluids. With proper constraints on initial data, we prove the existence of global attractors in generalized Sobolev spaces $H^{(1)}_{\delta}$ and $H^{(2)}_{\delta}$.


title: Femtosecond Photonic Viral Inactivation Probed Using Solid-State Nanopores

abstract: We report on the detection of inactivation of virus particles using femtosecond laser radiation by measuring the conductance of a solid state nanopore designed for detecting single virus particles. Conventional methods of assaying for viral inactivation based on plaque forming assays require 24-48 hours for bacterial growth. Nanopore conductance measurements provide information on morphological changes at a single virion level. We show that analysis of a time series of nanopore conductance can quantify the detection of inactivation, requiring only a few minutes from collection to analysis. Morphological changes were verified by Dynamic Light Scattering (DLS). Statistical analysis maximizing the information entropy provides a measure of the Log-reduction value. Taken together, our work provides a rapid method for assaying viral inactivation with femtosecond lasers using solid-state nanopores.


title: Hawking effect of AdS$_2$ black hole in the Jackiw-Teitelboim model

abstract: It might be tempting to consider that the two-dimensional anti-de Sitter black hole in the Jackiw-Teitelboim model is thermally hot by invoking the non-vanishing surface gravity. So, one might expect that the local temperature would also be blue-shifted near the horizon, while it would vanish at infinity because of the Tolman factor of the local temperature. In this paper, we will show that the local temperature vanishes everywhere, which respects the equivalence principle perfectly in contrast to the case of the asymptotically flat black holes recovering the equivalence principle just at the horizon. In other words, the origin of radiation for a fixed observer comes from the Unruh effect rather than the Hawking effect. If the inward and outward fluxes are defined non-trivially at the horizon, then the local temperature of the black hole takes the usual Tolman form; however, the divergent Tolman temperature at the horizon is incompatible with the finite local energy density there.


title: Superconductivity in the Hubbard model and its interplay with charge stripes and next-nearest hopping t'

abstract: We report a large-scale density-matrix renormalization group study of the lightly doped Hubbard model on 4-leg cylinders at hole doping concentration $\delta=12.5\%$. By keeping a large number of states for long system sizes, we are able to reveal a delicate interplay between superconductivity and charge and spin density wave orders tunable via next-nearest neighbor hopping t'. For finite t', the ground state is consistent with that of a Luther-Emery liquid, having `half-filled' charge stripes with power-law superconducting and charge-density-wave correlations of wave-length $\lambda=1/2\delta$, but short-range spin correlations. This is in direct contrast to the case with t'=0, where superconducting correlations fall off exponentially while charge- and spin-density modulations are dominant. Our results indicate that a route to robust long-range superconductivity involves destabilizing insulating charge stripes in the doped Hubbard model.


title: Informative Gene Selection for Microarray Classification via Adaptive Elastic Net with Conditional Mutual Information

abstract: Due to the advantage of achieving a better performance under weak regularization, elastic net has attracted wide attention in statistics, machine learning, bioinformatics, and other fields. In particular, a variation of the elastic net, adaptive elastic net (AEN), integrates the adaptive grouping effect. In this paper, we aim to develop a new algorithm: Adaptive Elastic Net with Conditional Mutual Information (AEN-CMI) that further improves AEN by incorporating conditional mutual information into the gene selection process. We apply this new algorithm to screen significant genes for two kinds of cancers: colon cancer and leukemia. Compared with other algorithms including Support Vector Machine, Classic Elastic Net and Adaptive Elastic Net, the proposed algorithm, AEN-CMI, obtains the best classification performance using the least number of genes.


title: Directed Non-Targeted Mass Spectrometry and Chemical Networking for Discovery of Eicosanoids

abstract: Eicosanoids and related species are critical, small bioactive mediators of human physiology and inflammation. While ~1100 distinct eicosanoids have been predicted to exist, to date, less than 150 of these molecules have been measured in humans, limiting our understanding of eicosanoids and their role in human biology. Using a directed non-targeted mass spectrometry approach in conjunction with computational chemical networking of spectral fragmentation patterns, we find over 500 discrete chemical signals highly consistent with known and putative eicosanoids in human plasma, including 46 putative novel molecules not previously described, thereby greatly expanding the breath of prior analytical strategies. In plasma samples from 1500 individuals, we find members of this expanded eicosanoid library hold close association with markers of inflammation, as well as clinical characteristics linked with inflammation, including advancing age and obesity. These experimental and computational approaches enable discovery of new chemical entities and will shed important insight into the role of bioactive molecules in human disease.


title: Understanding Regularized Spectral Clustering via Graph Conductance

abstract: This paper uses the relationship between graph conductance and spectral clustering to study (i) the failures of spectral clustering and (ii) the benefits of regularization. The explanation is simple. Sparse and stochastic graphs create a lot of small trees that are connected to the core of the graph by only one edge. Graph conductance is sensitive to these noisy `dangling sets'. Spectral clustering inherits this sensitivity. The second part of the paper starts from a previously proposed form of regularized spectral clustering and shows that it is related to the graph conductance on a `regularized graph'. We call the conductance on the regularized graph CoreCut. Based upon previous arguments that relate graph conductance to spectral clustering (e.g. Cheeger inequality), minimizing CoreCut relaxes to regularized spectral clustering. Simple inspection of CoreCut reveals why it is less sensitive to small cuts in the graph. Together, these results show that unbalanced partitions from spectral clustering can be understood as overfitting to noise in the periphery of a sparse and stochastic graph. Regularization fixes this overfitting. In addition to this statistical benefit, these results also demonstrate how regularization can improve the computational speed of spectral clustering. We provide simulations and data examples to illustrate these results.


title: Labeling Algorithm and Compact Routing Scheme for a Small World Network Model

abstract: This paper defines the toroidal small world labeling problem that asks for a labeling of the vertices of a network such that the labels possess information that allows a compact routing scheme in the network. We consider the problem over a small world network model we propose. Both the model and the compact routing scheme have applications in peer-to-peer networks. The proposed model is based on the model of Kleinberg (2000), and generates an undirected two-dimensional torus with one random long-range edge per vertex. These random edges create forbidden cycles that mimic the underlying torus topology, and this behavior confuses attempts for extracting the routing information from the network. We show that such forbidden cycles happen with small probability, allowing us to use a breadth-first search that finds the vertices positions on the torus. The positions are pairs of integer numbers that provides routing information to a greedy routing algorithm that finds small paths of the network. We present a linear time labeling algorithm that detects and removes the random edges, finds the underlying torus and labels almost all vertices through a breadth-first search. The labeling algorithm is then used by a compact routing scheme for the proposed small world model.


title: Experimental and Numerical Evaluations of the Viscoelastic Mechanical Behavior of Arterial Tissue

abstract: The viscoelastic properties of arterial tissue dictate vessel behavior in certain disease states, injury modalities, and during some endovascular procedures. In this study, we characterized the viscoelastic mechanical response of porcine abdominal aortic tissue via uniaxial mechanical experiments on vascular tissue samples coupled with full-field surface strain measurements using the 2-D digital image correlation (DIC) technique. The measured stress relaxation response and intimal surface strain field were used to identify material parameters for a proposed viscoelastic anisotropic (VA) constitutive model of the passive arterial wall. The obtained results show that the VA constitutive model is able to capture the viscoelastic mechanical behavior of the arterial tissue over clinically-relevant time scales. The identified material model and numerical simulations provide a comprehensive description of the passive viscoelastic tissue properties of the arterial wall, and a quantitative understanding of the spatial and directional variability underlying arterial tissue mechanical behavior.


title: PAC-learning in the presence of evasion adversaries

abstract: The existence of evasion attacks during the test phase of machine learning algorithms represents a significant challenge to both their deployment and understanding. These attacks can be carried out by adding imperceptible perturbations to inputs to generate adversarial examples and finding effective defenses and detectors has proven to be difficult. In this paper, we step away from the attack-defense arms race and seek to understand the limits of what can be learned in the presence of an evasion adversary. In particular, we extend the Probably Approximately Correct (PAC)-learning framework to account for the presence of an adversary. We first define corrupted hypothesis classes which arise from standard binary hypothesis classes in the presence of an evasion adversary and derive the Vapnik-Chervonenkis (VC)-dimension for these, denoted as the adversarial VC-dimension. We then show that sample complexity upper bounds from the Fundamental Theorem of Statistical learning can be extended to the case of evasion adversaries, where the sample complexity is controlled by the adversarial VC-dimension. We then explicitly derive the adversarial VC-dimension for halfspace classifiers in the presence of a sample-wise norm-constrained adversary of the type commonly studied for evasion attacks and show that it is the same as the standard VC-dimension, closing an open question. Finally, we prove that the adversarial VC-dimension can be either larger or smaller than the standard VC-dimension depending on the hypothesis class and adversary, making it an interesting object of study in its own right.


title: History state formalism for Dirac's theory

abstract: We propose a history state formalism for a Dirac particle. By introducing a reference quantum clock system it is first shown that Dirac's equation can be derived by enforcing a timeless Wheeler-DeWitt-like equation for a global state. The Hilbert space of the whole system constitutes a unitary representation of the Lorentz group with respect to a properly defined invariant product, and the proper normalization of global states directly ensures standard Dirac's norm. Moreover, by introducing a second quantum clock, the previous invariant product emerges naturally from a generalized continuity equation. The invariant parameter $\tau$ associated with this second clock labels history states for different particles, yielding an observable evolution in the case of an hypothetical superposition of different masses. Analytical expressions for both space-time density and electron-time entanglement are provided for two particular families of electron's states, the former including Pryce localized particles.


title: Vacuum Fluctuations and Boundary Conditions in a Global Monopole

abstract: We study the vacuum fluctuations of a massless scalar field $\hat{\Psi}$ on the background of a global monopole. Due to the nontrivial topology of the global monopole spacetime, characterized by a solid deficit angle parametrized by $\eta^2$, we expect that $\left<\hat{\Psi}^2\right>_{\text{ren}}$ and $\left<\hat{T}_{\mu\nu}\right>_{\text{ren}}$ are nonzero and proportional to $\eta^2$, so that they annul in the Minkowski limit $\eta\to0$. However, due to the naked singularity at the monopole core, the evolution of the scalar field is not unique. In fact, they are in one to one correspondence with the boundary conditions which turn into self-adjoint the spatial part of the wave operator. We show that only Dirichlet boundary condition corresponds to our expectations and gives zero contribution to the vacuum fluctuations in Minkowski limit. All other boundary conditions give nonzero contributions in this limit due to the nontrivial interaction between the field and the singularity.


title: Resonant scattering effect on the soft X-ray line emission from the hot interstellar medium: I. galactic bulges

abstract: Diffuse soft X-ray line emission is commonly used to trace the thermal and chemical properties of the hot interstellar medium, as well as its content, in nearby galaxies. Although resonant line scattering complicates the interpretation of the emission, it also offers an opportunity to measure the kinematics of the medium. We have implemented a direct Monte Carlo simulation scheme that enables us to account for resonant scattering effect in the medium, in principle, with arbitrary spatial, thermal, chemical, and kinematic distributions. Here we apply this scheme via dimensionless calculation to an isothermal, chemically uniform, and spherically symmetric medium with a radial density distribution characterized by a $\beta$-model. This application simultaneously account for both optical depth-dependent spatial distortion and intensity change of the resonant line emission due to the scattering, consistent with previous calculations. We further apply the modeling scheme to the OVII and OVIII emission line complex observed in the XMM-Newton RGS spectrum of the M31 bulge. This modeling, though with various limitations due to its simplicity, shows that the resonant scattering could indeed account for much of the spatial distortion of the emission, as well as the relative strengths of the lines, especially the large forbidden to resonant line ratio of the OVII He$\alpha$ triplet. We estimate the isotropic turbulence Mach number of the medium in M31 as $\sim0.17$ for the first time and the line-emitting gas temperature as $\sim2.3\times10^6$ K. We conclude that the resonant scattering may in general play an important role in shaping the soft X-ray spectra of diffuse hot gas in normal galaxies.


title: Anomalous Hall effect induced spin Hall magnetoresistance in an antiferromagnetic Cr2O3/Ta bilayer

abstract: The spin Hall magnetoresistance (SMR) and anomalous Hall effect (AHE) are observed in a Cr2O3/Ta structure. The structural and surface morphology of Cr2O3/Ta bilayers have been investigated. Temperature dependence of longitudinal and transverse resistances measurements confirm the relationship between SMR and AHE signals in Cr2O3/Ta structure. By means of temperature dependent magnetoresistance measurements, the physical origin of SMR in the Cr2O3/Ta structure is revealed, and the contribution to the SMR from the spin current generated by AHE has been proved. The so-called boundary magnetization due to the bulk antiferromagnetic order in Cr2O3 film may be responsible for the relationship of SMR and AHE in the Cr2O3/Ta bilayer.


title: Top quark decays with flavor violation in the B-LSSM

abstract: The decays of top quark $t\rightarrow c\gamma,\;t\rightarrow cg,\;t\rightarrow cZ,\;t\rightarrow ch$ are extremely rare processes in the standard model (SM). The predictions on the corresponding branching ratios in the SM are too small to be detected in the future, hence any measurable signal for the processes at the LHC is a smoking gun for new physics. In the extension of minimal supersymmetric standard model with an additional local $U(1)_{B-L}$ gauge symmetry (B-LSSM), new gauge interaction and new flavor changing interaction affect the theoretical evaluations on corresponding branching ratios of those processes. In this work, we analyze those processes in the B-LSSM, under a minimal flavor violating assumption for the soft breaking terms. Considering the constraints from updated experimental data, the numerical results imply $Br(t\rightarrow c\gamma)\sim5\times10^{-7}$, $Br(t\rightarrow cg)\sim2\times10^{-6}$, $Br(t\rightarrow cZ)\sim4\times10^{-7}$ and $Br(t\rightarrow ch)\sim3\times10^{-9}$ in our chosen parameter space. Simultaneously, new gauge coupling constants $g_{_B},\;g_{_{YB}}$ in the B-LSSM can also affect the numerical results of $Br(t\rightarrow c\gamma,\;cg,\;cZ,\;ch)$.


title: An Explainable Adversarial Robustness Metric for Deep Learning Neural Networks

abstract: Deep Neural Networks(DNN) have excessively advanced the field of computer vision by achieving state of the art performance in various vision tasks. These results are not limited to the field of vision but can also be seen in speech recognition and machine translation tasks. Recently, DNNs are found to poorly fail when tested with samples that are crafted by making imperceptible changes to the original input images. This causes a gap between the validation and adversarial performance of a DNN. An effective and generalizable robustness metric for evaluating the performance of DNN on these adversarial inputs is still missing from the literature. In this paper, we propose Noise Sensitivity Score (NSS), a metric that quantifies the performance of a DNN on a specific input under different forms of fix-directional attacks. An insightful mathematical explanation is provided for deeply understanding the proposed metric. By leveraging the NSS, we also proposed a skewness based dataset robustness metric for evaluating a DNN's adversarial performance on a given dataset. Extensive experiments using widely used state of the art architectures along with popular classification datasets, such as MNIST, CIFAR-10, CIFAR-100, and ImageNet, are used to validate the effectiveness and generalization of our proposed metrics. Instead of simply measuring a DNN's adversarial robustness in the input domain, as previous works, the proposed NSS is built on top of insightful mathematical understanding of the adversarial attack and gives a more explicit explanation of the robustness.


title: Sampling and Super-resolution of Sparse Signals Beyond the Fourier Domain

abstract: Recovering a sparse signal from its low-pass projections in the Fourier domain is a problem of broad interest in science and engineering and is commonly referred to as super-resolution. In many cases, however, Fourier domain may not be the natural choice. For example, in holography, low-pass projections of sparse signals are obtained in the Fresnel domain. Similarly, time-varying system identification relies on low-pass projections on the space of linear frequency modulated signals. In this paper, we study the recovery of sparse signals from low-pass projections in the Special Affine Fourier Transform domain (SAFT). The SAFT parametrically generalizes a number of well known unitary transformations that are used in signal processing and optics. In analogy to the Shannon's sampling framework, we specify sampling theorems for recovery of sparse signals considering three specific cases: (1) sampling with arbitrary, bandlimited kernels, (2) sampling with smooth, time-limited kernels and, (3) recovery from Gabor transform measurements linked with the SAFT domain. Our work offers a unifying perspective on the sparse sampling problem which is compatible with the Fourier, Fresnel and Fractional Fourier domain based results. In deriving our results, we introduce the SAFT series (analogous to the Fourier series) and the short time SAFT, and study convolution theorems that establish a convolution--multiplication property in the SAFT domain.


title: Exploiting wideband spectrum occupancy heterogeneity for weighted compressive spectrum sensing

abstract: Compressive sampling has shown great potential for making wideband spectrum sensing possible at sub-Nyquist sampling rates. As a result, there have recently been research efforts that aimed to develop techniques that leverage compressive sampling to enable compressed wideband spectrum sensing. These techniques consider homogeneous wideband spectrum where all bands are assumed to have similar PU traffic characteristics. In practice, however, wideband spectrum is not homogeneous, in that different spectrum bands could have different PU occupancy patterns. In fact, the nature of spectrum assignment, in which applications of similar types are often assigned bands within the same block, dictates that wideband spectrum is indeed heterogeneous, as different application types exhibit different behaviors. In this paper, we consider heterogeneous wideband spectrum, where we exploit this inherent, block-like structure of wideband spectrum to design efficient compressive spectrum sensing techniques that are well suited for heterogeneous wideband spectrum. We propose a weighted $\ell_1$-minimization sensing information recovery algorithm that achieves more stable recovery than that achieved by existing approaches while accounting for the variations of spectrum occupancy across both the time and frequency dimensions. Through intensive numerical simulations, we show that our approach achieves better performance when compared to the state-of-the-art approaches.


title: The different growth pathways of Brightest Cluster Galaxies and the Intra-Cluster Light

abstract: We study the growth pathways of Brightest Central Galaxies (BCGs) and Intra-Cluster Light (ICL) by means of a semi-analytic model. We assume that the ICL forms by stellar stripping of satellite galaxies and violent processes during mergers, and implement two independent models: (1) one considers both mergers and stellar stripping (named {\small STANDARD} model), and one considers only mergers (named {\small MERGERS} model). We find that BCGs and ICL form, grow and overall evolve at different times and with different timescales, but they show a clear co-evolution after redshift $z \sim 0.7-0.8$. Around 90\% of the ICL from stellar stripping is built-up in the innermost 150 Kpc from the halo centre and the dominant contribution comes from disk-like galaxies (B/T$<$0.4) through a large number of small/intermediate stripping events ($M_{strip}/M_{sat}<0.3$). The fractions of stellar mass in BCGs and in ICL over the total stellar mass within the virial radius of the halo evolve differently with time. At high redshift, the BCG accounts for the bulk of the mass, but its contribution gradually decreases with time and stays constant after $z\sim 0.4-0.5$. The ICL, instead, grows very fast and its contribution keeps increasing down to the present time. The {\small STANDARD} and the {\small MERGERS} models make very similar predictions in most of the cases, but predict different amounts of ICL associated to other galaxies within the virial radius of the group/cluster other than the BCG, at $z=0$. We then suggest that this quantity is a valid observable that can shed light on the relative importance of mergers and stellar stripping for the formation of the ICL.


title: $\mathcal{N}=1$ $\mathcal{D}=3$ Lifshitz-Wess-Zumino model: A paradigm of reconciliation between Lifshitz-like operators and supersymmetry

abstract: By imposing the weighted renormalization condition and the (super)symmetry requirements, we construct a Lifshitz-like extension of the three-dimensional Wess-Zumino model, with dynamical critical exponent z=2. In this context, the auxiliary field F plays a key role by introducing the appropriate Lifshitz operator in the bosonic sector of the theory, avoiding so undesirable time-space mixing derivatives and inconsistencies concerning the critical z exponent, as reported in the literature. The consistency of the proposed model is verified by building explicitly the susy algebra through the Noether method in the canonical formalism. This component-field Lifshitz-Wess-Zumino model is in addition rephrased in the Lifshitz superspace, a natural modification of the conventional one. The one-loop effective potential is computed to study the possibility of symmetry breaking. It is found that supersymmetry remains intact at one-loop order, while the U(1) phase symmetry suffers a spontaneous breakdown above the critical value of the renormalization point. By renormalizing the one-loop effective potential within the cutoff regularization scheme, it is observed an improvement of the UV behavior of the theory compared with the relativistic Wess-Zumino model.


title: SoPhie: An Attentive GAN for Predicting Paths Compliant to Social and Physical Constraints

abstract: This paper addresses the problem of path prediction for multiple interacting agents in a scene, which is a crucial step for many autonomous platforms such as self-driving cars and social robots. We present \textit{SoPhie}; an interpretable framework based on Generative Adversarial Network (GAN), which leverages two sources of information, the path history of all the agents in a scene, and the scene context information, using images of the scene. To predict a future path for an agent, both physical and social information must be leveraged. Previous work has not been successful to jointly model physical and social interactions. Our approach blends a social attention mechanism with a physical attention that helps the model to learn where to look in a large scene and extract the most salient parts of the image relevant to the path. Whereas, the social attention component aggregates information across the different agent interactions and extracts the most important trajectory information from the surrounding neighbors. SoPhie also takes advantage of GAN to generates more realistic samples and to capture the uncertain nature of the future paths by modeling its distribution. All these mechanisms enable our approach to predict socially and physically plausible paths for the agents and to achieve state-of-the-art performance on several different trajectory forecasting benchmarks.


title: JTAV: Jointly Learning Social Media Content Representation by Fusing Textual, Acoustic, and Visual Features

abstract: Learning social media content is the basis of many real-world applications, including information retrieval and recommendation systems, among others. In contrast with previous works that focus mainly on single modal or bi-modal learning, we propose to learn social media content by fusing jointly textual, acoustic, and visual information (JTAV). Effective strategies are proposed to extract fine-grained features of each modality, that is, attBiGRU and DCRNN. We also introduce cross-modal fusion and attentive pooling techniques to integrate multi-modal information comprehensively. Extensive experimental evaluation conducted on real-world datasets demonstrates our proposed model outperforms the state-of-the-art approaches by a large margin.


title: 3D Human Pose Estimation with 2D Marginal Heatmaps

abstract: Automatically determining three-dimensional human pose from monocular RGB image data is a challenging problem. The two-dimensional nature of the input results in intrinsic ambiguities which make inferring depth particularly difficult. Recently, researchers have demonstrated that the flexible statistical modelling capabilities of deep neural networks are sufficient to make such inferences with reasonable accuracy. However, many of these models use coordinate output techniques which are memory-intensive, not differentiable, and/or do not spatially generalise well. We propose improvements to 3D coordinate prediction which avoid the aforementioned undesirable traits by predicting 2D marginal heatmaps under an augmented soft-argmax scheme. Our resulting model, MargiPose, produces visually coherent heatmaps whilst maintaining differentiability. We are also able to achieve state-of-the-art accuracy on publicly available 3D human pose estimation data.


title: Modeling of drug diffusion in a solid tumor leading to tumor cell death

abstract: It has been shown recently that changing the fluidic properties of a drug can improve its efficacy in ablating solid tumors. We develop a modeling framework for tumor ablation, and present the simplest possible model for drug diffusion in a spherical tumor with leaky boundaries and assuming cell death eventually leads to ablation of that cell effectively making the two quantities numerically equivalent. The death of a cell after a given exposure time depends on both the concentration of the drug and the amount of oxygen available to the cell. Higher oxygen availability leads to cell death at lower drug concentrations. It can be assumed that a minimum concentration is required for a cell to die, effectively connecting diffusion with efficacy. The concentration threshold decreases as exposure time increases, which allows us to compute dose-response curves. Furthermore, these curves can be plotted at much finer time intervals compared to that of experiments, which is used to produce a dose-threshold-response surface giving an observer a complete picture of the drug's efficacy for an individual. In addition, since the diffusion, leak coefficients, and the availability of oxygen is different for different individuals and tumors, we produce artificial replication data through bootstrapping to simulate error. While the usual data-driven model with Sigmoidal curves use 12 free parameters, our mechanistic model only has two free parameters, allowing it to be open to scrutiny rather than forcing agreement with data. Even so, the simplest model in our framework, derived here, shows close agreement with the bootstrapped curves, and reproduces well established relations, such as Haber's rule.


title: Forecasting Crime with Deep Learning

abstract: The objective of this work is to take advantage of deep neural networks in order to make next day crime count predictions in a fine-grain city partition. We make predictions using Chicago and Portland crime data, which is augmented with additional datasets covering weather, census data, and public transportation. The crime counts are broken into 10 bins and our model predicts the most likely bin for a each spatial region at a daily level. We train this data using increasingly complex neural network structures, including variations that are suited to the spatial and temporal aspects of the crime prediction problem. With our best model we are able to predict the correct bin for overall crime count with 75.6% and 65.3% accuracy for Chicago and Portland, respectively. The results show the efficacy of neural networks for the prediction problem and the value of using external datasets in addition to standard crime data.


title: Berry-Esseen bound for the Parameter Estimation of Fractional Ornstein-Uhlenbeck Processes

abstract: For an Ornstein-Uhlenbeck process driven by fractional Brownian motion with Hurst index $H\in [\frac12,\frac34]$, we show the Berry-Ess\'een bound of the least squares estimator of the drift parameter. We use an approach based on Malliavin calculus given by Kim and Park \cite{kim 3}.


title: A Primer on Causal Analysis

abstract: We provide a conceptual map to navigate causal analysis problems. Focusing on the case of discrete random variables, we consider the case of causal effect estimation from observational data. The presented approaches apply also to continuous variables, but the issue of estimation becomes more complex. We then introduce the four schools of thought for causal analysis


title: Unfolding of a diblock chain and its anomalous diffusion induced by active particles

abstract: We study the structural and dynamical behaviors of a diblock copolymer chain in a bath of active Brownian particles (ABPs) by extensive Brownian dynamics simulation in a two-dimensional model system. Specifically, the A block of chain is self-attractive, while the B block is self-repulsive. We find, beyond a threshold, the A block unfolds with a pattern like extracting a woolen string from a ball. The critical force decreases with the increase of the B block length,NB, for short cases, then keeps a constant with further increase of NB. In addition, we find a power law exists between the unfolding time of chain and active force, Fa, as well as NB. Finally, we focus on the translational and rotational diffusion of chain, and find that both of them remain supper-diffusive at the long time limit for small active forces due to an asymmetry distribution of ABPs. Our results open new routes for manipulating polymer's behaviors with ABPs.


title: Design of 32-channel TDC Based on Single FPGA for {\mu}SR Spectrometer at CSNS

abstract: Muon Spin Rotation, Relaxation and Resonance ({\mu}SR) technology has an irreplaceable role in studying the microstructure and properties of materials, especially micro-magnetic properties. An experimental muon source is being built in China Spallation Neutron Source (CSNS) now. At the same time, a 128-channel {\mu}SR spectrometer as China's first {\mu}SR spectrometer is being developed. The time spectrum of {\mu}SR can be obtained by fitting the curve of positron count rate with time. This paper presents a 32-channel Time-to-Digital Converter (TDC) implemented in a Xilinx Virtex-6 Field Programmable Gate Array (FPGA) for measuring the positron's flight time of {\mu}SR Spectrometer. Signal of each channel is sampled by 16 equidistant shifted-phase 200 MHz sampling clocks, so the TDC bin size is 312.5ps. The measuring range is up to 327us. This TDC has the ability to store multiple hit signals in a short time with a deep hit-buffer up to 512. Time tag is added to each data to record the moment when the data was detected. Programmable time window and channel shielding give the flexibility to choose the time range and channels of interest. The delay of each channel can be calibrated. The data is transmitted to data acquisition system (DAQ) through Gigabit Ethernet. TDC and control logic are configured in real time by DAQ. The results of test show that the Full Width at Half Maximum (FWHM) precision of single channel is better than 273 ps with a low sensitivity to temperature and the linearity is pretty well.


title: THz Photodetector using sideband-modulated transport through surface states of a 3D Topological Insulator

abstract: The transport properties of the surface charge carriers of a three dimensional topological insulator under a terahertz (THz) field along with a resonant double barrier structure is theoretically analyzed within the framework of Floquet theory to explore the possibility of using such a device for photodetection purpose. We show that due to the contribution of elastic and inelastic scattering processes in the resulting transmission sidebands are formed in the conductance spectrum in somewhat similar way as in an optical cavity and this information can be used to detect the frequency of an unknown THz radiation. The dependence of the conductance on the bias voltage, the effect of THz radiation on resonances and the influence of zero energy points on the transmission spectrum are also discussed.


title: Near-Optimal Time and Sample Complexities for Solving Discounted Markov Decision Process with a Generative Model

abstract: In this paper we consider the problem of computing an $\epsilon$-optimal policy of a discounted Markov Decision Process (DMDP) provided we can only access its transition function through a generative sampling model that given any state-action pair samples from the transition function in $O(1)$ time. Given such a DMDP with states $S$, actions $A$, discount factor $\gamma\in(0,1)$, and rewards in range $[0, 1]$ we provide an algorithm which computes an $\epsilon$-optimal policy with probability $1 - \delta$ where \emph{both} the time spent and number of sample taken are upper bounded by \[ O\left[\frac{|S||A|}{(1-\gamma)^3 \epsilon^2} \log \left(\frac{|S||A|}{(1-\gamma)\delta \epsilon} \right) \log\left(\frac{1}{(1-\gamma)\epsilon}\right)\right] ~. \] For fixed values of $\epsilon$, this improves upon the previous best known bounds by a factor of $(1 - \gamma)^{-1}$ and matches the sample complexity lower bounds proved in Azar et al. (2013) up to logarithmic factors. We also extend our method to computing $\epsilon$-optimal policies for finite-horizon MDP with a generative model and provide a nearly matching sample complexity lower bound.


title: Newton-Kantorovitch method for decoupled forward-backward stochastic differential equations

abstract: We present and prove a Newton-Kantorovitch method for solving decoupled forward-backward stochastic differential equations (FBSDEs) involving smooth coefficients with uniformly bounded derivatives. As Newton's method is required a suitable initial condition to converge, we show that such initial conditions are solutions of a linear backward stochastic differential equation. In addition, we show that converges linearly to the solution.


title: Leave-out estimation of variance components

abstract: We propose leave-out estimators of quadratic forms designed for the study of linear models with unrestricted heteroscedasticity. Applications include analysis of variance and tests of linear restrictions in models with many regressors. An approximation algorithm is provided that enables accurate computation of the estimator in very large datasets. We study the large sample properties of our estimator allowing the number of regressors to grow in proportion to the number of observations. Consistency is established in a variety of settings where plug-in methods and estimators predicated on homoscedasticity exhibit first-order biases. For quadratic forms of increasing rank, the limiting distribution can be represented by a linear combination of normal and non-central $\chi^2$ random variables, with normality ensuing under strong identification. Standard error estimators are proposed that enable tests of linear restrictions and the construction of uniformly valid confidence intervals for quadratic forms of interest. We find in Italian social security records that leave-out estimates of a variance decomposition in a two-way fixed effects model of wage determination yield substantially different conclusions regarding the relative contribution of workers, firms, and worker-firm sorting to wage inequality than conventional methods. Monte Carlo exercises corroborate the accuracy of our asymptotic approximations, with clear evidence of non-normality emerging when worker mobility between blocks of firms is limited.


title: Dynamic optimal contract under parameter uncertainty with risk averse agent and principal

abstract: We consider a continuous time Principal-Agent model on a finite time horizon, where we look for the existence of an optimal contract both parties agreed on. Contrary to the main stream, where the principal is modelled as risk-neutral, we assume that both the principal and the agent have exponential utility, and are risk averse with same risk awareness level. Moreover, the agent's quality is unknown and modelled as a filtering term in the problem, which is revealed as time passes by. The principal can not observe the agent's real action, but can only recommend action levels to the agent. Hence, we have a \textit{moral hazard} problem. In this setting, we give an explicit solution to the optimal contract problem.


title: Deep Image Compression via End-to-End Learning

abstract: We present a lossy image compression method based on deep convolutional neural networks (CNNs), which outperforms the existing BPG, WebP, JPEG2000 and JPEG as measured via multi-scale structural similarity (MS-SSIM), at the same bit rate. Currently, most of the CNNs based approaches train the network using a L2 loss between the reconstructions and the ground-truths in the pixel domain, which leads to over-smoothing results and visual quality degradation especially at a very low bit rate. Therefore, we improve the subjective quality with the combination of a perception loss and an adversarial loss additionally. To achieve better rate-distortion optimization (RDO), we also introduce an easy-to-hard transfer learning when adding quantization error and rate constraint. Finally, we evaluate our method on public Kodak and the Test Dataset P/M released by the Computer Vision Lab of ETH Zurich, resulting in averaged 7.81% and 19.1% BD-rate reduction over BPG, respectively.


title: Decomposability and time consistency of risk averse multistage programs

abstract: Two approaches to time consistency of risk averse multistage stochastic problems were discussed in the recent literature. In one approach certain properties of the cor-responding risk measure are postulated which imply its decomposability. The other approach deals directly with conditional optimality of solutions of the considered problem. The aim of this paper is to discuss a relation between these two approaches.


title: On the Galerkin approximation and strong norm bounds for the stochastic Navier-Stokes equations with multiplicative noise

abstract: We investigate the convergence of the Galerkin approximation for the stochastic Navier-Stokes equations in an open bounded domain $\mathcal{O}$ with the non-slip boundary condition. We prove that \begin{equation*} \mathbb{E} \left[ \sup_{t \in [0,T]} \phi_1(\lVert (u(t)-u^n(t)) \rVert^2_V) \right] \rightarrow 0 \end{equation*} as $n \rightarrow \infty$ for any deterministic time $T > 0$ and for a specified moment function $\phi_1(x)$ where $u^n(t,x)$ denotes the Galerkin approximation of the solution $u(t,x)$. Also, we provide a result on uniform boundedness of the moment $\mathbb{E} [ \sup_{t \in [0,T]} \phi(\lVert u(t) \rVert^2_V) ] $ where $\phi$ grows as a single logarithm at infinity. Finally, we summarize results on convergence of the Galerkin approximation up to a deterministic time $T$ when the $V$-norm is replaced by the $H$-norm.


title: Facilitating Exploration with Interaction Snapshots under High Latency

abstract: Latency is, unfortunately, a reality when working with large datasets. Guaranteeing imperceptible latency for interactivity is often prohibitively expensive: the application developer may be forced to migrate data processing engines or deal with complex error bounds on samples, and to limit the application to users with high network bandwidth. Instead of relying on the backend, we propose a simple UX design---interaction snapshots. Responses of requests from the interactions are asynchronously loaded in "snapshots". With interaction snapshots, users can interact concurrently while the snapshots load. Our user study participants found it useful not to have to wait for each result and easily navigate to prior snapshots. For latency up to 5 seconds, participants were able to complete extrema, threshold, and trend identification tasks with little negative impact.


title: Six-dimensional heavenly equation. Dressing scheme and the hierarchy

abstract: We consider six-dimensional heavenly equation as a reduction in the framework of general six-dimensional linearly degenerate dispersionless hierarchy. We characterise the reduction in terms of wave functions, introduce generating relation, Lax-Sato equations and develop the dressing scheme for the reduced hierarchy. Using the dressing scheme, we construct a class of solutions for six-dimensional heavenly equation in terms of implicit functions.


title: Information Aggregation via Dynamic Routing for Sequence Encoding

abstract: While much progress has been made in how to encode a text sequence into a sequence of vectors, less attention has been paid to how to aggregate these preceding vectors (outputs of RNN/CNN) into fixed-size encoding vector. Usually, a simple max or average pooling is used, which is a bottom-up and passive way of aggregation and lack of guidance by task information. In this paper, we propose an aggregation mechanism to obtain a fixed-size encoding with a dynamic routing policy. The dynamic routing policy is dynamically deciding that what and how much information need be transferred from each word to the final encoding of the text sequence. Following the work of Capsule Network, we design two dynamic routing policies to aggregate the outputs of RNN/CNN encoding layer into a final encoding vector. Compared to the other aggregation methods, dynamic routing can refine the messages according to the state of final encoding vector. Experimental results on five text classification tasks show that our method outperforms other aggregating models by a significant margin. Related source code is released on our github page.


title: Boredom-driven curious learning by Homeo-Heterostatic Value Gradients

abstract: This paper presents the Homeo-Heterostatic Value Gradients (HHVG) algorithm as a formal account on the constructive interplay between boredom and curiosity which gives rise to effective exploration and superior forward model learning. We envisaged actions as instrumental in agent's own epistemic disclosure. This motivated two central algorithmic ingredients: devaluation and devaluation progress, both underpin agent's cognition concerning intrinsically generated rewards. The two serve as an instantiation of homeostatic and heterostatic intrinsic motivation. A key insight from our algorithm is that the two seemingly opposite motivations can be reconciled---without which exploration and information-gathering cannot be effectively carried out. We supported this claim with empirical evidence, showing that boredom-enabled agents consistently outperformed other curious or explorative agent variants in model building benchmarks based on self-assisted experience accumulation.


title: Influence of the ferroelectric polarization on the electronic structure of BaTiO3 thin films

abstract: Micrometric domains of precise ferroelectric polarization have been written into a 20 nm thick epitaxial thin film of BaTiO3(001) (BTO) on a Nb doped SrTiO3 (STO) substrate using PiezoForce Microscopy (PFM). The domain dependent electronic structure has been studied using fully energy-filtered PhotoEmission Microscopy (PEEM) and synchrotron radiation. Shifts, induced by ferroelectric polarization, of up to 300 meV are observed in the work function of the sample. The surface is Ba-O terminated. Polarization-induced distortion of the electronic structure is observed in the valence band and on the Ba 3d, Ti 2p and O 1s core levels of BTO. Polarization dependent surface adsorption is observed. A simple electrostatic model based on net surface charge is not sufficient to explain the observed modifications in the electronic levels.


title: Graph Compression Using Pattern Matching Techniques

abstract: Graphs can be used to represent a wide variety of data belonging to different domains. Graphs can capture the relationship among data in an efficient way, and have been widely used. In recent times, with the advent of Big Data, there has been a need to store and compute on large data sets efficiently. However, considering the size of the data sets in question, finding optimal methods to store and process the data has been a challenge. Therefore, in this paper, we study different graph compression techniques and propose novel algorithms to do the same. Specifically, given a graph G = (V, E), where V is the set of vertices and E is the set of edges, and |V| = n, we propose techniques to compress the adjacency matrix representation of the graph. Our algorithms are based on finding patterns within the adjacency matrix data, and replacing the common patterns with specific markers. All the techniques proposed here are lossless compression of graphs. Based on the experimental results, it is observed that our proposed techniques achieve almost 70% compression as compared to adjacency matrix representation. The results show that large graphs can be efficiently stored in smaller memory and exploit the parallel processing power of compute nodes as well as efficiently transfer data between resources.


title: On Dualities and Equivalences Between Physical Theories

abstract: The main aim of this paper is to make a remark about the relation between (i) dualities between theories, as `duality' is understood in physics and (ii) equivalence of theories, as `equivalence' is understood in logic and philosophy. The remark is that in physics, two theories can be dual, and accordingly get called `the same theory', though we interpret them as disagreeing---so that they are certainly not equivalent, as `equivalent' is normally understood. So the remark is simple: but, I shall argue, worth stressing---since often neglected. My argument for this is based on the account of duality developed by De Haro: which is illustrated here with several examples, from both elementary physics and string theory. Thus I argue that in some examples, including in string theory, two dual theories disagree in their claims about the world. I also spell out how this remark implies a limitation of proposals (both traditional and recent) to understand theoretical equivalence as either logical equivalence or a weakening of it.


title: Attention Based Fully Convolutional Network for Speech Emotion Recognition

abstract: Speech emotion recognition is a challenging task for three main reasons: 1) human emotion is abstract, which means it is hard to distinguish; 2) in general, human emotion can only be detected in some specific moments during a long utterance; 3) speech data with emotional labeling is usually limited. In this paper, we present a novel attention based fully convolutional network for speech emotion recognition. We employ fully convolutional network as it is able to handle variable-length speech, free of the demand of segmentation to keep critical information not lost. The proposed attention mechanism can make our model be aware of which time-frequency region of speech spectrogram is more emotion-relevant. Considering limited data, the transfer learning is also adapted to improve the accuracy. Especially, it's interesting to observe obvious improvement obtained with natural scene image based pre-trained model. Validated on the publicly available IEMOCAP corpus, the proposed model outperformed the state-of-the-art methods with a weighted accuracy of 70.4% and an unweighted accuracy of 63.9% respectively.


title: Using the AIDA Language to Formally Organize Scientific Claims

abstract: Scientific communication still mainly relies on natural language written in scientific papers, which makes the described knowledge very difficult to access with automatic means. We can therefore only make limited use of formal knowledge organization methods to support researchers and other interested parties with features such as automatic aggregations, fact checking, consistency checking, question answering, and powerful semantic search. Existing approaches to solve this problem by improving the scientific communication methods have either very restricted coverage, require formal logic skills on the side of the researchers, or depend on unreliable machine learning for the formalization of knowledge. Here, I propose an approach to this problem that is general, intuitive, and flexible. It is based on a unique kind of controlled natural language, called AIDA, consisting of English sentences that are atomic, independent, declarative, and absolute. Such sentences can then serve as nodes in a network of scientific claims linked to publications, researchers, and domain elements. I present here some small studies on preliminary applications of this language. The results indicate that it is well accepted by users and provides a good basis for the creation of a knowledge graph of scientific findings.


title: Cooperative emission mediated by cooperative energy transfer to a plasmonic antenna

abstract: We develop a theory of cooperative emission mediated by cooperative energy transfer (CET) from an ensemble of quantum emitters (QE) to plasmonic antenna at a rate equal to the sum of individual QE-plasmon energy transfer rates. If the antenna radiation efficiency is sufficiently high, the transferred energy is radiated away at approximately the same cooperative rate that scales with the ensemble size. We derive explicit expressions, in terms of local fields, for cooperative Purcell factor and enhancement factor for power spectrum valid for plasmonic structures of any shape with characteristic size smaller than the radiation wavelength. The radiated power spectrum retains the plasmon resonance lineshape with overall amplitude scaling with the ensemble size. If QEs are located in a region with nearly constant plasmon local density of states (LDOS), e.g., inside a plasmonic nanocavity, we demonstrate that the CET rate scales linearly with the number of excited QEs, consistent with the experiment, and can be tuned in a wide range by varying the excitation power. For QEs distributed in an extended region saturating the plasmon mode volume, we show that the cooperative Purcell factor has universal form independent of the system size. The CET mechanism incorporates the plasmon LDOS enhancement as well, giving rise to possibilities of controlling the emission rate beyond field enhancement limits.


title: Embeddedness, Convexity, and Rigidity of Hypersurfaces in Product Spaces

abstract: We establish the following Hadamard--Stoker type theorem: Let $f:M^n\rightarrow\mathscr{H}^n\times\mathbb R$ be a complete connected hypersurface with positive definite second fundamental form, where $\mathscr H^n$ is a Hadamard manifold. If the height function of $f$ has a critical point, then it is an embedding and $M$ is homeomorphic to $\mathbb S^n$ or $\mathbb R^n.$ Furthermore, $f(M)$ bounds a convex set in $\mathscr{H}^n\times\mathbb R.$ In addition, it is shown that, except for the assumption on convexity, this result is valid for hypersurfaces in $\mathbb S^n\times\mathbb R$ as well. We apply these theorems to show that a compact connected hypersurface in $\mathbb Q_\epsilon^n\times\mathbb R$ ($\epsilon=\pm 1$) is a rotational sphere, provided it has either constant mean curvature and positive-definite second fundamental form or constant sectional curvature greater than $(\epsilon +1)/2.$ We also prove that, for $\bar M=\mathscr H^n$ or $\mathbb S^n,$ any connected proper hypersurface $f:M^n\rightarrow\bar M^n \times\mathbb R$ with positive semi-definite second fundamental form and height function with no critical points is embedded and isometric to $\Sigma^{n-1}\times\mathbb R,$ where $\Sigma^{n-1}\subset\bar M^n$ is convex and homeomorphic to $\mathbb S^{n-1}$ (for $\bar M^n=\mathscr H^n$ we assume further that $f$ is cylindrically bounded). Analogous theorems for hypersurfaces in warped product spaces $\mathbb R\times_\rho\mathscr H^n$ and $\mathbb R\times_\rho\mathbb S^n$ are obtained. In all of these results, the manifold $M^n$ is assumed to have dimension $n\ge 3.$


title: Anharmonicity-Induced Criticality of Collective Excitation in a Trapped Bose-Einstein Condensate

abstract: We investigate the low energy excitations of a dilute atomic Bose gas confined in a anharmonic trap interacting with repulsive forces. The dispersion law of both surface and compression modes are derived and analyzed for large numbers of atoms in the trap, which show two branches of excitation and appear a two critical value. For a upper limit, BEC can be unstable with respect to some specific collective excitation, while for the lower limit, the frequency of collective excitation under anharmonic influence can be effectively lower than that without anharmonicity. Our work reveals the key role played by the anharmonicity and interatomic forces which introduce a rich structure in the dynamic behavior of these new many-body systems.


title: Multi-orbital nature of the spin fluctuations in Sr$_2$RuO$_4$

abstract: The spin susceptibility of strongly correlated Sr$_2$RuO$_4$ is known to display a rich structure in reciprocal space, with a prominent peak at ${\bf Q}_i$=$(0.3,0.3,0)$. It is still debated if the resulting incommensurate spin-density-wave fluctuations foster unconventional superconductivity at low temperature or compete therewith. By means of density functional theory combined with dynamical mean-field theory, we reveal the realistic multi-orbital signature of the (dynamic) spin susceptibility beyond existing weak-coupling approaches. The experimental fluctuation spectrum up to 80 meV is confirmed by theory. Furthermore, the peak at ${\bf Q}_i$ is shown to carry nearly equal contributions from each of the Ru$(4d)$-$t_\mathrm{2g}$ orbitals, pointing to a cooperative effect resulting in the dominant spin fluctuations.