title: Cube is a strict local maximizer for the illumination number

abstract: It was conjectured by Levi, Hadwiger, Gohberg and Markus that the boundary of any convex body in ${\mathbb R}^n$ can be illuminated by at most $2^n$ light sources, and, moreover, $2^n-1$ light sources suffice unless the body is a parallelotope. We show that if a convex body is close to the cube in the Banach-Mazur metric, and it is not a parallelotope, then indeed $2^n-1$ light sources suffice to illuminate its boundary. Equivalently, any convex body sufficiently close to the cube, but not isometric to it, can be covered by $2^n-1$ smaller homothetic copies of itself.


title: Invisible Tricorns in Real Slices of Rational Maps

abstract: One of the conspicuous features of real slices of bicritical rational maps is the existence of Tricorn-type hyperbolic components. Such a hyperbolic component is called invisible if the non-bifurcating sub-arcs on its boundary do not intersect the closure of any other hyperbolic component. Numerical evidence suggests an abundance of invisible Tricorn-type components in real slices of bicritical rational maps. In this paper, we study two different families of real bicritical maps and characterize invisible Tricorn-type components in terms of suitable topological properties in the dynamical planes of the representative maps. We use this criterion to prove the existence of infinitely many invisible Tricorn-type components in the corresponding parameter spaces. Although we write the proofs for two specific families, our methods apply to generic families of real bicritical maps.


title: Directional Design of Materials Based on the Multi-Objective Optimization: A Case Study of Two-Dimensional Thermoelectric SnSe

abstract: Directional design of functional materials with multi-objective constraints is a big challenge, whose performance and stability are determined by different physics factors entangled with each other complicatedly. In this work, we apply the multi-objective optimization based on the Pareto Efficiency and Particle-Swarm Optimization methods to design new functional materials directionally. As a demonstration, we achieve the thermoelectric design of 2D SnSe materials through the methods. We identify several novel metastable 2D SnSe structures with simultaneously lower free energy and better thermoelectric performance over the experimentally-reported monolayer structures. We hope our results about the multi-objective Pareto Optimization method can make a step towards the integrative design of multi-objective and multi-functional materials in the future.


title: Improving Shadow Suppression for Illumination Robust Face Recognition

abstract: 2D face analysis techniques, such as face landmarking, face recognition and face verification, are reasonably dependent on illumination conditions which are usually uncontrolled and unpredictable in the real world. An illumination robust preprocessing method thus remains a significant challenge in reliable face analysis. In this paper we propose a novel approach for improving lighting normalization through building the underlying reflectance model which characterizes interactions between skin surface, lighting source and camera sensor, and elaborates the formation of face color appearance. Specifically, the proposed illumination processing pipeline enables the generation of Chromaticity Intrinsic Image (CII) in a log chromaticity space which is robust to illumination variations. Moreover, as an advantage over most prevailing methods, a photo-realistic color face image is subsequently reconstructed which eliminates a wide variety of shadows whilst retaining the color information and identity details. Experimental results under different scenarios and using various face databases show the effectiveness of the proposed approach to deal with lighting variations, including both soft and hard shadows, in face recognition.


title: The Heating of Solar Coronal Loops by Alfven Wave Turbulence

abstract: In this paper we further develop a model for the heating of coronal loops by Alfven wave turbulence (AWT). The Alfven waves are assumed to be launched from a collection of kilogauss flux tubes in the photosphere at the two ends of the loop. Using a three-dimensional magneto-hydrodynamic (MHD) model for an active-region loop, we investigate how the waves from neighboring flux tubes interact in the chromosphere and corona. For a particular combination of model parameters we find that AWT can produce enough heat to maintain a peak temperature of about 2.5 MK, somewhat lower than the temperatures of 3-4 MK observed in the cores of active regions. The heating rates vary strongly in space and time, but the simulated heating events have durations less than 1 minute and are unlikely to reproduce the observed broad Differential Emission Measure distributions of active regions. The simulated spectral line non-thermal widths are predicted to be about 27 km/s, which is high compared to the observed values. Therefore, the present AWT model does not satisfy the observational constraints. An alternative "magnetic braiding" model is considered in which the coronal field lines are subject to slow random footpoint motions, but we find that such long period motions produce much less heating than the shorter period waves launched within the flux tubes. We discuss several possibilities for resolving the problem of producing sufficiently hot loops in active regions.


title: Sensitivity and discovery potential of the proposed nEXO experiment to neutrinoless double beta decay

abstract: The next-generation Enriched Xenon Observatory (nEXO) is a proposed experiment to search for neutrinoless double beta ($0\nu\beta\beta$) decay in $^{136}$Xe with a target half-life sensitivity of approximately $10^{28}$ years using $5\times10^3$ kg of isotopically enriched liquid-xenon in a time projection chamber. This improvement of two orders of magnitude in sensitivity over current limits is obtained by a significant increase of the $^{136}$Xe mass, the monolithic and homogeneous configuration of the active medium, and the multi-parameter measurements of the interactions enabled by the time projection chamber. The detector concept and anticipated performance are presented based upon demonstrated realizable background rates.


title: Limited Resource Semi-Quantum Key Distribution

abstract: A semi-quantum key distribution (SQKD) protocol allows a quantum user and a limited "classical" user to establish a shared secret key secure against an all-powerful adversary. In this work, we present a new SQKD protocol where the quantum user is also limited in her measurement capabilities. We describe the protocol, prove its security, and show its noise tolerance is as high as "fully quantum" QKD protocols.


title: The Anomalous Yield Behavior of Fused Silica Glass

abstract: We develop a critical-state model of fused silica plasticity on the basis of data mined from molecular dynamics (MD) calculations. The MD data is suggestive of an irreversible densification transition in volumetric compression resulting in permanent, or plastic, densification upon unloading. The MD data also reveals an evolution towards a critical state of constant volume under pressure-shear deformation. The trend towards constant volume is from above, when the glass is overconsolidated, or from below, when it is underconsolidated. We show that these characteristic behaviors are well-captured by a critical state model of plasticity, where the densification law for glass takes the place of the classical consolidation law of granular media and the locus of constant volume states denotes the critical-state line. A salient feature of the critical-state line of fused silica, as identified from the MD data, that renders its yield behavior anomalous is that it is strongly non-convex, owing to the existence of two well-differentiated phases at low and high pressures. We argue that this strong non-convexity of yield explains the patterning that is observed in molecular dynamics calculations of amorphous solids deforming in shear. We employ an explicit and exact rank-2 envelope construction to upscale the microscopic critical-state model to the macroscale. Remarkably, owing to the equilibrium constraint the resulting effective macroscopic behavior is still characterized by a non-convex critical-state line. Despite this lack of convexity, the effective macroscopic model is stable against microstructure formation and defines well-posed boundary-value problems.


title: Metric transforms yielding Gromov hyperbolic spaces

abstract: A real valued function $\varphi$ of one variable is called a metric transform if for every metric space $(X,d)$ the composition $d_\varphi = \varphi\circ d$ is also a metric on $X$. We give a complete characterization of the class of approximately nondecreasing, unbounded metric transforms $\varphi$ such that the transformed Euclidean half line $([0,\infty),|\cdot|_\varphi)$ is Gromov hyperbolic. As a consequence, we obtain metric transform rigidity for roughly geodesic Gromov hyperbolic spaces, that is, if $(X,d)$ is any metric space containing a rough geodesic ray and $\varphi$ is an approximately nondecreasing, unbounded metric transform such that the transformed space $(X,d_\varphi)$ is Gromov hyperbolic and roughly geodesic then $\varphi$ is an approximate dilation and the original space $(X,d)$ is Gromov hyperbolic and roughly geodesic.


title: A complexity problem for Borel graphs

abstract: We show that there is no simple (e.g. finite or countable) basis for Borel graphs with infinite Borel chromatic number. In fact, it is proved that the closed subgraphs of the shift graph on $[\mathbb{N}]^{<\mathbb{N}}$ with finite (or, equivalently, $\leq 3$) Borel chromatic number form a $\mathbf{\Sigma}^1_2$-complete set. This answers a question of Kechris and Marks and strengthens several earlier results.


title: DSCOVR: Randomized Primal-Dual Block Coordinate Algorithms for Asynchronous Distributed Optimization

abstract: Machine learning with big data often involves large optimization models. For distributed optimization over a cluster of machines, frequent communication and synchronization of all model parameters (optimization variables) can be very costly. A promising solution is to use parameter servers to store different subsets of the model parameters, and update them asynchronously at different machines using local datasets. In this paper, we focus on distributed optimization of large linear models with convex loss functions, and propose a family of randomized primal-dual block coordinate algorithms that are especially suitable for asynchronous distributed implementation with parameter servers. In particular, we work with the saddle-point formulation of such problems which allows simultaneous data and model partitioning, and exploit its structure by doubly stochastic coordinate optimization with variance reduction (DSCOVR). Compared with other first-order distributed algorithms, we show that DSCOVR may require less amount of overall computation and communication, and less or no synchronization. We discuss the implementation details of the DSCOVR algorithms, and present numerical experiments on an industrial distributed computing system.


title: On Graphene in the Interstellar Medium

abstract: The possible detection of C_{24}, a planar graphene, recently reported in several planetary nebulae by Garciaa-Hernandez et al. (2011, 2012) inspires us to explore whether and how much graphene could exist in the interstellar medium (ISM) and how it would reveal its presence through its ultraviolet (UV) extinction and infrared (IR) emission. In principle, interstellar graphene could arise from the photochemical processing of polycyclic aromatic hydrocarbon (PAH) molecules which are abundant in the ISM through a complete loss of their hydrogen atoms and/or from graphite which is thought to be a major dust species in the ISM through fragmentation caused by grain-grain collisional shattering. Both quantum-chemical computations and laboratory experiments have shown that the exciton-dominated electronic transitions in graphene cause a strong absorption band near 2755 Angstrom. We calculate the UV absorption of graphene and place an upper limit of ~5 ppm of C/H (i.e., ~1.9\% of the total interstellar C) on the interstellar graphene abundance. We also model the stochastic heating of graphene C_{24} excited by single starlight photons of the interstellar radiation field in the ISM and calculate its IR emission spectra. We also derive the abundance of graphene in the ISM to be <5 ppm of C/H by comparing the model emission spectra with that observed in the ISM.


