title: Decay of turbulence in a liquid metal duct flow with transverse magnetic field

abstract: Decay of honeycomb-generated turbulence in a duct with a static transverse magnetic field is studied via direct numerical simulations. The simulations follow the revealing experimental study of Sukoriansky et al. (1986), in particular the paradoxical observation of high-amplitude velocity fluctuations, which exist in the downstream portion of the flow when the strong transverse magnetic field is imposed in the entire duct including the honeycomb exit, but not in other configurations. It is shown that the fluctuations are caused by the large-scale quasi-two-dimensional structures forming in the flow at the initial stages of the decay and surviving the magnetic suppression. Statistical turbulence properties, such as the energy decay curves, two-point correlations and typical length scales are computed. The study demonstrates that turbulence decay in the presence of a magnetic field is a complex phenomenon critically depending on the state of the flow at the moment the field is introduced.


title: Performance of MPI sends of non-contiguous data

abstract: We present an experimental investigation of the performance of MPI derived datatypes. For messages up to the megabyte range most schemes perform comparably to each other and to manual copying into a regular send buffer. However, for large messages the internal buffering of MPI causes differences in efficiency. The optimal scheme is a combination of packing and derived types.


title: On the Reliability Roots of Simplicial Complexes and Matroids

abstract: Assume that the vertices of a graph $G$ are always operational, but the edges of $G$ fail independently with probability $q \in[0,1]$. The \emph{all-terminal reliability} of $G$ is the probability that the resulting subgraph is connected. The all-terminal reliability can be formulated into a polynomial in $q$, and it was conjectured \cite{BC1} that all the roots of (nonzero) reliability polynomials fall inside the closed unit disk. It has since been shown that there exist some connected graphs which have their reliability roots outside the closed unit disk, but these examples seem to be few and far between, and the roots are only barely outside the disk. In this paper we generalize the notion of reliability to simplicial complexes and matroids and investigate when, for small simplicial complexes and matroids, the roots fall inside the closed unit disk.


title: Morpho-MNIST: Quantitative Assessment and Diagnostics for Representation Learning

abstract: Revealing latent structure in data is an active field of research, having introduced exciting technologies such as variational autoencoders and adversarial networks, and is essential to push machine learning towards unsupervised knowledge discovery. However, a major challenge is the lack of suitable benchmarks for an objective and quantitative evaluation of learned representations. To address this issue we introduce Morpho-MNIST, a framework that aims to answer: "to what extent has my model learned to represent specific factors of variation in the data?" We extend the popular MNIST dataset by adding a morphometric analysis enabling quantitative comparison of trained models, identification of the roles of latent variables, and characterisation of sample diversity. We further propose a set of quantifiable perturbations to assess the performance of unsupervised and supervised methods on challenging tasks such as outlier detection and domain adaptation. Data and code are available at https://github.com/dccastro/Morpho-MNIST.


title: Methods and Concepts in Economic Complexity

abstract: Knowhow in societies accumulates as it gets transmitted from group to group, and from generation to generation. However, we lack of a unified quantitative formalism that takes into account the structured process for how this accumulation occurs, and this has precluded the development of a unified view of human development in the past and in the present. Here, we summarize a paradigm to understand and model this process. The paradigm goes under the general name of the Theory of Economic Complexity (TEC). Based on it, we present a combination of analytical, numerical and empirical results that illustrate how to characterize the process of development, providing measurable quantities that can be used to predict future developments. The emphasis is the quantification of the collective knowhow an economy has accumulated, and what are the directions in which it is likely to expand. As a case study we consider data on trade, which provides consistent data on the technological diversification of 200 countries across more than 50 years. The paradigm represented by TEC should be relevant for anthropologists, sociologists, and economists interested in the role of collective knowhow as the main determinant of the success and welfare of a society.


title: A User-based Visual Analytics Workflow for Exploratory Model Analysis

abstract: Many visual analytics systems allow users to interact with machine learning models towards the goals of data exploration and insight generation on a given dataset. However, in some situations, insights may be less important than the production of an accurate predictive model for future use. In that case, users are more interested in generating of diverse and robust predictive models, verifying their performance on holdout data, and selecting the most suitable model for their usage scenario. In this paper, we consider the concept of Exploratory Model Analysis (EMA), which is defined as the process of discovering and selecting relevant models that can be used to make predictions on a data source. We delineate the differences between EMA and the well-known term exploratory data analysis in terms of the desired outcome of the analytic process: insights into the data or a set of deployable models. The contributions of this work are a visual analytics system workflow for EMA, a user study, and two use cases validating the effectiveness of the workflow. We found that our system workflow enabled users to generate complex models, to assess them for various qualities, and to select the most relevant model for their task.


title: Dual Selection Games

abstract: Often, a given selection game studied in the literature has a known dual game. In dual games, a winning strategy for a player in either game may be used to create a winning strategy for the opponent in the dual. For example, the Rothberger selection game involving open covers is dual to the point-open game. This extends to a general theorem: if $\{\operatorname{ran}{f}:f\in\mathbf C(\mathcal R)\}$ is coinitial in $\mathcal A$ with respect to $\subseteq$, where $\mathbf C(\mathcal R)=\{f\in(\bigcup\mathcal R)^{\mathcal R}:R\in\mathcal R\Rightarrow f(R)\in R\}$ collects the choice functions on the set $\mathcal R$, then $G_1(\mathcal A,\mathcal B)$ and $G_1(\mathcal R,\neg\mathcal B)$ are dual selection games.


title: Adaptive Gaussian process surrogates for Bayesian inference

abstract: We present an adaptive approach to the construction of Gaussian process surrogates for Bayesian inference with expensive-to-evaluate forward models. Our method relies on the fully Bayesian approach to training Gaussian process models and utilizes the expected improvement idea from Bayesian global optimization. We adaptively construct training designs by maximizing the expected improvement in fit of the Gaussian process model to the noisy observational data. Numerical experiments on model problems with synthetic data demonstrate the effectiveness of the obtained adaptive designs compared to the fixed non-adaptive designs in terms of accurate posterior estimation at a fraction of the cost of inference with forward models.


title: A gentle introduction to anisotropic Banach spaces

abstract: The use of anisotropic Banach spaces has provided a wealth of new results in the study of hyperbolic dynamical systems in recent years, yet their application to specific systems is often technical and difficult to access. The purpose of this note is to provide an introduction to the use of these spaces in the study of hyperbolic maps and to highlight the important elements and how they work together. This is done via a concrete example of a family of dissipative Baker's transformations. Along the way, we prove an original result connecting such transformations with expanding maps via a continuous family of transfer operators acting on a single Banach space.


title: Fast SGL Fourier transforms for scattered data

abstract: Spherical Gauss-Laguerre (SGL) basis functions, i. e., normalized functions of the type $L_{n-l-1}^{(l + 1/2)}(r^2) r^l Y_{lm}(\vartheta,\varphi)$, $|m| \leq l < n \in \mathbb{N}$, $L_{n-l-1}^{(l + 1/2)}$ being a generalized Laguerre polynomial, $Y_{lm}$ a spherical harmonic, constitute an orthonormal polynomial basis of the space $L^2$ on $\mathbb{R}^3$ with radial Gaussian (multivariate Hermite) weight $\exp(-r^2)$. We have recently described fast Fourier transforms for the SGL basis functions based on an exact quadrature formula with certain grid points in $\mathbb{R}^3$. In this paper, we present fast SGL Fourier transforms for scattered data. The idea is to employ well-known basal fast algorithms to determine a three-dimensional trigonometric polynomial that coincides with the bandlimited function of interest where the latter is to be evaluated. This trigonometric polynomial can then be evaluated efficiently using the well-known non-equispaced FFT (NFFT). We proof an error estimate for our algorithms and validate their practical suitability in extensive numerical experiments.


title: Complexity of Training ReLU Neural Network

abstract: In this paper, we explore some basic questions on the complexity of training neural networks with ReLU activation function. We show that it is NP-hard to train a two-hidden layer feedforward ReLU neural network. If dimension of the input data and the network topology is fixed, then we show that there exists a polynomial time algorithm for the same training problem. We also show that if sufficient over-parameterization is provided in the first hidden layer of ReLU neural network, then there is a polynomial time algorithm which finds weights such that output of the over-parameterized ReLU neural network matches with the output of the given data.


title: Learning and Acting in Peripersonal Space: Moving, Reaching, and Grasping

abstract: The young infant explores its body, its sensorimotor system, and the immediately accessible parts of its environment, over the course of a few months creating a model of peripersonal space useful for reaching and grasping objects around it. Drawing on constraints from the empirical literature on infant behavior, we present a preliminary computational model of this learning process, implemented and evaluated on a physical robot. The learning agent explores the relationship between the configuration space of the arm, sensing joint angles through proprioception, and its visual perceptions of the hand and grippers. The resulting knowledge is represented as the peripersonal space (PPS) graph, where nodes represent states of the arm, edges represent safe movements, and paths represent safe trajectories from one pose to another. In our model, the learning process is driven by intrinsic motivation. When repeatedly performing an action, the agent learns the typical result, but also detects unusual outcomes, and is motivated to learn how to make those unusual results reliable. Arm motions typically leave the static background unchanged, but occasionally bump an object, changing its static position. The reach action is learned as a reliable way to bump and move an object in the environment. Similarly, once a reliable reach action is learned, it typically makes a quasi-static change in the environment, moving an object from one static position to another. The unusual outcome is that the object is accidentally grasped (thanks to the innate Palmar reflex), and thereafter moves dynamically with the hand. Learning to make grasps reliable is more complex than for reaches, but we demonstrate significant progress. Our current results are steps toward autonomous sensorimotor learning of motion, reaching, and grasping in peripersonal space, based on unguided exploration and intrinsic motivation.


title: An Empirical Comparison of Syllabuses for Curriculum Learning

