title: On the Throughput of Channels that Wear Out

abstract: This work investigates the fundamental limits of communication over a noisy discrete memoryless channel that wears out, in the sense of signal-dependent catastrophic failure. In particular, we consider a channel that starts as a memoryless binary-input channel and when the number of transmitted ones causes a sufficient amount of damage, the channel ceases to convey signals. Constant composition codes are adopted to obtain an achievability bound and the left-concave right-convex inequality is then refined to obtain a converse bound on the log-volume throughput for channels that wear out. Since infinite blocklength codes will always wear out the channel for any finite threshold of failure and therefore cannot convey information at positive rates, we analyze the performance of finite blocklength codes to determine the maximum expected transmission volume at a given level of average error probability. We show that this maximization problem has a recursive form and can be solved by dynamic programming. Numerical results demonstrate that a sequence of block codes is preferred to a single block code for streaming sources.


title: Searching for giga-Jansky fast radio bursts from the Milky Way with a global array of low-cost radio receivers

abstract: If fast radio bursts (FRBs) originate from galaxies at cosmological distances, then their all-sky rate implies that the Milky Way may host an FRB on average once every 30-1500 years. If FRBs repeat for decades or centuies, a local FRB could be active now. A typical Galactic FRB would produce a millisecond radio pulse with ~1 GHz flux density of ~3E10 Jy, comparable to the radio flux levels and frequencies of cellular communication devices (cell phones, Wi-Fi, GPS). We propose to search for Galactic FRBs using a global array of low-cost radio receivers. One possibility is to use the ~1GHz communication channel in cellular phones through a Citizens-Science downloadable application. Participating phones would continuously listen for and record candidate FRBs and would periodically upload information to a central data processing website, which correlates the incoming data from all participants, to identify the signature of a real, globe-encompassing, FRB from an astronomical distance. Triangulation of the GPS-based pulse arrival times reported from different locations will provide the FRB sky position, potentially to arc-second accuracy. Pulse arrival times from phones operating at diverse frequencies, or from an on-device de-dispersion search, will yield the dispersion measure (DM) which will indicate the FRB source distance within the Galaxy. A variant of this approach would be to use the built-in ~100 MHz FM-radio receivers present in cell phones for an FRB search at lower frequencies. Alternatively, numerous "software-defined radio" (SDR) devices, costing ~$10 US each, could be plugged into USB ports of personal computers around the world (particularly in radio quiet regions) to establish the global network of receivers.


title: Collisions of Dark Matter Axion Stars with Astrophysical Sources

abstract: If QCD axions form a large fraction of the total mass of dark matter, then axion stars could be very abundant in galaxies. As a result, collisions with each other, and with other astrophysical bodies, can occur. We calculate the rate and analyze the consequences of three classes of collisions, those occurring between a dilute axion star and: another dilute axion star, an ordinary star, or a neutron star. In all cases we attempt to quantify the most important astrophysical uncertainties; we also pay particular attention to scenarios in which collisions lead to collapse of otherwise stable axion stars, and possible subsequent decay through number changing interactions. Collisions between two axion stars can occur with a high total rate, but the low relative velocity required for collapse to occur leads to a very low total rate of collapses. On the other hand, collisions between an axion star and an ordinary star have a large rate, $\Gamma_\odot \sim 3000$ collisions/year/galaxy, and for sufficiently heavy axion stars, it is plausible that most or all such collisions lead to collapse. We identify in this case a parameter space which has a stable region and a region in which collision triggers collapse, which depend on the axion number ($N$) in the axion star, and a ratio of mass to radius cubed characterizing the ordinary star ($M_s/R_s^3$). Finally, we revisit the calculation of collision rates between axion stars and neutron stars, improving on previous estimates by taking cylindrical symmetry of the neutron star distribution into account. Collapse and subsequent decay through collision processes, if occurring with a significant rate, can affect dark matter phenomenology and the axion star mass distribution.


title: Unconstrained inverse quadratic programming problem

abstract: The paper covers a formulation of the inverse quadratic programming problem in terms of unconstrained optimization where it is required to find the unknown parameters (the matrix of the quadratic form and the vector of the quasi-linear part of the quadratic form) provided that approximate estimates of the optimal solution of the direct problem and those of the target function to be minimized in the form of pairs of values lying in the corresponding neighborhoods are only known. The formulation of the inverse problem and its solution are based on the least squares method. In the explicit form the inverse problem solution has been derived in the form a system of linear equations. The parameters obtained can be used for reconstruction of the direct quadratic programming problem and determination of the optimal solution and the extreme value of the target function, which were not known formerly. It is possible this approach opens new ways in over applications, for example, in neurocomputing and quadric surfaces fitting. Simple numerical examples have been demonstrated. A scenario in the Octave/MATLAB programming language has been proposed for practical implementation of the method.


title: New Multirectional Mean Value Inequality

abstract: We establish new and stronger inequality of Clarke-Ledyaev type by direct construction.


title: Parabolic problem with fractional time derivative with nonlocal and nonsingular Mittag-Leffler kernel

abstract: We prove H\"older regularity results for a class of nonlinear parabolic problem with fractional-time derivative with nonlocal and Mittag-Leffler nonsingular kernel. Existence of weak solutions via approximating solutions is proved. Moreover, the H\"{o}lder continuity of viscosity solutions is obtained. We get the similar results as those obtained by Allen (see {\url{https://arxiv.org/abs/1610.10073}})


title: Quantitative Analysis of Automatic Image Cropping Algorithms: A Dataset and Comparative Study

abstract: Automatic photo cropping is an important tool for improving visual quality of digital photos without resorting to tedious manual selection. Traditionally, photo cropping is accomplished by determining the best proposal window through visual quality assessment or saliency detection. In essence, the performance of an image cropper highly depends on the ability to correctly rank a number of visually similar proposal windows. Despite the ranking nature of automatic photo cropping, little attention has been paid to learning-to-rank algorithms in tackling such a problem. In this work, we conduct an extensive study on traditional approaches as well as ranking-based croppers trained on various image features. In addition, a new dataset consisting of high quality cropping and pairwise ranking annotations is presented to evaluate the performance of various baselines. The experimental results on the new dataset provide useful insights into the design of better photo cropping algorithms.


title: Low regularity Poincar\'e-Einstein metrics

abstract: We prove the existence of a $C^{1,1}$ conformally compact Einstein metric on the ball that has asymptotic sectional curvature decay to $-1$ plus terms of order $e^{-2r}$ where $r$ is the distance from any fixed compact set. This metric has no $C^2$ conformal compactification.


title: Accuracy of reconstruction of spike-trains with two near-colliding nodes

abstract: We consider a signal reconstruction problem for signals $F$ of the form $ F(x)=\sum_{j=1}^{d}a_{j}\delta\left(x-x_{j}\right),$ from their moments $m_k(F)=\int x^kF(x)dx.$ We assume $m_k(F)$ to be known for $k=0,1,\ldots,N,$ with an absolute error not exceeding $\epsilon > 0$. We study the "geometry of error amplification" in reconstruction of $F$ from $m_k(F),$ in situations where two neighboring nodes $x_i$ and $x_{i+1}$ near-collide, i.e $x_{i+1}-x_i=h \ll 1$. We show that the error amplification is governed by certain algebraic curves $S_{F,i},$ in the parameter space of signals $F$, along which the first three moments $m_0,m_1,m_2$ remain constant.


title: Noise Stability is computable and low dimensional

abstract: Questions of noise stability play an important role in hardness of approximation in computer science as well as in the theory of voting. In many applications, the goal is to find an optimizer of noise stability among all possible partitions of $\mathbb{R}^n$ for $n \geq 1$ to $k$ parts with given Gaussian measures $\mu_1,\ldots,\mu_k$. We call a partition $\epsilon$-optimal, if its noise stability is optimal up to an additive $\epsilon$. In this paper, we give an explicit, computable function $n(\epsilon)$ such that an $\epsilon$-optimal partition exists in $\mathbb{R}^{n(\epsilon)}$. This result has implications for the computability of certain problems in non-interactive simulation, which are addressed in a subsequent work.


title: Diluted banded random matrices: Scaling behavior of eigenfunction and spectral properties

abstract: We demonstrate that the normalised localization length $\beta$ of the eigenfunctions of diluted (sparse) banded random matrices follows the scaling law $\beta=x^*/(1+x^*)$. The scaling parameter of the model is defined as $x^*\propto(b_{eff}^2/N)^\delta$, where $b_{eff}$ is the average number of non-zero elements per matrix row, $N$ is the matrix size, and $\delta\sim 1$. Additionally, we show that $x^*$ also scales the spectral properties of the model (up to certain sparsity) characterized by the spacing distribution of eigenvalues.


title: Non interactive simulation of correlated distributions is decidable