title: Exact Exchange: a pathway for a Density Functional Theory of the Integer Quantum Hall Effect

abstract: It is shown here that the Exact Exchange (EE) formalism provides a natural and rigorous approach for a Density Functional Theory (DFT) of the Integer Quantum Hall Effect (IQHE). Application of a novel EE method to a quasi two-dimensional electron gas (q2DEG) subjected to a perpendicular magnetic field leads to the following main findings. \textit{i)} the microscopic exchange energy functional of the IQHE has been obtained, whose main feature being that it minimizes with a discontinuous derivative at every integer filling factor $\nu$; \textit{ii)} an analytical solution is found for the magnetic-field dependent EE potential, in the one-subband regime; \textit{iii)} as a consequence of \textit{i)}, the EE potential display sharp discontinuities at every integer $\nu$; and \textit{iv)} the widely used Local Spin Density Approximation (LSDA) is strongly violated for filling factors close to integer values.


title: Relieving the frustration through Mn$^{3+}$ substitution in Holmium Gallium Garnet

abstract: We present a study on the impact of Mn$^{3+}$ substitution in the geometrically frustrated Ising garnet Ho$_3$Ga$_5$O$_{12}$ using bulk magnetic measurements and low temperature powder neutron diffraction. We find that the transition temperature, $T_N$ = 5.8 K, for Ho$_3$MnGa$_4$O$_{12}$ is raised by almost 20 when compared to Ho$_3$Ga$_5$O$_{12}$. Powder neutron diffraction on Ho$_3$Mn$_x$Ga$_{5-x}$O$_{12}$ ($x$ = 0.5, 1) below $T_N$ shows the formation of a long range ordered ordered state with $\mathbf{k}$ = (0,0,0). Ho$^{3+}$ spins are aligned antiferromagnetically along the six crystallographic axes with no resultant moment while the Mn$^{3+}$ spins are oriented along the body diagonals, such that there is a net moment along [111]. The magnetic structure can be visualised as ten-membered rings of corner-sharing triangles of Ho$^{3+}$ spins with the Mn$^{3+}$ spins ferromagnetically coupled to each individual Ho$^{3+}$ spin in the triangle. Substitution of Mn$^{3+}$ completely relieves the magnetic frustration with $f = \theta_{CW}/T_N \approx 1.1$ for Ho$_3$MnGa$_4$O$_{12}$.


title: Finding paths in tree graphs with a quantum walk

abstract: In this paper, we analyze the potential for new types of searches using the formalism of scattering random walks on Quantum Computers. Given a particular type of graph consisting of nodes and connections, a "Tree Maze", we would like to find a selected final node as quickly as possible, faster than any classical search algorithm. We show that this can be done using a quantum random walk, both exactly through numerical calculations as well as analytically using eigenvectors and eigenvalues of the quantum system.


title: Solving the Insoluble: A New Rule for Quantization

abstract: The rules of canonical quantization normally offer good results, but sometimes they fail, e.g., leading to quantum triviality ($=$ free) for certain examples that are classically nontrivial ($\ne$ free). A new procedure, called Enhanced Quantization, relates classical models with their quantum partners differently and leads to satisfactory results for all systems. This paper features enhanced quantization procedures and provides highlights of two examples, a rotationally symmetric model and an ultralocal scalar model, for which canonical quantization fails while enhanced quantization succeeds.


title: A deep generative model for single-cell RNA sequencing with application to detecting differentially expressed genes

abstract: We propose a probabilistic model for interpreting gene expression levels that are observed through single-cell RNA sequencing. In the model, each cell has a low-dimensional latent representation. Additional latent variables account for technical effects that may erroneously set some observations of gene expression levels to zero. Conditional distributions are specified by neural networks, giving the proposed model enough flexibility to fit the data well. We use variational inference and stochastic optimization to approximate the posterior distribution. The inference procedure scales to over one million cells, whereas competing algorithms do not. Even for smaller datasets, for several tasks, the proposed procedure outperforms state-of-the-art methods like ZIFA and ZINB-WaVE. We also extend our framework to take into account batch effects and other confounding factors and propose a natural Bayesian hypothesis framework for differential expression that outperforms tradition DESeq2.


title: On Bi-free Multiplicative Convolution

abstract: In this paper, we study the partial bi-free $S$-transform of a pair $(a,b)$ of random variables, and the $S$-transform of the $2\times 2$ matrix-valued random variable $\left(\begin{matrix}a&0\\0&b\end{matrix}\right)$ associated with $(a,b)$ when restricted to upper triangular $2\times 2$ matrices. We first derive an explicit expression of bi-free multiplicative convolution (of probability measures on the bi-unit-sphere $\mathbb{T}^2$ of $\mathbb{C}^2$, or on $\mathbb{R}^2_+$ in $\mathbb{C}^2$) from a subordination equation for bi-free multiplicative convolution. We then show that, when $(a_1, b_1)$ and $(a_2,b_2)$ are bi-free, the $S$-transforms of $X_1=\left(\begin{matrix}a_1&0\\0&b_1\end{matrix}\right)$, $X_2=\left(\begin{matrix}a_2&0\\0&b_2\end{matrix}\right)$ satisfy Dykema's twisted multiplicative equation for free operator-valued random variables if and only if at least one of the two partial bi-free $S$-transforms of the pairs of random variables is the constant function 1 in a neighborhood of $(0,0)$. This is the case if and only if one of the two pairs, say $(a_1,b_1)$, has factoring two-band moments (that is, $\varphi(a_1^mb_1^n)=\varphi(a_1^m)\varphi(b_1^n)$, for all $m,n=1, 2, \cdots$). We thus find tons of bi-free pairs of random variables to which the $S$-transforms of the corresponding matrix-value random variables do not satisfy Dykema's twisted multiplicative formula. Finally, if both $(a_1,b_1)$ and $(a_2,b_2)$ have factoring two-band moments, we prove that the $\Psi$-transforms of $X_1$, $X_2$, and $X_1X_2$ satisfy a subordination equation.


title: The magnetic future of the Sun

abstract: We analyze space- and ground-based data for the old ($7.0\pm0.3$~Gyr) solar analogs 16 Cyg A and B. The stars were observed with the Cosmic Origins UV Spectrographs on the Hubble Space Telescope (HST) on 23 October 2015 and 3 February 2016 respectively, and with the Chandra X-ray Observatory on 7 February 2016. Time-series data in \ion{Ca}{2} data are used to place the UV data in context. The UV spectra of 18 Sco (3.7$\pm0.5$ Gyr), the Sun (4.6$\pm0.04$ Gyr) and $\alpha$ Cen A ($5.4_{-0.2}^{+1.2}$ Gyr), appear remarkably similar, pointing to a convergence of magnetic heating rates for G2 main-sequence stars older than $\approx 2-4$ Gyr. But the B component's X-ray (0.3-2.5 keV) flux lies 20$\times$ below a well-known minimum level reported by Schmitt. As reported for $\alpha$~Cen~A, the coronal temperature probably lies below that detectable in soft X-rays. No solar UV flux spectra of comparable resolution to stellar data exist, but they are badly needed for comparison with stellar data. Center-to-limb (C-L) variations are re-evaluated for lines such as \ion{Ca}{2} through to X-rays, with important consequences for observing activity cycles in such features. We also call into question work that has mixed solar intensity-intensity statistics with flux-flux relations of stars.


title: Variability Of The Magnetic Field Power Spectrum In The Solar Wind At Electron Scales

abstract: At the electron scales the power spectrum of solar-wind magnetic fluctuations can be highly variable and the dissipation mechanisms of the magnetic energy into the various particle species is under debate. In this paper we investigate data from the Cluster mission's STAFF Search Coil magnetometer when the level of turbulence is sufficiently high that the morphology of the power spectrum at electron scales can be investigated. The Cluster spacecraft sample a disturbed interval of plasma where two streams of solar wind interact. Meanwhile, several discontinuities (coherent structures) are seen in the large scale magnetic field, while at small scales several intermittent bursts of wave activity (whistler waves) are present. Several different morphologies of the power spectrum can be identified: (1) two power laws separated by a break (2) an exponential cutoff near the Taylor shifted electron scales and (3) strong spectral knees at the Taylor shifted electron scales. These different morphologies are investigated by using wavelet coherence, showing that in this interval a clear break and strong spectral knees are features which are associated with sporadic quasi parallel propagating whistler waves, even for short times. On the other hand, when no signatures of whistler waves at 0.1 - 0.2fce are present, a clear break is difficult to find and the spectrum is often more characteristic of a power law with an exponential cutoff.


title: Burn-In Demonstrations for Multi-Modal Imitation Learning

abstract: Recent work on imitation learning has generated policies that reproduce expert behavior from multi-modal data. However, past approaches have focused only on recreating a small number of distinct, expert maneuvers, or have relied on supervised learning techniques that produce unstable policies. This work extends InfoGAIL, an algorithm for multi-modal imitation learning, to reproduce behavior over an extended period of time. Our approach involves reformulating the typical imitation learning setting to include "burn-in demonstrations" upon which policies are conditioned at test time. We demonstrate that our approach outperforms standard InfoGAIL in maximizing the mutual information between predicted and unseen style labels in road scene simulations, and we show that our method leads to policies that imitate expert autonomous driving systems over long time horizons.


title: A simple data discretizer

abstract: Data discretization is an important step in the process of machine learning, since it is easier for classifiers to deal with discrete attributes rather than continuous attributes. Over the years, several methods of performing discretization such as Boolean Reasoning, Equal Frequency Binning, Entropy have been proposed, explored, and implemented. In this article, a simple supervised discretization approach is introduced. The prime goal of MIL is to maximize classification accuracy of classifier, minimizing loss of information while discretization of continuous attributes. The performance of the suggested approach is compared with the supervised discretization algorithm Minimum Information Loss (MIL), using the state-of-the-art rule inductive algorithms- J48 (Java implementation of C4.5 classifier). The presented approach is, indeed, the modified version of MIL. The empirical results show that the modified approach performs better in several cases in comparison to the original MIL algorithm and Minimum Description Length Principle (MDLP) .