abstract: Syllabuses for curriculum learning have been developed on an ad-hoc, per task basis and little is known about the relative performance of different syllabuses. We identify a number of syllabuses used in the literature. We compare the identified syllabuses based on their effect on the speed of learning and generalization ability of a LSTM network on three sequential learning tasks. We find that the choice of syllabus has limited effect on the generalization ability of a trained network. In terms of speed of learning our results demonstrate that the best syllabus is task dependent but that a recently proposed automated curriculum learning approach - Predictive Gain, performs very competitively against all identified hand-crafted syllabuses. The best performing hand-crafted syllabus which we term Look Back and Forward combines a syllabus which steps through tasks in the order of their difficulty with a uniform distribution over all tasks. Our experimental results provide an empirical basis for the choice of syllabus on a new problem that could benefit from curriculum learning. Additionally, insights derived from our results shed light on how to successfully design new syllabuses.


title: Deep Object Pose Estimation for Semantic Robotic Grasping of Household Objects

abstract: Using synthetic data for training deep neural networks for robotic manipulation holds the promise of an almost unlimited amount of pre-labeled training data, generated safely out of harm's way. One of the key challenges of synthetic data, to date, has been to bridge the so-called reality gap, so that networks trained on synthetic data operate correctly when exposed to real-world data. We explore the reality gap in the context of 6-DoF pose estimation of known objects from a single RGB image. We show that for this problem the reality gap can be successfully spanned by a simple combination of domain randomized and photorealistic data. Using synthetic data generated in this manner, we introduce a one-shot deep neural network that is able to perform competitively against a state-of-the-art network trained on a combination of real and synthetic data. To our knowledge, this is the first deep network trained only on synthetic data that is able to achieve state-of-the-art performance on 6-DoF object pose estimation. Our network also generalizes better to novel environments including extreme lighting conditions, for which we show qualitative results. Using this network we demonstrate a real-time system estimating object poses with sufficient accuracy for real-world semantic grasping of known household objects in clutter by a real robot.


title: Estimation of Personalized Effects Associated With Causal Pathways

abstract: The goal of personalized decision making is to map a unit's characteristics to an action tailored to maximize the expected outcome for that unit. Obtaining high-quality mappings of this type is the goal of the dynamic regime literature. In healthcare settings, optimizing policies with respect to a particular causal pathway may be of interest as well. For example, we may wish to maximize the chemical effect of a drug given data from an observational study where the chemical effect of the drug on the outcome is entangled with the indirect effect mediated by differential adherence. In such cases, we may wish to optimize the direct effect of a drug, while keeping the indirect effect to that of some reference treatment. [16] shows how to combine mediation analysis and dynamic treatment regime ideas to defines policies associated with causal pathways and counterfactual responses to these policies. In this paper, we derive a variety of methods for learning high quality policies of this type from data, in a causal model corresponding to a longitudinal setting of practical importance. We illustrate our methods via a dataset of HIV patients undergoing therapy, gathered in the Nigerian PEPFAR program.


title: Cursive Scene Text Analysis by Deep Convolutional Linear Pyramids

abstract: The camera captured images have various aspects to investigate. Generally, the emphasis of research depends on the interesting regions. Sometimes the focus could be on color segmentation, object detection or scene text analysis. The image analysis, visibility and layout analysis are the tasks easier for humans as suggested by behavioral trait of humans, but in contrast when these same tasks are supposed to perform by machines then it seems to be challenging. The learning machines always learn from the properties associated to provided samples. The numerous approaches are designed in recent years for scene text extraction and recognition and the efforts are underway to improve the accuracy. The convolutional approach provided reasonable results on non-cursive text analysis appeared in natural images. The work presented in this manuscript exploited the strength of linear pyramids by considering each pyramid as a feature of the provided sample. Each pyramid image process through various empirically selected kernels. The performance was investigated by considering Arabic text on each image pyramid of EASTR-42k dataset. The error rate of 0.17% was reported on Arabic scene text recognition.


title: Weak-Interaction Tests via Precision Superallowed $\beta$-Decay Studies at TRIUMF

abstract: Superallowed $\beta$-decay studies provide some of the best constraints on the possibility of additional quark generations, as well as limits on exotic currents in the weak interaction. The three experimental quantities that are required for performing these tests using $0^+\rightarrow0^+$ nuclear decays (branching ratio, half-life, and $Q$-value) can all be measured to high-precision with rare-isotope beams at the TRIUMF-ISAC facility in Vancouver, Canada. This proceeding presents a brief outline of the general experimental techniques used at TRIUMF over the past 15 years, as well as recent theoretical advances towards {\it ab-initio} isospin-symmetry-breaking corrections to superallowed nuclear decays.


title: Model-Preserving Sensitivity Analysis for Families of Gaussian Distributions

abstract: The accuracy of probability distributions inferred using machine-learning algorithms heavily depends on data availability and quality. In practical applications it is therefore fundamental to investigate the robustness of a statistical model to misspecification of some of its underlying probabilities. In the context of graphical models, investigations of robustness fall under the notion of sensitivity analyses. These analyses consist in varying some of the model's probabilities or parameters and then assessing how far apart the original and the varied distributions are. However, for Gaussian graphical models, such variations usually make the original graph an incoherent representation of the model's conditional independence structure. Here we develop an approach to sensitivity analysis which guarantees the original graph remains valid after any probability variation and we quantify the effect of such variations using different measures. To achieve this we take advantage of algebraic techniques to both concisely represent conditional independence and to provide a straightforward way of checking the validity of such relationships. Our methods are demonstrated to be robust and comparable to standard ones, which break the conditional independence structure of the model, using an artificial example and a medical real-world application.


title: A Hybrid Neural Network Framework and Application to Radar Automatic Target Recognition

abstract: Deep neural networks (DNNs) have found applications in diverse signal processing (SP) problems. Most efforts either directly adopt the DNN as a black-box approach to perform certain SP tasks without taking into account of any known properties of the signal models, or insert a pre-defined SP operator into a DNN as an add-on data processing stage. This paper presents a novel hybrid-NN framework in which one or more SP layers are inserted into the DNN architecture in a coherent manner to enhance the network capability and efficiency in feature extraction. These SP layers are properly designed to make good use of the available models and properties of the data. The network training algorithm of hybrid-NN is designed to actively involve the SP layers in the learning goal, by simultaneously optimizing both the weights of the DNN and the unknown tuning parameters of the SP operators. The proposed hybrid-NN is tested on a radar automatic target recognition (ATR) problem. It achieves high validation accuracy of 96\% with 5,000 training images in radar ATR. Compared with ordinary DNN, hybrid-NN can markedly reduce the required amount of training data and improve the learning performance.


title: FMIT: Feature Model Integration Techniques

abstract: Although feature models are widely used in practice, for example, representing variability in software product lines, their integration is still a challenge. Many integration techniques have been proposed, although none of these have proven to be fully effective. Integrating feature models becomes a difficult, costly, error-prone task. Since their transition occurs in a generalized and automated way, the techniques applied to compose the models end up giving rise to a final model, in many cases undesired, without taking into account the specific needs arising from the requirements determined by the analysts and developers. Therefore, this work proposes FMIT, a technique for integrating feature models. The FMIT is based on contemporary model integration strategies to increase the accuracy and quality of the integrated feature model. In this way, it will be possible to identify the degree of similarity between composite feature diagrams, to verify their accuracy, as well as to identify conflicts. In addition, this work proposes the development of a prototype based on the set of strategies, used to take decisions according to the requirements established during the integration of feature models, whether this is semi-automatic or automatic. To evaluate FMIT, experimental studies were conducted with 10 participants, including students and professionals. Participants performed 12 integration scenarios, 6 using the FMIT and 6 manually. The results suggest that FMIT improved accuracy by 43\% of the cases, as well as reduced the effort by 70\% to perform the integrations.


title: Collaborative Robot Learning from Demonstrations using Hidden Markov Model State Distribution

abstract: In robotics, there is need of an interactive and expedite learning method as experience is expensive. Robot Learning from Demonstration (RLfD) enables a robot to learn a policy from demonstrations performed by teacher. RLfD enables a human user to add new capabilities to a robot in an intuitive manner, without explicitly reprogramming it. In this work, we present a novel interactive framework, where a collaborative robot learns skills for trajectory based tasks from demonstrations performed by a human teacher. The robot extracts features from each demonstration called as key-points and learns a model of the demonstrated skill using Hidden Markov Model (HMM). Our experimental results show that the learned model can be used to produce a generalized trajectory based skill.


title: Large vacuum flux surfaces generated by tilted planar coils

abstract: Helical equilibria can be generated by arrangements of planar coils similar to tokamaks, but without a central solenoid and with the toroidal field (TF) coils tilted with respect to the vertical. This is known from earlier numerical works, e.g. P.E. Moroz, Phys.Plasmas 2, 4269 (1995). However, such concept tends to need large coils (of low aspect ratio) but form small plasmas (of large aspect ratio). Here it is numerically shown that larger, more attractive vacuum flux surfaces -- relative to the size of the device -- can be generated by carefully optimizing the inclination of the TF coils and currents in the various coil-sets. Vacuum configurations of aspect ratios as low as 4 are found for 6 tilted TF circular coils. Higher numbers of TF coils have advantages (smaller effective ripple) and disadvantages (lower rotational transform, smaller plasma). Finally, the aspect-ratio $A$ of the vacuum flux surfaces is quantified as a function of the ratio $A_c$ of the coil-radius to the radial location of the coil-center. It is found that, in order to minimize $A$, it is beneficial to interlink or marginally interlink the TF coils ($A_c \lesssim 1$).


title: FanStore: Enabling Efficient and Scalable I/O for Distributed Deep Learning