abstract: A basic problem in information theory is the following: Let $\mathbf{P} = (\mathbf{X}, \mathbf{Y})$ be an arbitrary distribution where the marginals $\mathbf{X}$ and $\mathbf{Y}$ are (potentially) correlated. Let Alice and Bob be two players where Alice gets samples $\{x_i\}_{i \ge 1}$ and Bob gets samples $\{y_i\}_{i \ge 1}$ and for all $i$, $(x_i, y_i) \sim \mathbf{P}$. What joint distributions $\mathbf{Q}$ can be simulated by Alice and Bob without any interaction? Classical works in information theory by G{\'a}cs-K{\"o}rner and Wyner answer this question when at least one of $\mathbf{P}$ or $\mathbf{Q}$ is the distribution on $\{0,1\} \times \{0,1\}$ where each marginal is unbiased and identical. However, other than this special case, the answer to this question is understood in very few cases. Recently, Ghazi, Kamath and Sudan showed that this problem is decidable for $\mathbf{Q}$ supported on $\{0,1\} \times \{0,1\}$. We extend their result to $\mathbf{Q}$ supported on any finite alphabet. We rely on recent results in Gaussian geometry (by the authors) as well as a new \emph{smoothing argument} inspired by the method of \emph{boosting} from learning theory and potential function arguments from complexity theory and additive combinatorics.


title: Motion Deblurring in the Wild

abstract: The task of image deblurring is a very ill-posed problem as both the image and the blur are unknown. Moreover, when pictures are taken in the wild, this task becomes even more challenging due to the blur varying spatially and the occlusions between the object. Due to the complexity of the general image model we propose a novel convolutional network architecture which directly generates the sharp image.This network is built in three stages, and exploits the benefits of pyramid schemes often used in blind deconvolution. One of the main difficulties in training such a network is to design a suitable dataset. While useful data can be obtained by synthetically blurring a collection of images, more realistic data must be collected in the wild. To obtain such data we use a high frame rate video camera and keep one frame as the sharp image and frame average as the corresponding blurred image. We show that this realistic dataset is key in achieving state-of-the-art performance and dealing with occlusions.


title: Designing a Safe Autonomous Artificial Intelligence Agent based on Human Self-Regulation

abstract: There is a growing focus on how to design safe artificial intelligent (AI) agents. As systems become more complex, poorly specified goals or control mechanisms may cause AI agents to engage in unwanted and harmful outcomes. Thus it is necessary to design AI agents that follow initial programming intentions as the program grows in complexity. How to specify these initial intentions has also been an obstacle to designing safe AI agents. Finally, there is a need for the AI agent to have redundant safety mechanisms to ensure that any programming errors do not cascade into major problems. Humans are autonomous intelligent agents that have avoided these problems and the present manuscript argues that by understanding human self-regulation and goal setting, we may be better able to design safe AI agents. Some general principles of human self-regulation are outlined and specific guidance for AI design is given.


title: Dynamic wrinkling and strengthening of a filament in a viscous fluid

abstract: We investigate the wrinkling dynamics of an elastic filament immersed in a viscous fluid submitted to compression at a finite rate with experiments and by combining geometric nonlinearities, elasticity, and slender body theory. The drag induces a dynamic lateral reinforcement of the filament leading to growth of wrinkles that coarsen over time. We discover a new dynamical regime characterized by a timescale with a non-trivial dependence on the loading rate, where the growth of the instability is super-exponential and the wavenumber is an increasing function of the loading rate. We find that this timescale can be interpreted as the characteristic time over which the filament transitions from the extensible to the inextensible regime. In contrast with our analysis with moving boundary conditions, Biot's analysis in the limit of infinitely fast loading leads to rate independent exponential growth and wavelength.


title: A New Class of Exact Solutions in Einstein-Maxwell-dilaton Theory

abstract: We find new solutions to the five-dimensional Einstein-Maxwell-dilaton theory with cosmological constant where the dilaton field couples to the electromagnetic field as well as to the cosmological term with two different coupling constants. The five-dimensional spacetime is non-stationary and is a conformally regular spacetime, everywhere. Both the dilaton field and the electromagnetic field depends on time and two spatial directions. The cosmological constant takes positive, negative or zero value, depending on the value of coupling constant. We study the physical properties of the spacetime and show that the solutions are unique in five dimensions and can't be uplifted to higher-dimensional Einstein-Maxwell theory or Einstein gravity in presence of cosmological constant. Moreover, we construct new solutions to the theory where both coupling constants are equal.


title: Long-term optical spectral monitoring of NGC 7469

abstract: We present the results of the long-term (20-year period, from 1996 to 2015) optical spectral monitoring of the Seyfert 1 galaxy NGC 7469. The variation in the light-curves of the broad He II {\lambda}4686A H{\beta} and H{\alpha} lines, and the continuum at 5100A and 6300A have been explored. The maximum of activity was in 1998, and the variability in the continuum and lines seems to have two periods of around 1200 and 2600 days, however these periodicities should be taken with caution because of the red-noise. Beside these periods, there are several short-term (1-5 days) flare-like events in the light-curves. There are good correlations between the continuum fluxes and H{\alpha} and H{\beta} line fluxes, but significantly smaller correlation between the He II and continuum. We found that the time-lags between the continuum and broad lines are different for H{\beta} (~20 l.d.) and H{\alpha} (~3 l.d.), and that He II also has a smaller lag (~2-3 l.d.). The H{\alpha} and H{\beta} line profiles show a slight red asymmetry, and the broad line profiles did not changed in the 20-year period. Using the lags and widths of H{\alpha} and H{\beta} we estimated the central black hole mass and found that it is ~(1-6)x$10^7 M_{\odot}$, which is in agreement with previous reverberation estimates.


title: MDS-Coded Distributed Caching for Low Delay Wireless Content Delivery

abstract: We investigate the use of maximum distance separable (MDS) codes to cache popular content to reduce the download delay of wireless content delivery. In particular, we consider a cellular system where devices roam in an out of a cell according to a Poisson random process. Popular content is cached in a limited number of the mobile devices using an MDS code and can be downloaded from the mobile devices using device-to-device communication. We derive an analytical expression for the delay incurred in downloading content from the wireless network and show that distributed caching using MDS codes can dramatically reduce the download delay with respect to the scenario where content is always downloaded from the base station and to the case of uncoded distributed caching.


title: Stochastic Least-Squares Petrov-Galerkin Method for Parameterized Linear Systems

abstract: We consider the numerical solution of parameterized linear systems where the system matrix, the solution, and the right-hand side are parameterized by a set of uncertain input parameters. We explore spectral methods in which the solutions are approximated in a chosen finite-dimensional subspace. It has been shown that the stochastic Galerkin projection technique fails to minimize any measure of the solution error [20]. As a remedy for this, we propose a novel stochastic least-squares Petrov--Galerkin (LSPG) method. The proposed method is optimal in the sense that it produces the solution that minimizes a weighted l2-norm of the residual over all solutions in a given finite-dimensional subspace. Moreover, the method can be adapted to minimize the solution error in different weighted l2-norms by simply applying a weighting function within the least-squares formulation. In addition, a goal-oriented semi-norm induced by an output quantity of interest can be minimized by defining a weighting function as a linear functional of the solution. We establish optimality and error bounds for the proposed method, and extensive numerical experiments show that the weighted LSPG methods outperforms other spectral methods in minimizing corresponding target weighted norms.


title: Average values of L-functions in even characteristic

abstract: Let $k = \mathbb{F}_{q}(T)$ be the rational function field over a finite field $\mathbb{F}_{q}$, where $q$ is a power of $2$. In this paper we solve the problem of averaging the quadratic $L$-functions $L(s, \chi_{u})$ over fundamental discriminants. Any separable quadratic extension $K$ of $k$ is of the form $K = k(x_{u})$, where $x_{u}$ is a zero of $X^2+X+u=0$ for some $u\in k$. We characterize the family $\mathcal I$ (resp. $\mathcal F$, $\mathcal F'$) of rational functions $u\in k$ such that any separable quadratic extension $K$ of $k$ in which the infinite prime $\infty = (1/T)$ of $k$ ramifies (resp. splits, is inert) can be written as $K = k(x_{u})$ with a unique $u\in\mathcal I$ (resp. $u\in\mathcal F$, $u\in\mathcal F'$). For almost all $s\in\mathbb C$ with ${\rm Re}(s)\ge \frac{1}2$, we obtain the asymptotic formulas for the summation of $L(s,\chi_{u})$ over all $k(x_{u})$ with $u\in \mathcal I$, all $k(x_{u})$ with $u\in \mathcal F$ or all $k(x_{u})$ with $u\in \mathcal F'$ of given genus. As applications, we obtain the asymptotic mean value formulas of $L$-functions at $s=\frac{1}2$ and $s=1$ and the asymptotic mean value formulas of the class number $h_{u}$ or the class number times regulator $h_{u} R_{u}$.


title: Bifurcations of multi-vortex configurations in rotating Bose--Einstein condensates

abstract: We analyze global bifurcations along the family of radially symmetric vortices in the Gross--Pitaevskii equation with a symmetric harmonic potential and a chemical potential $\mu$ under the steady rotation with frequency $\Omega$. The families are constructed in the small-amplitude limit when the chemical potential $\mu$ is close to an eigenvalue of the Schr\"{o}dinger operator for a quantum harmonic oscillator. We show that for $\Omega$ near $0$, the Hessian operator at the radially symmetric vortex of charge $m_{0}\in\mathbb{N}$ has $m_{0}(m_{0}+1)/2$ pairs of negative eigenvalues. When the parameter $\Omega$ is increased, $1+m_{0}(m_{0}-1)/2$ global bifurcations happen. Each bifurcation results in the disappearance of a pair of negative eigenvalues in the Hessian operator at the radially symmetric vortex. The distributions of vortices in the bifurcating families are analyzed by using symmetries of the Gross--Pitaevskii equation and the zeros of Hermite--Gauss eigenfunctions. The vortex configurations that can be found in the bifurcating families are the asymmetric vortex $(m_0 = 1)$, the asymmetric vortex pair $(m_0 = 2)$, and the vortex polygons $(m_0 \geq 2)$.


title: Membrane-Dependent Neuromorphic Learning Rule for Unsupervised Spike Pattern Detection

abstract: Several learning rules for synaptic plasticity, that depend on either spike timing or internal state variables, have been proposed in the past imparting varying computational capabilities to Spiking Neural Networks. Due to design complications these learning rules are typically not implemented on neuromorphic devices leaving the devices to be only capable of inference. In this work we propose a unidirectional post-synaptic potential dependent learning rule that is only triggered by pre-synaptic spikes, and easy to implement on hardware. We demonstrate that such a learning rule is functionally capable of replicating computational capabilities of pairwise STDP. Further more, we demonstrate that this learning rule can be used to learn and classify spatio-temporal spike patterns in an unsupervised manner using individual neurons. We argue that this learning rule is computationally powerful and also ideal for hardware implementations due to its unidirectional memory access.


title: Benchmarks for single-phase flow in fractured porous media

abstract: This paper presents several test cases intended to be benchmarks for numerical schemes for single-phase fluid flow in fractured porous media. A number of solution strategies are compared, including a vertex and a cell-centered finite volume method, a non-conforming embedded discrete fracture model, a primal and a dual extended finite element formulation, and a mortar discrete fracture model. The proposed benchmarks test the schemes by increasing the difficulties in terms of network geometry, e.g. intersecting fractures, and physical parameters, e.g. low and high fracture-matrix permeability ratio as well as heterogeneous fracture permeabilities. For each problem, the results presented by the participants are the number of unknowns, the approximation errors in the porous matrix and in the fractures with respect to a reference solution, and the sparsity and condition number of the discretized linear system. All data and meshes used in this study are publicly available for further comparisons.


title: Learning local trajectories for high precision robotic tasks : application to KUKA LBR iiwa Cartesian positioning

abstract: To ease the development of robot learning in industry, two conditions need to be fulfilled. Manipulators must be able to learn high accuracy and precision tasks while being safe for workers in the factory. In this paper, we extend previously submitted work which consists in rapid learning of local high accuracy behaviors. By exploration and regression, linear and quadratic models are learnt for respectively the dynamics and cost function. Iterative Linear Quadratic Gaussian Regulator combined with cost quadratic regression can converge rapidly in the final stages towards high accuracy behavior as the cost function is modelled quite precisely. In this paper, both a different cost function and a second order improvement method are implemented within this framework. We also propose an analysis of the algorithm parameters through simulation for a positioning task. Finally, an experimental validation on a KUKA LBR iiwa robot is carried out. This collaborative robot manipulator can be easily programmed into safety mode, which makes it qualified for the second industry constraint stated above.


title: Efficiency versus instability in plasma accelerators

abstract: Plasma wake-field acceleration is one of the main technologies being developed for future high-energy colliders. Potentially, it can create a cost-effective path to the highest possible energies for e+e- or {\gamma}-{\gamma} colliders and produce a profound effect on the developments for high-energy physics. Acceleration in a blowout regime, where all plasma electrons are swept away from the axis, is presently considered to be the primary choice for beam acceleration. In this paper, we derive a universal efficiency-instability relation, between the power efficiency and the key instability parameter of the trailing bunch for beam acceleration in the blowout regime. We also show that the suppression of instability in the trailing bunch can be achieved through BNS damping by the introduction of a beam energy variation along the bunch. Unfortunately, in the high efficiency regime, the required energy variation is quite high, and is not presently compatible with collider-quality beams. We would like to stress that the development of the instability imposes a fundamental limitation on the acceleration efficiency, and it is unclear how it could be overcome for high-luminosity linear colliders. With minor modifications, the considered limitation on the power efficiency is applicable to other types of acceleration.


title: Off-Diagonal Expansion Quantum Monte Carlo

abstract: We propose a Monte Carlo algorithm designed to simulate quantum as well as classical systems at equilibrium, bridging the algorithmic gap between quantum and classical thermal simulation algorithms. The method is based on a novel decomposition of the quantum partition function that can be viewed as a series expansion about its classical part. We argue that the algorithm is optimally suited to tackle quantum many-body systems that exhibit a range of behaviors from `fully-quantum' to `fully-classical', in contrast to many existing methods. We demonstrate the advantages of the technique by comparing it against existing schemes. We also illustrate how our method allows for the unification of quantum and classical thermal parallel tempering techniques into a single algorithm and discuss its practical significance.


title: VideoSet: A Large-Scale Compressed Video Quality Dataset Based on JND Measurement

abstract: A new methodology to measure coded image/video quality using the just-noticeable-difference (JND) idea was proposed. Several small JND-based image/video quality datasets were released by the Media Communications Lab at the University of Southern California. In this work, we present an effort to build a large-scale JND-based coded video quality dataset. The dataset consists of 220 5-second sequences in four resolutions (i.e., $1920 \times 1080$, $1280 \times 720$, $960 \times 540$ and $640 \times 360$). For each of the 880 video clips, we encode it using the H.264 codec with $QP=1, \cdots, 51$ and measure the first three JND points with 30+ subjects. The dataset is called the "VideoSet", which is an acronym for "Video Subject Evaluation Test (SET)". This work describes the subjective test procedure, detection and removal of outlying measured data, and the properties of collected JND data. Finally, the significance and implications of the VideoSet to future video coding research and standardization efforts are pointed out. All source/coded video clips as well as measured JND data included in the VideoSet are available to the public in the IEEE DataPort.


title: Probing stochastic inter-galactic magnetic fields using blazar-induced gamma ray halo morphology

abstract: Inter-galactic magnetic fields can imprint their structure on the morphology of blazar-induced gamma ray halos. We show that the halo morphology arises through the interplay of the source's jet and a two-dimensional surface dictated by the magnetic field. Through extensive numerical simulations, we generate mock halos created by stochastic magnetic fields with and without helicity, and study the dependence of the halo features on the properties of the magnetic field. We propose a sharper version of the Q-statistics and demonstrate its sensitivity to the magnetic field strength, the coherence scale, and the handedness of the helicity. We also identify and explain a new feature of the Q-statistics that can further enhance its power.


title: Nonuniqueness of nematic liquid crystal flows in dimension three

abstract: For suitable initial and boundary data, we construct infinitely many weak solutions to the nematic liquid crystal flows in dimension three. These solutions are in the axisymmetric class with bounded energy and backward bubbling at a large time.


title: Log-Linear Bayesian Additive Regression Trees for Multinomial Logistic and Count Regression Models

abstract: We introduce Bayesian additive regression trees (BART) for log-linear models including multinomial logistic regression and count regression with zero-inflation and overdispersion. BART has been applied to nonparametric mean regression and binary classification problems in a range of settings. However, existing applications of BART have been limited to models for Gaussian "data", either observed or latent. This is primarily because efficient MCMC algorithms are available for Gaussian likelihoods. But while many useful models are naturally cast in terms of latent Gaussian variables, many others are not -- including models considered in this paper. We develop new data augmentation strategies and carefully specified prior distributions for these new models. Like the original BART prior, the new prior distributions are carefully constructed and calibrated to be flexible while guarding against overfitting. Together the new priors and data augmentation schemes allow us to implement an efficient MCMC sampler outside the context of Gaussian models. The utility of these new methods is illustrated with examples and an application to a previously published dataset.


title: A detailed analysis of impact collision ion scattering spectroscopy of bismuth selenide

abstract: Impact collision ion scattering spectroscopy (ICISS), which is a variation of low energy ion scattering (LEIS) that employs large scattering angles, is performed on Bi2Se3 surfaces prepared by ion bombardment and annealing (IBA). ICISS angular scans are collected experimentally and simulated numerically along the [120] and [-1 -2 0] azimuths, and the match of the positions of the flux peaks shows that the top three atomic layers are bulk-terminated. A newly observed feature is identified as a minimum in the multiple scattering background when the ion beam incidence is along a low index direction. Calculated scans as a function of scattering angle are employed to identify the behavior of flux peaks to show whether they originate from shadowing, blocking or both. This new method for analysis of large-angle LEIS data is shown to be useful for accurately investigating complex surface structures.


title: Crime Topic Modeling

abstract: The classification of crime into discrete categories entails a massive loss of information. Crimes emerge out of a complex mix of behaviors and situations, yet most of these details cannot be captured by singular crime type labels. This information loss impacts our ability to not only understand the causes of crime, but also how to develop optimal crime prevention strategies. We apply machine learning methods to short narrative text descriptions accompanying crime records with the goal of discovering ecologically more meaningful latent crime classes. We term these latent classes "crime topics" in reference to text-based topic modeling methods that produce them. We use topic distributions to measure clustering among formally recognized crime types. Crime topics replicate broad distinctions between violent and property crime, but also reveal nuances linked to target characteristics, situational conditions and the tools and methods of attack. Formal crime types are not discrete in topic space. Rather, crime types are distributed across a range of crime topics. Similarly, individual crime topics are distributed across a range of formal crime types. Key ecological groups include identity theft, shoplifting, burglary and theft, car crimes and vandalism, criminal threats and confidence crimes, and violent crimes. Though not a replacement for formal legal crime classifications, crime topics provide a unique window into the heterogeneous causal processes underlying crime.


title: Background Scale Independence in Quantum Gravity

abstract: We study the background scale independence in single-metric approximation to the functional renormalization group equation (FRGE) for quantum gravity and show that it is possible to formulate it without using higher-derivative gauge fixing in arbitrary dimensions if we adopt the Landau gauge and suitable cutoff scheme. We discuss this problem for both the linear and exponential splits of the metric into background and fluctuations. The obtained modified Ward identity for the global rescalings of the background metric can be combined with the FRGE to give a manifestly scale-invariant solution. An explicit example of the FRGE is given for four-dimensional $f(R)$ gravity in this framework.


title: One-Loop Free Energy of Tensionless Type IIB String in AdS$_5\times$S$^5$

abstract: Considering the zero 't Hooft coupling limit of ${\cal N}=4$ super-Yang-Mills theory, the exact spectrum of all single-trace operators can be accessed in terms of the underlying $so(2,4)$ character. This makes it possible in turn to compute the one-loop free energy of the tensionless type IIB string theory in AdS$_5\times$S$^5$ background, with help of the recently developed method of character integral representation of zeta function (CIRZ). We calculate first the one-loop free energy of the string states in the $(p-1)$-th Regge trajectory and find the result to be $p$ times the free energy of a single ${\cal N}=4$ Maxwell multiplet. The full one-loop free energy is hence proportional to the divergent series $\sum_{p=2}^\infty p\,$. The divergence arises as a result of interrupting the regularization procedure in an intermediate stage. With a reorganization of states, we extract the finite part of free energy after summing over the Regge trajectories. This way gives us a finite result which is minus of the free energy of the ${\cal N}=4$ multiplet. Hence, this bulk one-loop result matches the -1 term in the $N^2-1$ factor of the boundary result.


title: Neumann boundary conditions with null external quasi-momenta in finite-systems

abstract: The order parameter of a critical system defined in a layered parallel plate geometry subject to Neumann boundary conditions at the limiting surfaces is studied. We utilize a one-particle irreducible vertex parts framework in order to study the critical behavior of such a system. The renormalized vertex parts are defined at zero external quasi-momenta, which makes the analysis particularly simple. The distance between the boundary plates $L$ characterizing the finite size system direction perpendicular to the hyperplanes plays a similar role here in comparison with our recent unified treatment for Neumann and Dirichlet boundary conditions. Critical exponents are computed using diagrammatic expansion at least up to two-loop order and are shown to be identical to those from the bulk theory (limit $L \rightarrow \infty$).

categories: cond-mat.stat-mech hep-th math-ph math.MP


title: Spin Analogues of Superconductivity and the Integer Quantum Hall Effect in an Array of Spin Chains

abstract: Motivated by the successful idea of using weakly-coupled quantum electronic wires to realize the quantum Hall effects and the quantum spin Hall effects, we theoretically construct two systems composed of weakly-coupled quantum spin chains, which can exhibit spin analogues of superconductivity and the integer quantum Hall effect. Specifically, a certain bilayer of two arrays of interacting spin chains is mapped, via the Jordan-Wigner transformation, to a negative-$U$ Hubbard model that exhibits superconductivity. In addition, an array of spin-orbit-coupled spin chains in the presence of an suitable external magnetic field is transformed to an array of quantum wires that exhibits the integer quantum Hall effect. The resultant spin superconductivity and spin integer quantum Hall effect can be characterized by their ability to transport spin without any resistance.


title: Curvature dimension inequalities on directed graphs

abstract: In this paper, we define the curvature dimension inequalities CD(m, K) on finite directed graphs modifying the case of undirected graphs. As a main result, we evaluate m and K on finite directed graphs.


title: Quantifying the information transduction of biochemical reaction cascades: transcription of mRNA from ligand stimulation

abstract: A cell has the ability to convert an environmental change into the expression of genetic information through a chain of intracellular signal transduction reactions. Here, we aimed to develop a method for quantifying this signal transduction. We showed that the channel capacities of individual steps in a given general model cascade were equivalent in an independent manner, and were given by the entropy production rate. Signal transduction was transmitted by fluctuation of the entropy production rate and quantified transduction was estimated by the work done in individual steps. If the individual step representing the modification to demodification of the signal molecules is considered to be a Szilard engine, the maximal work done is equivalent to the chemical potential change of the messenger that is consumed during the modification reaction. Our method was applicable to calculate the channel capacity of the MAPK cascade. In conclusion, our method is suitable for quantitative analyses of signal transduction.


title: A hot Jupiter around the very active weak-line T Tauri star TAP 26

abstract: We report the results of an extended spectropolarimetric and photometric monitoring of the weak-line T Tauri star TAP 26, carried out within the MaTYSSE programme with the ESPaDOnS spectropolarimeter at the 3.6 m Canada-France-Hawaii Telescope. Applying Zeeman-Doppler Imaging to our observations, concentrating in 2015 November and 2016 January and spanning 72 d in total, 16 d in 2015 November and 13 d in 2016 January, we reconstruct surface brightness and magnetic field maps for both epochs and demonstrate that both distributions exhibit temporal evolution not explained by differential rotation alone. We report the detection of a hot Jupiter (hJ) around TAP 26 using three different methods, two using Zeeman-Doppler Imaging (ZDI) and one Gaussian-Process Regression (GPR), with a false-alarm probability smaller than 6.10^-4. However, as a result of the aliasing related to the observing window, the orbital period cannot be uniquely determined; the orbital period with highest likelihood is 10.79 +/- 0.14 d followed by 8.99 +/- 0.09 d. Assuming the most likely period, and that the planet orbits in the stellar equatorial plane, we obtain that the planet has a minimum mass M.sin(i) of 1.66 +/- 0.31 M_Jup and orbits at 0.0968 +/- 0.0032 au from its host star. This new detection suggests that disc type II migration is efficient at generating newborn hJs, and that hJs may be more frequent around young T Tauri stars than around mature stars (or that the MaTYSSE sample is biased towards hJ-hosting stars).


title: iMapD: intrinsic Map Dynamics exploration for uncharted effective free energy landscapes

abstract: We describe and implement iMapD, a computer-assisted approach for accelerating the exploration of uncharted effective Free Energy Surfaces (FES), and more generally for the extraction of coarse-grained, macroscopic information from atomistic or stochastic (here Molecular Dynamics, MD) simulations. The approach functionally links the MD simulator with nonlinear manifold learning techniques. The added value comes from biasing the simulator towards new, unexplored phase space regions by exploiting the smoothness of the (gradually, as the exploration progresses) revealed intrinsic low-dimensional geometry of the FES.


title: On classification and invariants of second order non-parabolic linear partial differential equations in two variables

abstract: The paper deals with second order abstract linear partial differential equations (LPDE) over a partial differential field with two commuting differential operators. In terms of usual differential equations the main content can be presented as follows. The classification and invariants problems for second order LPDEs with respect to transformations \[x=x(\xi,\eta),\ y=y(\xi,\eta),\ u=h(x,y)v(\xi,\eta),\] where $x,y$ are independent and $u$ is the dependent variable of the LPDE, are considered. Solutions to these problems are given for different subclasses of non-parabolic LPDE which appear naturally in the equivalence problem of LPDE. A criterion for reducibility of such LPDE to LPDE with constant coefficients is offered as well.


title: Analytic properties of American option prices under a modified Black-Scholes equation with spatial fractional derivatives

abstract: This paper investigates analytic properties of American option prices under the finite moment log-stable (FMLS) model. Under this model the price of American options is characterised by the free boundary problem of a fractional partial differential equation (FPDE) system. Using the technique of approximation we prove that the American put price under the FMLS model is convex with respect the underlying price, and specify the impact of the tail index on option prices.


title: Asai cube L-functions and the local Langlands conjecture

abstract: Let $F$ be a non-archimedean locally compact field. We study a class of Langlands-Shahidi pairs $({\bf H},{\bf L})$, consisting of a quasi-split connected reductive group $\bf H$ over $F$ and a Levi subgroup $\bf L$ which is closely related to a product of restriction of scalars of ${\rm GL}_1$'s or ${\rm GL}_2$'s. We prove the compatibility of the resulting local factors with the Langlands correspondence. In particular, let $E$ be a cubic separable extension of $F$. We consider a simply connected quasi-split semisimple group $\bf H$ over $F$ of type $D_4$, with triality corresponding to $E$, and let $\bf L$ be its Levi subgroup with derived group ${\rm Res}_{E/F} {\rm SL}_2$. In this way we obtain Asai cube local factors attached to irreducible smooth representations of ${\rm GL}_2(E)$; we prove that they are Weil-Deligne factors obtained via the local Langlands correspondence for ${\rm GL}_2(E)$ and tensor induction from $E$ to $F$. A consequence is that Asai cube $\gamma$- and $\varepsilon$-factors become stable under twists by highly ramified characters.


title: Optical vortex knots - one photon at a time

abstract: Feynman described the double slit experiment as "a phenomenon which is impossible, absolutely impossible, to explain in any classical way and which has in it the heart of quantum mechanics". The double slit experiment, performed one photon at a time, dramatically demonstrates the particle-wave duality of quantum objects by generating a fringe pattern corresponding to the interference of light (a wave phenomenon) from two slits, even when there is only one photon (a particle) at a time passing through the apparatus. The particle-wave duality of light should also apply to complex three dimensional optical fields formed by multi-path interference, however, this has not been demonstrated. Here we observe particle-wave duality of a three dimensional field by generating a trefoil optical vortex knot - one photon at a time. This result demonstrates a fundamental physical principle, that particle-wave duality implies interference in both space (between spatially distinct modes) and time (through the complex evolution of the superposition of modes), and has implications for topologically entangled single photon states, orbital angular momentum multiplexing and topological quantum computing.


title: Statistical properties of electrochemical capacitance in disordered mesoscopic capacitors

abstract: We numerically investigate the statistical properties of electrochemical capacitance in disordered two-dimensional mesoscopic capacitors. Based on the tight-binding Hamiltonian, the Green's function formalism is adopted to study the average electrochemical capacitance, its fluctuation as well as the distribution of capacitance of the disordered mesoscopic capacitors for three different ensembles: orthogonal (symmetry index \beta=1), unitary (\beta=2), and symplectic (\beta=4). It is found that the electrochemical capacitance in the disordered systems exhibits universal behavior. In the case of single conducting channel, the electrochemical capacitance follows a symmetric Gaussian distribution at weak disorders as expected from the random matrix theory. In the strongly disordered regime, the distribution is found to be a sharply one-sided form with a nearly-constant tail in the large capacitance region. This behavior is due to the existence of the necklace states in disordered systems, which is characterized by the multi-resonance that gives rise to a large density of states. In addition, it is found that the necklace state also enhances the fluctuation of electrochemical capacitance in the case of single conducting channel. When the number of conducting channels increases, the influence of necklace states becomes less important. For large number of conducting channels, the electrochemical capacitance fluctuation develops a plateau region in the strongly disordered regime. The plateau value is identified as universal electrochemical capacitance fluctuation, which is independent of system parameters such as disorder strength, Fermi energy, geometric capacitance, and system size. Importantly, the universal electrochemical capacitance fluctuation is the same for all three ensembles, suggesting a super-universal behavior.


title: Numerical study of parametric pumping current in mesoscopic systems in the presence of magnetic field

abstract: We numerically study the parametric pumped current when magnetic field is applied both in the adiabatic and non-adiabatic regimes. In particular, we investigate the nature of pumped current for systems with resonance as well as anti-resonance. It is found that in the adiabatic regime, the pumped current changes sign across the sharp resonance with long lifetime while the non-adiabatic pumped current at finite frequency does not. When the lifetime of resonant level is short, the behaviors of adiabatic and non-adiabatic pumped current are similar with sign changes. Our results show that at the energy where complete transmission occurs the adiabatic pumped current is zero while non-adiabatic pumped current is non-zero. Different from the resonant case, both adiabatic and non-adiabatic pumped current are zero at anti-resonance with complete reflection. We also investigate the pumped current when the other system parameters such as magnetic field, pumped frequency, and pumping potentials. Interesting behaviors are revealed. Finally, we study the symmetry relation of pumped current for several systems with different spatial symmetry upon reversal of magnetic field. Different from the previous theoretical prediction, we find that a system with general inversion symmetry can pump out a finite current in the adiabatic regime. At small magnetic field, the pumped current has an approximate relation I(B) \approx I(-B) both in adiabatic and non-adiabatic regimes.


title: Photoinduced valley-polarized current of layered MoS2 by electric tuning

abstract: A photoinduced current of a layered MoS2-based transistor is studied from first-principles. Under the illumination of circular polarized light, a valley-polarized current is generated, which can be tuned by the gate voltage. For monolayer MoS2, the valley-polarized spin-up (down) electron current at K (K') points is induced by the right (left) circular polarized light. The valley polarization is found to reach +1.0 (-1.0) for the valley current that carried such a K (K') index. For bilayer MoS2, the spin-up (down) current can be induced at both K and K' valleys by the right (left) circular light. In contrast to monolayer MoS2, the photoinduced valley polarization shows asymmetric behavior upon reversal of the gate voltage. Our results show that the valley polarization of the photoinduced current can be modulated by the circular polarized light and the gate voltage. All the results can be well understood using a simple kp model.


title: Valley Seebeck effect in gate tunable zigzag graphene nanoribbons

abstract: We propose, for the first time, a valley Seebeck effect in gate tunable zigzag graphene nanoribbons as a result of the interplay between thermal gradient and valleytronics. A pure valley current is further generated by the thermal gradient as well as the external bias. In a broad temperature range, the pure valley current is found to be linearly dependent on the temperature gradient while it increases with the increasing temperature of one lead for a fixed thermal gradient. A valley field effect transistor (FET) driven by the temperature gradient is proposed that can turn on and off the pure valley current by gate voltage. The threshold gate voltage and on valley current are proportional to the temperature gradient. When the system switches on at positive gate voltage, the pure valley current is nearly independent of gate voltage. The valley transconductance is up to 30 {\mu}S if we take Ampere as the unit of the valley current. This valley FET may find potential application in future valleytronics and valley caloritronics.


title: Breaking through the bandwidth barrier in distributed fiber vibration sensing by sub-Nyquist randomized sampling

abstract: The round trip time of the light pulse limits the maximum detectable frequency response range of vibration in phase-sensitive optical time domain reflectometry ({\phi}-OTDR). We propose a method to break the frequency response range restriction of {\phi}-OTDR system by modulating the light pulse interval randomly which enables a random sampling for every vibration point in a long sensing fiber. This sub-Nyquist randomized sampling method is suits for detecting sparse-wideband-frequency vibration signals. Up to MHz resonance vibration signal with over dozens of frequency components and 1.153MHz single frequency vibration signal are clearly identified for a sensing range of 9.6km with 10kHz maximum sampling rate.


title: Improving Student Understanding of Magnetism

abstract: Formally investigating the sources of students' difficulties around specific subjects is crucial for developing appropriate strategies to help students. We have been studying difficulties in understanding magnetism encountered by students in the calculus-based introductory physics courses. A majority of these students are engineering majors. Student difficulties were assessed by administering written free-response questions and by developing and administering a research-based multiple-choice test to students. We also conducted individual interviews with a subset of students. Some of these interviews were lecture-demonstration based interviews in which students were asked to predict the outcomes of experiments, perform the experiments and reconcile the differences between their predictions and observations. Some of the common misconceptions found in magnetism are analogous to those found in electrostatics. Some additional difficulties are due to the non-intuitive three dimensional nature of the relation between magnetic field, magnetic force and velocity of the charged particles or direction of current. Another finding is that students often used their gut feeling and had more difficulty reconciling with the idea that the magnetic force and field are perpendicular to each other when they were shown actual lecture-demonstration setups and asked to predict outcomes of experiments than when they were asked to explain relation between magnetic force and field theoretically solely based upon an equation. We discuss the implications of this research to teaching and learning.


title: Advantages of Unfair Quantum Ground-State Sampling

abstract: The debate around the potential superiority of quantum annealers over their classical counterparts has been ongoing since the inception of the field by Kadowaki and Nishimori close to two decades ago. Recent technological breakthroughs in the field, which have led to the manufacture of experimental prototypes of quantum annealing optimizers with sizes approaching the practical regime, have reignited this discussion. However, the demonstration of quantum annealing speedups remains to this day an elusive albeit coveted goal. Here, we examine the power of quantum annealers to provide a different type of quantum enhancement of practical relevance, namely, their ability to serve as useful samplers from the ground-state manifolds of combinatorial optimization problems. We study, both numerically by simulating ideal stoquastic and non-stoquastic quantum annealing processes, and experimentally, using a commercially available quantum annealing processor, the ability of quantum annealers to sample the ground-states of spin glasses differently than classical thermal samplers. We demonstrate that i) quantum annealers in general sample the ground-state manifolds of spin glasses very differently than thermal optimizers, ii) the nature of the quantum fluctuations driving the annealing process has a decisive effect on the final distribution over ground-states, and iii) the experimental quantum annealer samples ground-state manifolds significantly differently than thermal and ideal quantum annealers. We illustrate how quantum annealers may serve as powerful tools when complementing standard sampling algorithms.