title: Dropout as a Low-Rank Regularizer for Matrix Factorization

abstract: Regularization for matrix factorization (MF) and approximation problems has been carried out in many different ways. Due to its popularity in deep learning, dropout has been applied also for this class of problems. Despite its solid empirical performance, the theoretical properties of dropout as a regularizer remain quite elusive for this class of problems. In this paper, we present a theoretical analysis of dropout for MF, where Bernoulli random variables are used to drop columns of the factors. We demonstrate the equivalence between dropout and a fully deterministic model for MF in which the factors are regularized by the sum of the product of squared Euclidean norms of the columns. Additionally, we inspect the case of a variable sized factorization and we prove that dropout achieves the global minimum of a convex approximation problem with (squared) nuclear norm regularization. As a result, we conclude that dropout can be used as a low-rank regularizer with data dependent singular-value thresholding.


title: Configuration spaces of products

abstract: We show that the configuration spaces of a product of parallelizable manifolds may be recovered from those of the factors as the Boardman-Vogt tensor product of right modules over the operads of little cubes of the appropriate dimension. We also discuss an analogue of this result for manifolds that are not necessarily parallelizable, which involves a new operad of skew little cubes.


title: Learning Phrase Embeddings from Paraphrases with GRUs

abstract: Learning phrase representations has been widely explored in many Natural Language Processing (NLP) tasks (e.g., Sentiment Analysis, Machine Translation) and has shown promising improvements. Previous studies either learn non-compositional phrase representations with general word embedding learning techniques or learn compositional phrase representations based on syntactic structures, which either require huge amounts of human annotations or cannot be easily generalized to all phrases. In this work, we propose to take advantage of large-scaled paraphrase database and present a pair-wise gated recurrent units (pairwise-GRU) framework to generate compositional phrase representations. Our framework can be re-used to generate representations for any phrases. Experimental results show that our framework achieves state-of-the-art results on several phrase similarity tasks.


title: Differentially Private Query Learning: from Data Publishing to Model Publishing

abstract: With the development of Big Data and cloud data sharing, privacy preserving data publishing becomes one of the most important topics in the past decade. As one of the most influential privacy definitions, differential privacy provides a rigorous and provable privacy guarantee for data publishing. Differentially private interactive publishing achieves good performance in many applications; however, the curator has to release a large number of queries in a batch or a synthetic dataset in the Big Data era. To provide accurate non-interactive publishing results in the constraint of differential privacy, two challenges need to be tackled: one is how to decrease the correlation between large sets of queries, while the other is how to predict on fresh queries. Neither is easy to solve by the traditional differential privacy mechanism. This paper transfers the data publishing problem to a machine learning problem, in which queries are considered as training samples and a prediction model will be released rather than query results or synthetic datasets. When the model is published, it can be used to answer current submitted queries and predict results for fresh queries from the public. Compared with the traditional method, the proposed prediction model enhances the accuracy of query results for non-interactive publishing. Experimental results show that the proposed solution outperforms traditional differential privacy in terms of Mean Absolute Value on a large group of queries. This also suggests the learning model can successfully retain the utility of published queries while preserving privacy.


title: On the Ontological Modeling of Trees

abstract: Trees -- i.e., the type of data structure known under this name -- are central to many aspects of knowledge organization. We investigate some central design choices concerning the ontological modeling of such trees. In particular, we consider the limits of what is expressible in the Web Ontology Language, and provide a reusable ontology design pattern for trees.


title: The electronic-structure origin of cation disorder in transition-metal oxides

abstract: Cation disorder is an important design criterion for technologically relevant transition-metal (TM) oxides, such as radiation-tolerant ceramics and Li-ion battery electrodes. In this letter, we use a combination of first-principles calculations, normal mode analysis, and band-structure arguments to pinpoint a specific electronic-structure effect that influences the stability of disordered phases. We find that the electronic configuration of a TM ion determines to which extent the structural energy is affected by site distortions. This mechanism explains the stability of disordered phases with large ionic radius differences and provides a concrete guideline for the discovery of novel disordered compositions.


title: Determination of Effective Synaptic Conductances Using Somatic Voltage Clamp

abstract: The interplay between excitatory and inhibitory neurons imparts rich functions of the brain. To understand the underlying synaptic mechanisms, a fundamental approach is to study the dynamics of excitatory and inhibitory conductances of each neuron. The traditional method of determining conductance employs the synaptic current-voltage (I-V) relation obtained via voltage clamp. Using theoretical analysis, electrophysiological experiments, and realistic simulations, here we demonstrate that the traditional method conceptually fails to measure the conductance due to the neglect of a nonlinear interaction between the clamp current and the synaptic current. Consequently, it incurs substantial measurement error, even giving rise to unphysically negative conductance as observed in experiments. To elucidate synaptic impact on neuronal information processing, we introduce the concept of effective conductance and propose a framework to determine it accurately. Our work suggests re-examination of previous studies involving conductance measurement and provides a reliable approach to assess synaptic influence on neuronal computation.


title: The second hyperpolarizability of systems described by the space-fractional Schrodinger equation

abstract: The static second hyperpolarizability is derived from the space-fractional Schr\"{o}dinger equation in the particle-centric view. The Thomas-Reiche-Kuhn sum rule matrix elements and the three-level ansatz determines the maximum second hyperpolarizability for a space-fractional quantum system. The total oscillator strength is shown to decrease as the space-fractional parameter $\alpha$ decreases, which reduces the optical response of a quantum system in the presence of an external field. This damped response is caused by the wavefunction dependent position and momentum commutation relation. Although the maximum response is damped, we show that the one-dimensional quantum harmonic oscillator is no longer a linear system for $\alpha \neq 1$, where the second hyperpolarizability becomes negative before ultimately damping to zero at the lower fractional limit of $\alpha \rightarrow 1/2$.


title: Neural network dose models for knowledge-based planning in pancreatic SBRT

abstract: Stereotactic body radiation therapy (SBRT) for pancreatic cancer requires a skillful approach to deliver ablative doses to the tumor while limiting dose to the highly sensitive duodenum, stomach, and small bowel. Here, we develop knowledge-based artificial neural network dose models (ANN-DMs) to predict dose distributions that would be approved by experienced physicians. Using dose distributions calculated by a commercial treatment planning system (TPS), physician-approved treatment plans were used to train ANN-DMs that could predict physician-approved dose distributions based on a set of geometric parameters (vary from voxel to voxel) and plan parameters (constant across all voxels for a given patient). Differences between TPS and ANN-DM dose distributions were used to evaluate model performance. ANN-DM design, including neural network structure and parameter choices, were evaluated to optimize dose model performance. Mean dose errors were less than 5% at all distances from the PTV, and mean absolute dose errors were on the order of 5%, but no more than 10%. Dose-volume histogram errors demonstrated good model performance above 25 Gy, but much larger errors were seen at lower doses. ANN-DM dose distributions showed excellent overall agreement with TPS dose distributions, and accuracy was substantially improved when each physician's treatment approach was taken into account by training their own dedicated models. In this manner, one could feasibly train ANN-DMs that could predict the dose distribution desired by a given physician for a given treatment site.


title: Unsupervised Real-Time Control through Variational Empowerment

abstract: We introduce a methodology for efficiently computing a lower bound to empowerment, allowing it to be used as an unsupervised cost function for policy learning in real-time control. Empowerment, being the channel capacity between actions and states, maximises the influence of an agent on its near future. It has been shown to be a good model of biological behaviour in the absence of an extrinsic goal. But empowerment is also prohibitively hard to compute, especially in nonlinear continuous spaces. We introduce an efficient, amortised method for learning empowerment-maximising policies. We demonstrate that our algorithm can reliably handle continuous dynamical systems using system dynamics learned from raw data. The resulting policies consistently drive the agents into states where they can use their full potential.


title: Darboux Transformation for the Nonisospectral and Variable-coefficient KdV Equation

abstract: With the nonuniform media taken into account, the nonisospectral and variable-coefficient Korteweg-de Vries equation, which describes various physical situations such as fluid dynamics and plasma, is under investigation in this paper. With appropriate selection of wave functions, the Darboux transformation is constructed, by which the multi-soliton solutions are derived and graphs are presented. The spectral parameters, coefficients and initial phase are discussed analytically and numerically to demonstrate their respective effect on the soliton dynamics, which plays a role in achieving the feasible soliton management with explicit conditions taken into account.


title: Exact and asymptotic enumeration of cyclic permutations according to descent set

abstract: Using a result of Gessel and Reutenauer, we find a simple formula for the number of cyclic permutations with a given descent set, by expressing it in terms of ordinary descent numbers (i.e., those counting all permutations with a given descent set). We then use this formula to show that, for almost all sets $I \subseteq [n-1]$, the fraction of size-$n$ permutations with descent set $I$ which are $n$-cycles is asymptotically $1/n$. As a special case, we recover a result of Stanley for alternating cycles. We also use our formula to count the cycles that do not have two consecutive descents.


title: An adaptive thresholding approach for automatic optic disk segmentation

abstract: Optic disk segmentation is a prerequisite step in automatic retinal screening systems. In this paper, we propose an algorithm for optic disk segmentation based on a local adaptive thresholding method. Location of the optic disk is validated by intensity and average vessel width of retinal images. Then an adaptive thresholding is applied on the temporal and nasal part of the optic disc separately. Adaptive thresholding, makes our algorithm robust to illumination variations and various image acquisition conditions. Moreover, experimental results on the DRIVE and KHATAM databases show promising results compared to the recent literature. In the DRIVE database, the optic disk in all images is correctly located and the mean overlap reached to 43.21%. The optic disk is correctly detected in 98% of the images with the mean overlap of 36.32% in the KHATAM database.


