title: An Ontological Learning Management System

abstract: The current learning systems typically lack the level of metacognitive awareness, self-directed learning, and time management skills. Most of the ontologically based learning management systems are in the proposed phase and those which are developed do not provide the necessary path guidance for proper learning. The systems available are not as adaptive from the viewpoint of the learner as required. Ontology engineering has become an important pillar for knowledge management and representation in recent years. The design, approach, and implementation of ontology in e-learning and m-learning systems have made them more effective. In this paper, we have proposed a system for the betterment of knowledge management and representation of associated data as compared to the previously available learning management systems. Here, we have presented the application and implementation of ontological engineering methodology in the Computer Science domain. For knowledge management, we have created a domain associated ontology which represents knowledge of a single domain. Subsequently, ontology has been created to manage a learner profile so that a learner may be aligned to a proper path of learning. The learner ontology will use the VARK learning model which classifies what kind of learning does the learner requires so that necessary resources could be provided.


title: Radon Daughter Plate-out Measurements at SNOLAB for Polyethylene and Copper

abstract: Polyethylene and copper samples were exposed to the underground air at SNOLAB for approximately three months while several environmental factors were monitored. Predictions of the radon-daughter plate-out rate are compared to the resulting surface activities, obtained from high-sensitivity measurements of alpha emissivity using the XIA UltraLo-1800 spectrometer at SMU. From these measurements, we determine an average $^{210}$Pb plate-out rate of 249 and 423 atoms/day/cm$^{2}$ for polyethylene and copper, respectively, when exposed to radon activity of 135 Bq/m$^{3}$ at SNOLAB. A time-dependent model of alpha activity is discussed for these materials placed in similar environmental conditions.


title: A Compressive Sensing Approach to Community Detection with Applications

abstract: The community detection problem for graphs asks one to partition the n vertices V of a graph G into k communities, or clusters, such that there are many intracluster edges and few intercluster edges. Of course this is equivalent to finding a permutation matrix P such that, if A denotes the adjacency matrix of G, then PAP^T is approximately block diagonal. As there are k^n possible partitions of n vertices into k subsets, directly determining the optimal clustering is clearly infeasible. Instead one seeks to solve a more tractable approximation to the clustering problem. In this paper we reformulate the community detection problem via sparse solution of a linear system associated with the Laplacian of a graph G and then develop a two-stage approach based on a thresholding technique and a compressive sensing algorithm to find a sparse solution which corresponds to the community containing a vertex of interest in G. Crucially, our approach results in an algorithm which is able to find a single cluster of size n_0 in O(nlog(n)n_0) operations and all k clusters in fewer than O(n^2ln(n)) operations. This is a marked improvement over the classic spectral clustering algorithm, which is unable to find a single cluster at a time and takes approximately O(n^3) operations to find all k clusters. Moreover, we are able to provide robust guarantees of success for the case where G is drawn at random from the Stochastic Block Model, a popular model for graphs with clusters. Extensive numerical results are also provided, showing the efficacy of our algorithm on both synthetic and real-world data sets.


title: Weighted real Egyptian numbers

abstract: Let $\mathcal A = (A_1,\ldots, A_n)$ be a sequence of nonempty finite sets of positive real numbers, and let $\mathcal{B} = (B_1,\ldots, B_n)$ be a sequence of infinite discrete sets of positive real numbers. A weighted real Egyptian number with numerators $\mathcal{A}$ and denominators $\mathcal{B}$ is a real number $c$ that can be represented in the form \[ c = \sum_{i=1}^n \frac{a_i}{b_i} \] with $a_i \in A_i$ and $b_i \in B_i$ for $i \in \{1,\ldots, n\}$. In this paper, classical results of Sierpinski for Egyptian fractions are extended to the set of weighted real Egyptian numbers.


title: Graphical Lasso and Thresholding: Equivalence and Closed-form Solutions

abstract: Graphical Lasso (GL) is a popular method for learning the structure of an undirected graphical model, which is based on an $l_1$ regularization technique. The objective of this paper is to compare the computationally-heavy GL technique with a numerically-cheap heuristic method that is based on simply thresholding the sample covariance matrix. To this end, two notions of sign-consistent and inverse-consistent matrices are developed, and then it is shown that the thresholding and GL methods are equivalent if: (i) the thresholded sample covariance matrix is both sign-consistent and inverse-consistent, and (ii) the gap between the largest thresholded and the smallest un-thresholded entries of the sample covariance matrix is not too small. By building upon this result, it is proved that the GL method---as a conic optimization problem---has an explicit closed-form solution if the thresholded sample covariance matrix has an acyclic structure. This result is then generalized to arbitrary sparse support graphs, where a formula is found to obtain an approximate solution of GL. Furthermore, it is shown that the approximation error of the derived explicit formula decreases exponentially fast with respect to the length of the minimum-length cycle of the sparsity graph. The developed results are demonstrated on synthetic data, functional MRI data, traffic flows for transportation networks, and massive randomly generated data sets. We show that the proposed method can obtain an accurate approximation of the GL for instances with the sizes as large as $80,000\times 80,000$ (more than 3.2 billion variables) in less than 30 minutes on a standard laptop computer running MATLAB, while other state-of-the-art methods do not converge within 4 hours.


title: The white dwarf binary pathways survey - II. Radial velocities of 1453 FGK stars with white dwarf companions from LAMOST DR4

abstract: We present the second paper of a series of publications aiming at obtaining a better understanding regarding the nature of type Ia supernovae (SNIa) progenitors by studying a large sample of detached F, G and K main sequence stars in close orbits with white dwarf companions (i.e. WD+FGK binaries). We employ the LAMOST (Large Sky Area Multi-Object Fibre Spectroscopic Telescope) data release 4 spectroscopic data base together with GALEX (Galaxy Evolution Explorer) ultraviolet fluxes to identify 1,549 WD+FGK binary candidates (1,057 of which are new), thus doubling the number of known sources. We measure the radial velocities of 1,453 of these binaries from the available LAMOST spectra and/or from spectra obtained by us at a wide variety of different telescopes around the globe. The analysis of the radial velocity data allows us to identify 24 systems displaying more than 3sigma radial velocity variation that we classify as close binaries. We also discuss the fraction of close binaries among WD+FGK systems, which we find to be ~10 per cent, and demonstrate that high-resolution spectroscopy is required to efficiently identify double-degenerate SNIa progenitor candidates.


title: Dynamic Bayesian Influenza Forecasting in the United States with Hierarchical Discrepancy

abstract: Timely and accurate forecasts of seasonal influenza would assist public health decision-makers in planning intervention strategies, efficiently allocating resources, and possibly saving lives. For these reasons, influenza forecasts are consequential. Producing timely and accurate influenza forecasts, however, have proven challenging due to noisy and limited data, an incomplete understanding of the disease transmission process, and the mismatch between the disease transmission process and the data-generating process. In this paper, we introduce a dynamic Bayesian (DB) flu forecasting model that exploits model discrepancy through a hierarchical model. The DB model allows forecasts of partially observed flu seasons to borrow discrepancy information from previously observed flu seasons. We compare the DB model to all models that competed in the CDC's 2015--2016 flu forecasting challenge. The DB model outperformed all models, indicating the DB model is a leading influenza forecasting model.


title: Exact fluctuations of nonequilibrium steady states from approximate auxiliary dynamics

abstract: We describe a framework to significantly reduce the computational effort to evaluate large deviation functions of time integrated observables within nonequilibrium steady states. We do this by incorporating an auxiliary dynamics into trajectory based Monte Carlo calculations, through a transformation of the system's propagator using an approximate guiding function. This procedure importance samples the trajectories that most contribute to the large deviation function, mitigating the exponentially complexity of such calculations. We illustrate the method by studying driven diffusions and interacting lattice models in one and two dimensions. Our work offers an avenue to calculate large deviation functions for high dimensional systems driven far from equilibrium.

categories: cond-mat.stat-mech cond-mat.soft physics.chem-ph physics.comp-ph


title: A question proposed by K. Mahler on exceptional sets of transcendental functions with integer coefficients: solution of a Mahler's problem

abstract: In this paper, we shall prove that any subset of $\overline{\mathbb Q}\cap B(0,1)$, which is closed under complex conjugation and which contains the element $0$, is the exceptional set of uncountably many transcendental functions, analytic in the unit ball, with integer coefficients. This solves a strong version of an old question proposed by K. Mahler (1976).


title: Temporal evolution of the high-energy irradiation and water content of TRAPPIST-1 exoplanets

abstract: The ultracool dwarf star TRAPPIST-1 hosts seven Earth-size transiting planets, some of which could harbour liquid water on their surfaces. UV observations are essential to measure their high-energy irradiation, and to search for photodissociated water escaping from their putative atmospheres. Our new observations of TRAPPIST-1 Ly-$\alpha$ line during the transit of TRAPPIST-1c show an evolution of the star emission over three months, preventing us from assessing the presence of an extended hydrogen exosphere. Based on the current knowledge of the stellar irradiation, we investigated the likely history of water loss in the system. Planets b to d might still be in a runaway phase, and planets within the orbit of TRAPPIST-1g could have lost more than 20 Earth oceans after 8 Gyr of hydrodynamic escape. However, TRAPPIST-1e to h might have lost less than 3 Earth oceans if hydrodynamic escape stopped once they entered the habitable zone. We caution that these estimates remain limited by the large uncertainty on the planet masses. They likely represent upper limits on the actual water loss because our assumptions maximize the XUV-driven escape, while photodissociation in the upper atmospheres should be the limiting process. Late-stage outgassing could also have contributed significant amounts of water for the outer, more massive planets after they entered the habitable zone. While our results suggest that the outer planets are the best candidates to search for water with the JWST, they also highlight the need for theoretical studies and complementary observations in all wavelength domains to determine the nature of the TRAPPIST-1 planets, and their potential habitability.


title: Learning Invariant Riemannian Geometric Representations Using Deep Nets