title: Dynamic Base Station Repositioning to Improve Performance of Drone Small Cells

abstract: Thanks to the recent advancements in drone technology, it has become viable and cost-effective to quickly deploy small cells in areas of urgent needs by using a drone as a cellular base station. In this paper, we explore the benefit of dynamically repositioning the drone base station in the air to reduce the distance between the BS and the mobile user equipment, thereby improving the spectral efficiency of the small cell. In particular, we propose algorithms to autonomously control the repositioning of the drone in response to users activities and movements. We demonstrate that, compared to a drone hovering at a fixed location, dynamic repositioning of the drone moving with high speed can increase spectral efficiency by 15\%. However, considering the tradeoff between spectral efficiency and energy efficiency of the drone, we show that 10.5\% spectral efficiency gain can be obtained without negatively affecting energy consumption of the drone.


title: Magnetic separatrix as the source region of the plasma supply for an active-region filament

abstract: Solar filaments can be formed via chromospheric evaporation followed by condensation in the corona or by the direct injection of cool plasma from the chromosphere to the corona. In this paper, with high-resolution H$\alpha$ data observed by the 1.6 m New Solar Telescope of the Big Bear Solar Observatory on 2015 August 21, we confirmed that an active-region filament is maintained by the continuous injection of cold chromospheric plasma. We find that the filament is rooted along a bright ridge in H$\alpha$, which corresponds to the intersection of a magnetic quasi-separatrix layer with the solar surface. This bright ridge consists of many small patches and the sizes of these patches are comparable to the width of the filament threads. It is found that upflows originate from the brighter patches of the ridge, whereas the downflows move toward the weaker patches of the ridge. The whole filament is composed of two opposite directional streams, implying that longitudinal oscillations are not the only cause of the counterstreamings, and unidirectional siphon flows with alternative directions are another possibility.