abstract: Emerging Deep Learning (DL) applications introduce heavy I/O workloads on computer clusters. The inherent long lasting, repeated, and random file access pattern can easily saturate the metadata and data service and negatively impact other users. In this paper, we present FanStore, a transient runtime file system that optimizes DL I/O on existing hardware/software stacks. FanStore distributes datasets to the local storage of compute nodes, and maintains a global namespace. With the techniques of system call interception, distributed metadata management, and generic data compression, FanStore provides a POSIX-compliant interface with native hardware throughput in an efficient and scalable manner. Users do not have to make intrusive code changes to use FanStore and take advantage of the optimized I/O. Our experiments with benchmarks and real applications show that FanStore can scale DL training to 512 compute nodes with over 90\% scaling efficiency.


title: On two-weight norm inequalities for positive dyadic operators

abstract: Let $\sigma$ and $\omega$ be locally finite Borel measures on $\mathbb{R}^d$, and let $p\in(1,\infty)$ and $q\in(0,\infty)$. We study the two-weight norm inequality $$ \lVert T(f\sigma) \rVert_{L^q(\omega)}\leq C \lVert f \rVert_{L^p(\sigma)}, \quad \text{for all} \, \, f \in L^p(\sigma), $$ for both the positive summation operators $T=T_\lambda(\cdot \sigma)$ and positive maximal operators $T=M_\lambda(\cdot \sigma)$. Here, for a family $\{\lambda_Q\}$ of non-negative reals indexed by the dyadic cubes $Q$, these operators are defined by $$ T_\lambda(f\sigma):=\sum_Q \lambda_Q \langle f\rangle^\sigma_Q 1_Q \quad\text{ and } \quad M_\lambda(f\sigma):=\sup_Q \lambda_Q \langle f\rangle^\sigma_Q 1_Q, $$ where $\langle f\rangle^\sigma_Q:=\frac{1}{\sigma(Q)} \int_Q |f| d \sigma.$ We obtain new characterizations of the two-weight norm inequalities in the following cases: 1. For $T=T_\lambda(\cdot\sigma)$ in the subrange $q<p$. Under the additional assumption that $\sigma$ satisfies the $A_\infty$ condition with respect to $\omega$, we characterize the inequality in terms of a simple integral condition. The proof is based on characterizing the multipliers between certain classes of Carleson measures. 2. For $T=M_\lambda(\cdot \sigma)$ in the subrange $q<p$. We introduce a scale of simple conditions that depends on an integrability parameter and show that, on this scale, the sufficiency and necessity are separated only by an arbitrarily small integrability gap. 3. For the summation operators $T=T_\lambda(\cdot\sigma)$ in the subrange $1<q<p$. We characterize the inequality for summation operators by means of related inequalities for maximal operators $T=M_\lambda(\cdot \sigma)$. This maximal-type characterization is an alternative to the known potential-type characterization.


title: Effective Cloud Detection and Segmentation using a Gradient-Based Algorithm for Satellite Imagery; Application to improve PERSIANN-CCS

abstract: Being able to effectively identify clouds and monitor their evolution is one important step toward more accurate quantitative precipitation estimation and forecast. In this study, a new gradient-based cloud-image segmentation technique is developed using tools from image processing techniques. This method integrates morphological image gradient magnitudes to separable cloud systems and patches boundaries. A varying scale-kernel is implemented to reduce the sensitivity of image segmentation to noise and capture objects with various finenesses of the edges in remote-sensing images. The proposed method is flexible and extendable from single- to multi-spectral imagery. Case studies were carried out to validate the algorithm by applying the proposed segmentation algorithm to synthetic radiances for channels of the Geostationary Operational Environmental Satellites (GOES-R) simulated by a high-resolution weather prediction model. The proposed method compares favorably with the existing cloud-patch-based segmentation technique implemented in the PERSIANN-CCS (Precipitation Estimation from Remotely Sensed Information using Artificial Neural Network - Cloud Classification System) rainfall retrieval algorithm. Evaluation of event-based images indicates that the proposed algorithm has potential to improve rain detection and estimation skills with an average of more than 45% gain comparing to the segmentation technique used in PERSIANN-CCS and identifying cloud regions as objects with accuracy rates up to 98%.


title: Signature of inflation in the stochastic gravitational wave background generated by cosmic string networks

abstract: A cosmic string network created during an inflationary stage in the early Universe - here defined as i-string network - is expected to enter a transient stretching regime during inflation, in which its characteristic length is stretched to scales much larger than the Hubble radius, before attaining a standard evolution once the network re-enters the Hubble volume after inflation. During the stretching regime, the production of cosmic string loops and consequent emission of gravitational radiation are significantly suppressed. Here, we compute the power spectrum of the stochastic gravitational wave background generated by i-string networks using the velocity-dependent one scale model to describe the network dynamics, and demonstrate that this regime introduces a high-frequency signature on an otherwise standard spectrum of the stochastic gravitational wave background generated by cosmic strings. We argue that, if observed by current or forthcoming experiments, this signature would provide strong evidence for i-strings and, therefore, for (primordial) inflation. We also develop a simple single-parameter algorithm for the computation of the stochastic gravitational wave background generated by i-strings from that of a standard cosmic string network, which may be useful in the determination of the observational constraints to be obtained by current and forthcoming gravitational wave experiments.


title: Distant White Dwarfs in the US Naval Observatory Flagstaff Station Parallax Sample

abstract: This paper presents new trigonometric parallaxes and proper motions for 214 stars. The measurements were made at the US Naval Observatory Flagstaff Station (NOFS) between 1989 and 2017, and the average uncertainty in the parallax values is 0.6 mas. We find good agreement with Gaia Data Release 2 measurements for the stars in common, although there may be a small systematic offset similar to what has been found by other investigators. The sample is matched to catalogs and the literature to create a photometric dataset which spans the ultraviolet to the mid-infrared. New mid-infrared photometry is obtained for nineteen stars from archived Spitzer mosaics. New optical spectroscopy is presented for seven systems and additional spectra were obtained from the literature. We identify a sub-sample of 179 white dwarfs (WDs) at distances of 25 - 200 pc. Their spectral energy distributions (SEDs) are analyzed using model atmospheres. The models reproduce the entire flux-calibrated SED very well and provide the atmospheric chemical composition, temperature, surface gravity, mass and cooling age of each WD. Twenty-six WDs are newly classified and twelve systems are presented as candidate unresolved binaries. We confirm one WD+red dwarf system and identify two WDs as candidate dust disk systems. Twelve old and high-velocity systems are identified as candidate thick disk or halo objects. The WDs in the sample generally have Galactic disk-like ages of < 8 Gyr and masses close to the canonical 0.6 M_Sun.


title: Patient Risk Assessment and Warning Symptom Detection Using Deep Attention-Based Neural Networks

abstract: We present an operational component of a real-world patient triage system. Given a specific patient presentation, the system is able to assess the level of medical urgency and issue the most appropriate recommendation in terms of best point of care and time to treat. We use an attention-based convolutional neural network architecture trained on 600,000 doctor notes in German. We compare two approaches, one that uses the full text of the medical notes and one that uses only a selected list of medical entities extracted from the text. These approaches achieve 79% and 66% precision, respectively, but on a confidence threshold of 0.6, precision increases to 85% and 75%, respectively. In addition, a method to detect warning symptoms is implemented to render the classification task transparent from a medical perspective. The method is based on the learning of attention scores and a method of automatic validation using the same data.


title: The first ultraviolet detection of the Large Magellanic Cloud pulsar PSR B0540-69 and its multi-wavelength properties

abstract: We observed the young ($\sim 1700$ yrs) pulsar PSR B0540-69 in the near-ultraviolet (UV) for the first time with the Space Telescope Imaging Spectrograph (STIS) aboard the {\em Hubble Space Telescope}. Imaging observations with the NUV- and FUV-MAMA detectors in TIME-TAG mode allowed us to clearly detect the pulsar in two bands around 2350\AA\ and 1590\AA, with magnitudes $m_{\rm NUV} =21.449 \pm 0.019$ and $m_{\rm FUV} =21.832 \pm 0.103$. We also detected the pulsar-wind nebula (PWN) in the NUV-MAMA image, with a morphology similar to that observed in the optical and near-infrared (IR). The extinction-corrected NUV and FUV pulsar fluxes are compatible with a very steep power law spectrum $F_{\nu} \propto \nu^{-\alpha}$ with spectral index $\alpha_{\rm UV} \sim 3$, non compatible with a Rayleigh Jeans spectrum, indicating a non-thermal origin of the emission. The comparison with the optical/near-IR power-law spectrum (spectral index $\alpha_{\rm O,nIR} \sim 0.7$), indicates an abrupt turn-off at wavelengths below 2500 \AA, not observed in other pulsars. We detected pulsations in both the NUV and FUV data at the 50 ms pulsar period. In both cases, the pulse profile features two peaks closely spaced in phase, as observed in the optical and X-ray light curves. The NUV/FUV peaks are also aligned in phase with those observed in the radio (1.4 GHz), optical, X, and $\gamma$-ray light curves, like in the Crab pulsar, implying a similar beaming geometry across all wavelengths. PSR B0540-69 is now the fifth isolated pulsar, together with Crab, Vela, PSR\, B0656+14, and the radio-quiet Geminga, detected in the optical, near-UV, near-IR, X-rays and $\gamma$-rays, and seen to pulsate in at least four of these energy bands.


title: Direct observation of electron capture & emission processes by the time domain charge pumping measurement of MoS2 FET

abstract: Understanding interface properties in MoS2 field effect transistors with a high-k gate insulator is critical for improving the performance of the device. Here, by applying the time domain charge pumping method, the elementary process for capture and emission of electrons to the interface states is monitored directly using a fast acquisition system. The main outcome is the detection of the clear difference in the capture and emission process of electrons to the interface states. In addition to the transient current response for gate capacitance, the current peak is observed during electron capture, while the broad tail is detected during electron emission. This different behavior is associated with the fact that the time constant for electron capture is much shorter than that for electron emission. Moreover, Dit is evaluated to be in the range of 1012 - 1013 cm-2 eV-1, which is comparable with that estimated from subthreshold swing.