title: Diagonalization of indefinite saddle point forms

abstract: We obtain sufficient conditions that ensure block diagonalization (by a direct rotation) of sign-indefinite symmetric sesquilinear forms as well as the associated operators that are semi-bounded neither from below nor from above. In the semi-bounded case, we refine the obtained results and, as an example, revisit the block Stokes Operator from fluid dynamics.


title: CM-GANs: Cross-modal Generative Adversarial Networks for Common Representation Learning

abstract: It is known that the inconsistent distribution and representation of different modalities, such as image and text, cause the heterogeneity gap that makes it challenging to correlate such heterogeneous data. Generative adversarial networks (GANs) have shown its strong ability of modeling data distribution and learning discriminative representation, existing GANs-based works mainly focus on generative problem to generate new data. We have different goal, aim to correlate heterogeneous data, by utilizing the power of GANs to model cross-modal joint distribution. Thus, we propose Cross-modal GANs to learn discriminative common representation for bridging heterogeneity gap. The main contributions are: (1) Cross-modal GANs architecture is proposed to model joint distribution over data of different modalities. The inter-modality and intra-modality correlation can be explored simultaneously in generative and discriminative models. Both of them beat each other to promote cross-modal correlation learning. (2) Cross-modal convolutional autoencoders with weight-sharing constraint are proposed to form generative model. They can not only exploit cross-modal correlation for learning common representation, but also preserve reconstruction information for capturing semantic consistency within each modality. (3) Cross-modal adversarial mechanism is proposed, which utilizes two kinds of discriminative models to simultaneously conduct intra-modality and inter-modality discrimination. They can mutually boost to make common representation more discriminative by adversarial training process. To the best of our knowledge, our proposed CM-GANs approach is the first to utilize GANs to perform cross-modal common representation learning. Experiments are conducted to verify the performance of our proposed approach on cross-modal retrieval paradigm, compared with 10 methods on 3 cross-modal datasets.


title: Linear Progress with Exponential Decay in Weakly Hyperbolic Groups

abstract: A random walk $w_n$ on a separable, geodesic hyperbolic metric space $X$ converges to the boundary $\partial X$ with probability one when the step distribution supports two independent loxodromics. In particular, the random walk makes positive linear progress. Progress is known to be linear with exponential decay when (1) the step distribution has exponential tail and (2) the action on $X$ is acylindrical. We extend exponential decay to the non-acylindrical case.


title: Solitons and breathers for nonisospectral mKdV equation with Darboux transformation

abstract: Under investigation in this paper is the nonisospectral and variable coefficients modified Kortweg-de Vries (vc-mKdV) equation, which manifests in diverse areas of physics such as fluid dynamics, ion acoustic solitons and plasma mechanics. With the degrees of restriction reduced, a simplified constraint is introduced, under which the vc-mKdV equation is an integrable system and the spectral flow is time-varying. The Darboux transformation for such equation is constructed, which gives rise to the generation of variable kinds of solutions including the double-breather coherent structure, periodical soliton-breather and localized solitons and breathers. In addition, the effect of variable coefficients and initial phases is discussed in terms of the soliton amplitude, polarity, velocity and width, which might provide feasible soliton management with certain conditions taken into account.


title: Characterization of an Ionization Readout Tile for nEXO

abstract: A new design for the anode of a time projection chamber, consisting of a charge-detecting "tile", is investigated for use in large scale liquid xenon detectors. The tile is produced by depositing 60 orthogonal metal charge-collecting strips, 3~mm wide, on a 10~\si{\cm} $\times$ 10~\si{\cm} fused-silica wafer. These charge tiles may be employed by large detectors, such as the proposed tonne-scale nEXO experiment to search for neutrinoless double-beta decay. Modular by design, an array of tiles can cover a sizable area. The width of each strip is small compared to the size of the tile, so a Frisch grid is not required. A grid-less, tiled anode design is beneficial for an experiment such as nEXO, where a wire tensioning support structure and Frisch grid might contribute radioactive backgrounds and would have to be designed to accommodate cycling to cryogenic temperatures. The segmented anode also reduces some degeneracies in signal reconstruction that arise in large-area crossed-wire time projection chambers. A prototype tile was tested in a cell containing liquid xenon. Very good agreement is achieved between the measured ionization spectrum of a $^{207}$Bi source and simulations that include the microphysics of recombination in xenon and a detailed modeling of the electrostatic field of the detector. An energy resolution $\sigma/E$=5.5\% is observed at 570~\si{keV}, comparable to the best intrinsic ionization-only resolution reported in literature for liquid xenon at 936~V/\si{cm}.


title: Subsampling for Ridge Regression via Regularized Volume Sampling

abstract: Given $n$ vectors $\mathbf{x}_i\in \mathbb{R}^d$, we want to fit a linear regression model for noisy labels $y_i\in\mathbb{R}$. The ridge estimator is a classical solution to this problem. However, when labels are expensive, we are forced to select only a small subset of vectors $\mathbf{x}_i$ for which we obtain the labels $y_i$. We propose a new procedure for selecting the subset of vectors, such that the ridge estimator obtained from that subset offers strong statistical guarantees in terms of the mean squared prediction error over the entire dataset of $n$ labeled vectors. The number of labels needed is proportional to the statistical dimension of the problem which is often much smaller than $d$. Our method is an extension of a joint subsampling procedure called volume sampling. A second major contribution is that we speed up volume sampling so that it is essentially as efficient as leverage scores, which is the main i.i.d. subsampling procedure for this task. Finally, we show theoretically and experimentally that volume sampling has a clear advantage over any i.i.d. sampling when labels are expensive.


title: Reconfigurable Antennas in mmWave MIMO Systems

abstract: The key obstacle to achieving the full potential of the millimeter wave (mmWave) band has been the poor propagation characteristics of wireless signals in this band. One approach to overcome this issue is to use antennas that can support higher gains while providing beam adaptability and diversity, i.e., reconfigurable antennas. In this article, we present a new architecture for mmWave multiple-input multiple-output (MIMO) communications that uses a new class of reconfigurable antennas. More specifically, the proposed lens-based antennas can support multiple radiation patterns while using a single radio frequency chain. Moreover, by using a beam selection network, each antenna beam can be steered in the desired direction. Further, using the proposed reconfigurable antenna in a MIMO architecture, we propose a new signal processing algorithm that uses the additional degrees of freedom provided by the antennas to overcome propagation issues at mmWave frequencies. Our simulation results show that the proposed reconfigurable antenna MIMO architecture significantly enhances the performance of mmWave communication systems.


title: Video Classification With CNNs: Using The Codec As A Spatio-Temporal Activity Sensor

abstract: We investigate video classification via a two-stream convolutional neural network (CNN) design that directly ingests information extracted from compressed video bitstreams. Our approach begins with the observation that all modern video codecs divide the input frames into macroblocks (MBs). We demonstrate that selective access to MB motion vector (MV) information within compressed video bitstreams can also provide for selective, motion-adaptive, MB pixel decoding (a.k.a., MB texture decoding). This in turn allows for the derivation of spatio-temporal video activity regions at extremely high speed in comparison to conventional full-frame decoding followed by optical flow estimation. In order to evaluate the accuracy of a video classification framework based on such activity data, we independently train two CNN architectures on MB texture and MV correspondences and then fuse their scores to derive the final classification of each test video. Evaluation on two standard datasets shows that the proposed approach is competitive to the best two-stream video classification approaches found in the literature. At the same time: (i) a CPU-based realization of our MV extraction is over 977 times faster than GPU-based optical flow methods; (ii) selective decoding is up to 12 times faster than full-frame decoding; (iii) our proposed spatial and temporal CNNs perform inference at 5 to 49 times lower cloud computing cost than the fastest methods from the literature.


title: A New Framework for Determination of Excitatory and Inhibitory Conductances Using Somatic Clamp

abstract: The interaction between excitation and inhibition is crucial for brain computation. To understand synaptic mechanisms underlying brain function, it is important to separate excitatory and inhibitory inputs to a target neuron. In the traditional method, after applying somatic current or voltage clamp, the excitatory and inhibitory conductances are determined from the synaptic current-voltage (I-V) relation --- the slope corresponds to the total conductance and the intercept corresponds to the reversal current. Because of the space clamp effect, the measured conductance in general deviates substantially from the local conductance on the dendrite. Therefore, the interpretation of the conductance measured by the traditional method remains to be clarified. In this work, based on the investigation of an idealized ball-and-stick neuron model and a biologically realistic pyramidal neuron model, we first demonstrate both analytically and numerically that the conductance determined by the traditional method has no clear biological interpretation due to the neglect of a nonlinear interaction between the clamp current and the synaptic current across the spatial dendrites. As a consequence, the traditional method can induce an arbitrarily large error of conductance measurement, sometimes even leads to unphysically negative conductance. To circumvent the difficulty of elucidating synaptic impact on neuronal computation using the traditional method, we then propose a framework to determine the effective conductance that reflects directly the functional impact of synaptic inputs on action potential initiation and thereby neuronal information processing. Our framework has been further verified in realistic neuron simulations, thus greatly improves upon the traditional approach by providing a reliable and accurate assessment of the role of synaptic activity in neuronal computation.


title: Deep Learning in a Generalized HJM-type Framework Through Arbitrage-Free Regularization

abstract: We introduce a regularization approach to arbitrage-free factor-model selection. The considered model selection problem seeks to learn the closest arbitrage-free HJM-type model to any prespecified factor-model. An asymptotic solution to this, a priori computationally intractable, problem is represented as the limit of a 1-parameter family of optimizers to computationally tractable model selection tasks. Each of these simplified model-selection tasks seeks to learn the most similar model, to the prescribed factor-model, subject to a penalty detecting when the reference measure is a local martingale-measure for the entire underlying financial market. A simple expression for the penalty terms is obtained in the bond market withing the affine-term structure setting, and it is used to formulate a deep-learning approach to arbitrage-free affine term-structure modelling. Numerical implementations are also performed to evaluate the performance in the bond market.