title: Coordinated Autonomous Vehicle Parking for Vehicle-to-Grid Services: Formulation and Distributed Algorithm

abstract: Autonomous vehicles (AVs) will revolutionarize ground transport and take a substantial role in the future transportation system. Most AVs are likely to be electric vehicles (EVs) and they can participate in the vehicle-to-grid (V2G) system to support various V2G services. Although it is generally infeasible for EVs to dictate their routes, we can design AV travel plans to fulfill certain system-wide objectives. In this paper, we focus on the AVs looking for parking and study how they can be led to appropriate parking facilities to support V2G services. We formulate the Coordinated Parking Problem (CPP), which can be solved by a standard integer linear program solver but requires long computational time. To make it more practical, we develop a distributed algorithm to address CPP based on dual decomposition. We carry out a series of simulations to evaluate the proposed solution methods. Our results show that the distributed algorithm can produce nearly optimal solutions with substantially less computational time. A coarser time scale can improve computational time but degrade the solution quality resulting in possible infeasible solution. Even with communication loss, the distributed algorithm can still perform well and converge with only little degradation in speed.


title: Examples of descent up to nilpotence

abstract: We give a survey of the ideas of descent and nilpotence. We focus on examples arising from chromatic homotopy theory and from group actions, as well as a few examples in algebra.