title: Accumulation-mode two-dimensional field-effect transistor: Operation mechanism and thickness scaling rule

abstract: Understanding the operation mode of a two-dimensional (2D) material-based field-effect transistor (FET) is one of the most essential issues in the study of electronics and physics. The existing Schottky barrier-FET model for devices with global back gate and metallic contacts overemphasizes the metal-2D contact effect, and the widely observed residual conductance cannot be explained by this model. Here, an accumulation-mode FET model, which directly reveals 2D channel transport properties, is developed based on a partial top-gate MoS2 FET with metallic contacts and a channel thickness of 0.65~118 nm. The operation mechanism of an accumulation-mode FET is validated and clarified by carefully performed capacitance measurements. A depletion capacitance-quantum capacitance transition is observed. After the analysis of the MoS2 accumulation-mode FET, we have confirmed that most 2D-FETs show accumulation-mode behavior. The universal thickness scaling rule of 2D-FETs is then proposed, which provides guidance for future research on 2D materials.


title: Game-Theoretic Model and Experimental Investigation of Cyber Wargaming

abstract: We demonstrate that game-theoretic calculations serve as a useful tool for assisting cyber wargaming teams in identifying useful strategies. We note a significant similarity between formulating cyber wargaming strategies and the methodology known in military practice as Course of Action (COA) generation. For scenarios in which the attacker must penetrate multiple layers in a defense-in-depth security configuration, an accounting of attacker and defender costs and penetration probabilities provides cost-utility payoff matrices and penetration probability matrices. These can be used as decision tools by both the defender and attacker. Inspection of the matrices allows players to deduce preferred strategies (or COAs) based on game-theoretical equilibrium solutions. The matrices also help in analyzing anticipated effects of potential human-based choices of wargame strategies and counter-strategies. We describe a mathematical game-theoretic formalism and offer detailed analysis of a table-top cyber wargame executed at the US Army Research Laboratory. Our analysis shows how game-theoretical calculations can provide an effective tool for decision-making during cyber wargames.


title: Kinetic simulations of mildly relativistic shocks I: particle acceleration in high Mach number shocks

abstract: We use fully kinetic particle-in-cell simulations with unprecedentedly large transverse box sizes to study particle acceleration in weakly-magnetized mildly relativistic shocks traveling at a velocity $\approx 0.75c$ and a Mach number of 15. We examine both subluminal (quasi-parallel) and superluminal (quasi-perpendicular) magnetic field orientations. We find that quasi-parallel shocks are mediated by a filamentary non-resonant (Bell) instability driven by non-thermal ions, producing magnetic fluctuations on scales comparable to the ion gyro-radius. In quasi-parallel shocks, both electrons and ions are accelerated into non-thermal power-laws whose maximum energy grows linearly with time. The upstream heating of electrons is small, and the two species enter the shock front in rough thermal equilibrium. The shock's structure is complex; the current of reflected non-thermal ions evacuates cavities in the upstream which form filaments of amplified magnetic fields once advected downstream. At late times, $10\%$ of the shock's energy goes into non-thermal protons and $\gtrsim10\%$ into magnetic fields. We find that properly capturing the magnetic turbulence driven by the non-thermal ions is important for properly measuring the energy fraction of non-thermal electrons, $\epsilon_e$. We find $\epsilon_e\sim 5\times10^{-4}$ for quasi-parallel shocks with $v=0.75c$, slightly larger than what was measured in simulations of non-relativistic shocks. In quasi-perpendicular shocks, no non-thermal power-law develops in ions or electrons. The ion acceleration efficiency in quasi-parallel shocks suggests that astrophysical objects that could host mildly relativistic quasi-parallel shocks -- for example, the jets of active galactic nuclei or microquasars -- may be important sources of cosmic rays and their secondaries, such as gamma-rays and neutrinos.


title: Controlling decoherence speed limit of a single impurity atom in a Bose-Einstein-condensate reservoir

abstract: We study the decoherence speed limit (DSL) of a single impurity atom immersed in a Bose-Einsteincondensed (BEC) reservoir when the impurity atom is in a double-well potential. We demonstrate how the DSL of the impurity atom can be manipulated by engineering the BEC reservoir and the impurity potential within experimentally realistic limits. We show that the DSL can be controlled by changing key parameters such as the condensate scattering length, the effective dimension of the BEC reservoir, and the spatial configuration of the double-well potential imposed on the impurity. We uncover the physical mechanisms of controlling the DSL at root of the spectral density of the BEC reservoir.


title: Using Deep Reinforcement Learning to Learn High-Level Policies on the ATRIAS Biped

abstract: Learning controllers for bipedal robots is a challenging problem, often requiring expert knowledge and extensive tuning of parameters that vary in different situations. Recently, deep reinforcement learning has shown promise at automatically learning controllers for complex systems in simulation. This has been followed by a push towards learning controllers that can be transferred between simulation and hardware, primarily with the use of domain randomization. However, domain randomization can make the problem of finding stable controllers even more challenging, especially for underactuated bipedal robots. In this work, we explore whether policies learned in simulation can be transferred to hardware with the use of high-fidelity simulators and structured controllers. We learn a neural network policy which is a part of a more structured controller. While the neural network is learned in simulation, the rest of the controller stays fixed, and can be tuned by the expert as needed. We show that using this approach can greatly speed up the rate of learning in simulation, as well as enable transfer of policies between simulation and hardware. We present our results on an ATRIAS robot and explore the effect of action spaces and cost functions on the rate of transfer between simulation and hardware. Our results show that structured policies can indeed be learned in simulation and implemented on hardware successfully. This has several advantages, as the structure preserves the intuitive nature of the policy, and the neural network improves the performance of the hand-designed policy. In this way, we propose a way of using neural networks to improve expert designed controllers, while maintaining ease of understanding.


title: Constraining small scale magnetic fields through plasma lensing: Application to the Black widow eclipsing pulsar binary

abstract: In regions with strongly varying electron density, radio emission can be magnified significantly by plasma lensing. In the presence of magnetic fields, magnification in time and frequency will be different for two circular polarizations. We show how these effects can be used to measure or constrain the magnetic field parallel to the line of sight, $B_\parallel$, as well as its spatial structure, $\sigma_{B_\parallel}$, in the lensing region. In addition, we discuss how generalized Faraday rotation can constrain the strength of the perpendicular field, $B_\perp$. We attempt to make such measurements for the Black Widow pulsar, PSR~B1957+20, in which plasma lensing was recently discovered. For this system, pressure equilibrium suggests $B\gtrsim 20\,$G at the interface between the pulsar and companion winds, where the radio eclipse starts and ends, and where most lensing occurs. We find no evidence for large-scale magnetic fields, with, on average, $B_\parallel=0.02\pm0.09\,$G over the egress lensing region. From individual lensing events, we strongly constrain small scale magnetic structure to $\sigma_B<10\,$mG, thus excluding scenarios with a strong but rapidly varying field. Finally, from the lack of reduction of average circular polarization in the same region, we rule out a strong, quasi-transverse field. We cannot identify any plausible scenario in which a large magnetic field in this system is concealed, leaving the nature of the interface between the pulsar and companion winds an enigma. Our method can be applied to other sources showing plasma lensing, including other eclipsing pulsars and fast radio bursts, to study the local properties of the magnetic field.


title: Robin's inequality for 20-free integers

abstract: In 1984, Robin showed that the Riemann Hypothesis for $\zeta$ is equivalent to demonstrating $\sigma(n) < e^\gamma n \log \log n$ for all $n > 5040$. Robin's inequality has since been proven for various infinite families of power-free integers: $5$-free integers, $7$-free integers, and $11$-free integers. We extend these results to cover $20$-free integers.


title: Harmonic maps and biharmonic Riemanian submersions

abstract: Characterizations for Riemannian submersions to be harmonic or biharmonic are shown. Examples of biharmonic but not harmonic Riemannian submersions are shown.


title: Asymptotic behavior of the principal eigenvalue of a linear second order elliptic operator with small/large diffusion coefficient and its application

abstract: In this article, we are concerned with the following eigenvalue problem of a linear second order elliptic operator: \begin{equation} \nonumber -D\Delta \phi -2\alpha\nabla m(x)\cdot \nabla\phi+V(x)\phi=\lambda\phi\ \ \hbox{ in }\Omega, \end{equation} complemented by a general boundary condition including Dirichlet boundary condition and Robin boundary condition: $$ \frac{\partial\phi}{\partial n}+\beta(x)\phi=0 \ \ \hbox{ on }\partial\Omega, $$ where $\beta\in C(\partial\Omega)$ allows to be positive, sign-changing or negative, and $n(x)$ is the unit exterior normal to $\partial\Omega$ at $x$. The domain $\Omega\subset\mathbb{R}^N$ is bounded and smooth, the constants $D>0$ and $\alpha>0$ are, respectively, the diffusive and advection coefficients, and $m\in C^2(\bar\Omega),\,V\in C(\bar\Omega)$ are given functions. We aim to investigate the asymptotic behavior of the principal eigenvalue of the above eigenvalue problem as the diffusive coefficient $D\to0$ or $D\to\infty$. Our results, together with those of \cite{CL2,DF,Fr} where the Nuemann boundary case (i.e., $\beta=0$ on $\partial\Omega$) and Dirichlet boundary case were studied, reveal the important effect of advection and boundary conditions on the asymptotic behavior of the principal eigenvalue. We also apply our results to a reaction-diffusion-advection equation which is used to describe the evolution of a single species living in a heterogeneous stream environment and show some interesting behaviors of the species persistence and extinction caused by the buffer zone and small/large diffusion rate.


title: Generative Adversarial Active Learning for Unsupervised Outlier Detection