abstract: Non-Euclidean constraints are inherent in many kinds of data in computer vision and machine learning, typically as a result of specific invariance requirements that need to be respected during high-level inference. Often, these geometric constraints can be expressed in the language of Riemannian geometry, where conventional vector space machine learning does not apply directly. The central question this paper deals with is: How does one train deep neural nets whose final outputs are elements on a Riemannian manifold? To answer this, we propose a general framework for manifold-aware training of deep neural networks -- we utilize tangent spaces and exponential maps in order to convert the proposed problem into a form that allows us to bring current advances in deep learning to bear upon this problem. We describe two specific applications to demonstrate this approach: prediction of probability distributions for multi-class image classification, and prediction of illumination-invariant subspaces from a single face-image via regression on the Grassmannian. These applications show the generality of the proposed framework, and result in improved performance over baselines that ignore the geometry of the output space. In addition to solving this specific problem, we believe this paper opens new lines of enquiry centered on the implications of Riemannian geometry on deep architectures.


title: Topological Hochschild homology of X(n)

abstract: We show that Ravenel's spectrum $X(2)$ is the versal $E_1$-$S$-algebra of characteristic $\eta$. This implies that every $E_1$-$S$-algebra $R$ of characteristic $\eta$ admits an $E_1$-ring map $X(2)\to R$, i.e. an $\mathbb{A}_\infty$ complex orientation of degree 2. This implies that $R^\ast(\mathbb{C}P^2)\cong R_\ast[x]/x^3$. Additionally, if $R$ is an $\mathbb{E}_2$-ring Thom spectrum admitting a map (of homotopy ring spectra) from $X(2)$, e.g. $X(n)$, its topological Hochschild homology has a simple description.


title: Integrated nanoplasmonic quantum interfaces for room temperature single photon sources

abstract: We describe a general analytical framework of a nanoplasmonic cavity-emitter system interacting with a dielectric photonic waveguide. Taking into account emitter quenching and dephasing, our model directly reveals the single photon extraction efficiency, $\eta$, as well as the indistinguishability, $I$, of photons coupled into the waveguide mode. Rather than minimizing the cavity modal volume, our analysis predicts an optimum modal volume to maximize $\eta$ that balances waveguide coupling and spontaneous emission rate enhancement. Surprisingly, our model predicts that near-unity indistinguishability is possible, but this requires a much smaller modal volume, implying a fundamental performance trade-off between high $\eta$ and $I$ at room temperature. Finally, we show that maximizing $\eta I$ requires that the system has to be driven in the weak coupling regime because quenching effects and decreased waveguide coupling drastically reduce $\eta$ in the strong coupling regime.


title: On Scheduling a Photolithography Process Containing Cluster Tools

abstract: Photolithography is typically the bottleneck process in semiconductor manufacturing. In this paper, we present a model for optimizing the scheduling of the photolithography process in the presence of both individual and cluster tools. The combination of these individual and cluster tools that process various layers (stages) of the semiconductor manufacturing process flow is a special type of flexible flowshop. We seek separately to minimize total weighted completion time and maximize on-time delivery performance. Experimental results suggest that our solution algorithms show promise for real world implementation as they can help to improve resource utilization, reduce job completion times, and decrease unnecessary delays in a wafer fab.


title: Thermal Molecular Focusing: Tunable Cross Effect of Phoresis and Advection

abstract: The control of solute fluxes through either microscopic phoresis or hydrodynamic advection is a fundamental way to transport molecules, which are ubiquitously present in nature and technology. We study the transport of large solute such as DNA driven by a time-dependent thermal field in a polymer solution. Heat propagation of a single heat spot moving back and forth gives rise to the molecular focusing of DNA with frequency-tunable control. We developed a theoretical model, where heat conduction, viscoelastic expansion of walls, and the viscosity gradient of a smaller solute are coupled, and that can explain the underlying hydrodynamic focusing and its interplay with phoretic transports. This cross effect may allow one to design a unique miniaturized pump in microfluidics.


title: On Kalman's functor for bounded hemi-implicative semilattices and hemi-implicative lattices

abstract: Hemi-implicative semilattices (lattices), originally defined under the name of weak implicative semilattices (lattices), were introduced by the second author of the present paper. A hemi-implicative semilattice is an algebra $(H,\wedge,\rightarrow,1)$ of type $(2,2,0)$ such that $(H,\wedge)$ is a meet semilattice, $1$ is the greatest element with respect to the order, $a\rightarrow a = 1$ for every $a\in H$ and for every $a$, $b$, $c\in H$, if $a\leq b\rightarrow c$ then $a\wedge b \leq c$. A bounded hemi-implicative semilattice is an algebra $(H,\wedge,\rightarrow,0,1)$ of type $(2,2,0,0)$ such that $(H,\wedge,\rightarrow,1)$ is a hemi-implicative semilattice and $0$ is the first element with respect to the order. A hemi-implicative lattice is an algebra $(H,\wedge,\vee,\rightarrow,0,1)$ of type $(2,2,2,0,0)$ such that $(H,\wedge,\vee,0,1)$ is a bounded distributive lattice and the reduct algebra $(H,\wedge,\rightarrow,1)$ is a hemi-implicative semilattice. In this paper we introduce an equivalence for the categories of bounded hemi-implicative semilattices and hemi-implicative lattices, respectively, which is motivated by an old construction due J. Kalman that relates bounded distributive lattices and Kleene algebras.


title: On-and-off chip cooling of a Coulomb blockade thermometer down to 2.8 mK

abstract: Cooling nanoelectronic devices below 10 mK is a great challenge since thermal conductivities become very small, thus creating a pronounced sensitivity to heat leaks. Here, we overcome these difficulties by using adiabatic demagnetization of \emph{both} the electronic leads \emph{and} the large metallic islands of a Coulomb blockade thermometer. This reduces the external heat leak through the leads and also provides on-chip refrigeration, together cooling the thermometer down to 2.8$\pm$0.1 mK. We present a thermal model which gives a good qualitative account and suggests that the main limitation is heating due to pulse tube vibrations. With better decoupling, temperatures below 1 mK should be within reach, thus opening the door for microkelvin nanoelectronics.


title: Automatically Generating Commit Messages from Diffs using Neural Machine Translation

abstract: Commit messages are a valuable resource in comprehension of software evolution, since they provide a record of changes such as feature additions and bug repairs. Unfortunately, programmers often neglect to write good commit messages. Different techniques have been proposed to help programmers by automatically writing these messages. These techniques are effective at describing what changed, but are often verbose and lack context for understanding the rationale behind a change. In contrast, humans write messages that are short and summarize the high level rationale. In this paper, we adapt Neural Machine Translation (NMT) to automatically "translate" diffs into commit messages. We trained an NMT algorithm using a corpus of diffs and human-written commit messages from the top 1k Github projects. We designed a filter to help ensure that we only trained the algorithm on higher-quality commit messages. Our evaluation uncovered a pattern in which the messages we generate tend to be either very high or very low quality. Therefore, we created a quality-assurance filter to detect cases in which we are unable to produce good messages, and return a warning instead.


title: Fast scrambling in holographic Einstein-Podolsky-Rosen pair

abstract: We demonstrate that a holographic model of the Einstein-Podolsky-Rosen pair exhibits fast scrambling. Strongly entangled quark and antiquark in $\mathcal{N}=4$ super Yang-Mills theory are considered. Their gravity dual is a fundamental string whose endpoints are uniformly accelerated in opposite direction. We slightly increase the acceleration of the endpoint and show that it quickly destroys the correlation between the quark and antiquark. The proper time scale of the destruction is $\tau_\ast\sim \beta \ln S$ where $\beta$ is the inverse Unruh temperature and $S$ is the entropy of the accelerating quark. We also evaluate the Lyapunov exponent from correlation function as $\lambda_L=2\pi/\beta$, which saturates the Lyapunov bound. Our results suggest that the fast scrambling or saturation of the Lyapunov bound do not directly imply the existence of an Einstein dual. When we slightly decrease the acceleration, the quark and antiquark are causally connected and an "one-way traversable wormhole" is created on the worldsheet. It causes the divergence of the correlation function between the quark and antiquark.

categories: hep-th cond-mat.dis-nn gr-qc nlin.CD quant-ph


title: An Exploratory Study of Field Failures

abstract: Field failures, that is, failures caused by faults that escape the testing phase leading to failures in the field, are unavoidable. Improving verification and validation activities before deployment can identify and timely remove many but not all faults, and users may still experience a number of annoying problems while using their software systems. This paper investigates the nature of field failures, to understand to what extent further improving in-house verification and validation activities can reduce the number of failures in the field, and frames the need of new approaches that operate in the field. We report the results of the analysis of the bug reports of five applications belonging to three different ecosystems, propose a taxonomy of field failures, and discuss the reasons why failures belonging to the identified classes cannot be detected at design time but shall be addressed at runtime. We observe that many faults (70%) are intrinsically hard to detect at design-time.


title: Integer sorting on multicores: some (experiments and) observations

abstract: There have been many proposals for sorting integers on multicores/GPUs that include radix-sort and its variants or other approaches that exploit specialized hardware features of a particular multicore architecture. Comparison-based algorithms have also been used. Network-based algorithms have also been used with primary example Batcher's bitonic sorting algorithm. Although such a latter approach is theoretically "inefficient", if there are few keys to sort, it can lead to better running times as it has low overhead and is simple to implement. In this work we perform an experimental study of integer sorting on multicore processors using not only multithreading but also multiprocessing parallel programming approaches. Our implementations work under Open MPI, MulticoreBSP, and BSPlib. We have implemented serial and parallel radix-sort for various radixes and also some previously little explored or unexplored variants of bitonic-sort and odd-even transposition sort. We offer our observations on a performance evaluation using the MBSP model of such algorithm implementations on multiple platforms and architectures and multiple programming libraries. If we can conclude anything is that modeling their performance by taking into consideration architecture dependent features such as the structure and characteristics of multiple memory hierarchies is difficult and more often than not unsuccessful or unreliable. However we can still draw some very simple conclusions using traditional architecture independent parallel modeling.


title: Inferring Narrative Causality between Event Pairs in Films