title: Abstract Wiener measure using abelian Yang-Mills action on $\mathbb{R}^4$

abstract: Let $\mathfrak{g}$ be the Lie algebra of a compact Lie group. For a $\mathfrak{g}$-valued 1-form $A$, consider the Yang-Mills action \begin{equation} S_{{\rm YM}}(A) = \int_{\mathbb{R}^4} \left|dA + A \wedge A \right|^2\ d\omega \nonumber \end{equation} using the standard metric on $T\mathbb{R}^4$. When we consider the Lie group ${\rm U}(1)$, the Lie algebra $\mathfrak{g}$ is isomorphic to $\mathbb{R} \otimes i$, thus $A \wedge A = 0$. For some simple closed loop $C$, we want to make sense of the following path integral, \begin{equation} \frac{1}{Z}\ \int_{A \in \mathcal{A} /\mathcal{G}} \exp \left[ \int_{C} A\right] e^{-\frac{1}{2}\int_{\mathbb{R}^4}|dA|^2\ d\omega}\ DA, \nonumber \end{equation} whereby $DA$ is some Lebesgue type of measure on the space of $\mathfrak{g}$-valued 1-forms, modulo gauge transformations, $\mathcal{A} /\mathcal{G}$, and $Z$ is some partition function. We will construct an Abstract Wiener space for which we can define the above Yang-Mills path integral rigorously, using renormalization techniques found in lattice gauge theory. We will further show that the Area Law formula do not hold in the abelian Yang-Mills theory.


title: On the Reliability Function of the Common-Message Broadcast Channel with Variable-Length Feedback

abstract: We derive upper and lower bounds on the reliability function for the common-message discrete memoryless broadcast channel with variable-length feedback. We show that the bounds are tight when the broadcast channel is stochastically degraded. For the achievability part, we adapt Yamamoto and Itoh's coding scheme by controlling the expectation of the maximum of a set of stopping times. For the converse part, we adapt Burnashev's proof techniques for establishing the reliability functions for (point-to-point) discrete memoryless channels with variable-length feedback and sequential hypothesis testing.