abstract: Outlier detection is an important topic in machine learning and has been used in a wide range of applications. In this paper, we approach outlier detection as a binary-classification issue by sampling potential outliers from a uniform reference distribution. However, due to the sparsity of data in high-dimensional space, a limited number of potential outliers may fail to provide sufficient information to assist the classifier in describing a boundary that can separate outliers from normal data effectively. To address this, we propose a novel Single-Objective Generative Adversarial Active Learning (SO-GAAL) method for outlier detection, which can directly generate informative potential outliers based on the mini-max game between a generator and a discriminator. Moreover, to prevent the generator from falling into the mode collapsing problem, the stop node of training should be determined when SO-GAAL is able to provide sufficient information. But without any prior information, it is extremely difficult for SO-GAAL. Therefore, we expand the network structure of SO-GAAL from a single generator to multiple generators with different objectives (MO-GAAL), which can generate a reasonable reference distribution for the whole dataset. We empirically compare the proposed approach with several state-of-the-art outlier detection methods on both synthetic and real-world datasets. The results show that MO-GAAL outperforms its competitors in the majority of cases, especially for datasets with various cluster types or high irrelevant variable ratio.


title: Harmonic morphisms of Allof-Wallach spaces of positive curvature

abstract: An infinite family of distinct harmonic morphisms with minimal circle fibers from the 7-dimensional homogeneous Allof-Wallach spaces of positive curvature onto the 6-dimensional flag manifolds is given.


title: Learning Confidence Sets using Support Vector Machines

abstract: The goal of confidence-set learning in the binary classification setting is to construct two sets, each with a specific probability guarantee to cover a class. An observation outside the overlap of the two sets is deemed to be from one of the two classes, while the overlap is an ambiguity region which could belong to either class. Instead of plug-in approaches, we propose a support vector classifier to construct confidence sets in a flexible manner. Theoretically, we show that the proposed learner can control the non-coverage rates and minimize the ambiguity with high probability. Efficient algorithms are developed and numerical studies illustrate the effectiveness of the proposed method.


title: Novel Multi Agent Models for Chemical Self-assembly

abstract: The chemical self-assembly has been considered as one of most important scientific problems in the 21th Century; however, since the process of self-assembly is very complex, there is few mathematic theory for it currently. This paper provides a novel multi-agent model for chemical self-assembly, where the interaction between agents adopts the classic Lennard-Jones potential. Under this model, we propose an optimal problem by taking the temperature as the control input, and choosing the internal energy as the optimal object. A numerical solution for our optimal problem is also developed. Simulations show that our control scheme can improve the product of self-assembly. Further more, we give a strict analysis for the self-assembly model without noise, which corresponds to an attraction-repulsion multi-agent system, and prove it converges to a stable configuration eventually.


title: Inverse Transport Networks

abstract: We introduce inverse transport networks as a learning architecture for inverse rendering problems where, given input image measurements, we seek to infer physical scene parameters such as shape, material, and illumination. During training, these networks are evaluated not only in terms of how close they can predict groundtruth parameters, but also in terms of whether the parameters they produce can be used, together with physically-accurate graphics renderers, to reproduce the input image measurements. To en- able training of inverse transport networks using stochastic gradient descent, we additionally create a general-purpose, physically-accurate differentiable renderer, which can be used to estimate derivatives of images with respect to arbitrary physical scene parameters. Our experiments demonstrate that inverse transport networks can be trained efficiently using differentiable rendering, and that they generalize to scenes with completely unseen geometry and illumination better than networks trained without appearance- matching regularization.


title: Boundary-guided Feature Aggregation Network for Salient Object Detection

abstract: Fully convolutional networks (FCN) has significantly improved the performance of many pixel-labeling tasks, such as semantic segmentation and depth estimation. However, it still remains non-trivial to thoroughly utilize the multi-level convolutional feature maps and boundary information for salient object detection. In this paper, we propose a novel FCN framework to integrate multi-level convolutional features recurrently with the guidance of object boundary information. First, a deep convolutional network is used to extract multi-level feature maps and separately aggregate them into multiple resolutions, which can be used to generate coarse saliency maps. Meanwhile, another boundary information extraction branch is proposed to generate boundary features. Finally, an attention-based feature fusion module is designed to fuse boundary information into salient regions to achieve accurate boundary inference and semantic enhancement. The final saliency maps are the combination of the predicted boundary maps and integrated saliency maps, which are more closer to the ground truths. Experiments and analysis on four large-scale benchmarks verify that our framework achieves new state-of-the-art results.


title: Linear and Nonlinear Rheology and Structural Relaxation in Dense Glassy and Jammed Soft Repulsive Microgel Suspensions

abstract: We present an integrated experimental and quantitative theoretical study of the mechanics of self-crosslinked, neutral, repulsive pNIPAM microgel suspensions over concentration (c) range spanning the fluid, glassy and putative "soft jammed" regimes. In the glassy regime we measure a linear elastic dynamic shear modulus over 3 decades which follows an apparent power law concentration dependence G'~$c^{5.64}$, followed by a sharp crossover to a nearly linear growth at high concentrations. We formulate a theoretical approach to address all three regimes within a single conceptual Brownian dynamics framework. A minimalist single particle description is constructed that allows microgel size to vary with concentration due to steric de-swelling effects. Using a Hertzian repulsion interparticle potential and a suite of statistical mechanical theories, quantitative predictions under quiescent conditions of microgel collective structure, dynamic localization length, elastic modulus, and the structural relaxation time are made. Based on a constant inter-particle repulsion strength parameter which is determined by requiring the theory to reproduce the linear elastic shear modulus over the entire concentration regime, we demonstrate good agreement between theory and experiment. Theoretical predictions of how quiescent structural relaxation time changes under deformation, and how the yield stress and strain change as a function of concentration has been made. Reasonable agreement with our observations is obtained. To the best of our knowledge, this is the first attempt to quantitatively understand structure, quiescent relaxation and shear elasticity, and nonlinear yielding of dense microgel suspensions using microscopic force based theoretical methods that include activated hopping processes. We expect our approach will be useful for other soft polymeric particle suspensions in the core-shell family.


title: Form factors for pions couplings to constituent quarks under weak magnetic field

abstract: By considering a one loop background field method for a nonperturbative one gluon exchange between quarks, leading electromagnetic form factors of pion couplings to constituent quarks are derived by means of from a large quark and gluon effective masses expansion. With this calculation the leading anisotropic corrections induced by a weak magnetic field for the pion-constituent quark form factors are obtained. Besides that, few effective couplings which emerge only due to the weak magnetic field that break chiral and isospin symmetry are also found. Numerical estimations for these magnetic field corrections are presented for two different nonperturbative gluon propagators. All of these corrections are ultraviolet finite and their relative values are found to be of order of $(eB_0/{M^*}^2)^n$ being $n=1$ for the vector and axial pion couplings and $n=2$ for the pseudoscalar and scalar ones. The corresponding anisotropic corrections to the Strong averaged quadratic radii in the plane perpendicular to the magnetic field are also calculated as functions of the quark effective mass.


title: Helical Topological Edge States in a Quadrupole Phase

abstract: Topological electric quadrupole is a recently proposed concept that extends the theory of electric polarization of crystals to higher orders. Such a quadrupole phase supports topological states localized on both edges and corners. In this work, we show that in a quadrupole phase of honeycomb lattice, topological helical edge states and pseudo-spin-polarized corner states appear by making use of a pseudo-spin degree of freedom related to point group symmetry. Furthermore, we argue that a general condition for emergence of helical edge states in a (pseudo-)spinful quadrupole phase is mirror or time-reversal symmetry. Our results offers a way of generating topological helical edge states without spin-orbital couplings.


title: Systems of Polynomials with At Least One Positive Real Zero

abstract: In this paper, we prove several theorems on systems of polynomials with at least one positive real zero based on the theory of conceive polynomials. These theorems provide sufficient conditions for systems of multivariate polynomials admitting at least one positive real zero in terms of their Newton polytopes and combinatorial structure. Moreover, a class of polynomials attaining their global minimums in the first quadrant are given, which is useful in polynomial optimization.


title: Nonlinear nonlocal multicontinua upscaling framework and its applications

abstract: In this paper, we discuss multiscale methods for nonlinear problems. The main idea of these approaches is to use local constraints and solve problems in oversampled regions for constructing macroscopic equations. These techniques are intended for problems without scale separation and high contrast, which often occur in applications. For linear problems, the local solutions with constraints are used as basis functions. This technique is called Constraint Energy Minimizing Generalized Multiscale Finite Element Method (CEM-GMsFEM). GMsFEM identifies macroscopic quantities based on rigorous analysis. In corresponding upscaling methods, the multiscale basis functions are selected such that the degrees of freedom have physical meanings, such as averages of the solution on each continuum. This paper extends the linear concepts to nonlinear problems, where the local problems are nonlinear. The main concept consists of: (1) identifying macroscopic quantities; (2) constructing appropriate oversampled local problems with coarse-grid constraints; (3) formulating macroscopic equations. We consider two types of approaches. In the first approach, the solutions of local problems are used as basis functions (in a linear fashion) to solve nonlinear problems. This approach is simple to implement; however, it lacks the nonlinear interpolation, which we present in our second approach. In this approach, the local solutions are used as a nonlinear forward map from local averages (constraints) of the solution in oversampling region. This local fine-grid solution is further used to formulate the coarse-grid problem. Both approaches are discussed on several examples and applied to single-phase and two-phase flow problems, which are challenging because of convection-dominated nature of the concentration equation.


title: Weak detection in the spiked Wigner model

abstract: We consider the weak detection problem in a rank-one spiked Wigner data matrix where the signal-to-noise ratio is small so that reliable detection is impossible. We propose a hypothesis test on the presence of the signal by utilizing the linear spectral statistics of the data matrix. The test is data-driven and does not require prior knowledge about the distribution of the signal or the noise. When the noise is Gaussian, the proposed test is optimal in the sense that its error matches that of the likelihood ratio test, which minimizes the sum of the Type-I and Type-II errors. If the density of the noise is known and non-Gaussian, the error of the test can be lowered by applying an entrywise transformation to the data matrix. We establish a central limit theorem for the linear spectral statistics of general rank-one spiked Wigner matrices as an intermediate step.