abstract: To understand narrative, humans draw inferences about the underlying relations between narrative events. Cognitive theories of narrative understanding define these inferences as four different types of causality, that include pairs of events A, B where A physically causes B (X drop, X break), to pairs of events where A causes emotional state B (Y saw X, Y felt fear). Previous work on learning narrative relations from text has either focused on "strict" physical causality, or has been vague about what relation is being learned. This paper learns pairs of causal events from a corpus of film scene descriptions which are action rich and tend to be told in chronological order. We show that event pairs induced using our methods are of high quality and are judged to have a stronger causal relation than event pairs from Rel-grams.


title: Unsupervised Induction of Contingent Event Pairs from Film Scenes

abstract: Human engagement in narrative is partially driven by reasoning about discourse relations between narrative events, and the expectations about what is likely to happen next that results from such reasoning. Researchers in NLP have tackled modeling such expectations from a range of perspectives, including treating it as the inference of the contingent discourse relation, or as a type of common-sense causal reasoning. Our approach is to model likelihood between events by drawing on several of these lines of previous work. We implement and evaluate different unsupervised methods for learning event pairs that are likely to be contingent on one another. We refine event pairs that we learn from a corpus of film scene descriptions utilizing web search counts, and evaluate our results by collecting human judgments of contingency. Our results indicate that the use of web search counts increases the average accuracy of our best method to 85.64% over a baseline of 50%, as compared to an average accuracy of 75.15% without web search.


title: Multi-state effective Hamiltonian and size-consistency corrections in stochastic configuration interactions

abstract: Model space quantum Monte Carlo (MSQMC) is an extension of full configuration interaction QMC (FCIQMC) that allows us to calculate quasi-degenerate and excited electronic states by sampling the effective Hamiltonian in the model space. We introduce a novel algorithm based on the state-selective partitioning for the effective Hamiltonian using left eigenvectors to calculate several electronic states simultaneously at much less computational cost than the original MSQMC with the energy dependent partitioning. The sampling of walkers in MSQMC is analyzed in the single reference limit using a stochastic algorithm for higher-order perturbation energies by the analogy of the deterministic case utilizing a full configuration interaction program. We further develop size-consistency corrections of the initiator adaptation (i-MSQMC) in three different ways, i.e. the coupled electron pair approximation, a posteriori, and second-order pertrubative corrections. It is clearly demonstrated that most of the initiator error is originating from the deficiency of proper scaling of correlation energy due to its truncated CI nature of the initiator approximation, and that the greater part of the error can be recovered by the size-consistency corrections developed in this work.


title: Crystallization in melts and poor-solvent solutions of semiflexible polymers: extensive DPD study

abstract: In the present work, crystallization in melts and poor-solvent solutions of semiflexible polymers with different concentration was studied by means of dissipative particle dynamics simulation technique. We use a coarse-grained polymer model trying to catch general principles of crystallization in such systems on large time and length scales. We observe the crystallization process starting from an initial randomly prepared system with different polymer volume fractions in a poor solvent. Because the solvent is very poor, the macrophase polymer-solvent separation takes place very fast and is accompanied by partial polymer crystallization. We have found that the overall crystalline fraction at the end of crystallization process decreases upon increasing the polymer volume fraction in the initial randomly prepared system, while the steady-state crystallization speed is almost the same at polymer volume fractions larger than 50\%. At the same time, the average crystallite size differs considerably and has a maximum value in the systems with 90\% polymer volume fraction. We assume that this polymer concentration is an optimal value in a sense of a balance between the amount of polymer material available for increasing crystallite size and chain entanglements preventing crystallites growth and merging.


title: Detection of vertical muons with the HAWC water Cherenkov detectors and its application to gamma/hadron discrimination

abstract: The HAWC observatory reconstructs atmospheric showers induced by very high-energy primary cosmic- or gamma-rays. Cosmic-rays with hadronic nature are several orders of magnitude more frequent than primary gammas and therefore it is necessary to perform the best gamma-hadron discrimination for each event to be able to survey the gamma-ray sky. Single vertical muons cause very characteristic signals in the 300 HAWC detectors. Particularly, those events that produce signals in the four PMTs of a tank corresponding to muons traversing the tank in the central region. In such cases exist the correlation that the largest signal is also the earliest to arrive, this due to the geometry of the Cherenkov light emission and the fact that relativistic muons travel faster than the light in the medium. This pattern is different from that of electrons, positrons and low energy gammas showering in the first meter of the tank water, which is typical of electromagnetic showers. In this work, we explored whether in air-showers it is possible to identify muon signals and whether this can help in the gamma/hadron discrimination of the recorded events.


title: Efficient computation of two-dimensional steady free-surface flows

abstract: We consider a family of steady free-surface flow problems in two dimensions, concentrating on the effect of nonlinearity on the train of gravity waves that appear downstream of a disturbance. By exploiting standard complex variable techniques, these problems are formulated in terms of a coupled system of Bernoulli's equation and an integral equation. When applying a numerical collocation scheme, the Jacobian for the system is dense, as the integral equation forces each of the algebraic equations to depend on each of the unknowns. We present here a strategy for overcoming this challenge, which leads to a numerical scheme that is much more efficient than what is normally employed for these types of problems, allowing for many more grid points over the free surface. In particular, we provide a simple recipe for constructing a sparse approximation to the Jacobian that is used as a preconditioner in a Jacobian-free Newton-Krylov method for solving the nonlinear system. We use this approach to compute numerical results for a variety of prototype problems including flows past pressure distributions, a surface-piercing object and bottom topographies.


title: Finite Sample Inference for Targeted Learning

abstract: The Highly-Adaptive-Lasso(HAL)-TMLE is an efficient estimator of a pathwise differentiable parameter in a statistical model that at minimal (and possibly only) assumes that the sectional variation norm of the true nuisance parameters are finite. It relies on an initial estimator (HAL-MLE) of the nuisance parameters by minimizing the empirical risk over the parameter space under the constraint that sectional variation norm is bounded by a constant, where this constant can be selected with cross-validation. In the formulation of the HALMLE this sectional variation norm corresponds with the sum of absolute value of coefficients for an indicator basis. Due to its reliance on machine learning, statistical inference for the TMLE has been based on its normal limit distribution, thereby potentially ignoring a large second order remainder in finite samples. In this article, we present four methods for construction of a finite sample 0.95-confidence interval that use the nonparametric bootstrap to estimate the finite sample distribution of the HAL-TMLE or a conservative distribution dominating the true finite sample distribution. We prove that it consistently estimates the optimal normal limit distribution, while its approximation error is driven by the performance of the bootstrap for a well behaved empirical process. We demonstrate our general inferential methods for 1) nonparametric estimation of the average treatment effect based on observing on each unit a covariate vector, binary treatment, and outcome, and for 2) nonparametric estimation of the integral of the square of the multivariate density of the data distribution.


title: Outflows of hot molecular gas in ultra-luminous infra-red galaxies mapped with VLT-SINFONI

abstract: We present the detection and morphological characterization of hot molecular gas outflows in nearby ultra-luminous infrared galaxies, using the near-IR integral-field spectrograph SINFONI on the VLT. We detect outflows observed in the 2.12 micron H$_{2}$ 1-0 S(1) line for three out of four ULIRGs analyzed; IRAS 12112+0305, 14348-1447, and 22491-1808. The outflows are mapped on scales of 0.7-1.6 kpc, show typical outflow velocities of 300-500 km/s, and appear to originate from the nuclear region. The outflows comprise hot molecular gas masses of ~6-8x10$^3$ M(sun). Assuming a hot-to-cold molecular gas mass ratio of 6x10$^{-5}$, as found in nearby luminous IR galaxies, the total (hot+cold) molecular gas mass in these outflows is expected to be ~1x10$^{8}$ M(sun). This translates into molecular mass outflow rates of ~30-85 M(sun)/yr, which is a factor of a few lower than the star formation rate in these ULIRGs. In addition, most of the outflowing molecular gas does not reach the escape velocity of these merger systems, which implies that the bulk of the outflowing molecular gas is re-distributed within the system and thus remains available for future star formation. The fastest H$_{2}$ outflow is seen in the Compton-thick AGN of IRAS 14348-1447, reaching a maximum outflow velocity of ~900 km/s. Another ULIRG, IRAS 17208-0014, shows asymmetric H$_{2}$ line profiles different from the outflows seen in the other three ULIRGs. We discuss several alternative explanations for its line asymmetries, including a very gentle galactic wind, internal gas dynamics, low-velocity gas outside the disk, or two superposed gas disks. We do not detect the hot molecular counterpart to the outflow previously detected in CO(2-1) in IRAS 17208-0014, but we note that our SINFONI data are not sensitive enough to detect this outflow if it has a small hot-to-cold molecular gas mass ratio of < 9x10$^{-6}$.


title: NEOWISE Reactivation Mission Year Three: Asteroid Diameters and Albedos

abstract: The Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE) reactivation mission has completed its third year of surveying the sky in the thermal infrared for near-Earth asteroids and comets. NEOWISE collects simultaneous observations at 3.4 um and 4.6 um of solar system objects passing through its field of regard. These data allow for the determination of total thermal emission from bodies in the inner solar system, and thus the sizes of these objects. In this paper we present thermal model fits of asteroid diameters for 170 NEOs and 6110 MBAs detected during the third year of the survey, as well as the associated optical geometric albedos. We compare our results with previous thermal model results from NEOWISE for overlapping sample sets, as well as diameters determined through other independent methods, and find that our diameter measurements for NEOs agree to within 26% (1-sigma) of previously measured values. Diameters for the MBAs are within 17% (1-sigma). This brings the total number of unique near-Earth objects characterized by the NEOWISE survey to 541, surpassing the number observed during the fully cryogenic mission in 2010.


title: Closed Timelike Curves and "Effective" Superluminal Travel with Naked Line Singularities