title: Control of Dynamics in Brain Networks

abstract: The ability to effectively control brain dynamics holds great promise for the enhancement of cognitive function in humans, and the betterment of their quality of life. Yet, successfully controlling dynamics in neural systems is challenging, in part due to the immense complexity of the brain and the large set of interactions that can drive any single change. While we have gained some understanding of the control of single neurons, the control of large-scale neural systems -- networks of multiply interacting components -- remains poorly understood. Efforts to address this gap include the construction of tools for the control of brain networks, mostly adapted from control and dynamical systems theory. Informed by current opportunities for practical intervention, these theoretical contributions provide models that draw from a wide array of mathematical approaches. We present intriguing recent developments for effective strategies of control in dynamic brain networks, and we also describe potential mechanisms that underlie such processes. We review efforts in the control of general neurophysiological processes with implications for brain development and cognitive function, as well as the control of altered neurophysiological processes in medical contexts such as anesthesia administration, seizure suppression, and deep-brain stimulation for Parkinson's disease. We conclude with a forward-looking discussion regarding how emerging results from network control -- especially approaches that deal with nonlinear dynamics or more realistic trajectories for control transitions -- could be used to directly address pressing questions in neuroscience.

categories: q-bio.QM cond-mat.dis-nn cond-mat.soft physics.bio-ph q-bio.NC


title: Suboptimum Low Complexity Joint Multi-target Detection and Localization for Noncoherent MIMO Radar with Widely Separated Antennas

abstract: In this paper, the problems of simultaneously detecting and localizing multiple targets are considered for noncoherent multiple-input multiple-output (MIMO) radar with widely separated antennas. By assuming a prior knowledge of target number, an optimal solution to this problem is presented first. It is essentially a maximum-likelihood (ML) estimator searching parameters of interest in a high dimensional space. However, the complexity of this method increases exponentially with the number G of targets.Besides, without the prior information of the number of targets, a multi-hypothesis testing strategy to determine the number of targets is required, which further complicates this method. Therefore, we split the joint maximization into G disjoint optimization problems by clearing the interference from previously declared targets. In this way, we derive two fast and robust suboptimal solutions which allow trading performance for a much lower implementation complexity which is almost independent of the number of targets. In addition, the multi-hypothesis testing is no longer required when target number is unknown. Simulation results show the proposed algorithms can correctly detect and accurately localize multiple targets even when targets share common range bins in some paths.


title: CENTURION: Incentivizing Multi-Requester Mobile Crowd Sensing

abstract: The recent proliferation of increasingly capable mobile devices has given rise to mobile crowd sensing (MCS) systems that outsource the collection of sensory data to a crowd of participating workers that carry various mobile devices. Aware of the paramount importance of effectively incentivizing participation in such systems, the research community has proposed a wide variety of incentive mechanisms. However, different from most of these existing mechanisms which assume the existence of only one data requester, we consider MCS systems with multiple data requesters, which are actually more common in practice. Specifically, our incentive mechanism is based on double auction, and is able to stimulate the participation of both data requesters and workers. In real practice, the incentive mechanism is typically not an isolated module, but interacts with the data aggregation mechanism that aggregates workers' data. For this reason, we propose CENTURION, a novel integrated framework for multi-requester MCS systems, consisting of the aforementioned incentive and data aggregation mechanism. CENTURION's incentive mechanism satisfies truthfulness, individual rationality, computational efficiency, as well as guaranteeing non-negative social welfare, and its data aggregation mechanism generates highly accurate aggregated results. The desirable properties of CENTURION are validated through both theoretical analysis and extensive simulations.


title: On approximation of Ginzburg-Landau minimizers by $\mathbb S^1$-valued maps in domains with vanishingly small holes

abstract: We consider a two-dimensional Ginzburg-Landau problem on an arbitrary domain with a finite number of vanishingly small circular holes. A special choice of scaling relation between the material and geometric parameters (Ginzburg-Landau parameter vs hole radius) is motivated by a recently dsicovered phenomenon of vortex phase separation in superconducting composites. We show that, for each hole, the degrees of minimizers of the Ginzburg-Landau problems in the classes of $\mathbb S^1$-valued and $\mathbb C$-valued maps, respectively, are the same. The presence of two parameters that are widely separated on a logarithmic scale constitutes the principal difficulty of the analysis that is based on energy decomposition techniques.


title: High-Assurance Separation Kernels: A Survey on Formal Methods

abstract: Separation kernels provide temporal/spatial separation and controlled information flow to their hosted applications. They are introduced to decouple the analysis of applications in partitions from the analysis of the kernel itself. More than 20 implementations of separation kernels have been developed and widely applied in critical domains, e.g., avionics/aerospace, military/defense, and medical devices. Formal methods are mandated by the security/safety certification of separation kernels and have been carried out since this concept emerged. However, this field lacks a survey to systematically study, compare, and analyze related work. On the other hand, high-assurance separation kernels by formal methods still face big challenges. In this paper, an analytical framework is first proposed to clarify the functionalities, implementations, properties and standards, and formal methods application of separation kernels. Based on the proposed analytical framework, a taxonomy is designed according to formal methods application, functionalities, and properties of separation kernels. Research works in the literature are then categorized and overviewed by the taxonomy. In accordance with the analytical framework, a comprehensive analysis and discussion of related work are presented. Finally, four challenges and their possible technical directions for future research are identified, e.g. specification bottleneck, multicore and concurrency, and automation of full formal verification.


title: Quark orbital dynamics in the proton from Lattice QCD -- from Ji to Jaffe-Manohar orbital angular momentum

abstract: Given a Wigner distribution simultaneously characterizing quark transverse positions and momenta in a proton, one can directly evaluate their cross-product, i.e., quark orbital angular momentum. The aforementioned distribution can be obtained by generalizing the proton matrix elements of quark bilocal operators which define transverse momentum-dependent parton distributions (TMDs); the transverse momentum information is supplemented with transverse position information by introducing an additional nonzero momentum transfer. A gauge connection between the quarks must be specified in the quark bilocal operators; the staple-shaped gauge link path used in TMD calculations yields the Jaffe-Manohar definition of orbital angular momentum, whereas a straight path yields the Ji definition. An exploratory lattice calculation, performed at the pion mass m_pi = 518 MeV, is presented which quasi-continuously interpolates between the two definitions and demonstrates that their difference can be clearly resolved. The resulting Ji orbital angular momentum is confronted with traditional evaluations based on Ji's sum rule. Jaffe-Manohar orbital angular momentum is enhanced in magnitude compared to its Ji counterpart.


title: A Novel Quantum Image Compression Method Based on JPEG

abstract: Quantum image processing has been a hot topic. The first step of it is to store an image into qubits, which is called quantum image preparation. Different quantum image representations may have different preparation methods. In this paper, we use GQIR (the generalized quantum image representation) to represent an image, and try to decrease the operations used in preparation, which is also known as quantum image compression. Our compression scheme is based on JPEG (named from its inventor: the Joint Photographic Experts Group) --- the most widely used method for still image compression in classical computers. We input the quantized JPEG coefficients into qubits and then convert them into pixel values. Theoretical analysis and experimental results show that the compression ratio of our scheme is obviously higher than that of the previous compression method.


title: An Explicit Determination of the Springer Morphism

abstract: Let $G$ be a simply connected semisimple algebraic group over $\mathbb{C}$ and let $\rho :G\rightarrow GL(V_\lambda)$ be an irreducible representation of highest weight $\lambda$. Suppose that $\rho$ has finite kernel. Springer defined adjoint-invariant regular map with Zariski dense image from the group to its Lie algebra, $\theta_\lambda:G\rightarrow\mathfrak{g}$, which depends on $\lambda$ [Kumar]. By a lemma in Kumar's recent paper, $\theta_\lambda$ takes the maximal torus to its Lie algebra $\mathfrak{t}$. Thus, for a given simple group $G$ and an irreducible representation $V_\lambda$, one may write $\theta_\lambda (t)=\sum\limits_{i=1}^n c_i(t)\check{\alpha_i}$, where the simple co-roots $\{\check{\alpha_i}\}$ are a basis for $\mathfrak{t}$. We give a complete determination of these coefficients $c_i(t)$ for any simple group $G$ as a sum over the weights of the torus action on $V_\lambda$.


title: Algorithms for Optimal Replica Placement Under Correlated Failure in Hierarchical Failure Domains

abstract: In data centers, data replication is the primary method used to ensure availability of customer data. To avoid correlated failure, cloud storage infrastructure providers model hierarchical failure domains using a tree, and avoid placing a large number of data replicas within the same failure domain (i.e. on the same branch of the tree). Typical best practices ensure that replicas are distributed across failure domains, but relatively little is known concerning optimization algorithms for distributing data replicas. Using a hierarchical model, we answer how to distribute replicas across failure domains optimally. We formulate a novel optimization problem for replica placement in data centers. As part of our problem, we formalize and explain a new criterion for optimizing a replica placement. Our overall goal is to choose placements in which correlated failures disable as few replicas as possible. We provide two optimization algorithms for dependency models represented by trees. We first present an $O(n + \rho \log \rho)$ time dynamic programming algorithm for placing $\rho$ replicas of a single file on the leaves (representing servers) of a tree with $n$ vertices. We next consider the problem of placing replicas of $m$ blocks of data, where each block may have different replication factors. For this problem, we give an exact algorithm which runs in polynomial time when the skew, the difference in the number of replicas between the largest and smallest blocks of data, is constant.


title: Exact Computation of Influence Spread by Binary Decision Diagrams