title: Standing lattice solitons in the discrete NLS equation with saturation

abstract: We consider standing lattice solitons for discrete nonlinear Schrodinger equation with saturation (NLSS), where so-called transparent points were recently discovered. These transparent points are the values of the governing parameter (e.g., the lattice spacing) for which the Peierls-Nabarro barrier vanishes. In order to explain the existence of transparent points, we study a solitary wave solution in the continuous NLSS and analyse the singularities of its analytic continuation in the complex plane. The existence of a quadruplet of logarithmic singularities nearest to the real axis is proven and applied to two settings: (i) the fourth-order differential equation arising as the next-order continuum approximation of the discrete NLSS and (ii) the advance-delay version of the discrete NLSS.


title: A theoretical framework for deep locally connected ReLU network

abstract: Understanding theoretical properties of deep and locally connected nonlinear network, such as deep convolutional neural network (DCNN), is still a hard problem despite its empirical success. In this paper, we propose a novel theoretical framework for such networks with ReLU nonlinearity. The framework explicitly formulates data distribution, favors disentangled representations and is compatible with common regularization techniques such as Batch Norm. The framework is built upon teacher-student setting, by expanding the student forward/backward propagation onto the teacher's computational graph. The resulting model does not impose unrealistic assumptions (e.g., Gaussian inputs, independence of activation, etc). Our framework could help facilitate theoretical analysis of many practical issues, e.g. overfitting, generalization, disentangled representations in deep networks.


title: Throughput Optimization in FDD MU-MISO Wireless Powered Communication Networks

abstract: In this paper, we consider a frequency-division duplexing (FDD) multiple-user multiple-input-single-output (MU-MISO) wireless-powered communication network (WPCN) consisting of one hybrid data-and-energy access point (HAP) with multiple antennas which coordinates energy/information transfer to/from several single-antenna wireless devices (WD). Typically, in such a system, wireless energy transfer (WET) requires such techniques as energy beamforming (EB) for efficient transfer of energy to the WDs. Yet, efficient EB can only be accomplished if channel state information (CSI) is available to the transmitter, which, in FDD systems is only achieved through uplink (UL) feedback. Therefore, while in our scheme we use the downlink (DL) channels for WET only, the UL channel frames are split into two phases: the CSI feedback phase during which the WDs feed CSI back to the HAP and the WIT phase where the HAP performs wireless information transmission (WIT) via space-division-multiple-access (SDMA). To ensure rate fairness among the WDs, this paper maximizes the minimum WIT data rate among the WDs. Using an iterative solution, the original optimization problem can be relaxed into two sub-problems whose convexity conditions are derived. Finally, the behavior of this system when the number of HAP antennas increases is analyzed. Simulation results verify the truthfulness of our analysis.


title: OGLE-2017-BLG-0039: Microlensing Event with Light from the Lens Identified from Mass Measurement

abstract: We present the analysis of the caustic-crossing binary microlensing event OGLE-2017-BLG-0039. Thanks to the very long duration of the event, with an event time scale $t_{\rm E}\sim 130$ days, the microlens parallax is precisely measured despite its small value of $\pie\sim 0.06$. The analysis of the well-resolved caustic crossings during both the source star's entrance and exit of the caustic yields the angular Einstein radius $\thetae\sim 0.6$~mas. The measured $\pie$ and $\thetae$ indicate that the lens is a binary composed of two stars with masses $\sim 1.0~M_\odot$ and $\sim 0.15~M_\odot$, and it is located at a distance of $\sim 6$ kpc. From the color and brightness of the lens estimated from the determined lens mass and distance, it is expected that $\sim 2/3$ of the $I$-band blended flux comes from the lens. Therefore, the event is a rare case of a bright lens event for which high-resolution follow-up observations can confirm the nature of the lens.


title: Nilmanifolds and their associated non local fields

abstract: For six dimensional nilmanifolds we build a module $\mathcal{H}$ of an affine Kac Moody vertex algebras. Then, we associate some logarithmic fields for the module $\mathcal{H}$ and we study their singularities. We also presented a physics motivation behind this construction. We study a particular case, we show that when the nilmanifold $N$ is a $k$ degree $S^1$--fibration over the two torus and a choice of $l \in \mathbb{Z} \simeq H^3(N, \mathbb{Z})$ the fields associated to the space $\mathcal{H}$ have tri-logarithm singularities whenever $kl \neq 0$.


title: Limits on Electroweak Instanton-Induced Processes with Multiple Boson Production

abstract: Recently, the CMS collaboration has reported their search for electroweak instanton-like processes with anomalous $B+L$ violation assuming multi-fermion but zero-boson final states. On the other hand, many theoretical studies suggest that anomalous $B+L$ processes may have an observably large production rate only if their final state contains a large number of electroweak gauge bosons. In this paper, we compare collider signatures of zero- and multi-boson events of anomalous $B+L$ violation at the LHC and derive an upper limit on the cross-section for the multi-boson process by recasting the CMS analysis.


title: Spatial-Temporal Inference of Urban Traffic Emissions Based on Taxi Trajectories and Multi-Source Urban Data

abstract: Vehicle trajectory data collected via GPS-enabled devices have played increasingly important roles in estimating network-wide traffic, given their broad spatial-temporal coverage and representativeness of traffic dynamics. This paper exploits taxi GPS data, license plate recognition (LPR) data and geographical information for reconstructing the spatial and temporal patterns of urban traffic emissions. Vehicle emission factor models are employed to estimate emissions based on taxi trajectories. The estimated emissions are then mapped to spatial grids of urban areas to account for spatial heterogeneity. To extrapolate emissions from the taxi fleet to the whole vehicle population, we use Gaussian process regression models supported by geographical features to estimate the spatially heterogeneous traffic volume and fleet composition. Unlike previous studies, this paper utilizes the taxi GPS data and LPR data to disaggregate vehicle and emission characteristics through space and time in a large-scale urban network. The results of a case study in Hangzhou, China, reveal high-resolution spatiotemporal patterns of traffic flows and emissions and identify emission hotspots at different locations. This study provides an accessible means of inferring the environmental impact of urban traffic with multi-source data that are now widely available in urban areas.


title: Embedded-State Latent Conditional Random Fields for Sequence Labeling

abstract: Complex textual information extraction tasks are often posed as sequence labeling or \emph{shallow parsing}, where fields are extracted using local labels made consistent through probabilistic inference in a graphical model with constrained transitions. Recently, it has become common to locally parametrize these models using rich features extracted by recurrent neural networks (such as LSTM), while enforcing consistent outputs through a simple linear-chain model, representing Markovian dependencies between successive labels. However, the simple graphical model structure belies the often complex non-local constraints between output labels. For example, many fields, such as a first name, can only occur a fixed number of times, or in the presence of other fields. While RNNs have provided increasingly powerful context-aware local features for sequence tagging, they have yet to be integrated with a global graphical model of similar expressivity in the output distribution. Our model goes beyond the linear chain CRF to incorporate multiple hidden states per output label, but parametrizes their transitions parsimoniously with low-rank log-potential scoring matrices, effectively learning an embedding space for hidden states. This augmented latent space of inference variables complements the rich feature representation of the RNN, and allows exact global inference obeying complex, learned non-local output constraints. We experiment with several datasets and show that the model outperforms baseline CRF+RNN models when global output constraints are necessary at inference-time, and explore the interpretable latent structure.


title: Positivities of knots and links and the defect of Bennequin inequality

abstract: We discuss relations among various positivities of knots and links, such as strong quasipositivity and quasipositivity. We give several pieces of supporting evidence for conjectural statements concerning these positivities and the defect of Bennequin inequality. Finally, we determine strong quasipositivity and quasipositivity for knots up to 12 crossings (with two exceptions for quasipositivity).


title: Morphology of soft and slippery contact via fluid drainage

abstract: The dynamic of contact formation between soft materials immersed in a fluid is accompanied by fluid drainage and elastic deformation. As a result, controlling the coupling between lubrication pressure and elasticity provides strategies to design materials with reversible and dynamic adhesion to wet or flooded surfaces. We characterize the elastic deformation of a soft coating with nanometer-scale roughness as it approaches and contacts a rigid surface in a fluid environment. The lubrication pressure during the approach causes elastic deformation and prevents contact formation. We observe deformation profiles that are drastically different from those observed for elastic half-space when the thickness of the soft coating is comparable to the hydrodynamic radius. In contrast, we show that surface roughness favors fluid drainage without altering the elastic deformation. As a result, the coupling between elasticity and slip (caused by surface roughness) can lead to trapped fluid pockets in the contact region.


title: Model confidence sets and forecast combination: An application to age-specific mortality

abstract: Model averaging combines forecasts obtained from a range of models, and it often produces more accurate forecasts than a forecast from a single model. The crucial part of forecast accuracy improvement in using the model averaging lies in the determination of optimal weights from a finite sample. If the weights are selected sub-optimally, this can affect the accuracy of the model-averaged forecasts. Instead of choosing the optimal weights, we consider trimming a set of models before equally averaging forecasts from the selected superior models. Motivated by Hansen, Lunde and Nason (2011), we apply and evaluate the model confidence set procedure when combining mortality forecasts. The proposed model averaging procedure is motivated by Samuels and Sekkel (2017) based on the concept of model confidence sets as proposed by Hansen et al. (2011) that incorporates the statistical significance of the forecasting performance. As the model confidence level increases, the set of superior models generally decreases. The proposed model averaging procedure is demonstrated via national and sub-national Japanese mortality for retirement ages between 60 and 100+. Illustrated by national and sub-national Japanese mortality for ages between 60 and 100+, the proposed model-average procedure gives the smallest interval forecast errors, especially for males. We find that robust out-of-sample point and interval forecasts may be obtained from the trimming method. By robust, we mean robustness against model misspecification.