title: Benefits from Superposed Hawkes Processes

abstract: The superposition of temporal point processes has been studied for many years, although the usefulness of such models for practical applications has not be fully developed. We investigate superposed Hawkes process as an important class of such models, with properties studied in the framework of least squares estimation. The superposition of Hawkes processes is demonstrated to be beneficial for tightening the upper bound of excess risk under certain conditions, and we show the feasibility of the benefit in typical situations. The usefulness of superposed Hawkes processes is verified on synthetic data, and its potential to solve the cold-start problem of recommendation systems is demonstrated on real-world data.


title: Photothermal effects during nanodiamond synthesis from a carbon aerogel in a laser-heated diamond anvil cell

abstract: Nanodiamonds have emerged as promising materials for quantum computing, biolabeling, and sensing due to their ability to host color centers with remarkable photostability and long spin-coherence times at room temperature. Recently, a bottom-up, high-pressure, high-temperature (HPHT) approach was demonstrated for growing nanodiamonds with color centers from amorphous carbon precursors in a laser-heated diamond anvil cell (LH-DAC) that was supported by a near-hydrostatic noble gas pressure medium. However, a detailed understanding of the photothermal heating and its effect on diamond growth, including the phase conversion conditions and the temperature-dependence of color center formation, has not been reported. In this work, we measure blackbody radiation during LH-DAC synthesis of nanodiamond from carbon aerogel to examine these temperature-dependent effects. Blackbody temperature measurements suggest that nanodiamond growth can occur at 16.3 GPa and 1800 K. We use Mie theory and analytical heat transport to develop a predictive photothermal heating model. This model demonstrates that melting the noble gas pressure medium during laser heating decreases the local thermal conductivity to drive a high spatial resolution of phase conversion to diamond. Finally, we observe a temperature-dependent formation of nitrogen vacancy centers and interpret this phenomenon in the context of HPHT carbon vacancy diffusion using CB{\Omega} theory.


title: Computing the maximum matching width is NP-hard

abstract: The maximum matching width is a graph width parameter that is defined on a branch-decomposition over the vertex set of a graph. In this short paper, we prove that the problem of computing the maximum matching width is NP-hard.


title: Cutting a part from many measures

abstract: Holmsen, Kyn\v{c}l and Valculescu recently conjectured that if a finite set $X$ with $\ell n$ points in $\mathbb{R}^d$ that is colored by $m$ different colors can be partitioned into $n$ subsets of $\ell$ points each, such that each subset contains points of at least $d$ different colors, then there exists such a partition of $X$ with the additional property that the convex hulls of the $n$ subsets are pairwise disjoint. We prove a continuous analogue of this conjecture, generalized so that each subset contains points of at least $c$ different colors, where we also allow $c$ to be greater than $d$. Furthermore, we give lower bounds on the fraction of the points each of the subsets contains from $c$ different colors. For example, when $n\geq 2$, $d\geq 2$, $c\geq d$ with $m\geq n(c-d)+d$ are integers, and $\mu_1, \dots, \mu_m$ are $m$ positive finite absolutely continuous measures on $\mathbb{R}^d$, we prove that there exists a partition of $\mathbb{R}^d$ into $n$ convex pieces which equiparts the measures $\mu_1, \dots, \mu_{d-1}$, and in addition every piece of the partition has positive measure with respect to at least $c$ of the measures $\mu_1, \dots, \mu_m$.


title: Stability of nonlocal Dirichlet integrals and implications for peridynamic correspondence material modeling

abstract: Nonlocal gradient operators are basic elements of nonlocal vector calculus that play important roles in nonlocal modeling and analysis. In this work, we extend earlier analysis on nonlocal gradient operators. In particular, we study a nonlocal Dirichlet integral that is given by a quadratic energy functional based on nonlocal gradients. Our main finding, which differs from claims made in previous studies, is that the coercivity and stability of this nonlocal continuum energy functional may hold for some properly chosen nonlocal interaction kernels but may fail for some other ones. This can be significant for possible applications of nonlocal gradient operators in various nonlocal models. In particular, we discuss some important implications for the peridynamic correspondence material models.


title: Quantum phase transition in a realistic double-quantum-dot system

abstract: Observing quantum phase transitions in mesoscopic systems is a daunting task, thwarted by the difficulty of experimentally varying the magnetic interactions, the typical driving force behind these phase transitions. Here we demonstrate that in realistic coupled double-dot systems, the level energy difference between the two dots, which can be easily tuned experimentally, can drive the system through a phase transition, when its value crosses the difference between the intra- and inter-dot Coulomb repulsion. Using the numerical renormalization group and the semi-analytic slave-boson mean-field theory, we study the nature of this phase transition, and demonstrate, by mapping the Hamiltonian into an even-odd basis, that indeed the competition between the dot level energy difference and the difference in repulsion energies governs the sign and magnitude of the effective magnetic interaction. The observational consequences of this transition are discussed.


title: Netrunner Mate-in-1 or -2 is Weakly NP-Hard

abstract: We prove that deciding whether the Runner can win this turn (mate-in-1) in the Netrunner card game generalized to allow decks to contain an arbitrary number of copies of a card is weakly NP-hard. We also prove that deciding whether the Corp can win within two turns (mate-in-2) in this generalized Netrunner is weakly NP-hard.


title: Entangled state fusion with Rydberg atoms

abstract: We propose a scheme for preparation of large-scale entangled $GHZ$ states and $W$ states with neutral Rydberg atoms. The scheme mainly depends on Rydberg antiblockade effect, i.e., as the Rydberg-Rydberg-interaction (RRI) strength and the detuning between the atom transition frequency and the classical laser frequency satisfies some certain conditions, the effective Rabi oscillation between the two ground states and the two excitation Rydberg states would be generated. The prominent advantage is that both two-multiparticle $GHZ$ states and two-multiparticle $W$ states can be fused in this model, especially the success probability for fusion of $GHZ$ states can reach unit. In addition, the imperfections induced by the spontaneous emission is also discussed through numerical simulation.


title: Hom and Ext, Revisited

abstract: Let $R$ be a commutative Noetherian local ring and $M,N$ be finitely generated $R$-modules. We prove a number of results of the form: if $\mbox{Hom}_R(M,N)$ has some nice properties and $\mbox{Ext}^{1 \leq i \leq n}_R(M,N)=0$ for some $n$, then $M$ (and sometimes $N$) must be be close to free. Our methods are quite elementary, yet they suffice to give a unified treatment, simplify, and sometimes extend a number of results in the literature.


title: Monomial ideals with large projective dimension

abstract: Let S be a polynomial ring in n variables, over an arbitrary field. We give the total, graded, and multigraded Betti numbers of S/M, for every monomial ideal M in S. We also give an explicit characterization of all monomial ideals M in S for which the quotient S/M has projective dimension n.


title: Condensate $\langle A_{\mu}^+A_{\mu}^-\rangle$ and massive magnetic potential in Euclidean gauge theories

abstract: Euclidean $SU(2)$ gauge theory is studied in a nonlinear gauge. In this theory, ghost condensation happens and gauge fields acquire tachyonic masses. It is shown that these tachyonic masses are removed by a gauge field condensate $\langle A_{\mu}^+A_{\mu}^-\rangle$. Because of the ghost condensation, monopole solutions are included naturally. We find the condensate $\langle A_{\mu}^+A_{\mu}^-\rangle$ makes the magnetic potential massive.


title: Hierarchical semantic segmentation using modular convolutional neural networks

abstract: Image recognition tasks that involve identifying parts of an object or the contents of a vessel can be viewed as a hierarchical problem, which can be solved by initial recognition of the main object, followed by recognition of its parts or contents. To achieve such modular recognition, it is necessary to use the output of one recognition method (which identifies the general object) as the input for a second method (which identifies the parts or contents). In recent years, convolutional neural networks have emerged as the dominant method for segmentation and classification of images. This work examines a method for serially connecting convolutional neural networks for semantic segmentation of materials inside transparent vessels. It applies one fully convolutional neural net to segment the image into vessel and background, and the vessel region is used as an input for a second net which recognizes the contents of the glass vessel. Transferring the segmentation map generated by the first nets to the second net was performed using the valve filter attention method that involves using different filters on different segments of the image. This modular semantic segmentation method outperforms a single step method in which both the vessel and its contents are identified using a single net. An advantage of the modular neural net is that it allows networks to be built from existing trained modules, as well the transfer and reuse of trained net modules without the need for any retraining of the assembled net.


title: Splitting theorem for sheaves of holomorphic $k$-vectors on complex contact manifolds

abstract: A complex contact structure $\gamma$ is defined by a system of holomorphic local $1$-forms satisfying the completely non-integrability condition. The contact structure induces a subbundle ${\rm Ker}\, \gamma$ of the tangent bundle and a line bundle $L$. In this paper, we prove that the sheaf of holomorphic $k$-vectors on a complex contact manifold splits into the sum of $\mathcal{O}(\bigwedge^{k}{\rm Ker}\, \gamma)$ and $\mathcal{O}(L\otimes \bigwedge^{k-1} {\rm Ker}\, \gamma)$ as sheaves of {\it $\mathbb{C}$-module}. The theorem induces the short exact sequence of cohomology of holomorphic $k$-vectors, and we obtain vanishing theorems for the cohomology of $\mathcal{O}(\bigwedge^{k} \ker \gamma)$.


title: Parametric t-Distributed Stochastic Exemplar-centered Embedding

abstract: Parametric embedding methods such as parametric t-SNE (pt-SNE) have been widely adopted for data visualization and out-of-sample data embedding without further computationally expensive optimization or approximation. However, the performance of pt-SNE is highly sensitive to the hyper-parameter batch size due to conflicting optimization goals, and often produces dramatically different embeddings with different choices of user-defined perplexities. To effectively solve these issues, we present parametric t-distributed stochastic exemplar-centered embedding methods. Our strategy learns embedding parameters by comparing given data only with precomputed exemplars, resulting in a cost function with linear computational and memory complexity, which is further reduced by noise contrastive samples. Moreover, we propose a shallow embedding network with high-order feature interactions for data visualization, which is much easier to tune but produces comparable performance in contrast to a deep neural network employed by pt-SNE. We empirically demonstrate, using several benchmark datasets, that our proposed methods significantly outperform pt-SNE in terms of robustness, visual effects, and quantitative evaluations.