abstract: Evaluating influence spread in social networks is a fundamental procedure to estimate the word-of-mouth effect in viral marketing. There are enormous studies about this topic; however, under the standard stochastic cascade models, the exact computation of influence spread is known to be #P-hard. Thus, the existing studies have used Monte-Carlo simulation-based approximations to avoid exact computation. We propose the first algorithm to compute influence spread exactly under the independent cascade model. The algorithm first constructs binary decision diagrams (BDDs) for all possible realizations of influence spread, then computes influence spread by dynamic programming on the constructed BDDs. To construct the BDDs efficiently, we designed a new frontier-based search-type procedure. The constructed BDDs can also be used to solve other influence-spread related problems, such as random sampling without rejection, conditional influence spread evaluation, dynamic probability update, and gradient computation for probability optimization problems. We conducted computational experiments to evaluate the proposed algorithm. The algorithm successfully computed influence spread on real-world networks with a hundred edges in a reasonable time, which is quite impossible by the naive algorithm. We also conducted an experiment to evaluate the accuracy of the Monte-Carlo simulation-based approximation by comparing exact influence spread obtained by the proposed algorithm.


title: Simple way to calculate UV-finite one-loop quantum energy in Randall-Sundrum model

abstract: The surprising simplicity of Barvinsky-Nesterov or equivalently Gelfand-Yaglom methods of calculation of quantum determinants permits to obtain compact expressions for UV-finite difference of one-loop quantum energies for two arbitrary values of parameter of the double-trace asymptotic boundary conditions. This result generalizes Gubser and Mitra calculation for particular case of difference of "regular" and "irregular" one-loop energies in one-brane RS-model. Approach developed in the paper also allows to get "in one line" the one-loop quantum energies in two-brane RS-model. The relationship between "one-loop" expressions corresponding to mixed Robin and to double-trace asymptotic boundary conditions is traced.


title: Quantization of vector bundles on Lagrangian subvarieties

abstract: We consider a smooth Lagrangian subvariety Y in a smooth algebraic variety X with an algebraic symplectic from. For a vector bundle E on Y and a choice Oh of deformation quantization of the structure sheaf of X, we establish when E admits a deformation quantization to a module over Oh. If the necessary conditions hold, we describe the set of equivalence classes of such quantizations.


title: Gap engineering in strained fold-like armchair graphene nanoribbons

abstract: Strain fold-like deformations on armchair graphene nanoribbons (AGNRs) can be properly engineered in experimental setups, and could lead to a new controlling tool for gaps and transport properties. Here, we analyze the electronic properties of folded AGNRs relating the electronic responses and the mechanical deformation. An important and universal parameter for the gap engineering is the ribbon percent width variation, i.e., the difference between the deformed and undeformed ribbon widths. AGNRs bandgap can be tuned mechanically in a well defined bounded range of energy values, eventually leading to a metallic system. This characteristic provides a new controllable degree of freedom that allows manipulation of electronic currents. We show that the numerical results are analytically predicted by solving the Dirac equation for the strained system.


title: Dense Small Cell Networks: From Noise-Limited to Dense Interference-Limited

abstract: Considering both non-line-of-sight (NLoS) and line-of-sight (LoS) transmissions, the transitional behaviors from noise-limited regime to dense interference-limited regime have been investigated for the fifth generation (5G) small cell networks (SCNs). Besides, we identify four performance regimes based on base station (BS) density, i.e., (i) the noise-limited regime, (ii) the signal-dominated regime, (iii) the interference-dominated regime, and (iv) the interference-limited regime. To characterize the performance regime, we propose a unified framework analyzing the future 5G wireless networks over generalized shadowing/fading channels, in which the user association schemes based on the strongest instantaneous received power (SIRP) and the strongest average received power (SARP) can be studied, while NLoS/LoS transmissions and multi-slop path loss model are considered. Simulation results indicate that different factors, i.e., noise, desired signal, and interference, successively and separately dominate the network performance with the increase of BS density. Hence, our results shed new light on the design and management of SCNs in urban and rural areas with different BS deployment densities.


title: Continuous body dynamics and the Mathisson-Papapetrou-Dixon equations

abstract: We show that an effective particle Lagrangian yields the Mathisson-Papapetrou-Dixon (MPD) equations. The spin of the effective particle is defined without any reference to a fixed body frame or angular velocity variable. We then demonstrate that a continuous body, defined by a congruence of world lines and described by a general action, can be rewritten as an effective particle. We analyze the gauge freedom of the body and show that a natural center of mass condition is related to a spin supplementary condition.


title: Abnormal Event Detection in Videos using Spatiotemporal Autoencoder

abstract: We present an efficient method for detecting anomalies in videos. Recent applications of convolutional neural networks have shown promises of convolutional layers for object detection and recognition, especially in images. However, convolutional neural networks are supervised and require labels as learning signals. We propose a spatiotemporal architecture for anomaly detection in videos including crowded scenes. Our architecture includes two main components, one for spatial feature representation, and one for learning the temporal evolution of the spatial features. Experimental results on Avenue, Subway and UCSD benchmarks confirm that the detection accuracy of our method is comparable to state-of-the-art methods at a considerable speed of up to 140 fps.


title: Stochastic Optimal Control for Modeling Reaching Movements in the Presence of Obstacles: Theory and Simulation

abstract: In many human-in-the-loop robotic applications such as robot-assisted surgery and remote teleoperation, predicting the intended motion of the human operator may be useful for successful implementation of shared control, guidance virtual fixtures, and predictive control. Developing computational models of human movements is a critical foundation for such motion prediction frameworks. With this motivation, we present a computational framework for modeling reaching movements in the presence of obstacles. We propose a stochastic optimal control framework that consists of probabilistic collision avoidance constraints and a cost function that trades-off between effort and end-state variance in the presence of a signal-dependent noise. First, we present a series of reformulations to convert the original non-linear and non-convex optimal control into a parametric quadratic programming problem. We show that the parameters can be tuned to model various collision avoidance strategies, thereby capturing the quintessential variability associated with human motion. Then, we present a simulation study that demonstrates the complex interaction between avoidance strategies, control cost, and the probability of collision avoidance. The proposed framework can benefit a variety of applications that require teleoperation in cluttered spaces, including robot-assisted surgery. In addition, it can also be viewed as a new optimizer which produces smooth and probabilistically-safe trajectories under signal dependent noise.


title: Gravitational radiation driven capture in unequal mass black hole encounters

abstract: The gravitational radiation driven capture (GR capture) between unequal mass black holes without spins has been investigated with numerical relativistic simulations. We adopt the parabolic approximation which assumes that the gravitational wave radiation from a weakly hyperbolic orbit is the same as that from the parabolic orbit having the same pericenter distance. Using the radiated energies from the parabolic orbit simulations, we have obtained the critical impact parameter ($b_{\rm crit}$) for the GR capture for weakly hyperbolic orbit as a function of initial energy. The most energetic encounters occur around the boundary between the direct merging and the fly-by orbits, and can emit several percent of initial total ADM energy at the peak. When the total mass is fixed, energy and angular momentum radiated in the case of unequal mass black holes are smaller than those of equal mass black holes having the same initial orbital angular momentum for the fly-by orbits. We have compared our results with two different Post-Newtonian (PN) approximations, the exact parabolic orbit (EPO) and PN corrected orbit (PNCO). We find that the agreement between the EPO and the numerical relativity breaks down for very close encounters ($\it{e.g.}$, $b_{\rm crit} \lesssim 100$ M), and it becomes worse for higher mass ratios. For instance, the critical impact parameters can differ by more than $50\%$ from those obtained in EPO if the relative velocity at infinity $v_{\infty}$ is larger than 0.1 for the mass ratio of $m_{1}/m_{2}=16$. The PNCO gives more consistent results than EPO, but it also underestimates the critical impact parameter for the GR capture at $b_{\rm crit} \lesssim40$ M.


title: Tight Relaxations for Polynomial Optimization and Lagrange Multiplier Expressions

abstract: This paper proposes tight semidefinite relaxations for polynomial optimization. The optimality conditions are investigated. We show that generally Lagrange multipliers can be expressed as polynomial functions in decision variables over the set of critical points. The polynomial expressions can be determined by linear equations. Based on these expressions, new Lasserre type semidefinite relaxations are constructed for solving polynomial optimization. We show that the hierarchy of new relaxations has finite convergence, or equivalently, the new relaxations are tight for a finite relaxation order.


title: Asymptotically Locally Optimal Weight Vector Design for a Tighter Correlation Lower Bound of Quasi-Complementary Sequence Sets

abstract: A quasi-complementary sequence set (QCSS) refers to a set of two-dimensional matrices with low non-trivial aperiodic auto- and cross- correlation sums. For multicarrier code-division multiple-access applications, the availability of large QCSSs with low correlation sums is desirable. The generalized Levenshtein bound (GLB) is a lower bound on the maximum aperiodic correlation sum of QCSSs. The bounding expression of GLB is a fractional quadratic function of a weight vector $\mathbf{w}$ and is expressed in terms of three additional parameters associated with QCSS: the set size $K$, the number of channels $M$, and the sequence length $N$. It is known that a tighter GLB (compared to the Welch bound) is possible only if the condition $M\geq2$ and $K\geq \overline{K}+1$, where $\overline{K}$ is a certain function of $M$ and $N$, is satisfied. A challenging research problem is to determine if there exists a weight vector which gives rise to a tighter GLB for \textit{all} (not just \textit{some}) $K\geq \overline{K}+1$ and $M\geq2$, especially for large $N$, i.e., the condition is {asymptotically} both necessary and sufficient. To achieve this, we \textit{analytically} optimize the GLB which is (in general) non-convex as the numerator term is an indefinite quadratic function of the weight vector. Our key idea is to apply the frequency domain decomposition of the circulant matrix (in the numerator term) to convert the non-convex problem into a convex one. Following this optimization approach, we derive a new weight vector meeting the aforementioned objective and prove that it is a local minimizer of the GLB under certain conditions.


title: The hidden-charm multiquark states

abstract: Since 2003 many charmonium-like states were observed experimentally. Especially those charged charmonium-like $Z_c$ states and bottomonium-like $Z_b$ states cannot be accommodated within the naive quark model, which are good candidates of either the hidden-charm tetraquark states or molecules composed of a pair of charmed mesons. In 2015, the LHCb Collaboration discovered two hidden-charm pentaquark states, which are also beyond the quark model. In this talk, we review the current experimental progress and investigate various theoretical interpretations of these candidates of the multiquark states. We list the puzzles and theoretical challenges of these models when confronted with the experimental data. We also discuss possible future measurements which may distinguish the theoretical schemes on the underlying structures of the hidden-charm multiquark states.