title: Cocktail BPSK: Cross Power Utilization for High Data Rates

abstract: This paper proposes a novel transmission strategy, referred to as cocktail BPSK, where two independent BPSK symbols are layered with various weights to be transmitted and, capitalizing on the cross power utilization, will be demodulated at the receiver without interference from each other. To evaluate the performance of the proposed scheme, its achievable data rate is formulated over additive white Gaussian noise (AWGN) channels. Based on the theoretical analysis, numerical results are provided for the performance comparisons between the proposed scheme and conventional transmission schemes, which substantiate the validity of the proposed scheme. Specifically, when the signal-to-noise power ratio (SNR) is small, the achievable data rate of the proposed scheme outperforms the channel capacity achieved by Gaussian-distributed inputs.


title: Characterization and Performance of PADME's Cherenkov-Based Small-Angle Calorimeter

abstract: The PADME experiment, at the Laboratori Nazionali di Frascati (LNF), in Italy, will search for invisible decays of the hypothetical dark photon via the process $e^+e^-\rightarrow \gamma A'$, where the $A'$ escapes detection. The dark photon mass range sensitivity in a first phase will be 1 to 24 MeV. We report here on measurement and simulation studies of the performance of the Small-Angle Calorimeter, a component of PADME's detector dedicated to rejecting 2- and 3-gamma backgrounds. The crucial requirement is a timing resolution of less than 200 ps, which is satisfied by the choice of PbF$_2$ crystals and the newly released Hamamatsu R13478UV photomultiplier tubes (PMTs). We find a timing resolution of 81 ps (with double-peak separation resolution of 1.8 ns) and a single-crystal energy resolution of 5.7%/$\sqrt{E}$ with light yield of 2.07 photo-electrons per MeV, using 100 to 400 MeV electrons at the Beam Test Facility of LNF. We also propose the investigation of a two-PMT solution coupled to a single PbF$_2$ crystal for higher-energy applications, which has potentially attractive features.


title: What privacy concerns do parents have about children's mobile apps, and how can they stay SHARP?