title: Shortcuts to adiabaticity in non-Hermitian quantum systems without rotating-wave approximation

abstract: The technique of shortcuts to adiabaticity (STA) has attracted broad attention due to their possible applications in quantum information processing and quantum control. However, most studies published so far have been only focused on Hermitian systems under the rotating-wave approximation (RWA). In this paper, we propose a modified STA technique to realize population transfer for a non-Hermitian system without RWA. We work out an exact expression for the control function and present examples consisting of two- and three-level systems with decay to show the theory. The results suggest that the STA technique presented here is robust for fast passages. We also find that the decay has small effect on the population transfer in the three-level system. To shed more light on the physics behind this result, we reduce the quantum three-level system to an effective two-level one with large detunings. The STA technique of effective two-level system is studied. Thereby the high-fidelity population transfer can be implemented in non-Hermitian systems by our method, and it works even without RWA.


title: MinDelay: Low-latency Forwarding and Caching Algorithms for Information-Centric Networks

abstract: We present a new unified framework for minimizing congestion-dependent network cost in information-centric networks by jointly optimizing forwarding and caching strategies. As caching variables are integer-constrained, the resulting optimization problem is NP-hard. To make progress, we focus on a relaxed version of the optimization problem, where caching variables are allowed to be real-valued. We develop necessary optimality conditions for the relaxed problem, and leverage this result to design MinDelay, an adaptive and distributed joint forwarding and caching algorithm, based on the conditional gradient algorithm. The MinDelay algorithm elegantly yields feasible routing variables and integer caching variables at each iteration, and can be implemented in a distributed manner with low complexity and overhead. Over a wide range of network topologies, simulation results show that MinDelay typically has significantly better delay performance in the low to moderate request rate regions. Furthermore, the MinDelay and VIP algorithms complement each other in delivering superior delay performance across the entire range of request arrival rates.


title: Mean Field Game Approach to Production and Exploration of Exhaustible Commodities

abstract: In a game theoretic framework, we study energy markets with a continuum of homogenous producers who produce energy from an exhaustible resource such as oil. Each producer simultaneously optimizes production rate that drives her revenues, as well as exploration effort to replenish her reserves. This exploration activity is modeled through a controlled point process that leads to stochastic increments to reserves level. The producers interact with each other through the market price that depends on the aggregate production. We employ a mean field game approach to solve for a Markov Nash equilibrium and develop numerical schemes to solve the resulting system of non-local HJB and transport equations with non-local coupling. A time-stationary formulation is also explored, as well as the fluid limit where exploration becomes deterministic.


title: On pseudospectral bound for non-selfadjoint operators and its application to stability of Kolmogorov flows

abstract: We study the stability of the Kolmogorov flows which are stationary solutions to the two-dimensional Navier-Stokes equations in the presence of the shear external force. We establish the linear stability estimate when the viscosity coefficient $\nu$ is sufficiently small, where the enhanced dissipation is rigorously verified in the time scale $O(\nu^{-\frac12})$ for solutions to the linearized problem, which has been numerically conjectured and is much shorter than the usual viscous time scale $O(\nu^{-1})$. Our approach is based on the detailed analysis for the resolvent problem. We also provide the abstract framework which is applicable to the resolvent estimate for the Kolmogorov flows.


title: Tracking Moving Agents via Inexact Online Gradient Descent Algorithm

abstract: Multi-agent systems are being increasingly deployed in challenging environments for performing complex tasks such as multi-target tracking, search-and-rescue, and intrusion detection. Notwithstanding the computational limitations of individual robots, such systems rely on collaboration to sense and react to the environment. This paper formulates the generic target tracking problem as a time-varying optimization problem and puts forth an inexact online gradient descent method for solving it sequentially. The performance of the proposed algorithm is studied by characterizing its dynamic regret, a notion common to the online learning literature. Building upon the existing results, we provide improved regret rates that not only allow non-strongly convex costs but also explicating the role of the cumulative gradient error. Two distinct classes of problems are considered: one in which the objective function adheres to a quadratic growth condition, and another where the objective function is convex but the variable belongs to a compact domain. For both cases, results are developed while allowing the error to be either adversarial or arising from a white noise process. Further, the generality of the proposed framework is demonstrated by developing online variants of existing stochastic gradient algorithms and interpreting them as special cases of the proposed inexact gradient method. The efficacy of the proposed inexact gradient framework is established on a multi-agent multi-target tracking problem, while its flexibility is exemplified by generating online movie recommendations for Movielens $10$M dataset.


title: Solution of the $k$-th eigenvalue problem in large-scale electronic structure calculations

abstract: We consider computing the $k$-th eigenvalue and its corresponding eigenvector of a generalized Hermitian eigenvalue problem of $n\times n$ large sparse matrices. In electronic structure calculations, several properties of materials, such as those of optoelectronic device materials, are governed by the eigenpair with a material-specific index $k.$ We present a three-stage algorithm for computing the $k$-th eigenpair with validation of its index. In the first stage of the algorithm, we propose an efficient way of finding an interval containing the $k$-th eigenvalue $(1 \ll k \ll n)$ with a non-standard application of the Lanczos method. In the second stage, spectral bisection for large-scale problems is realized using a sparse direct linear solver to narrow down the interval of the $k$-th eigenvalue. In the third stage, we switch to a modified shift-and-invert Lanczos method to reduce bisection iterations and compute the $k$-th eigenpair with validation. Numerical results with problem sizes up to 1.5 million are reported, and the results demonstrate the accuracy and efficiency of the three-stage algorithm.


title: When Point Process Meets RNNs: Predicting Fine-Grained User Interests with Mutual Behavioral Infectivity

abstract: Predicting fine-grained interests of users with temporal behavior is important to personalization and information filtering applications. However, existing interest prediction methods are incapable of capturing the subtle degreed user interests towards particular items, and the internal time-varying drifting attention of individuals is not studied yet. Moreover, the prediction process can also be affected by inter-personal influence, known as behavioral mutual infectivity. Inspired by point process in modeling temporal point process, in this paper we present a deep prediction method based on two recurrent neural networks (RNNs) to jointly model each user's continuous browsing history and asynchronous event sequences in the context of inter-user behavioral mutual infectivity. Our model is able to predict the fine-grained interest from a user regarding a particular item and corresponding timestamps when an occurrence of event takes place. The proposed approach is more flexible to capture the dynamic characteristic of event sequences by using the temporal point process to model event data and timely update its intensity function by RNNs. Furthermore, to improve the interpretability of the model, the attention mechanism is introduced to emphasize both intra-personal and inter-personal behavior influence over time. Experiments on real datasets demonstrate that our model outperforms the state-of-the-art methods in fine-grained user interest prediction.


title: Probabilistic representation of weak solutions to a parabolic boundary value problem on a non-smooth domain

abstract: The probabilistic representation of weak solutions to a parabolic boundary value problem is established in the following framework. The boundary value problem consists of a second order parabolic equation defined on a time-varying Lipschitz domain in a Euclidean space and of a mixed boundary condition composed of a Robin and the homogeneous Dirichlet conditions. It is assumed that the time-varying domain is included in a fixed smooth domain and that a certain part of the boundary of the time-varying domain is also included in the boundary of the fixed domain, say the fixed boundary. The Robin condition is imposed on a part of the boundary included in the fixed one and the Dirichlet condition on the rest of the boundary. Such a parabolic boundary value problem always has a unique weak solution for given data; however it does not possess a classical or strong solution in general, even in the case of equations with constant coefficients. The stochastic solution to the boundary value problem is also considered and, by showing the equality between both the solutions, it is obtained the probabilistic representation for the weak solution. Furthermore, it is ensured that, for the weak solution, the stochastic solution gives a version which is continuous up to the lateral boundary of the domain except the border of the adjoining place imposed each of the boundary conditions. As an application, it is shown the continuity property of a functional (cost function) related to an optimal stopping problem motivated by an inverse problem determining the shape of a domain.


title: TASI Lectures on Indirect Detection of Dark Matter

abstract: These lectures, presented at TASI 2016: Anticipating the Next Discoveries in Particle Physics, provide an introduction to some key methods and tools of indirect dark matter searches. Topics covered include estimation of dark matter signals, thermal freezeout and related scenarios, potential effects of dark matter annihilation on the early universe, modeling photon signals from annihilation or decay, and a brief and qualitative introduction to diffusive propagation of cosmic rays. The second half of the notes gives a status report (circa summer 2016) on selected experimental searches, the resulting constraints and some potential signal candidates. These notes are intended as an introduction to indirect dark matter searches for graduate students, focusing on back-of-the-envelope estimates and useful concepts rather than detailed quantitative computations.


title: Homogeneous 3-dimensional permutation structures

abstract: We provide a classification of the homogeneous 3-dimensional permutation structures, i.e. homogeneous structures in a language of 3 linear orders, partially answering a question of Cameron. We also arrive at a natural description of all known homogeneous finite-dimensional permutation structures by modifying the language used in the construction from a previous paper, completing the "census" begun there.


title: Light nuclei production in Pb+Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV

abstract: Using the coalescence model based on the phase-space distribution of nucleons from an extended blast-wave model that includes the space-momentum correlation of high momentum nucleons, we study the transverse momentum spectra and elliptic flows of deuteron and helium-3 in Pb+Pb collisions at the energy $\sqrt{s_{NN}}=2.76$ TeV. We find that the measured elliptic flows of deuteron can be satisfactorily described if nucleons of large transverse momenta are more spread in space when their momenta are more aligned along the reaction plane.


title: On complexity of mutlidistance graph recognition in $\mathbb{R}^1$