title: Topological Spin Liquid with Symmetry-Protected Edge States

abstract: Topological spin liquids are robust quantum states of matter with long-range entanglement and possess many exotic properties such as the fractional statistics of the elementary excitations. Yet these states, short of local parameters like all topological states, are elusive for conventional experimental probes. In this work, we combine theoretical analysis and quantum Monte Carlo numerics on a frustrated spin model which hosts a $\mathbb Z_2$ topological spin liquid ground state, and demonstrate that the presence of symmetry-protected gapless edge modes is a characteristic feature of the state, originating from the nontrivial symmetry fractionalization of the elementary excitations. Experimental observation of these modes on the edge would directly indicate the existence of the topological spin liquids in the bulk, analogous to the fact that the observation of Dirac edge states confirmed the existence of topological insulators.


title: Universality and stability of the edge states of chiral-symmetric topological semimetals and surface states of the Luttinger semimetal

abstract: We theoretically demonstrate that the chiral structure of the nodes of nodal semimetals is responsible for the existence and universal local properties of the edge states in the vicinity of the nodes. We perform a general analysis of the edge states for an isolated node of a 2D semimetal, protected by {\em chiral symmetry} and characterized by the topological winding number $N$. We derive the asymptotic chiral-symmetric boundary conditions and find that there are $N+1$ universal classes of them. The class determines the numbers of flat-band edge states on either side off the node in the 1D spectrum and the winding number $N$ gives the {\em total} number of edge states. We then show that the edge states of chiral nodal semimetals are {\em robust}: they persist in a finite-size {\em stability region} of parameters of chiral-asymmetric terms. This significantly extends the notion of 2D and 3D topological nodal semimetals. We demonstrate that the Luttinger model with a quadratic node for $j=\frac32$ electrons is a 3D topological semimetal in this new sense and predict that $\alpha$-Sn, HgTe, possibly Pr$_2$Ir$_2$O$_7$, and many other semimetals described by it are topological and exhibit surface states.


title: Basics of thermal field theory -- a tutorial on perturbative computations

abstract: These lecture notes, suitable for a two-semester introductory course or self-study, offer an elementary and self-contained exposition of the basic tools and concepts that are encountered in practical computations in perturbative thermal field theory. Selected applications to heavy ion collision physics and cosmology are outlined in the last chapter.


title: The second flight of the SUNRISE balloon-borne solar observatory: overview of instrument updates, the flight, the data and first results

abstract: The SUNRISE balloon-borne solar observatory, consisting of a 1~m aperture telescope that provided a stabilized image to a UV filter imager and an imaging vector polarimeter, carried out its second science flight in June 2013. It provided observations of parts of active regions at high spatial resolution, including the first high-resolution images in the Mg~{\sc ii}~k line. The obtained data are of very high quality, with the best UV images reaching the diffraction limit of the telescope at 3000~\AA\ after Multi-Frame Blind Deconvolution reconstruction accounting for phase-diversity information. Here a brief update is given of the instruments and the data reduction techniques, which includes an inversion of the polarimetric data. Mainly those aspects that evolved compared with the first flight are described. A tabular overview of the observations is given. In addition, an example time series of a part of the emerging active region NOAA AR~11768 observed relatively close to disk centre is described and discussed in some detail. The observations cover the pores in the trailing polarity of the active region, as well as the polarity inversion line where flux emergence was ongoing and a small flare-like brightening occurred in the course of the time series. The pores are found to contain magnetic field strengths ranging up to 2500~G and, while large pores are clearly darker and cooler than the quiet Sun in all layers of the photosphere, the temperature and brightness of small pores approach or even exceed those of the quiet Sun in the upper photosphere.


title: Self-Trapping of G-Mode Oscillations in Relativistic Thin Disks, Revisited

abstract: We examine by a perturbation method how the self-trapping of g-mode oscillations in geometrically thin relativistic disks is affected by uniform vertical magnetic fields. Disks which we consider are isothermal in the vertical direction, but are truncated at a certain height by presence of hot coronae. We find that the characteristics of self-trapping of axisymmetric g-mode oscillations in non-magnetized disks is kept unchanged in magnetized disks at least till a strength of the fields, depending on vertical thickness of disks. These magnetic fields become stronger as the disk becomes thinner. This result suggests that trapped g-mode oscillations still remain as one of possible candidates of quasi-periodic oscillations observed in black-hole and neutron-star X-ray binaries in the cases where vertical magnetic fields in disks are weak.


title: A time dependent approach to model X-ray and $\gamma$--ray light curves of Mrk 421 observed during the flare in February 2010

abstract: Strong X-ray and $\gamma$--ray flares have been detected in February 2010 from the high synchrotron peaked blazar Mrk 421 (z=0.031). With the motivation of understanding the physics involved in this flaring activity, we study the variability of the source in X-ray and $\gamma$--ray energy bands during the period February 10-23, 2010 (MJD 55237-55250). We use near simultaneous X-ray data collected by \emph{MAXI}, \emph{Swift}-XRT and $\gamma$--ray data collected by \emph{Fermi}-LAT and \emph{TACTIC} along with the optical V-band observations by \emph{SPOL} at Steward Observatory. We observe that the variation in the one day averaged flux from the source during the flare is characterized by fast rise and slow decay. Besides, the TeV $\gamma$--ray flux shows a strong correlation with the X-ray flux, suggesting the former to be an outcome of synchrotron self Compton emission process. To model the observed X-ray and $\gamma$--ray light curves, we numerically solve the kinetic equation describing the evolution of particle distribution in the emission region. The injection of particle distribution into the emission region, from the putative acceleration region, is assumed to be a time dependent power law. The synchrotron and synchrotron self Compton emission from the evolving particle distribution in the emission region are used to reproduce the X-ray and $\gamma$--ray flares successfully. Our study suggests that the flaring activity of Mrk 421 can be an outcome of an efficient acceleration process associated with the increase in underlying non-thermal particle distribution.


title: Bayesian Semiparametric Estimation of Cancer-specific Age-at-onset Penetrance with Application to Li-Fraumeni Syndrome

abstract: Penetrance, which plays a key role in genetic research, is defined as the proportion of individuals with the genetic variants (i.e., {genotype}) that cause a particular trait and who have clinical symptoms of the trait (i.e., {phenotype}). We propose a Bayesian semiparametric approach to estimate the cancer-specific age-at-onset penetrance in the presence of the competing risk of multiple cancers. We employ a Bayesian semiparametric competing risk model to model the duration until individuals in a high-risk group develop different cancers, and accommodate family data using family-wise likelihoods. We tackle the ascertainment bias arising when family data are collected through probands in a high-risk population in which disease cases are more likely to be observed. We apply the proposed method to a cohort of 186 families with Li-Fraumeni syndrome identified through probands with sarcoma treated at MD Anderson Cancer Center from 1944 to 1982.


title: Emergence and scaling of spin turbulence in quenched antiferromagnetic spinor Bose-Einstein condensates

abstract: We investigate the phase transition dynamics of a quasi-2D antiferromagnetic spin-1 Bose-Einstein condensate from the easy-axis polar phase to the easy-plane polar phase, which is initiated by suddenly changing the sign of the quadratic Zeeman energy $q$. We observe the emergence and decay of spin turbulence and the formation of half-quantum vortices (HQVs) in the quenched condensate. The characteristic time and length scales of the turbulence generation dynamics are proportional to $|q|^{-1/2}$ as inherited from the dynamic instability of the initial state. In the evolution of the spin turbulence, spin wave excitations develop from large to small length scales, suggesting a direct energy cascade, and the spin population for the axial polar domains exhibit a nonexponential decay. The final equilibrated condensate contains HQVs, and the number is found to increase and saturate with increasing $|q|$. Our results demonstrate the time-space scaling properties of the phase transition dynamics near the critical point and the peculiarities of the spin turbulence state of the antiferromagnetic spinor condensate.


title: The cumulative Kolmogorov filter for model-free screening in ultrahigh dimensional data

abstract: We propose a cumulative Kolmogorov filter to improve the fused Kolmogorov filter proposed by Zou (2015) via cumulative slicing. We establish an improved asymptotic result under relaxed assumptions and numerically demonstrate its enhanced finite sample performance.


title: Calabi-Yau geometry and electrons on 2d lattices

abstract: The B-model approach of topological string theory leads to difference equations by quantizing algebraic mirror curves. It is known that these quantum mechanical systems are solved by the refined topological strings. Recently, it was pointed out that the quantum eigenvalue problem for a particular Calabi--Yau manifold, known as local $\mathbb{F}_0$, is closely related to the Hofstadter problem for electrons on a two-dimensional square lattice. In this paper, we generalize this idea to a more complicated Calabi--Yau manifold. We find that the local $\mathcal{B}_3$ geometry, which is a three-point blow-up of local $\mathbb{P}^2$, is associated with electrons on a triangular lattice. This correspondence allows us to use known results in condensed matter physics to investigate the quantum geometry of the toric Calabi--Yau manifold.


title: Modulations in martensitic Heusler alloys originate from nanotwin ordering

abstract: Heusler alloys exhibiting magnetic and martensitic transitions enable applications like magnetocaloric refrigeration and actuation based on the magnetic shape memory effect. Their outstanding functional properties depend on low hysteresis losses and low actuation fields. These are only achieved if the atomic positions deviate from a tetragonal lattice by periodic displacements. The origin of the so-called modulated structures is the subject of much controversy: They are either explained by phonon softening or adaptive nanotwinning. Here we used large-scale density functional theory calculations on the Ni2MnGa prototype system to demonstrate interaction energy between twin boundaries. Minimizing the interaction energy resulted in the experimentally observed ordered modulations at the atomic scale, it explained that a/b twin boundaries are stacking faults at the mesoscale, and contributed substantially to the macroscopic hysteresis losses. Furthermore, we found that phonon softening paves the transformation path towards the nanotwinned martensite state. This unified both opposing concepts to explain modulated martensite.