abstract: Tablet computers are widely used by young children. A report in 2016 shows that children aged 5 to 15 years are spending more time online than watching TV. A 2017 update of the same report shows that parents are becoming more concerned about their children's online risks compared to the previous year. Parents are working hard to protect their children's online safety. An increasing number of parents are setting up content filtering at home or having regular discussions with their children regarding online risks. However, although risks related to Social Media platforms or social video sharing sites (like YouTube) are widely known, risks posed by mobile applications or games (i.e. `apps') are less known. Behind the cute characters, apps used by children can not only have the possibility of exposing them to age-inappropriate content or excessive in-app promotions, but may also make a large amount of their personal information accessible to third-party online marketing and advertising industry. Such practices are not unique to children's apps, but young children are probably less capable of resisting the resulting personalised advertisements and game promotions. In this report, we present findings from our online survey of 220 parents with children aged 6-10, mainly from the U.K. and other western countries, regarding their privacy concerns and expectations of their children's use of mobile apps. Parents play a key role in children's use of digital technology, especially for children under 10 years old. Recent reports have highlighted parents' lack of sufficient support for choosing appropriate digital content for their children. Our report sheds some initial light on parents' key struggles and points to immediate steps and possible areas of future development.


title: Learning and Planning with a Semantic Model

abstract: Building deep reinforcement learning agents that can generalize and adapt to unseen environments remains a fundamental challenge for AI. This paper describes progresses on this challenge in the context of man-made environments, which are visually diverse but contain intrinsic semantic regularities. We propose a hybrid model-based and model-free approach, LEArning and Planning with Semantics (LEAPS), consisting of a multi-target sub-policy that acts on visual inputs, and a Bayesian model over semantic structures. When placed in an unseen environment, the agent plans with the semantic model to make high-level decisions, proposes the next sub-target for the sub-policy to execute, and updates the semantic model based on new observations. We perform experiments in visual navigation tasks using House3D, a 3D environment that contains diverse human-designed indoor scenes with real-world objects. LEAPS outperforms strong baselines that do not explicitly plan using the semantic content.


title: Heegaard Floer invariants of contact structures on links of surface singularities

abstract: Let a contact 3-manifold $(Y, \xi_0)$ be the link of a normal surface singularity equipped with its canonical contact structure $\xi_0$. We prove a special property of such contact 3-manifolds of "algebraic" origin: the Heegaard Floer invariant $c^+(\xi_0)\in HF^+(-Y)$ cannot lie in the image of the $U$-action on $HF^+(-Y)$. It follows that Karakurt's "height of $U$-tower" invariants are always 0 for canonical contact structures on singularity links, which contrasts the fact that the height of $U$-tower can be arbitrary for general fillable contact structures. Our proof uses the interplay between the Heegaard Floer homology and N\'emethi's lattice cohomology.


title: Proximal Recursion for Solving the Fokker-Planck Equation

abstract: We develop a new method to solve the Fokker-Planck or Kolmogorov's forward equation that governs the time evolution of the joint probability density function of a continuous-time stochastic nonlinear system. Numerical solution of this equation is fundamental for propagating the effect of initial condition, parametric and forcing uncertainties through a nonlinear dynamical system, and has applications encompassing but not limited to forecasting, risk assessment, nonlinear filtering and stochastic control. Our methodology breaks away from the traditional approach of spatial discretization for solving this second-order partial differential equation (PDE), which in general, suffers from the "curse-of-dimensionality". Instead, we numerically solve an infinite dimensional proximal recursion in the space of probability density functions, which is theoretically equivalent to solving the Fokker-Planck-Kolmogorov PDE. We show that the dual formulation along with the introduction of an entropic regularization, leads to a smooth convex optimization problem that can be implemented via suitable block co-ordinate iteration and has fast convergence due to certain contraction property that we establish. This approach enables meshless implementation leading to remarkably fast computation.


title: UVM Based Reusable Verification IP for Wishbone Compliant SPI Master Core

abstract: The System on Chip design industry relies heavily on functional verification to ensure that the designs are bug-free. As design engineers are coming up with increasingly dense chips with much functionality, the functional verification field has advanced to provide modern verification techniques. In this paper, we present verification of a wishbone compliant Serial Peripheral Interface (SPI) Master core using a System Verilog based standard verification methodology, the Universal Verification Methodology (UVM). The reason for using UVM factory pattern with parameterized classes is to develop a robust and reusable verification IP. SPI is a full duplex communication protocol used to interface components most likely in embedded systems. We have verified an SPI Master IP core design that is wishbone compliant and compatible with SPI protocol and bus and furnished the results of our verification. We have used QuestaSim for simulation and analysis of waveforms, Integrated Metrics Center, Cadence for coverage analysis. We also propose interesting future directions for this work in developing reliable systems.


title: Comparison of Self-Aware and Organic Computing Systems

abstract: With increasing complexity and heterogeneity of computing devices, it has become crucial for system to be autonomous, adaptive to dynamic environment, robust, flexible, and having so called self-*properties. These autonomous systems are called organic computing(OC) systems. OC system was proposed as a solution to tackle complex systems. Design time decisions have been shifted to run time in highly complex and interconnected systems as it is very hard to consider all scenarios and their appropriate actions in advance. Consequently, Self-awareness becomes crucial for these adaptive autonomous systems. To cope with evolving environment and changing user needs, system need to have knowledge about itself and its surroundings. Literature review shows that for autonomous and intelligent systems, researchers are concerned about knowledge acquisition, representation and learning which is necessary for a system to adapt. This paper is written to compare self-awareness and organic computing by discussing their definitions, properties and architecture.


title: Using Multi-task and Transfer Learning to Solve Working Memory Tasks

abstract: We propose a new architecture called Memory-Augmented Encoder-Solver (MAES) that enables transfer learning to solve complex working memory tasks adapted from cognitive psychology. It uses dual recurrent neural network controllers, inside the encoder and solver, respectively, that interface with a shared memory module and is completely differentiable. We study different types of encoders in a systematic manner and demonstrate a unique advantage of multi-task learning in obtaining the best possible encoder. We show by extensive experimentation that the trained MAES models achieve task-size generalization, i.e., they are capable of handling sequential inputs 50 times longer than seen during training, with appropriately large memory modules. We demonstrate that the performance achieved by MAES far outperforms existing and well-known models such as the LSTM, NTM and DNC on the entire suite of tasks.


title: A New Sparse SOS Decomposition Algorithm Based on Term Sparsity

abstract: A new sparse SOS decomposition algorithm is proposed based on a new sparsity pattern, called cross sparsity patterns. The new sparsity pattern focuses on the sparsity of terms and thus is different from the well-known correlative sparsity pattern which focuses on the sparsity of variables though the sparse SOS decomposition algorithms based on these two sparsity patterns both take use of chordal extensions/chordal decompositions. Moreover, it is proved that the SOS decomposition obtained by the new sparsity pattern is always a refinement of the block-diagonalization obtained by the sign-symmetry method. %Because the new sparsity pattern covers more sparse polynomials than correlative sparsity pattern, Various experiments show that the new algorithm dramatically saves the computational cost compared to existing tools and can handle some really huge polynomials.


title: Global or Local Pure-Condensible Atmospheres: Importance of Horizontal Latent Heat Transport

abstract: The distribution of a pure condensible planetary atmosphere in equilibrium with a surface reservoir is revisited employing the energy budget of the climate system, emphasizing the atmospheric horizontal latent heat transport. This configuration is applicable to icy Solar System bodies such as Triton as well as a range of possible exoplanet atmospheres, including water or $\mathrm{CO_2}$ iceballs or ocean worlds, and lava planets with mineral vapor atmospheres. Climate regimes for slowly rotating planets with the hot-spot near the substellar point, and for rapidly rotating planets with a warm equatorial belt, are both treated. A non-dimensional parameter controlling the fractional variation of the surface pressure is derived; it measures whether the pure condensible atmosphere is global or localized. The global pure condensible atmosphere with the non-dimensional parameter much less than order of unity is maintained by the strong horizontal latent heat transport associated with an "evaporation/sublimation-driven flow" from warm to cold places that compensates for the incoming differential radiative forcing. We show that the variation of surface temperature can be estimated in terms of this non-dimensional parameter if it is not too large. In the case of a pure water-vapor atmosphere with an ice or liquid surface, we show that the atmosphere is thick enough to maintain nearly isothermal surface conditions even when the substellar surface temperature is around the freezing point. Finally, it is proposed that the evaporation/sublimation-driven flow regime for global atmospheres could be detected via its effect on the inhomogeneous distribution of minor non-condensible components in the atmosphere.


title: Elastic deformation during dynamic force measurements in viscous fluids

abstract: Understanding and harnessing the coupling between lubrication pressure and elasticity provides materials design strategies for applications such as adhesives, coatings, microsensors, and biomaterials. Elastic deformation of compliant solids caused by viscous forces can also occur during dynamic force measurements in instruments such as the surface forces apparatus (SFA) or the atomic force microscope (AFM). We briefly review hydrodynamic interactions in the presence of soft, deformable interfaces in the lubrication limit. More specifically, we consider the scenario of two surfaces approaching each other in a viscous fluid where one or both surfaces is deformable, which is also relevant to many force measurement systems. In this article the basic theoretical background of the elastohydrodynamic problem is detailed, followed by a discussion of experimental validation and considerations, especially for the role of elastic deformation on surface forces measurements. Finally, current challenges to our understanding of soft hydrodynamic interactions, such as the consideration of substrate layering, poroelasticity, viscoelasticity, surface heterogeneity, as well as their implications are discussed.


title: Encoding Robust Representation for Graph Generation

abstract: Generative networks have made it possible to generate meaningful signals such as images and texts from simple noise. Recently, generative methods based on GAN and VAE were developed for graphs and graph signals. However, the mathematical properties of these methods are unclear, and training good generative models is difficult. This work proposes a graph generation model that uses a recent adaptation of Mallat's scattering transform to graphs. The proposed model is naturally composed of an encoder and a decoder. The encoder is a Gaussianized graph scattering transform, which is robust to signal and graph manipulation. The decoder is a simple fully connected network that is adapted to specific tasks, such as link prediction, signal generation on graphs and full graph and signal generation. The training of our proposed system is efficient since it is only applied to the decoder and the hardware requirements are moderate. Numerical results demonstrate state-of-the-art performance of the proposed system for both link prediction and graph and signal generation.


title: Spin-orbit crossed susceptibility in topological Dirac semimetals

abstract: We theoretically study the spin-orbit crossed susceptibility of topological Dirac semimetals. Because of strong spin-orbit coupling, the orbital motion of electrons is modulated by Zeeman coupling, which contributes to orbital magnetization. We find that the spin-orbit crossed susceptibility is proportional to the separation of the Dirac points and it is highly anisotropic. The orbital magnetization is induced only along the rotational symmetry axis. We also study the conventional spin susceptibility. The spin susceptibility exhibits anisotropy and the spin magnetization is induced only along the perpendicular to the rotational symmetry axis in contrast to the spin-orbit crossed susceptibility. We quantitatively compare the two susceptibilities and find that they can be comparable.


title: Adaptive Input Representations for Neural Language Modeling

abstract: We introduce adaptive input representations for neural language modeling which extend the adaptive softmax of Grave et al. (2017) to input representations of variable capacity. There are several choices on how to factorize the input and output layers, and whether to model words, characters or sub-word units. We perform a systematic comparison of popular choices for a self-attentional architecture. Our experiments show that models equipped with adaptive embeddings are more than twice as fast to train than the popular character input CNN while having a lower number of parameters. On the WikiText-103 benchmark we achieve 18.7 perplexity, an improvement of 10.5 perplexity compared to the previously best published result and on the Billion Word benchmark, we achieve 23.02 perplexity.


title: Semiclassical Lattice effects on interband tunneling of a two-state system

abstract: Previously, we have shown that the transition probability of the Landau-Zener problem in periodic lattice systems becomes large by taking into account the nonlinearity of the energy spectra, compared with the probability by the conventional Landau-Zener formula. The enhancement comes from the nonlinearity peculiar to the periodic lattice system, and this effect from the lattice on transition action cannot be neglected in the transition process. In the present paper, we first give a brief review of the previous work, and construct the transfer matrix of the Landau-Zener problem by the semiclassical description for lattice systems. Next, we study a ladder lattice system and show that the transition action obtains a phase due to the nonlinearity. Then, we consider the double-passage problem of the ladder system within the semiclassical description. We find the oscillation of the probability by the transition phase by the lattice effect. This phase comes from the semiclassical analysis unlike the Stokes phase, and we show that the oscillation is mainly contributed by the transition phase by the lattice effect, when the hybridization of the ladder is strong.


title: Minimax Lower Bounds for $\mathcal{H}_\infty$-Norm Estimation

abstract: The problem of estimating the $\mathcal{H}_\infty$-norm of an LTI system from noisy input/output measurements has attracted recent attention as an alternative to parameter identification for bounding unmodeled dynamics in robust control. In this paper, we study lower bounds for $\mathcal{H}_\infty$-norm estimation under a query model where at each iteration the algorithm chooses a bounded input signal and receives the response of the chosen signal corrupted by white noise. We prove that when the underlying system is an FIR filter, $\mathcal{H}_\infty$-norm estimation is no more efficient than model identification for passive sampling. For active sampling, we show that norm estimation is at most a factor of $\log{r}$ more sample efficient than model identification, where $r$ is the length of the filter. We complement our theoretical results with experiments which demonstrate that a simple non-adaptive estimator of the norm is competitive with state-of-the-art adaptive norm estimation algorithms.


title: Answering Analytical Queries on Text Data with Temporal Term Histograms

abstract: Temporal text, i.e., time-stamped text data are found abundantly in a variety of data sources like newspapers, blogs and social media posts. While today's data management systems provide facilities for searching full-text data, they do not provide any simple primitives for performing analytical operations with text. This paper proposes the temporal term histograms (TTH) as an intermediate primitive that can be used for analytical tasks. We propose an algebra, with operators and equivalence rules for TTH and present a reference implementation on a relational database system.


title: Current control of systems with a Peierls distortion by magnetic field

abstract: We study the tunneling phenomenon of a ladder system with a Peierls distortion in a magnetic flux, and the response of electrons the insulator is investigated when the tunneling current flows on one-dimensional gapped chains along the external electric field. Without the magnetic field, the ladder system is insulated by the charge density wave order. Then, by the increase of the magnetic field, it becomes metallic with the disappearance of the distortion of the lattice. Finally, the gap appears, and it becomes a insulator. At the metallic state, the topological transition also occurs. To show this phenomenon, we consider the distortion by the phonon in the ladder model, and calculate the distortion gap and the transition probability by using both Landau-Zener formula and the instanton method. The transition to the metallic states will be applied to the current control by the magnetic field.


title: Efficiently testing local optimality and escaping saddles for ReLU networks

abstract: We provide a theoretical algorithm for checking local optimality and escaping saddles at nondifferentiable points of empirical risks of two-layer ReLU networks. Our algorithm receives any parameter value and returns: local minimum, second-order stationary point, or a strict descent direction. The presence of $M$ data points on the nondifferentiability of the ReLU divides the parameter space into at most $2^M$ regions, which makes analysis difficult. By exploiting polyhedral geometry, we reduce the total computation down to one convex quadratic program (QP) for each hidden node, $O(M)$ (in)equality tests, and one (or a few) nonconvex QP. For the last QP, we show that our specific problem can be solved efficiently, in spite of nonconvexity. In the benign case, we solve one equality constrained QP, and we prove that projected gradient descent solves it exponentially fast. In the bad case, we have to solve a few more inequality constrained QPs, but we prove that the time complexity is exponential only in the number of inequality constraints. Our experiments show that either benign case or bad case with very few inequality constraints occurs, implying that our algorithm is efficient in most cases.


title: Engineering long spin coherence times of spin-orbit systems

abstract: Spin-orbit coupling fundamentally alters spin qubits, opening pathways to improve the scalability of quantum computers via long distance coupling mediated by electric fields, photons, or phonons. It also allows for new engineered hybrid and topological quantum systems. However, spin qubits with intrinsic spin-orbit coupling are not yet viable for quantum technologies due to their short ($\sim1~\mu$s) coherence times $T_2$, while qubits with long $T_2$ have weak spin-orbit coupling making qubit coupling short-ranged and challenging for scale-up. Here we show that an intrinsic spin-orbit coupled "generalised spin" with total angular momentum $J=\tfrac{3}{2}$, which is defined by holes bound to boron dopant atoms in strained $^{28}\mathrm{Si}$, has $T_2$ rivalling the electron spins of donors and quantum dots in $^{28}\mathrm{Si}$. Using pulsed electron paramagnetic resonance, we obtain $0.9~\mathrm{ms}$ Hahn-echo and $9~\mathrm{ms}$ dynamical decoupling $T_2$ times, where strain plays a key role to reduce spin-lattice relaxation and the longitudinal electric coupling responsible for decoherence induced by electric field noise. Our analysis shows that transverse electric dipole can be exploited for electric manipulation and qubit coupling while maintaining a weak longitudinal coupling, a feature of $J=\tfrac{3}{2}$ atomic systems with a strain engineered quadrupole degree of freedom. These results establish single-atom hole spins in silicon with quantised total angular momentum, not spin, as a highly coherent platform with tuneable intrinsic spin-orbit coupling advantageous to build artificial quantum systems and couple qubits over long distances.