abstract: Let $\mathcal{A}$ be a set of positive numbers. A graph $G$ is called an $\mathcal{A}$-embeddable graph in $\mathbb{R}^d$ if the vertices of $G$ can be positioned in $\mathbb{R}^d$ so that the distance between endpoints of any edge is an element of $\mathcal{A}$. We consider the computational problem of recognizing $\mathcal{A}$-embeddable graphs in $\mathbb{R}^1$ and classify all finite sets $\mathcal{A}$ by complexity of this problem in several natural variations.


title: Digital Currency Design for Sustainable Active Debris Removal in Space

abstract: Orbital debris remains as an obstacle to further space development. While efforts are ongoing to avoid newly launched objects becoming debris, the number of debris would still continue to grow because of collisions. ADR (Active Debris Removal) is an effective measure, but building a sustainable economic model for ADR remains as a difficult problem. We propose that the cost of removal can be paid by circulating digital currency tokens on a blockchain platform whose values may decrease and/or increase over time, issued by global cooperation (a consortium) of parties interested in space development, in exchange with proofs of ADR. The tokens pay their cost by themselves through contributions by the token holders, who are likely to be benefited by removal of debris. This scheme imposes virtually no cost to the consortium. We have generalized this concept as POD (Proof of Disposal), which, we believe, provides a more accountable foundation for solving social problems with digital currency than many ICO (Initial Coin or Cryptoasset Offering) in practice today. We evaluated the feasibility of our proposal through a simulation. We conclude that dynamic estimation of the economic values of each ADR and automated pricing of tokens that represent the orbital debris being removed are indeed possible. Actual prototyping of the proposed digital currency system is ongoing.


title: On the Origin of Banded Structure in Dusty Protoplanetary Discs: HL Tau and TW Hya

abstract: Recent observations of HL Tau revealed remarkably detailed structure within the system's circumstellar disc. A range of hypotheses have been proposed to explain the morphology, including, e.g., planet-disc interactions, condensation fronts, and secular gravitational instabilities. While embedded planets seem to be able to explain some of the major structure in the disc through interactions with gas and dust, the substructure, such as low-contrast rings and bands, are not so easily reproduced. Here, we show that dynamical interactions between three planets (only two of which are modelled) and an initial population of large planetesimals can potentially explain both the major and minor banded features within the system. In this context, the small grains, which are coupled to the gas and reveal the disc morphology, are produced by the collisional evolution of the newly-formed planetesimals, which are ubiquitous in the system and are decoupled from the gas.


title: A Complementary Inequality to the Information Monotonicity for Tsallis Relative Operator Entropy

abstract: We establish a reverse inequality for Tsallis relative operator entropy involving a positive linear map. In addition, we present converse of Ando's inequality, for each parameter. We give examples to compare our results with the known results by Furuta and Seo. In particular, we establish an extension and a reverse of the L\"owner-Heinz inequality under certain condition. Some interesting consequences of inner product spaces and norm inequalities are also presented.


title: Open Type Nodal line Topological Semimetal in Two Dimensional B2C

abstract: Topological semimetals, including Dirac semimetals, Weyl semimetals, and nodal line semimetals, receive enormous research interest due to their intrinsic topological nature and fascinating properties. In present work, with the help of density functional theory (DFT), we find that nodal line semimetals and Dirac semimetals can coexist in the low energy electron of two dimensional (2D) monolayer B2C. Intriguingly, besides type-I and type-II Dirac fermions, a kind of open nodal line appears around the Fermi level. The low energy electronic nature of B2C sheet can be described by a tight-binding (TB) model relied on the basis of B-py, pz and C-py, pz states. One of the most merit of the system is that the energy windows of these four types topological semimetals are different that can be easily distinguished in experiments. Moreover, the system provides an excellent platform for studying the interplay between different type semimetals.


title: Lorentz-boosted evanescent waves

abstract: Polarization, spin, and helicity are important properties of electromagnetic waves. It is commonly believed that helicity is invariant under the Lorentz transformations. This is indeed so for plane waves and their localized superpositions. However, this is not the case for evanescent waves, which are well-defined only in a half-space, and are characterized by complex wave vectors. Here we describe transformations of evanescent electromagnetic waves and their polarization/spin/helicity properties under the Lorentz boosts along the three spatial directions.


title: Two phase helium cooling characteristics in Cable-in Conduit Conductors

abstract: Cable-in-Conduit Conductors (CICCs) are used in the fabrication of superconducting fusion grade magnets. It acts as a narrow cryostat to provide cryo-stability with direct contact of coolant fluid to conductor. The superconducting magnets are cooled using forced flow (FF), supercritical helium or two phase (TP) cooling through void space in the CICC. Thermo-hydraulics using supercritical helium single phase flow is well-known and established. Research topic of behavior of forced flow, two phase (TP) helium cooling in CICC involves perceived risks of the CICC running into flow chocking and possible thermo-acoustic oscillations leading to flow instabilities. This research work involves study of forced flow two phase helium cooling in CICC wound superconducting magnets. The TP flow provides cryo-stability by the latent heat of helium not by enthalpy as in case of CICC being cooled with supercritical helium. Study reveals some attractive regimes in the case of TP cooling, at a given mass flow rate of single phase helium at the inlet and a heat flux acting on the CICC. Analysis carried out predicts significant gains with TP cooling on a prototype CICC, which is circular in cross section and appropriate for fusion devices for high magnetic field applications. These general formalisms may be extended to specific magnets wound with CICC. This paper describes analysis of TP cooling of a CICC.


title: Cyclic Isogenies for Abelian Varieties with Real Multiplication

abstract: We study quotients of principally polarized abelian varieties with real multiplication by Galois-stable finite subgroups and describe when these quotients are principally polarizable. We use this characterization to provide an algorithm to compute explicit cyclic isogenies from kernel for abelian varieties with real multiplication over finite fields. Our algorithm is polynomial in the size of the finite field as well as the degree of the isogeny and is based on Mumford's theory of theta functions and theta embeddings. Recently, the algorithm has been successfully applied to obtain new results on the discrete logarithm problem in genus 2 as well as to study the discrete logarithm problem in genus 3.


title: Tungsten disulfide-gold nanohole hybrid metasurfaces for nonlinear metalens in the visible region

abstract: Recently, nonlinear hybrid metasurface comes into an attractive new concept in the research of nanophotonics and nanotechnology. It is composed of semiconductors with an intrinsically large nonlinear susceptibility and traditional plasmonic metasurfaces, offering opportunities for efficiently generating and manipulating nonlinear optical responses. A high second-harmonic generation (SHG) conversion efficiency has been demonstrated in the mid-infrared region by using multi-quantum-well (MQW) based plasmonic metasurfaces. However, it has yet to be demonstrated in the visible region. Here we present a new type of nonlinear hybrid metasurfaces for the visible region, which consists of a single layer of tungsten disulfide (WS2) and a phased gold nanohole array. The results indicate that a large SHG susceptibility of ~0.1 nm/V at 810 nm is achieved, which is 2~3 orders of magnitude larger than that of typical plasmonic metasurfaces. Nonlinear metalenses with the focal lengths of 30 {\mu}m, 50 {\mu}m and 100 {\mu}m are demonstrated experimentally, providing a direct evidence for both generating and manipulating SH signals based on the nonlinear hybrid metasurfaces. It shows great potential applications in designing of integrated, ultra-thin, compacted and efficient nonlinear optical devices, such as frequency converters, nonlinear holography and generation of nonlinear optical vortex beam.


title: Resource Allocation for D2D Communications Underlaying a NOMA-Based Cellular Network

abstract: This letter investigates the power control and channel assignment problem in device-to-device (D2D) communications underlaying a non-orthogonal multiple access (NOMA) cellular network. With the successive interference cancellation decoding order constraints, our target is to maximize the sum rate of D2D pairs while guaranteeing the minimum rate requirements of NOMA-based cellular users. Specifically, the optimal conditions for power control of cellular users on each subchannel are derived first. Then, based on these results, we propose a dual-based iterative algorithm to solve the resource allocation problem. Simulation results validate the superiority of proposed resource allocation algorithm over the existing orthogonal multiple access scheme.


title: Inversion of the decay of the cyclotron line energy in Her X-1

abstract: Recent observations of Her X-1 with NuSTAR and INTEGRAL in 2016 have provided evidence that the 20-year decay of the cyclotron line energy found between 1996 and 2015 has ended and that an inversion with a new increase, possibly similar to the one observed around 1990-1993, has started. We consider this a strong motivation for further observations and for enhanced efforts to significantly improve our theoretical understanding of the accretion process in binary X-ray pulsars. We speculate about the physics behind the long-term decay and its inversion, a possible cyclic behavior, and correlations with other variable observables.


title: Notes on toric varieties from Mori theoretic viewpoint, II

abstract: We give new estimates of lengths of extremal rays of birational type for toric varieties. We can see that our new estimates are the best by constructing some examples explicitly. As applications, we discuss the nefness and pseudo-effectivity of adjoint bundles of projective toric varieties. We also treat some generalizations of Fujita's freeness and very ampleness for toric varieties.


title: GHCLNet: A Generalized Hierarchically tuned Contact Lens detection Network

abstract: Iris serves as one of the best biometric modality owing to its complex, unique and stable structure. However, it can still be spoofed using fabricated eyeballs and contact lens. Accurate identification of contact lens is must for reliable performance of any biometric authentication system based on this modality. In this paper, we present a novel approach for detecting contact lens using a Generalized Hierarchically tuned Contact Lens detection Network (GHCLNet) . We have proposed hierarchical architecture for three class oculus classification namely: no lens, soft lens and cosmetic lens. Our network architecture is inspired by ResNet-50 model. This network works on raw input iris images without any pre-processing and segmentation requirement and this is one of its prodigious strength. We have performed extensive experimentation on two publicly available data-sets namely: 1)IIIT-D 2)ND and on IIT-K data-set (not publicly available) to ensure the generalizability of our network. The proposed architecture results are quite promising and outperforms the available state-of-the-art lens detection algorithms.