abstract: We examine closed timelike curves (CTCs) and "effective" superluminal travel in a spacetime containing naked line singularities, which we call "wires". Each wire may be straight-line singularity or a ring singularity. The Weak Energy Condition (WEC) is preserved in all well-defined regions of the spacetime. (The singularities themselves are not well-defined, so the WEC is undefined there, but it is never explicitly violated.) Parallel to the wire, "effective" superluminal travel is possible, in that the wire may be used as a shortcut between distant regions of spacetime. Our purpose in presenting the superluminal aspects of the wire is to dispel the commonly held view that explicit WEC violation is necessary for effective superluminal travel, whereas in truth the strictures against superluminal travel are more complicated. We also demonstrate how the existence of such "wires" could create CTCs. We present a model spacetime which contains two wires which are free to move relative to each other. This spacetime is asymptotically flat: It becomes a Minkowski spacetime a finite distance away from each of the wires. The CTCs under investigation do not need to enter the wires' singularities, and can be confined to regions that are weak-field: This means that if these wires were physically possible, they would present causality problems even in nonsingular, energetically realistic regions of the spacetime. We conclude that the Weak Energy Condition alone is not sufficient to prevent superluminal travel in asymptotically flat spacetimes.


title: Affine map equivalence versus critical set equivalence for quadratic maps of the plane

abstract: In recent work [Nien et al. 2016], the authors enumerated a classification of quadratic maps of the plane according to their critical sets and images. It is straightforward to show that quadratic maps which are affinely map equivalent are also equivalent in the critical set classification. The question remained whether maps that are equivalent in the critical set classification are also affinely map equivalent. This paper establishes a complete enumeration of the affine map equivalence classes. As a consequence, the relationship between affine map equivalence and critical set equivalence is established. In short, there are eighteen affine map equivalence classes. Three pairs of those classes have critical sets and images that match, but each pair has some other geometric property, preserved by affine map equivalence, that does not match. The other twelve affine map equivalence classes and the critical set equivalences are in one-to-one correspondence.


title: Estimation in Semiparametric Quantile Factor Models

abstract: We propose an estimation methodology for a semiparametric quantile factor panel model. We provide tools for inference that are robust to the existence of moments and to the form of weak cross-sectional dependence in the idiosyncratic error term. We apply our method to daily stock return data.


title: Frequency Vectoralization and Frequency Birefringence

abstract: In view of momentum continuity at a temporal slab, it is shown that instantaneous switching of an isotropic medium to an anisotropic medium offers the incident frequency a directional property- a counterintuitive process which is called frequency vectoralization. By expressing the dispersion diagrams before and after the temporal transition, a general analogy between spatial and temporal interfaces is given and the concept of frequency birefringence, i.e. double frequency jump, will be explained. Furthermore, it will be shown that an anisotropic temporal interface diffracts a monochromatic beam in both the spatial and spectral domains.


title: Resistive MHD Modelling of Quasi-Single Helicity State in the KTX Regimes

abstract: The potential formation of the quasi-single-helicity (QSH) state in the Keda Torus eXperiment (KTX) is investigated in resistive MHD simulations using the NIMROD code. We focus on the effects of finite resistivity on the mode structure and characteristics of the dominant linear and nonlinear resistive tearing-mode in a finite $\beta$, cylindrical configuration of reversed field pinch model for KTX. In the typical resistive regimes of KTX where Lundquist number $S=5 \times 10^4$, the plasma transitions to a steady QSH state after evolving through an initial transient phase with multiple helicities. The dominant mode of the QSH state develops from the dominant linear tearing mode instability. In lower $\beta$ regime, the QSH state are intermittent and short in duration; in higher $\beta$ regime, the QSH state persists for a longer time and should be more observable in experiment.


title: Desorption of hydrocarbon chains by association with ionic and nonionic surfactants under flow as a mechanism for enhanced oil recovery

abstract: The need to extract oil from wells where it is embedded on the surfaces of rocks has led to the development of new and improved enhanced oil recovery techniques. One of those is the injection of surfactants with water vapor, which promotes desorption of oil that can then be extracted using pumps, as the surfactants encapsulate the oil in foams. However, the mechanisms that lead to the optimal desorption of oil and the best type of surfactants to carry out desorption are not well known yet, which warrants the need to carry out basic research on this topic. In this work, we report non equilibrium dissipative particle dynamics simulations of model surfactants and oil molecules adsorbed on surfaces, with the purpose of studying the efficiency of the surfactants to desorb hydrocarbon chains, that are found adsorbed over flat surfaces. The model surfactants studied correspond to nonionic and cationic surfactants, and the hydrocarbon desorption is studied as a function of surfactant concentration under increasing Poiseuille flow. We obtain various hydrocarbon desorption isotherms for every model of surfactant proposed, under flow. Nonionic surfactants are found to be the most effective to desorb oil and the mechanisms that lead to this phenomenon are presented and discussed.


title: Temperature effects on nuclear pseudospin symmetry in the Dirac-Hartree-Bogoliubov formalism

abstract: We present finite temperature Dirac-Hartree-Bogoliubov (FTDHB) calculations for the tin isotope chain to study the dependence of pseudospin on the nuclear temperature. In the FTDHB calculation, the density dependence of the self-consistent relativistic mean fields, the pairing, and the vapor phase that takes into account the unbound nucleon states are considered self-consistently. The mean field potentials obtained in the FTDHB calculations are fit by Woods-Saxon (WS) potentials to examine how the WS parameters are related to the energy splitting of the pseudospin pairs as the temperature increases. We find that the nuclear potential surface diffuseness is the main driver for the pseudospin splittings and that it increases as the temperature grows. We conclude that pseudospin symmetry is better realized when the nuclear temperature increases. The results confirm the findings of previous works using RMF theory at $T=0$, namely that the correlation between the pseudospin splitting and the parameters of the Woods-Saxon potentials implies that pseudospin symmetry is a dynamical symmetry in nuclei. We show that the dynamical nature of the pseudospin symmetry remains when the temperature is considered in a realistic calculation of the tin isotopes, such as that of the Dirac-Hartree-Bogoliubov formalism.


title: On the error rate of conditional quasi-Monte Carlo for discontinuous functions

abstract: This paper studies the rate of convergence for conditional quasi-Monte Carlo (QMC), which is a counterpart of conditional Monte Carlo. We focus on discontinuous integrands defined on the whole of $R^d$, which can be unbounded. Under suitable conditions, we show that conditional QMC not only has the smoothing effect (up to infinitely times differentiable), but also can bring orders of magnitude reduction in integration error compared to plain QMC. Particularly, for some typical problems in options pricing and Greeks estimation, conditional randomized QMC that uses $n$ samples yields a mean error of $O(n^{-1+\epsilon})$ for arbitrarily small $\epsilon>0$. As a by-product, we find that this rate also applies to randomized QMC integration with all terms of the ANOVA decomposition of the discontinuous integrand, except the one of highest order.


title: Applying Generic Coding with Help to Uniformizations