title: $W'$ Pair Production in the Light of CMS Searches

abstract: For the first time, the pair production of the heavy charged gauge bosons, known as $W'$ bosons is considered, when both decay to $\tau$ leptons. The reported detailed efficiency of object/event selection by the CMS experiment is used to find the lower limit on the mass of $W'$ boson. Various assumptions for the coupling of the new gauge boson are examined and the results are reported. In the case of a SM-like $W'$ boson, masses below 290 GeV are excluded at 95\% confidence level. The method can be used to constrain the new models with the similar final state.


title: High Throughput 2D Spatial Image Filters on FPGAs

abstract: FPGAs are well established in the signal processing domain, where their fine-grained programmable nature allows the inherent parallelism in these applications to be exploited for enhanced performance. As architectures have evolved, FPGA vendors have added more heterogeneous resources to allow often-used functions to be implemented with higher performance, at lower power and using less area. DSP blocks, for example, have evolved from basic multipliers to support the multiply-accumulate operations that are the core of many signal processing tasks. While more features were added to DSP blocks, their structure and connectivity has been optimised primarily for one-dimensional signal processing. Basic operations in image processing are similar, but performed in a two-dimensional structure, and hence, many of the optimisations in newer DSP blocks are not exploited when mapping image processing algorithms to them. We present a detailed study of two-dimensional spatial filter implementation on FPGAs, showing how to maximise performance through exploitation of DSP block capabilities, while also presenting a lean border pixel management policy.


title: A Multispin Algorithm for the Kob-Andersen Stochastic Dynamics on Regular Lattices

abstract: The aim of the paper is to propose an algorithm based on the Multispin Coding technique for the Kob-Andersen glassy dynamics. We first give motivations to speed up the numerical simulation in the context of spin glass models [M. Mezard, G. Parisi, M. Virasoro, Spin Glass Theory and Beyond (World Scientific, Singapore, 1987)], after defining the Markovian dynamics as in [W. Kob, H.C. Andersen, Phys. Rev. E 48, 4364 (1993)] as well as the related interesting observables, we extend it to the more general framework of random regular graphs, listing at the same time some known analytical results [C. Toninelli, G. Biroli, D.S. Fisher, J. Stat. Phys. 120, 167 (2005)]. The purpose of this work is a dual one, firstly, we describe how bitwise operators can be used to build up the algorithm by carefully exploiting the way data are stored on a computer. Since it was first introduced [M. Creutz, L. Jacobs, C. Rebbi, Phys. Rev. D 20, 1915 (1979), C. Rebbi, R.H. Swendsen, Phys. Rev. D 21, 4094 (1980)], this technique has been widely used to perform Monte Carlo simulations for Ising and Potts spin systems, however, it can be successfully adapted to more complex systems in which microscopic parameters may assume boolean values. Secondly, we introduce a random graph in which a characteristic para- meter allows to tune the possible transition point. A consistent part is devoted to listing the numerical results obtained by running numerical simulations.


title: Asymptotic enumeration of perfect matchings in $m$-barrel fullerene graphs

abstract: A connected planar cubic graph is called an $m$-barrel fullerene and denoted by $F(m,k)$, if it has the following structure: The first circle is an $m$-gon. Then $m$-gon is bounded by $m$ pentagons. After that we have additional k layers of hexagons. At the last circle $m$-pentagons connected to the second $m$-gon. In this paper we asymptotically count by two different methods the number of perfect matchings in $m$-barrel fullerene graphs, as the number of hexagonal layers is large, and show that the results are equal.


title: Metric systolicity and two-dimensional Artin groups

abstract: We introduce the notion of metrically systolic simplicial complexes. We study geometric and large-scale properties of such complexes and of groups acting on them geometrically. We show that all two-dimensional Artin groups act geometrically on metrically systolic complexes. As direct corollaries we obtain new results on two-dimensional Artin groups and all their finitely presented subgroups: we prove that the Conjugacy Problem is solvable, and that the Dehn function is quadratic. We also show several large-scale features of finitely presented subgroups of two-dimensional Artin groups, lying background for further studies concerning their quasi-isometric rigidity.


title: BrainSegNet : A Segmentation Network for Human Brain Fiber Tractography Data into Anatomically Meaningful Clusters

abstract: The segregation of brain fiber tractography data into distinct and anatomically meaningful clusters can help to comprehend the complex brain structure and early investigation and management of various neural disorders. We propose a novel stacked bidirectional long short-term memory(LSTM) based segmentation network, (BrainSegNet) for human brain fiber tractography data classification. We perform a two-level hierarchical classification a) White vs Grey matter (Macro) and b) White matter clusters (Micro). BrainSegNet is trained over three brain tractography data having over 250,000 fibers each. Our experimental evaluation shows that our model achieves state-of-the-art results. We have performed inter as well as intra class testing over three patient's brain tractography data and achieved a high classification accuracy for both macro and micro levels both under intra as well as inter brain testing scenario.


title: Evolution of initial discontinuities in the DNLS equation theory

abstract: We present the full classification of wave patterns evolving from an initial step-like discontinuity for arbitrary choice of boundary conditions at the discontinuity location in the DNLS equation theory. In this non-convex dispersive hydrodynamics problem, solutions of the Whitham modulation equations are mapped to parameters of a modulated wave by two-valued functions what makes situation much richer than that for a convex case of the NLS equation type. In particular, new types of simple-wave-like structures appear as building elements of the whole wave pattern. The developed here theory can find applications to propagation of light pulses in fibers and to the theory of Alfv\'en dispersive shock waves.


title: Hyper-reduction over nonlinear manifolds for large nonlinear mechanical systems

abstract: Common trends in model order reduction of large nonlinear finite-element-discretized systems involve the introduction of a linear mapping into a reduced set of unknowns, followed by Galerkin projection of the governing equations onto a constant reduction basis. Though this reduces the number of unknowns in the system, the computational cost for obtaining the solution could still be high due to the prohibitive computational costs involved in the evaluation of nonlinear terms. Hyper-reduction methods are then seen as a fast way of approximating the nonlinearity in the system of equations. In the finite element context, the energy conserving sampling and weighing (ECSW) method has emerged as a stability and structure-preserving method for hyper-reduction. Classical hyper-reduction techniques, however, are applicable only in the context of linear mappings into the reduction subspace. In this work, we extend the concept of hyper-reduction using ECSW to general nonlinear mappings, while retaining its desirable stability and structure-preserving properties. As a proof of concept, the proposed hyper-reduction technique is demonstrated over models of a flat plate and a realistic wing structure, whose dynamics has been shown to evolve over a nonlinear (quadratic) manifold. An online speed-up of over one thousand times relative to the full system has been obtained for the wing structure using the proposed method, which is higher than its linear counterpart using the ECSW.


title: Rota-Baxter bisystems and covariant bialgebras

abstract: The aim of this paper is first to introduce and study Rota-Baxter cosystems and bisystems as generalization of Rota-Baxter coalgebras and bialgebras, respectively, with various examples. The second purpose is to provide an alternative definition of covariant bialgebras via coderivations. Furthermore, we consider coquasitriangular covariant bialgebras which are generalization of coquasitriangular infinitesimal bialgebras, coassociative Yang-Baxter pairs, coassociative Yang-Baxter equation, double coalgebras and dendriform coalgebras. We study some of their properties and relationships with Rota-Baxter cosystems and bisystems.


title: Density split statistics: joint model of counts and lensing in cells

abstract: We present density split statistics, a framework that studies lensing and counts-in-cells as a function of foreground galaxy density, thereby providing a large-scale measurement of both 2-point and 3-point statistics. Our method extends our earlier work on trough lensing and is summarized as follows: given a foreground (low redshift) population of galaxies, we divide the sky into subareas of equal size but distinct galaxy density. We then measure lensing around uniformly spaced points separately in each of these subareas, as well as counts-in-cells statistics (CiC). The lensing signals trace the matter density contrast around regions of fixed galaxy density. Through the CiC measurements this can be related to the density profile around regions of fixed matter density. Together, these measurements constitute a powerful probe of cosmology, the skewness of the density field and the connection of galaxies and matter. In this paper we show how to model both the density split lensing signal and CiC from basic ingredients: a non-linear power spectrum, clustering hierarchy coefficients from perturbation theory and a parametric model for galaxy bias and shot-noise. Using N-body simulations, we demonstrate that this model is sufficiently accurate for a cosmological analysis on year 1 data from the Dark Energy Survey.


title: An Improved Modified Cholesky Decomposition Method for Precision Matrix Estimation

abstract: The modified Cholesky decomposition is commonly used for precision matrix estimation given a specified order of random variables. However, the order of variables is often not available or cannot be pre-determined. In this work, we propose to address the variable order issue in the modified Cholesky decomposition for sparse precision matrix estimation. The key idea is to effectively combine a set of estimates obtained from multiple permutations of variable orders, and to efficiently encourage the sparse structure for the resultant estimate by the thresholding technique on the ensemble Cholesky factor matrix. The consistent property of the proposed estimate is established under some weak regularity conditions. Simulation studies are conducted to evaluate the performance of the proposed method in comparison with several existing approaches. The proposed method is also applied into linear discriminant analysis of real data for classification.


title: Contraction of Dirac matrices via chord diagrams

abstract: Chord diagrams and combinatorics of word algebras are used to model products of Dirac matrices, their traces, and contractions. A simple formula for the result of arbitrary contractions is derived, simplifying and extending an old contraction algorithm due to Kahane. This formula is then used to express the Schwinger parametric integrand of a QED Feynman integral in a much simplified form, with the entire internal tensor structure eliminated. Possible next steps for further simplification, including a specific conjecture, are discussed.


title: Irreducible polynomials of bounded height

abstract: The goal of this paper is to prove that a random polynomial with i.i.d. random coefficients taking values uniformly in $\{1,\ldots, 210\}$ is irreducible with probability tending to $1$ as the degree tends to infinity. Moreover, we prove that the Galois group of the random polynomial contains the alternating group, again with probability tending to $1$.