abstract: This is a follow up to a paper by the author where the disjointness relation for (the graphs of) definable functions from ${^\omega \omega}$ to ${^\omega \omega}$ is analyzed. In that paper, for each $a \in {^\omega \omega}$ we defined a Baire class one function $f_a^{GC} : {^\omega \omega} \to {^\omega \omega}$ which encoded $a$ in a certain sense. Given $g : {^\omega \omega} \to {^\omega \omega}$, let $\Psi(g)$ be the statement that $g$ is disjoint from at most countably many of the functions $f_a^{GC}$. We show the consistency strength of $(\forall g)\, \Psi(g)$ is at most one inaccessible cardinal. We show that $\mbox{AD}^+$ implies $(\forall g)\, \Psi(g)$. Finally, we show that assuming large cardinals, $(\forall g)\, \Psi(g)$ holds in models of the form $L(\mathbb{R} [\mathcal{U}]$ where $\mathcal{U}$ is a selective ultrafilter on $\omega$.


title: A hypothetical effect of the Maxwell-Proca electromagnetic stresses on galaxy rotation curves

abstract: The Maxwell-Proca electrodynamics corresponding to a finite photon mass causes a substantial change of the Maxwell stress tensor and, under certain circumstances, may cause the electromagnetic stresses to act effectively as "negative pressure." The paper describes a model where this negative pressure imitates gravitational pull and may produce forces comparable to gravity and even become dominant. The effect is associated with the random magnetic fields in the galactic disk with a scale exceeding the photon Compton wavelength. The presence of a weaker regular field does not affect the forces under consideration. The stresses act predominantly on the interstellar gas and cause an additional force pulling the gas towards the center and towards the galactic plane. The stars do not experience any significant direct force but get involved in this process via a "recycling loop" where rapidly evolving massive stars are formed from the gas undergoing galactic rotation and then lose their masses back to the gas within a time shorter than roughly 1/6 of the rotation period. This makes their dynamics inseparable from that of the rotating gas. The lighter, slowly evolving stars, as soon as they are formed, lose connection to the gas and are confined within the galaxy only gravitationally. Numerical examples based on the parameters of our galaxy reveal both opportunities and challenges of this model and motivate further analysis. The critical issue is the plausibility of formation of the irregular magnetic field that would be force free. Another challenge is developing a predictive model of the evolution of the gaseous and stellar population of the galaxy under the aforementioned scenario. It may be interesting to also explore possible broader cosmological implications of the negative-pressure model.


title: On Upward Drawings of Trees on a Given Grid

abstract: Computing a minimum-area planar straight-line drawing of a graph is known to be NP-hard for planar graphs, even when restricted to outerplanar graphs. However, the complexity question is open for trees. Only a few hardness results are known for straight-line drawings of trees under various restrictions such as edge length or slope constraints. On the other hand, there exist polynomial-time algorithms for computing minimum-width (resp., minimum-height) upward drawings of trees, where the height (resp., width) is unbounded. In this paper we take a major step in understanding the complexity of the area minimization problem for strictly-upward drawings of trees, which is one of the most common styles for drawing rooted trees. We prove that given a rooted tree $T$ and a $W\times H$ grid, it is NP-hard to decide whether $T$ admits a strictly-upward (unordered) drawing in the given grid.


title: Leveraging Deep Neural Network Activation Entropy to cope with Unseen Data in Speech Recognition

abstract: Unseen data conditions can inflict serious performance degradation on systems relying on supervised machine learning algorithms. Because data can often be unseen, and because traditional machine learning algorithms are trained in a supervised manner, unsupervised adaptation techniques must be used to adapt the model to the unseen data conditions. However, unsupervised adaptation is often challenging, as one must generate some hypothesis given a model and then use that hypothesis to bootstrap the model to the unseen data conditions. Unfortunately, reliability of such hypotheses is often poor, given the mismatch between the training and testing datasets. In such cases, a model hypothesis confidence measure enables performing data selection for the model adaptation. Underlying this approach is the fact that for unseen data conditions, data variability is introduced to the model, which the model propagates to its output decision, impacting decision reliability. In a fully connected network, this data variability is propagated as distortions from one layer to the next. This work aims to estimate the propagation of such distortion in the form of network activation entropy, which is measured over a short- time running window on the activation from each neuron of a given hidden layer, and these measurements are then used to compute summary entropy. This work demonstrates that such an entropy measure can help to select data for unsupervised model adaptation, resulting in performance gains in speech recognition tasks. Results from standard benchmark speech recognition tasks show that the proposed approach can alleviate the performance degradation experienced under unseen data conditions by iteratively adapting the model to the unseen datas acoustic condition.


title: Upper and Lower Bounds on the Capacity of Amplitude-Constrained MIMO Channels

abstract: In this work, novel upper and lower bounds for the capacity of channels with arbitrary constraints on the support of the channel input symbols are derived. As an immediate practical application, the case of multiple-input multiple-output channels with amplitude constraints is considered. The bounds are shown to be within a constant gap if the channel matrix is invertible and are tight in the high amplitude regime for arbitrary channel matrices. Moreover, in the high amplitude regime, it is shown that the capacity scales linearly with the minimum between the number of transmit and receive antennas, similarly to the case of average power-constrained inputs.


title: The reactor antineutrino anomaly and low energy threshold neutrino experiments

abstract: Short distance reactor antineutrino experiments measure an antineutrino spectrum a few percent lower than expected from theoretical predictions. In this work we study the potential of low energy threshold reactor experiments in the context of a light sterile neutrino signal. We discuss the perspectives of the recently detected coherent elastic neutrino-nucleus scattering in future reactor antineutrino experiments. We find that the expectations to improve the current constraints on the mixing with sterile neutrinos are promising. We also analyse the measurements of antineutrino scattering off electrons from short distance reactor experiments. In this case, the statistics is not competitive with inverse beta decay experiments, although future experiments might play a role when compare it with the Gallium anomaly.


title: On Dissipative Nonlinear Evolutional Pseudo-Differential Equations

abstract: First, using the uniform decomposition in both physical and frequency spaces, we obtain an equivalent norm on modulation spaces. Secondly, we consider the Cauchy problem for the dissipative evolutionary pseudo-differential equation \partial_t u + A(x,D) u = F\big((\partial^\alpha_x u)_{|\alpha|\leq \kappa}\big), \ \ u(0,x)= u_0(x), where $A(x,D)$ is a dissipative pseudo-differential operator and $F(z)$ is a multi-polynomial. We will develop the uniform decomposition techniques in both physical and frequency spaces to study its local well posedness in modulation spaces $M^s_{p,q}$ and in Sobolev spaces $H^s$. Moreover, the local solution can be extended to a global one in $L^2$ and in $H^s$ ($s>\kappa+d/2$) for certain nonlinearities.


title: High-Frequency Jump Tests: Which Test Should We Use?

abstract: We conduct an extensive evaluation of price jump tests based on high-frequency financial data. After providing a concise review of multiple alternative tests, we document the size and power of all tests in a range of empirically relevant scenarios. Particular focus is given to the robustness of test performance to the presence of jumps in volatility and microstructure noise, and to the impact of sampling frequency. The paper concludes by providing guidelines for empirical researchers about which test to choose in any given setting.


title: Helical Majorana fermions at the interface of Weyl semimetal and d-wave superconductor: Application to Iridates and high-Tc Cuprates

abstract: Majorana bound states exist inside a core of half-quantum vortex in spin-triplet superconductors. Despite intense efforts, they have not been discovered in any spin-triplet superconductor candidate materials. After the success of topological insulator, another route to achieve Majorana fermions has been suggested: heterostructures of s-wave superconductors and topological insulator or semimetal with strong spin-orbit coupling provide an effective spinless $p+ip$ pairing which supports Majorana bound states in a single vortex. This theoretical observation has stimulated both experimental and theoretical communities to search for Majorana fermions, and recently a localized Majorana state at the end of one-dimensional wire has been reported. Here we study the two-dimensional interface of Weyl semimetal and d-wave superconductor which promotes a pair of helical Majorana fermions propagating along the edge of the interface. We suggest that Iridium oxide layer IrO$_2$ classified as two-dimensional Weyl semimetal in close proximity to high temperature Cuprates would be a best example to explore these helical Majorana modes.


title: Action Classification and Highlighting in Videos

abstract: Inspired by recent advances in neural machine translation, that jointly align and translate using encoder-decoder networks equipped with attention, we propose an attentionbased LSTM model for human activity recognition. Our model jointly learns to classify actions and highlight frames associated with the action, by attending to salient visual information through a jointly learned soft-attention networks. We explore attention informed by various forms of visual semantic features, including those encoding actions, objects and scenes. We qualitatively show that soft-attention can learn to effectively attend to important objects and scene information correlated with specific human actions. Further, we show that, quantitatively, our attention-based LSTM outperforms the vanilla LSTM and CNN models used by stateof-the-art methods. On a large-scale youtube video dataset, ActivityNet, our model outperforms competing methods in action classification.


title: Period mappings and properties of the augmented Hodge line bundle

abstract: Let $P$ be the image of a period map. We discuss progress towards a conjectural Hodge theoretic completion $\overline{P}$, an analogue of the Satake-Baily-Borel compactification in the classical case. The set $\overline{P}$ is defined and given the structure of a compact Hausdorff topological space. We conjecture that it admits the structure of a compact complex analytic variety. We verify this conjecture when $\mathrm{dim} P \le 2$. In general, $\overline{P}$ admits a finite cover $\overline{S}$ (also a compact Hausdorff space, and constructed from Stein factorizations of period maps). Assuming that $\overline{S}$ is a compact complex analytic variety, we show that a lift of the augmented Hodge line bundle $\Lambda$ extends to an ample line bundle, giving $\overline{P}$ the structure of a projective normal variety. Our arguments rely on refined positivity properties of Chern forms associated to various Hodge bundles; properties that might be of independent interest.


title: Thermodynamics of Radiation Pressure and Photon Momentum

abstract: Theoretical analyses of radiation pressure and photon momentum in the past 150 years have focused almost exclusively on classical and/or quantum theories of electrodynamics. In these analyses, Maxwell's equations, the properties of polarizable and/or magnetizable material media, and the stress tensors of Maxwell, Abraham, Minkowski, Chu, and Einstein-Laub have typically played prominent roles. Each stress tensor has subsequently been manipulated to yield its own expressions for the electromagnetic (EM) force, torque, energy, and linear as well as angular momentum densities of the EM field. This paper presents an alternative view of radiation pressure from the perspective of thermal physics, invoking the properties of blackbody radiation in conjunction with empty as well as gas-filled cavities that contain EM energy in thermal equilibrium with the container's walls. In this type of analysis, Planck's quantum hypothesis, the spectral distribution of the trapped radiation, the entropy of the photon gas, and Einstein's A and B coefficients play central roles.


title: Decompositions of amplituhedra

abstract: The (tree) amplituhedron A(n,k,m) is the image in the Grassmannian Gr(k,k+m) of the totally nonnegative part of Gr(k,n), under a (map induced by a) linear map which is totally positive. It was introduced by Arkani-Hamed and Trnka in 2013 in order to give a geometric basis for the computation of scattering amplitudes in N=4 supersymmetric Yang-Mills theory. In the case relevant to physics (m=4), there is a collection of recursively-defined 4k-dimensional BCFW cells in the totally nonnegative part of Gr(k,n), whose images conjecturally "triangulate" the amplituhedron--that is, their images are disjoint and cover a dense subset of A(n,k,4). In this paper, we approach this problem by first giving an explicit (as opposed to recursive) description of the BCFW cells. We then develop sign-variational tools which we use to prove that when k=2, the images of these cells are disjoint in A(n,k,4). We also conjecture that for arbitrary even m, there is a decomposition of the amplituhedron A(n,k,m) involving precisely M(k, n-k-m, m/2) top-dimensional cells (of dimension km), where M(a,b,c) is the number of plane partitions contained in an a x b x c box. This agrees with the fact that when m=4, the number of BCFW cells is the Narayana number N(n-3, k+1).


title: Non-invasive acquisition of fetal ECG from the maternal xyphoid process: a feasibility study in pregnant sheep and a call for open data sets

abstract: Objective: The utility of fetal heart rate (FHR) monitoring can only be achieved with an acquisition sampling rate that preserves the underlying physiological information on the millisecond time scale (1000 Hz rather than 4 Hz). For such acquisition, fetal ECG (fECG) is required, rather than the ultrasound to derive FHR. We tested one recently developed algorithm, SAVER, and two widely applied algorithms to extract fECG from a single channel maternal ECG signal recorded over the xyphoid process rather than the routine abdominal signal. Approach: At 126dG, ECG was attached to near-term ewe and fetal shoulders, manubrium and xyphoid processes (n=12). FECG served as the ground-truth to which the fetal ECG signal extracted from the simultaneously-acquired maternal ECG was compared. All fetuses were in good health during surgery (pH 7.29+/-0.03, pO2 33.2+/-8.4, pCO2 56.0+/-7.8, O2Sat 78.3+/-7.6, lactate 2.8+/-0.6, BE -0.3+/-2.4). Main result: In all animals, single lead fECG extraction algorithm could not extract fECG from the maternal ECG signal over the xyphoid process with the F1 less than 50%. Significance: The applied fECG extraction algorithms might be unsuitable for the maternal ECG signal over the xyphoid process, or the latter does not contain strong enough fECG signal, although the lead is near the mother's abdomen. Fetal sheep model is widely used to mimic various fetal conditions, yet ECG recordings in a public data set form are not available to test the predictive ability of fECG and FHR. We are making this data set openly available to other researchers to foster non-invasive fECG acquisition in this animal model.


title: Ap\'ery sets of shifted numerical monoids

abstract: A numerical monoid is an additive submonoid of the non-negative integers. Given a numerical monoid $S$, consider the family of "shifted" monoids $M_n$ obtained by adding $n$ to each generator of $S$. In this paper, we characterize the Ap\'ery set of $M_n$ in terms of the Ap\'ery set of the base monoid $S$ when $n$ is sufficiently large. We give a highly efficient algorithm for computing the Ap\'ery set of $M_n$ in this case, and prove that several numerical monoid invariants, such as the genus and Frobenius number, are eventually quasipolynomial as a function of $n$.


title: Optical and Near-Infrared Polarimetry of Non-Periodic Comet C/2013 US10 (Catalina)

abstract: We present an optical and near-infrared (hereafter NIR) polarimetric study of a comet C/2013 US10 (Catalina) observed on UT 2015 December 17-18 at phase angles of $\alpha$=52.1 deg - 53.1 deg. Additionally, we obtained an optical spectrum and multi-band images to examine the influence of gas emission. We find that the observed optical signals are significantly influenced by gas emission, that is, the gas-to-total intensity ratio varies from 5 to 30 % in the $R_{\rm C}$ and 3 to 18 % in the $I_{\rm C}$ bands, depending on the position in the coma. We derive the `gas-free dust polarization degrees' of 13.8$\pm$1.0 % in the $R_{\rm C}$ and 12.5$\pm$1.1 % in the $I_{\rm C}$ bands and a gray polarimetric color, i.e., -8.7$\pm$9.9 % $\mu \mathrm{m}$$^{-1}$ in optical and 1.6$\pm$0.9 % $\mu \mathrm{m}$$^{-1}$ in NIR. The increments of polarization obtained from the gas correction show that the polarimetric properties of the dust in this low-polarization comet are not different from those in high-polarization comets. In this process, the cometocentric distance dependence of polarization has disappeared. We also find that the $R_{\rm C}$-band polarization degree of the southeast dust tail, which consists of large dust particles (100 $\mu \mathrm{m}$ - 1 mm), is similar to that in the outer coma where small and large ones are mixed. Our study confirms that the dichotomy of cometary polarization does not result from the difference of dust properties, but from depolarizing gas contamination. This conclusion can provide a strong support for similarity in origin of comets.


title: The global mass functions of 35 Galactic globular clusters: I. Observational data and correlations with cluster parameters

abstract: We have derived the global mass functions of a sample of 35 Galactic globular clusters by comparing deep Hubble Space Telescope photometry with suitable multimass dynamical models. For a subset of 29 clusters with available radial velocity information we were also able to determine dynamical parameters, mass-to-light ratios and the mass fraction of dark remnants. The derived global mass functions are well described by single power-laws in the mass range $0.2 < m/M_\odot < 0.8$ with mass function slopes $\alpha>-1$. Less evolved clusters show deviations from a single-power law, indicating that the original shape of their mass distribution was not a power-law. We find a tight anticorrelation between the present-day mass function slopes and the half-mass relaxation times, which can be understood if clusters started from the same universal IMF and internal dynamical evolution is the main driver in shaping the present-day mass functions. Alternatively, IMF differences correlated with the present-day half-mass relaxation time are needed to explain the observed correlation. The large range of mass function slopes seen for our clusters implies that most globular clusters are dynamically highly evolved, a fact that seems difficult to reconcile with standard estimates for the dynamical evolution of clusters. The mass function slopes also correlate with the dark remnant fractions indicating a preferential retention of massive remnants in clusters subject to high mass-loss rates.


title: The Global Mass Functions of 35 Galactic globular clusters: II. Clues on the Initial Mass Function and Black Hole Retention Fraction

abstract: In this paper we compare the mass function slopes of Galactic globular clusters recently determined by Sollima & Baumgardt (2017) with a set of dedicated N-body simulations of star clusters containing between 65,000 to 200,000 stars. We study clusters starting with a range of initial mass functions (IMFs), black hole retention fractions and orbital parameters in the parent galaxy. We find that the present-day mass functions of globular clusters agree well with those expected for star clusters starting with Kroupa or Chabrier IMFs, and are incompatible with clusters starting with single power-law mass functions for the low-mass stars. The amount of mass segregation seen in the globular clusters studied by Sollima & Baumgardt (2017) can be fully explained by two-body relaxation driven mass segregation from initially unsegregated star clusters. Based on the present-day global mass functions, we expect that a typical globular cluster in our sample has lost about 75% of its mass since formation, while the most evolved clusters have already lost more than 90% of their initial mass and should dissolve within the next 1 to 2 Gyr. Most clusters studied by Sollima & Baumgardt also show a large difference between their central and global MF slopes, implying that the majority of Galactic globular clusters is either near or already past core collapse. The strong mass segregation seen in most clusters also implies that only a small fraction of all black holes formed in globular clusters still reside in them.


title: Influence of subgrid scale models on the buffer sublayer in channel flow

abstract: Subgrid-scale (SGS) models are critical in large-eddy simulations (LES) of turbulent flows. In this paper we conduct a comparative study on different SGS models, including one-k-equation, wall-adapting local eddy-viscosity (WALE), Sigma and shear-constrained model. Wall-resolved LES simulations of channel flows are performed with a finite volume code at shear Reynolds number Re_\tau = 395. In the simulations, the buffer sublayer turns out to be the most sensitive to the SGS model. Through the analysis of the mean velocity, the second-order moments, the SGS viscosity and the fluctuating vorticities, it is shown that the WALE and Sigma model outperform others significantly in terms of fluctuations while the constrained model improves slightly the mean velocity. The results also indicate that the SGS dissipation influences strongly the velocity fluctuations but not the mean flow nor the log-layer mismatch.


title: Leveraging local h-index to identify and rank influential spreaders in networks

abstract: Identifying influential nodes in complex networks has received increasing attention for its great theoretical and practical applications in many fields. Traditional methods, such as degree centrality, betweenness centrality, closeness centrality, and coreness centrality, have more or less disadvantages in detecting influential nodes, which have been illustrated in related literatures. Recently, the h-index, which is utilized to measure both the productivity and citation impact of the publications of a scientist or scholar, has been introduced to the network world to evaluate a node's spreading ability. However, this method assigns too many nodes with the same value, which leads to a resolution limit problem in distinguishing the real influence of these nodes. In this paper, we propose a local h-index centrality (LH-index) method for identifying and ranking influential nodes in networks. The LH-index method simultaneously takes into account of h-index values of the node itself and its neighbors, which is based on the idea that a node connects to more influential nodes will also be influential. According to the simulation results with the stochastic Susceptible-Infected-Recovered (SIR) model in four real world networks and several simulated networks, we demonstrate the effectivity of the LH-index method in identifying influential nodes in networks.


title: Improved ArtGAN for Conditional Synthesis of Natural Image and Artwork

abstract: This paper proposes a series of new approaches to improve Generative Adversarial Network (GAN) for conditional image synthesis and we name the proposed model as ArtGAN. One of the key innovation of ArtGAN is that, the gradient of the loss function w.r.t. the label (randomly assigned to each generated image) is back-propagated from the categorical discriminator to the generator. With the feedback from the label information, the generator is able to learn more efficiently and generate image with better quality. Inspired by recent works, an autoencoder is incorporated into the categorical discriminator for additional complementary information. Last but not least, we introduce a novel strategy to improve the image quality. In the experiments, we evaluate ArtGAN on CIFAR-10 and STL-10 via ablation studies. The empirical results showed that our proposed model outperforms the state-of-the-art results on CIFAR-10 in terms of Inception score. Qualitatively, we demonstrate that ArtGAN is able to generate plausible-looking images on Oxford-102 and CUB-200, as well as able to draw realistic artworks based on style, artist, and genre. The source code and models are available at: https://github.com/cs-chan/ArtGAN


title: Towards On-Chip Optical FFTs for Convolutional Neural Networks

abstract: Convolutional neural networks have become an essential element of spatial deep learning systems. In the prevailing architecture, the convolution operation is performed with Fast Fourier Transforms (FFT) electronically in GPUs. The parallelism of GPUs provides an efficiency over CPUs, however both approaches being electronic are bound by the speed and power limits of the interconnect delay inside the circuits. Here we present a silicon photonics based architecture for convolutional neural networks that harnesses the phase property of light to perform FFTs efficiently. Our all-optical FFT is based on nested Mach-Zender Interferometers, directional couplers, and phase shifters, with backend electro-optic modulators for sampling. The FFT delay depends only on the propagation delay of the optical signal through the silicon photonics structures. Designing and analyzing the performance of a convolutional neural network deployed with our on-chip optical FFT, we find dramatic improvements by up to 10^4 when compared to state-of-the-art GPUs when exploring a compounded figure-of-merit given by power per convolution over area. At a high level, this performance is enabled by mapping the desired mathematical function, an FFT, synergistically onto hardware, in this case optical delay interferometers.


title: A secure blind watermarking scheme based on DCT domain of the scrambled image

abstract: This paper investigates a secure blind watermarking scheme. The main idea of the scheme not only protects the watermark information but also the embedding positions. To achieve a higher level of security, we propose a sub key generation mechanism based on the singular value decomposition and hash function, where sub keys depend on both the main key and the feature codes of the original image. The different sub keys ensure that the embedding positions are randomly selected for different original images. Watermark is embedded in the Discrete Cosine Transform (DCT) coefficients of the scrambled original image. Simulation results show that such embedded method resolves well the contradiction of imperceptibility and robustness. Based on good correlation properties of chaotic sequences, we design a detection method, which can accurately compute geometric transformation (rotation and translation transformations) parameters. The security analysis, including key space analysis, key sensitivity analysis, cryptanalysis, and the comparison results demonstrate that the proposed watermarking scheme also achieves high security.


title: Localisation property of Battle-Lemari\'{e} wavelets' sums

abstract: Explicit formulae are given for a type of Battle-Lemari\'{e} scaling functions and related wavelets. Compactly supported sums of their translations are established and applied to alternative norm characterization of sequence spaces isometrically isomorphic to Nikolskii-Besov spaces on R.


title: DolphinAtack: Inaudible Voice Commands

abstract: Speech recognition (SR) systems such as Siri or Google Now have become an increasingly popular human-computer interaction method, and have turned various systems into voice controllable systems(VCS). Prior work on attacking VCS shows that the hidden voice commands that are incomprehensible to people can control the systems. Hidden voice commands, though hidden, are nonetheless audible. In this work, we design a completely inaudible attack, DolphinAttack, that modulates voice commands on ultrasonic carriers (e.g., f > 20 kHz) to achieve inaudibility. By leveraging the nonlinearity of the microphone circuits, the modulated low frequency audio commands can be successfully demodulated, recovered, and more importantly interpreted by the speech recognition systems. We validate DolphinAttack on popular speech recognition systems, including Siri, Google Now, Samsung S Voice, Huawei HiVoice, Cortana and Alexa. By injecting a sequence of inaudible voice commands, we show a few proof-of-concept attacks, which include activating Siri to initiate a FaceTime call on iPhone, activating Google Now to switch the phone to the airplane mode, and even manipulating the navigation system in an Audi automobile. We propose hardware and software defense solutions. We validate that it is feasible to detect DolphinAttack by classifying the audios using supported vector machine (SVM), and suggest to re-design voice controllable systems to be resilient to inaudible voice command attacks.


title: A Novel Scheduling Framework Leveraging Hardware Cache Partitioning for Cache-Side-Channel Elimination in Clouds

abstract: While there exist many isolation mechanisms that are available to cloud service providers, including virtual machines, containers, etc., the problem of side-channel increases in importance as a remaining security vulnerability, particularly in the presence of shared caches and multicore processors. In this paper we present a hardware-software mechanism that improves the isolation of cloud processes in the presence of shared caches on multicore chips. Combining the Intel CAT architecture that enables cache partitioning on the fly with novel scheduling techniques and state cleansing mechanisms, we enable cache-side-channel free computing for Linux-based containers and virtual machines, in particular, those managed by KVM. We do a preliminary evaluation of our system using a CPU bound workload. Our system allows Simultaneous Multithreading (SMT) to remain enabled and does not require application level changes.


title: Coherent Atomically-Thin Superlattices with Engineered Strain

abstract: Epitaxy forms the basis of modern electronics and optoelectronics. We report coherent atomically-thin superlattices, in which different transition metal dichalcogenide monolayers--despite large lattice mismatches--are repeated and integrated without dislocations. Grown by a novel omnidirectional epitaxy, these superlattices display fully-matched lattice constants across heterointerfaces while maintaining a surprisingly isotropic lattice structure and triangular symmetry. This strong epitaxial strain is precisely engineered via the nanoscale supercell dimensions, thereby enabling broad tuning of the optical properties and producing photoluminescence peak shifts as large as 250 meV. We present theoretical models to explain this coherent growth as well as the energetic interplay governing the flat-rippled configuration space in these strained monolayers. Such coherent superlattices provide novel building blocks with targeted functionalities at the atomically-thin monolayer limit.


title: Virtual Reality: Blessings and Risk Assessment

abstract: Objectives: This paper presents an up-to-date overview of research performed in the Virtual Reality (VR) environment ranging from definitions, its presence in the various fields, and existing market players and their projects in the VR technology. Further an attempt is made to gain an insight on the psychological mechanism underlying experience in using VR device. Methods: Our literature survey is based on the research articles, analysis of the projects of various companies and their findings for different areas of interest. Findings: In our literature survey we observed that the recent advances in virtual reality enabling technologies have led to variety of virtual devices that facilitate people to interact with the digital world. In fact in the past two decades researchers have tried to integrate reality and VR in the form of intuitive computer interface. Improvements: This has led to variety of potential benefits of VR in many applications such as News, Healthcare, Entertainment, Tourism, Military and Defence etc. However despite the extensive research efforts in creating virtual system environments it is yet to become apparent in normal daily life.


title: Searching for new young stars in the northern hemisphere: The Pisces Moving Group

abstract: Using the kinematically unbiased technique described in Binks, Jeffries & Maxted (2015), we present optical spectra for a further 122 rapidly-rotating (rotation periods < 6 days), X-ray active FGK stars, selected from the SuperWASP survey. We identify 17 new examples of young, probably single stars with ages of < 200 Myr and provide additional evidence for a new northern hemisphere kinematic association: the Pisces Moving Group (MG). The group consists of 14 lithium-rich G- and K-type stars, that have a dispersion of only $\sim $3 kms$^{-1}$ in each Galactic space velocity coordinate. The group members are approximately co-eval in the colour-magnitude diagram, with an age of 30-50 Myr, and have similar, though not identical, kinematics to the Octans-Near MG.


title: Hopf bifurcation for a delayed diffusive logistic population model in the advective heterogeneous environment

abstract: In this paper, we investigate a delayed reaction-diffusion-advection equation, which models the population dynamics in the advective heterogeneous environment. The existence of the nonconstant positive steady state and associated Hopf bifurcation are obtained. A weighted inner product associated with the advection rate is introduced to compute the normal forms, which is the main difference between Hopf bifurcation for delayed reaction-diffusion-advection model and that for delayed reaction-diffusion model. Moreover, we find that the spatial scale and advection can affect Hopf bifurcation in the heterogenous environment.


title: Confidence Intervals that Utilize Uncertain Prior Information about Exogeneity in Panel Data

abstract: Consider panel data modelled by a linear random intercept model that includes a time-varying covariate. Suppose that we have uncertain prior information that this covariate is exogenous. We present a new confidence interval for the slope parameter that utilizes this uncertain prior information. This interval has minimum coverage probability very close to its nominal coverage. Let the scaled expected length of this new confidence interval be its expected length divided by the expected length of the confidence interval, with the same minimum coverage, constructed using the fixed effects model. This new interval has scaled expected length that (a) is substantially less than 1 when the prior information is correct, (b) has a maximum value that is not too much larger than 1 and (c) is close to 1 when the data strongly contradict the prior information. We illustrate the properties of this new interval using an airfare data set.


title: Quantum measurement in two-dimensional conformal field theories: Application to quantum energy teleportation

abstract: We construct a set of quasi-local measurement operators in 2D CFT, and then use them to proceed the quantum energy teleportation (QET) protocol and show it is viable. These measurement operators are constructed out of the projectors constructed from shadow operators, but further acting on the product of two spatially separated primary fields. They are equivalently the OPE blocks in the large central charge limit up to some UV-cutoff dependent normalization but the associated probabilities of outcomes are UV-cutoff independent. We then adopt these quantum measurement operators to show that the QET protocol is viable in general. We also check the CHSH inequality a la OPE blocks.


title: Video Summarization with Attention-Based Encoder-Decoder Networks

abstract: This paper addresses the problem of supervised video summarization by formulating it as a sequence-to-sequence learning problem, where the input is a sequence of original video frames, the output is a keyshot sequence. Our key idea is to learn a deep summarization network with attention mechanism to mimic the way of selecting the keyshots of human. To this end, we propose a novel video summarization framework named Attentive encoder-decoder networks for Video Summarization (AVS), in which the encoder uses a Bidirectional Long Short-Term Memory (BiLSTM) to encode the contextual information among the input video frames. As for the decoder, two attention-based LSTM networks are explored by using additive and multiplicative objective functions, respectively. Extensive experiments are conducted on three video summarization benchmark datasets, i.e., SumMe, and TVSum. The results demonstrate the superiority of the proposed AVS-based approaches against the state-of-the-art approaches,with remarkable improvements from 0.8% to 3% on two datasets,respectively..


title: Differentiable cellular automata

abstract: We describe a class of cellular automata (CAs) that are end-to-end differentiable. DCAs interpolate the behavior of ordinary CAs through rules that act on distributions of states. The gradient of a DCA with respect to its parameters can be computed with an iterative propagation scheme that uses previously-computed gradients and values. Gradient-based optimization over DCAs could be used to find ordinary CAs with desired properties.


title: A Note on the Spectral Transfer Morphisms for Affine Hecke Algebras

abstract: E. Opdam introduced the tool of spectral transfer morphism (STM) of affine Hecke algebras to study the formal degrees of unipotent discrete series representations. He established a uniqueness property of STM for the affine Hecke algebras associated of unipotent discrete series representations. Based on this result, Opdam gave an explanation for Lusztig's arithmetic/geometric correspondence (in Lusztig's classification of unipotent representations of $p$-adic adjoint simple groups) in terms of harmonic analysis, and partitioned the unipotent discrete series representations into $L$-packets based on the Lusztig-Langlands parameters. The present paper provides some omitted details for the argument of the uniqueness property of STM. In the last section, we prove that three finite morphisms of algebraic tori are spectral transfer morphisms, and hence complete the proof of the uniqueness property.


title: Neutrinos in Large Extra Dimensions and Short-Baseline $\nu_e$ Appearance

abstract: We show that, in the presence of bulk masses, sterile neutrinos propagating in large extra dimensions (LED) can induce electron-neutrino appearance effects. This is in contrast to what happens in the standard LED scenario and hence LED models with explicit bulk masses have the potential to address the MiniBooNE and LSND appearance results, as well as the reactor and Gallium anomalies. A special feature in our scenario is that the mixing of the first KK modes to active neutrinos can be suppressed, making the contribution of heavier sterile neutrinos to oscillations relatively more important. We study the implications of this neutrino mass generation mechanism for current and future neutrino oscillation experiments, and show that the Short-Baseline Neutrino Program at Fermilab will be able to efficiently probe such a scenario. In addition, this framework leads to massive Dirac neutrinos and thus precludes any signal in neutrinoless double beta decay experiments.


title: Necessary and sufficient conditions for boundedness of commutators of bilinear Hardy-Littlewood maximal function

abstract: Let $\mathcal{M}$ be the bilinear Hardy-Littlewood maximal function and $\vec{b}=(b,b)$ be a collection of locally integrable functions. In this paper, the authors establish characterizations of the weighted {\rm BMO} space in terms of several different commutators of bilinear Hardy-Littlewood maximal function, respectively; these commutators include the maximal iterated commutator $\mathcal{M}_{\Pi \vec{b}}$, the maximal linear commutator $\mathcal{M}_{\Sigma\vec{b}}$, the iterated commutator $[\Pi \vec{b},\mathcal{M}]$ and the linear commutator $[\Sigma \vec{b},\mathcal{M}]$.


title: Generalised contact geometry as reduced generalised complex geometry

abstract: Generalised contact structures are studied from the point of view of reduced generalised complex structures, naturally incorporating non-coorientable structures as non-trivial fibering. The infinitesimal symmetries are described in detail, with a geometric description given in terms of gerbes. As an application of the reduction procedure, generalised coK\"ahler structures are defined in a way which extends the K\"ahler/coK\"ahler correspondence.


title: Thresholds Optimization for One-Bit Feedback Multi-User Scheduling

abstract: We propose a new one-bit feedback scheme with scheduling decision based on the maximum expected weighted rate. We show the concavity of the $2$-user case and provide the optimal solution which achieves the maximum weighted rate of the users. For the general asymmetric M-user case, we provide a heuristic method to achieve the maximum expected weighted rate. We show that the sum rate of our proposed scheme is very close to the sum rate of the full channel state information case, which is the upper bound performance.


title: Cutting sequences, regular polygons, and the Veech group

abstract: We describe the cutting sequences associated to geodesic flow on regular polygons, in terms of a combinatorial process called "derivation." This work is an extension of some of the ideas and results in Smillie and Ulcigrai's recent paper, where the analysis was made for the regular octagon. It turns out that the main structural properties of the octagon generalize in a natural way.


title: Phase separated magnetic ground state in Mn$_3$Ga$_{0.45}$Sn$_{0.55}$C

abstract: Existence of non-ergodic ground states is considered as a precursor to a first order long range magnetostructural transformation. Mn$_3$Ga$_{0.45}$Sn$_{0.55}$C lies compositionally between two compounds, Mn$_3$GaC and Mn$_3$SnC, undergoing first order magnetic transformation. Mn$_3$Ga$_{0.45}$Sn$_{0.55}$C though crystallizes in single phase cubic structure, exhibits more than one long range magnetic transitions. Using a combination of magnetization, ac susceptibility, neutron diffraction and XAFS techniques it is shown that, though Mn$_3$Ga$_{0.45}$Sn$_{0.55}$C exhibits long range magnetic order, it presents a cluster glassy ground state due to formation of magnetically ordered Ga rich and Sn rich clusters. The clusters are big enough to present signatures of long range magnetic order but are distributed in such way that it limits interactions between two clusters of the same type leading to a frozen magnetic state at low temperatures. The main reason for such a cluster glass state is the difference in local structure of Mn atoms that find themselves in Ga rich and Sn rich clusters.


title: Mechanism of magnetostructural transformation in multifunctional Mn$_3$GaC

abstract: Mn$_3$GaC undergoes a ferromagnetic to antiferromagnetic, volume discontinuous cubic-cubic phase transition as a function of temperature, pressure and magnetic field. Through a series of temperature dependent x-ray absorption fine structure spectroscopy experiments at the Mn K and Ga K edge, it is shown that the first order magnetic transformation in Mn$_3$GaC is entirely due to distortions in Mn sub-lattice and with a very little role for Mn-C interactions. The distortion in Mn sub-lattice results in long and short Mn-Mn bonds with the longer Mn-Mn bonds favoring ferromagnetic interactions and the shorter Mn-Mn bonds favoring antiferromagnetic interactions. At the first order transition, the shorter Mn-Mn bonds exhibit an abrupt decrease in their length resulting in an antiferromagnetic ground state and a strained lattice.


title: Translations in the exponential Orlicz space with Gaussian weight

abstract: We study the continuity of space translations on non-parametric exponential families based on the exponential Orlicz space with Gaussian reference density.


title: Finite-error metrological bounds on the multi-parameter Hamiltonian estimation

abstract: Estimation of multiple parameters in an unknown Hamiltonian is investigated. We present upper and lower bounds on the time required to complete the estimation within a prescribed tolerance $\delta$. The lower bound is given on the basis of the Cram\'er-Rao inequality, where the quantum Fisher information is bounded by the squared evolution time. The upper bound is obtained by an explicit construction of estimation procedures. By comparing the cases with different numbers of Hamiltonian channels, we also find that the few-channel procedure with adaptive feedback and the many-channel procedure with entanglement are equivalent in that they require the same amount of time resource up to a constant factor.


title: Definite Determinantal Representations of Multivariate Polynomials

abstract: In this paper, we consider the problem of representing a multivariate polynomial as the determinant of a definite (monic) symmetric/Hermitian linear matrix polynomial (LMP). Such a polynomial is known as determinantal polynomial. Determinantal polynomials can characterize the feasible sets of semidefinite programming problems that motivates us to deal with this problem. We introduce the notion of generalized mixed discriminant of matrices which translates the determinantal representation problem into computing a point of a real variety of a specified ideal. We develop an algorithm to determine such a determinantal representation of a bivariate polynomial of degree $d$. Then we propose a heuristic method to obtain a monic symmetric determinantal representation of a multivariate polynomial of degree $d$.


title: Cech Closure Spaces: A Unified Framework for Discrete and Continuous Homotopy

abstract: Motivated by constructions in topological data analysis and algebraic combinatorics, we study homotopy theory on the category of Cech closure spaces $\mathbf{Cl}$, the category whose objects are sets endowed with a Cech closure operator and whose morphisms are the continuous maps between them. We introduce new classes of Cech closure structures on metric spaces, graphs, and simplicial complexes, and we show how each of these cases gives rise to an interesting homotopy theory. In particular, we show that there exists a natural family of Cech closure structures on metric spaces which produces a non-trivial homotopy theory for finite metric spaces, i.e. point clouds, the spaces of interest in topological data analysis. We then give a Cech closure structure to graphs and simplicial complexes which may be used to construct a new combinatorial (as opposed to topological) homotopy theory for each skeleton of those spaces. We further show that there is a Seifert-van Kampen theorem for closure spaces, a well-defined notion of persistent homotopy, and an associated interleaving distance. As an illustration of the difference with the topological setting, we calculate the fundamental group for the circle, `circular graphs', and the wedge of circles endowed with different closure structures. Finally, we produce a continuous map from the topological circle to `circular graphs' which, given the appropriate closure structures, induces an isomorphism on the fundamental groups.


title: Definite Determinantal Representations via Orthostochastic Matrices

abstract: Determinantal polynomials play a crucial role in semidefinite programming problems. Helton-Vinnikov proved that real zero (RZ) bivariate polynomials are determinantal. However, it leads to a challenging problem to compute such a determinantal representation. We provide a necessary and sufficient condition for the existence of definite determinantal representation of a bivariate polynomial by identifying its coefficients as scalar products of two vectors where the scalar products are defined by orthostochastic matrices. This alternative condition enables us to develop a method to compute a monic symmetric/Hermitian determinantal representations for a bivariate polynomial of degree $d$. In addition, we propose a computational relaxation to the determinantal problem which turns into a problem of expressing the vector of coefficients of the given polynomial as convex combinations of some specified points. We also characterize the range set of vector coefficients of a certain type of determinantal bivariate polynomials.


title: A Note on Plus-Contacts, Rectangular Duals, and Box-Orthogonal Drawings

abstract: A plus-contact representation of a planar graph $G$ is called $c$-balanced if for every plus shape $+_v$, the number of other plus shapes incident to each arm of $+_v$ is at most $ c \Delta +O(1)$, where $\Delta$ is the maximum degree of $G$. Although small values of $c$ have been achieved for a few subclasses of planar graphs (e.g., $2$- and $3$-trees), it is unknown whether $c$-balanced representations with $c<1$ exist for arbitrary planar graphs. In this paper we compute $(1/2)$-balanced plus-contact representations for all planar graphs that admit a rectangular dual. Our result implies that any graph with a rectangular dual has a 1-bend box-orthogonal drawings such that for each vertex $v$, the box representing $v$ is a square of side length $\frac{deg(v)}{2}+ O(1)$.


title: Optimal Dynamic Cloud Network Control

abstract: Distributed cloud networking enables the deployment of a wide range of services in the form of interconnected software functions instantiated over general purpose hardware at multiple cloud locations distributed throughout the network. We consider the problem of optimal service delivery over a distributed cloud network, in which nodes are equipped with both communication and computation resources. We address the design of distributed online solutions that drive flow processing and routing decisions, along with the associated allocation of cloud and network resources. For a given set of services, each described by a chain of service functions, we characterize the cloud network capacity region and design a family of dynamic cloud network control (DCNC) algorithms that stabilize the underlying queuing system, while achieving arbitrarily close to minimum cost with a tradeoff in network delay. The proposed DCNC algorithms make local decisions based on the online minimization of linear and quadratic metrics obtained from an upper bound on the Lyapunov drift-plus-penalty of the cloud network queuing system. Minimizing a quadratic vs. a linear metric is shown to improve the cost-delay tradeoff at the expense of increased computational complexity. Our algorithms are further enhanced with a shortest transmission-plus-processing distance bias that improves delay performance without compromising throughput or overall cloud network cost. We provide throughput and cost optimality guarantees, convergence time analysis, and extensive simulations in representative cloud network scenarios.


title: Matched disturbance rejection for energy-shaping controlled underactuated mechanical systems

abstract: In this paper, we present a method of applying integral action to enhance the robustness of energy shaping controllers for underactuated mechanical systems with matched disturbances. Previous works on this problem have required a number of technical assumptions to be satisfied, restricting the class of systems for which the proposed solution applies. The design proposed in this paper relaxes some of these technical assumptions.


title: Controlling of ZnO nanostructures morphology for guiding the structural and optical properties

abstract: Most of the applications of the nano structures dependent on the morphology which affects the opto electronics properties. This research article provides a pathway of guiding optical properties like band-gap and fluorescence properties by controlled growth of nano-rods, -flowers, -needles or- tubes without external chemical doping, by simple hydro thermal method by controlling over synthesis parameter, temperature.