title: Double inclusive small-x gluon production and their azimuthal correlations in a biased ensemble

abstract: We consider double $gg\rightarrow g$ production in the presence of a bias on the unintegrated gluon distribution of the colliding hadron or nuclei. Such bias could be due to the selection of configurations with a greater number of gluons or higher mean transverse momentum squared or, more generally, due to a modified spectral shape of the gluon distribution in the hadrons. Hence, we consider reweighted functional averages over the stochastic ensemble of small-x gluons. We evaluate explicitly the double inclusive gluon transverse momentum spectrum in high-energy collisions, and their azimuthal correlations, for a few simple examples of biases.


title: Large-scale dynamo action of magnetized Taylor-Couette flows

abstract: A conducting Taylor-Couette flow with quasi-Keplerian rotation law containing a toroidal magnetic field serves as a mean-field dynamo model of the Tayler-Spruit-type. The flows are unstable against nonaxisymmetric perturbations which form electromotive forces defining $\alpha$ effect and eddy diffusivity. If both degenerated modes with $m=\pm 1$ are excited with the same power then the global $\alpha$ effect vanishes and a dynamo cannot work. It is shown, however, that the Tayler instability produces finite $\alpha$ effects if only an isolated mode is considered but this intrinsic helicity of the single-mode is too low for an $\alpha^2$ dynamo. Moreover, an $\alpha\Om$ dynamo model with quasi-Keplerian rotation requires a minimum magnetic Reynolds number of rotation of ${\rm Rm}\simeq 2.000$ to work. Whether it really works depends on assumptions about the turbulence energy. For a steeper-than-quadratic dependence of the turbulence intensity on the magnetic field, however, dynamos are only excited if the resulting magnetic eddy diffusivity approximates its microscopic value, $\eta_{\rm T}\simeq \eta$. By basically lower or larger eddy diffusivities the dynamo instability is suppressed.


title: Experimental investigation of water emulsion fuel stability

abstract: The combustion of liquid fuels emulsified with water have long generated interest in the internal combustion engine research community. Typically, these fuels consist of small quantities of water emulsified with ultrasonification or other mechanical methods into a pure or multicomponent hydrocarbon fuel. These emulsion fuels promise significant advantages over base liquid fuels, such as better fuel economy, colder combustion temperatures, less NOx emissions, and so on. However, a significant practical disadvantage of these fuels is that they are prone to phase separation after they have been prepared. Till date, an objective but economical method of identifying the various degrees of phase separation has not been identified. Present research presents such a method and shows its utilization in analyzing the stability of water and hydrocarbon fuel emulsions over time without the addition of chemical stabilizers. It is expected that present research will pave the way in establishing this method to study the stability of other specialized multicomponent fluids.


title: Unsupervised star, galaxy, qso classification: Application of HDBSCAN

abstract: Classification will be an important first step for upcoming surveys that will detect billions of new sources such as LSST and Euclid, as well as DESI, 4MOST and MOONS. The application of traditional methods of model fitting and colour-colour selections will face significant computational constraints, while machine-learning (ML) methods offer a viable approach to tackle datasets of that volume. While supervised learning methods can perform very well for classification tasks, the creation of representative and accurate training sets is a resource and time consuming task. We present a viable alternative using an unsupervised ML method to separate stars, galaxies and QSOs using photometric data. The heart of our work uses HDBSCAN to find the star, galaxy and QSO clusters in a multidimensional colour space. We optimized the hyperparameters and input attributes of three separate HDBSCAN runs, each to select a particular object class, and thus treat the output of each separate run as a binary classifier. We subsequently consolidate the output to give our final classifications, optimized on their F1 scores. We explore the use of Random Forest and PCA as part of the pre-processing stage for feature selection and dimensionality reduction. Using our dataset of 50,000 spectroscopically labelled objects we obtain an F1 score of 98.9, 98.9 and 93.13 respectively for star, galaxy and QSO selection using our unsupervised learning method. We find that careful attribute selection is a vital part of accurate classification with HDBSCAN. We applied our classification to a subset of the SDSS spectroscopic catalogue and demonstrate the potential of our approach in correcting misclassified spectra useful for DESI and 4MOST. Finally, we create a multiwavelength catalogue of 2.7 million sources using the KiDS, VIKING and ALLWISE surveys and publish corresponding classifications and photometric redshifts.


title: Exploring the Grain Properties in the Disk of HL Tau with an Evolutionary Model

abstract: We model the ALMA and VLA millimeter radial profiles of the disk around HL Tau to constrain the properties of the dust grains. We adopt the disk evolutionary models of Lynden-Bell \& Pringle and calculate their temperature and density structure and emission. These disks are heated by the internal viscosity and irradiated by the central star and a warm envelope. We consider a dust size distribution $n(a) da \propto a^{-3.5} da $, and vary the maximum grain size in the atmosphere and the midplane, $a_{\rm max}=100\ \mu$m, 1 mm, and 1cm. We also include dust settling and vary the dust-to-gas mass ratio from 1 to 9 times the ISM value. We find that the models that can fit the observed level of emission along the profiles at all wavelengths have an atmosphere with a maximum grain size $a_{\rm max} = 100 \ \mu$m, and a midplane with $a_{\rm max}=1$ cm. The disk substructure, with a deficit of emission in the gaps, can be due to dust properties in these regions that are different from those in the rings. We test an opacity effect (different $a_{\rm max}$) and a dust mass deficit (smaller dust-to-gas mass ratio) in the gaps. We find that the emission profiles are better reproduced by models with a dust deficit in the gaps, although a combined effect is also possible. These models have a global dust-to-gas mass ratio twice the ISM value, needed to reach the level of emission of the 7.8 mm VLA profile.


title: Harmonic Mean Point Processes: Proportional Rate Error Minimization for Obtundation Prediction

abstract: In healthcare, the highest risk individuals for morbidity and mortality are rarely those with the greatest modifiable risk. By contrast, many machine learning formulations implicitly attend to the highest risk individuals. We focus on this problem in point processes, a popular modeling technique for the analysis of the temporal event sequences in electronic health records (EHR) data with applications in risk stratification and risk score systems. We show that optimization of the log-likelihood function also gives disproportionate attention to high risk individuals and leads to poor prediction results for low risk individuals compared to ones at high risk. We characterize the problem and propose an adjusted log-likelihood formulation as a new objective for point processes. We demonstrate the benefits of our method in simulations and in EHR data of patients admitted to the critical care unit for intracerebral hemorrhage.


title: Second Order Threshold Dynamics Schemes for Two Phase Motion by Mean Curvature

abstract: The threshold dynamics algorithm of Merriman, Bence, and Osher is only first order accurate in the two-phase setting. Its accuracy degrades further to half order in the multi-phase setting, a shortcoming it has in common with other related, more recent algorithms such as the equal surface tension version of the Voronoi implicit interface method. As a first, rigorous step in addressing this shortcoming, we present two different second order accurate versions of two-phase threshold dynamics. Unlike in previous efforts in this direction, we present careful consistency calculations for both of our algorithms. The first algorithm is consistent with its limit (motion by mean curvature) up to second order in any space dimension. The second achieves second order accuracy only in dimension two, but comes with a rigorous stability guarantee (unconditional energy stability) in any dimension -- a first for high order schemes of its type.


title: Dynamics of Bose-Einstein Recondensation in Higher Bands

abstract: Motivated by recent experiments, we explore the kinetics of Bose-Einstein condensation in the upper band of a double well optical lattice. These experiments engineer a non-equilibrium situation in which the highest energy state in the band is macroscopically occupied. The system subsequently relaxes and the condensate moves to the lowest energy state. We model this process, finding that the kinetics occurs in three phases: The condensate first evaporates, forming a highly non-equilibrium gas with no phase coherence. Energy is then redistributed among the noncondensed atoms. Finally the atoms recondense. We calculate the time-scales for each of these phases, and explain how this scenario can be verified through future experiments.


title: Scattering in quantum dots via noncommutative rational functions

abstract: In the customary random matrix model for transport in quantum dots with $M$ internal degrees of freedom coupled to a chaotic environment via $N\ll M$ channels, the density $\rho$ of transmission eigenvalues is computed from a specific invariant ensemble for which explicit formula for the joint probability density of all eigenvalues is available. We revisit this problem in the large $N$ regime allowing for (i) arbitrary ratio $\phi := N/M \le 1$; and (ii) general distributions for the matrix elements of the Hamiltonian of the quantum dot. In the limit $\phi \to 0$ we recover the formula for the density $\rho$ that Beenakker (Rev. Mod. Phys., 69:731-808, 1997) has derived for a special matrix ensemble. We also prove that the inverse square root singularity of the density at zero and full transmission in Beenakker's formula persists for any $\phi <1$ but in the borderline case $\phi=1$ an anomalous $\lambda^{-2/3}$ singularity arises at zero. To access this level of generality we develop the theory of global and local laws on the spectral density of a large class of noncommutative rational expressions in large random matrices with i.i.d. entries.


title: Semi-Supervised Multi-Organ Segmentation through Quality Assurance Supervision

abstract: Human in-the-loop quality assurance (QA) is typically performed after medical image segmentation to ensure that the systems are performing as intended, as well as identifying and excluding outliers. By performing QA on large-scale, previously unlabeled testing data, categorical QA scores can be generatedIn this paper, we propose a semi-supervised multi-organ segmentation deep neural network consisting of a traditional segmentation model generator and a QA involved discriminator. A large-scale dataset of 2027 volumes are used to train the generator, whose 2-D montage images and segmentation mask with QA scores are used to train the discriminator. To generate the QA scores, the 2-D montage images were reviewed manually and coded 0 (success), 1 (errors consistent with published performance), and 2 (gross failure). Then, the ResNet-18 network was trained with 1623 montage images in equal distribution of all three code labels and achieved an accuracy 94% for classification predictions with 404 montage images withheld for the test cohort. To assess the performance of using the QA supervision, the discriminator was used as a loss function in a multi-organ segmentation pipeline. The inclusion of QA-loss function boosted performance on the unlabeled test dataset from 714 patients to 951 patients over the baseline model. Additionally, the number of failures decreased from 606 (29.90%) to 402 (19.83%). The contributions of the proposed method are threefold: We show that (1) the QA scores can be used as a loss function to perform semi-supervised learning for unlabeled data, (2) the well trained discriminator is learnt by QA score rather than traditional true/false, and (3) the performance of multi-organ segmentation on unlabeled datasets can be fine-tuned with more robust and higher accuracy than the original baseline method.


title: Coordinated Management of Processor Configuration and Cache Partitioning to Optimize Energy under QoS Constraints

abstract: An effective way to improve energy efficiency is to throttle hardware resources to meet a certain performance target, specified as a QoS constraint, associated with all applications running on a multicore system. Prior art has proposed resource management (RM) frameworks in which the share of the last-level cache (LLC) assigned to each processor and the voltage-frequency (VF) setting for each processor is managed in a coordinated fashion to reduce energy. A drawback of such a scheme is that, while one core gives up LLC resources for another core, the performance drop must be compensated by a higher VF setting which leads to a quadratic increase in energy consumption. By allowing each core to be adapted to exploit instruction and memory-level parallelism (ILP/MLP), substantially higher energy savings are enabled. This paper proposes a coordinated RM for LLC partitioning, processor adaptation, and per-core VF scaling. A first contribution is a systematic study of the resource trade-offs enabled when trading between the three classes of resources in a coordinated fashion. A second contribution is a new RM framework that utilizes these trade-offs to save more energy. Finally, a challenge to accurately model the impact of resource throttling on performance is to predict the amount of MLP with high accuracy. To this end, the paper contributes with a mechanism that estimates the effect of MLP over different processor configurations and LLC allocations. Overall, we show that up to 18% of energy, and on average 10%, can be saved using the proposed scheme.


title: Deep Multi-task Prediction of Lung Cancer and Cancer-free Progression from Censored Heterogenous Clinical Imaging

abstract: Annual low dose computed tomography (CT) lung screening is currently advised for individuals at high risk of lung cancer (e.g., heavy smokers between 55 and 80 years old). The recommended screening practice significantly reduces all-cause mortality, but the vast majority of screening results are negative for cancer. If patients at very low risk could be identified based on individualized, image-based biomarkers, the health care resources could be more efficiently allocated to higher risk patients and reduce overall exposure to ionizing radiation. In this work, we propose a multi-task (diagnosis and prognosis) deep convolutional neural network to improve the diagnostic accuracy over a baseline model while simultaneously estimating a personalized cancer-free progression time (CFPT). A novel Censored Regression Loss (CRL) is proposed to perform weakly supervised regression so that even single negative screening scans can provide small incremental value. Herein, we study 2287 scans from 1433 de-identified patients from the Vanderbilt Lung Screening Program (VLSP) and Molecular Characterization Laboratories (MCL) cohorts. Using five-fold cross-validation, we train a 3D attention-based network under two scenarios: (1) single-task learning with only classification, and (2) multi-task learning with both classification and regression. The single-task learning leads to a higher AUC compared with the Kaggle challenge winner pre-trained model (0.878 v. 0.856), and multi-task learning significantly improves the single-task one (AUC 0.895, p<0.01, McNemar test). In summary, the image-based predicted CFPT can be used in follow-up year lung cancer prediction and data assessment.


title: An Unethical Optimization Principle

abstract: If an artificial intelligence aims to maximise risk-adjusted return, then under mild conditions it is disproportionately likely to pick an unethical strategy unless the objective function allows sufficiently for this risk. Even if the proportion ${\eta}$ of available unethical strategies is small, the probability ${p_U}$ of picking an unethical strategy can become large; indeed unless returns are fat-tailed ${p_U}$ tends to unity as the strategy space becomes large. We define an Unethical Odds Ratio Upsilon (${\Upsilon}$) that allows us to calculate ${p_U}$ from ${\eta}$, and we derive a simple formula for the limit of ${\Upsilon}$ as the strategy space becomes large. We give an algorithm for estimating ${\Upsilon}$ and ${p_U}$ in finite cases and discuss how to deal with infinite strategy spaces. We show how this principle can be used to help detect unethical strategies and to estimate ${\eta}$. Finally we sketch some policy implications of this work.


title: The Barwise-Schlipf Theorem

abstract: In 1975 Barwise and Schlipf published a landmark paper whose main theorem asserts that a nonstandard model $\mathcal{M}$ of PA (Peano arithmetic) is recursively saturated iff $\mathcal{M}$ has an expansion that satisfies the subsystem $\Delta_1^1$-${\sf CA}_0$ of second order arithmetic. In this paper we identify a crucial error in the Barwise-Schlipf proof of the right-to-left direction of the theorem, and additionally, we offer a correct proof of the problematic direction.


title: On the generic family of Cayley graphs of a finite group

abstract: Let $G$ be a finite group. For each $m>1$ we define the symmetric canonical subset $S=S(m)$ of the Cartesian power $G^m$ and we consider the family of Cayley graphs $\mathscr{G}_m(G)=Cay(G^m,S)$. We describe properties of these graphs and show that for a fixed $m>1$ and groups $G$ and $H$ the graphs $\mathscr{G}_m(G)$ and $\mathscr{G}_m(H)$ are isomorphic if and only if the groups $G$ and $H$ are isomorphic. We describe also the groups of automorphisms $\mathbf{Aut}(\mathscr{G}_m(G))$. It is shown that if $G$ is a non-abelian group, then $\mathbf{Aut}(\mathscr{G}_m(G))\simeq \big(G^m \rtimes \mathbf{Aut}(G)\big)\rtimes D_{m+1}$, where $D_{m+1}$ is the dihedral group of order $2m+2$. If $G$ is an abelian group (with some exceptions for $m=3$), then $\mathbf{Aut}(\mathscr{G}_m(G))\simeq G^m\rtimes \big(\mathbf{Aut}(G)\times S_{m+1}\big)$, where $S_{m+1}$ is the symmetric group of degree $m+1$. As an example of application we discuss relations between Cayley graphs $\mathscr{G}_m(G)$ and Bergman-Isaacs Theorem on rings with fixed-point-free group actions.


title: Rigidity of Kac-Schwarz operators

abstract: In his work on the mathematical formulation of 2d quantum gravity Schwarz established a rigidity result for Kac-Schwarz operators for the n-KdV hierarchies. Later on, Adler and van Moerbeke as well as Fastr\'{e} obtained different proofs of this result. We give yet another proof of the rigidity, one that in fact holds for all Drinfeld-Sokolov hierarchies.


title: Numerical solutions for a class of singular boundary value problems arising in the theory of epitaxial growth

abstract: The existence of numerical solutions to a fourth order singular boundary value problem arising in the theory of epitaxial growth is studied. An iterative numerical method is applied on a second order nonlinear singular boundary value problem which is the exact result of the reduction of this fourth order singular boundary value problem. It turns out that the existence or nonexistence of numerical solutions fully depends on the value of a parameter. We show that numerical solutions exist for small positive values of this parameter. For large positive values of the parameter, we find nonexistence of solutions. We also observe existence of solutions for negative values of the parameter and determine the range of parameter values which separates existence and nonexistence of solutions. This parameter has a clear physical meaning as it describes the rate at which new material is deposited onto the system. This fact allows us to interpret the physical significance of our results.


title: Detecting Patterns of Physiological Response to Hemodynamic Stress via Unsupervised Deep Learning

abstract: Monitoring physiological responses to hemodynamic stress can help in determining appropriate treatment and ensuring good patient outcomes. Physicians' intuition suggests that the human body has a number of physiological response patterns to hemorrhage which escalate as blood loss continues, however the exact etiology and phenotypes of such responses are not well known or understood only at a coarse level. Although previous research has shown that machine learning models can perform well in hemorrhage detection and survival prediction, it is unclear whether machine learning could help to identify and characterize the underlying physiological responses in raw vital sign data. We approach this problem by first transforming the high-dimensional vital sign time series into a tractable, lower-dimensional latent space using a dilated, causal convolutional encoder model trained purely unsupervised. Second, we identify informative clusters in the embeddings. By analyzing the clusters of latent embeddings and visualizing them over time, we hypothesize that the clusters correspond to the physiological response patterns that match physicians' intuition. Furthermore, we attempt to evaluate the latent embeddings using a variety of methods, such as predicting the cluster labels using explainable features.


title: A two-player portfolio tracking game

abstract: We study the competition of two strategic agents for liquidity in the benchmark portfolio tracking setup of Bank, Soner, Vo{\ss} (2017). Specifically, both agents track their own stochastic running trading targets while interacting through common aggregated temporary and permanent price impact \`a la Almgren and Chriss (2001). The resulting stochastic linear quadratic differential game with terminal state constraints allows for a unique and explicitly available open-loop Nash equilibrium. Our results reveal how the equilibrium strategies of the two players take into account the other agent's trading targets: either in an exploitative intent or by providing liquidity to the competitor, depending on the relation between temporary and permanent price impact. As a consequence, different behavioral patterns can emerge as optimal in equilibrium. These insights complement and extend existing studies in the literature on predatory trading models examined in the context of optimal portfolio liquidation games.


title: On the Axiomatic Systems of Singular Cohomology Theory

abstract: On the category of pairs of topological spaces having a homotopy type of $CW$ complexes the singular (co)homology theory was axiomatically studied by J.Milnor. In particular, Milnor gave additivity axiom for a (co)homology theory and proved that any additive (co)homology theory on the given category is isomorphic to the singular (co)homology. On the other hand, the singular homology is a homology with compact support \cite{3}. In the paper \cite{6}, L. Mdzinarishvili proposed {\it partially compact support property} for a cohomology theory and gave another axiomatic characterization of the singular cohomology theory \cite{6}. In this paper, we will give additional different axiomatic characterizations of the singular cohomology theory. Moreover, we will study connections of the mentioned axiomatic systems.


title: Imperfect bifurcations in opinion dynamics under external fields

abstract: We investigate, through a kinetic-exchange model, the impact that an external field, like advertising and propaganda, has on opinion dynamics. We address the situations where two opposite alternatives can be selected but the possibility of indecision also exists. In this model, individuals influence each other through pairwise interactions, which can be of agreement or disagreement, and there are also external fields that can skew decision making. Two parameters are used to model the interactions with the field: one measures the sensitivity of the individuals to be influenced, another quantifies in which direction. We study this model in a fully connected social network scenario, by means of numerical simulations of the kinetic exchange dynamics and analytical results derived from the mean-field rate equations. We show how the external bias gives rise to imperfect bifurcations, and cusp catastrophes, allowing abrupt changes and hysteresis depending on the level of disagreement in interpersonal interactions and on the strength of the external influence.


title: Robust Fitting for Generalized Additive Models for Location, Scale and Shape

abstract: The validity of estimation and smoothing parameter selection for the wide class of generalized additive models for location, scale and shape (GAMLSS) relies on the correct specification of a likelihood function. Deviations from such assumption are known to mislead any likelihood-based inference and can hinder penalization schemes meant to ensure some degree of smoothness for non-linear effects. We propose a general approach to achieve robustness in fitting GAMLSSs by limiting the contribution of observations with low log-likelihood values. Robust selection of the smoothing parameters can be carried out either by minimizing information criteria that naturally arise from the robustified likelihood or via an extended Fellner-Schall method. The latter allows for automatic smoothing parameter selection and is particularly advantageous in applications with multiple smoothing parameters. We also address the challenge of tuning robust estimators for models with non-linear effects by proposing a novel median downweighting proportion criterion. This enables a fair comparison with existing robust estimators for the special case of generalized additive models, where our estimator competes favorably. The overall good performance of our proposal is illustrated by further simulations in the GAMLSS setting and by an application to functional magnetic resonance brain imaging using bivariate smoothing splines.


title: KPsec: Secure End-to-End Communications for Multi-Hop Wireless Networks

abstract: The security of cyber-physical systems, from self-driving cars to medical devices, depends on their underlying multi-hop wireless networks. Yet, the lack of trusted central infrastructures and limited nodes' resources make securing these networks challenging. Recent works on key pre-distribution schemes, where nodes communicate over encrypted overlay paths, provide an appealing solution because of their distributed, computationally light-weight nature. Alas, these schemes share a glaring security vulnerability: the two ends of every overlay link can decrypt---and potentially modify and alter---the message. Plus, the longer overlay paths impose traffic overhead and increase latency. We present a novel routing mechanism, KPsec, to address these issues. KPsec deploys multiple disjoint paths and an initial key-exchange phase to secure end-to-end communications. After the initial key-exchange phase, traffic in KPsec follows the shortest paths and, in contrast to key pre-distribution schemes, intermediate nodes cannot decrypt it. We measure the security and performance of KPsec as well as three state-of-the-art key pre-distribution schemes using a real 10-node testbed and large-scale simulations. Our experiments show that, in addition to its security benefits, KPsec results in $5-15\%$ improvement in network throughput, up to $75\%$ reduction in latency, and an order of magnitude reduction in energy consumption.


title: Distributed Online Convex Optimization with Improved Dynamic Regret

abstract: In this paper, we consider the problem of distributed online convex optimization, where a group of agents collaborate to track the global minimizers of a sum of time-varying objective functions in an online manner. Specifically, we propose a novel distributed online gradient descent algorithm that relies on an online adaptation of the gradient tracking technique used in static optimization. We show that the dynamic regret bound of this algorithm has no explicit dependence on the time horizon and, therefore, can be tighter than existing bounds especially for problems with long horizons. Our bound depends on a new regularity measure that quantifies the total change in the gradients at the optimal points at each time instant. Furthermore, when the optimizer is approximatly subject to linear dynamics, we show that the dynamic regret bound can be further tightened by replacing the regularity measure that captures the path length of the optimizer with the accumulated prediction errors, which can be much lower in this special case. We present numerical experiments to corroborate our theoretical results.


title: Effective field theories as a novel probe of fine-tuning of cosmic inflation

abstract: The leading account of several salient observable features of our universe today is provided by the theory of cosmic inflation. But an important and thus far intractable question is whether inflation is generic, or whether it is finely tuned---requiring very precisely specified initial conditions. In this paper I argue that a recent, model-independent characterization of inflation, known as the 'effective field theory (EFT) of inflation', promises to address this question in a thoroughly modern and significantly more comprehensive way than in the existing literature. To motivate and provide context for this claim, I distill three core problems with the theory of inflation, which I dub the permissiveness problem, the initial conditions problem, and the multiverse problem. I argue that the initial conditions problem lies within the scope of EFTs of inflation as they are currently conceived, whereas the other two problems remain largely intractable: their solution must await a more complete description of the very early universe. I highlight recent work that addresses the initial conditions problem within the context of a dynamical systems analysis of a specific (state-of-the-art) EFT of inflation, and conclude with a roadmap for how such work might be extended to realize the promise claimed above.


title: Managing access to primary care clinics using robust scheduling templates

abstract: An important challenge confronting healthcare is the effective management of access to primary care. Robust appointment scheduling policies/templates can help strike an effective balance between the lead-time to an appointment (a.k.a. indirect waiting time, measuring the difference between a patient's desired and actual appointment dates) and waiting times at the clinic on the day of the appointment (a.k.a. direct waiting time). We propose methods for identifying effective appointment scheduling templates using a two-stage stochastic mixed-integer linear program model. The model embeds simulation for accurate evaluation of direct waiting times and uses sample average approximation method for computational efficiency. The model accounts for patients' no-show behaviors, provider availability, overbooking, demand uncertainty, and overtime constraints. The model allows the scheduling templates to be potentially updated at regular intervals while minimizing the patient expected waiting times and balancing provider utilization. Proposed methods are validated using data from the U.S. Department of Veterans Affairs (VA) primary care clinics.


title: Turnaround physics beyond spherical symmetry

abstract: The concept of turnaround surface in an accelerating universe is generalized to arbitrarily large deviations from spherical symmetry, to close the gap between the idealized theoretical literature and the real world observed by astronomers. As an analytical application, the characterization of turnaround surface is applied to small deviations from spherical symmetry, recovering a previous result while extending it to scalar-tensor gravity.


title: Modelling the distributions of white dwarf atmospheric pollution: a low Mg abundance for accreted planetesimals?

abstract: The accretion of planetesimals onto white dwarf atmospheres allows determination of the composition of this polluting material. This composition is usually inferred from observed pollution levels by assuming it originated from a single body. This paper instead uses a stochastic model wherein polluting planetesimals are chosen randomly from a mass distribution, finding that the single body assumption is invalid in >20% of cases. Planetesimal compositions are modelled assuming parent bodies that differentiated into core, mantle and crust components. Atmospheric levels of Ca, Mg and Fe in the model are compared to a sample of 230 DZ white dwarfs for which such pollution is measured. A good fit is obtained when each planetesimal has its core, mantle and crust fractions chosen independently from logit-normal distributions which lead to average mass fractions of $f_\text{Cru}=0.15$, $f_\text{Man}=0.49$ and $f_\text{Cor}=0.36$. However, achieving this fit requires a factor 4 depletion of Mg relative to stellar material. This depletion is unlikely to originate in planetesimal formation processes, but might occur from heating while the star is on the giant branch. Alternatively the accreted material has stellar abundance, and either the inferred low Mg abundance was caused by an incorrect assumption that Mg sinks slower than Ca and Fe, or there are unmodelled biases in the observed sample. Finally, the model makes predictions for the timescale on which the observed pollutant composition varies, which should be the longer of the sinking and disc timescales, implying variability on decadal timescales for DA white dwarfs.


title: Heterogeneities Shape Passive Intracellular Transport

abstract: A living cell's interior is one of the most complex and intrinsically dynamic systems, providing an elaborate interplay between cytosolic crowding and ATP-driven motion, which controls cellular functionality. Here, we investigated two distinct fundamental features of the merely passive, not-bio-motor shuttled material transport within the cytoplasm of Dictyostelium discoideum cells: the anomalous non-linear scaling of the mean-squared displacement of a 150nm-diameter particle and non-Gaussian distribution of increments. Relying on single-particle tracking data of 320,000 data points, we performed a systematic analysis of four possible origins for non-Gaussian transport: (1) sample-based variability, (2) rare occurring strong motion events, (3) ergodicity breaking/ageing, and (4) spatio-temporal heterogeneities of the intracellular medium. After excluding the first three reasons, we investigated the remaining hypothesis of a heterogeneous cytoplasm as cause for non-Gaussian transport. A novel fit model with randomly distributed diffusivities implementing medium heterogeneities suits the experimental data. Strikingly, the non-Gaussian feature is independent of the cytoskeleton condition and lag time. This reveals that efficiency and consistency of passive intracellular transport and the related anomalous scaling of the mean-squared displacement are regulated by cytoskeleton components, while cytoplasmic heterogeneities are responsible for the generic, non-Gaussian distribution of increments.


title: Energy-Efficient Beamforming and Cooperative Jamming in IRS-Assisted MISO Networks

abstract: The energy-efficient (EE) design problem is very important to design a good performance network. Besides this, the security of communication is also important. Based on the new IRS technology, the EE maximization problem is formed in this paper for secure communication. We proposed a solution that outperformed other methods.


title: Mapping standing-wave cavity modes with a commercial scanning near-field microscope tip

abstract: We describe a method to map the standing-wave pattern inside a Fabry-Perot optical cavity with sub-wavelength resolution by perturbing it with a commercially available scanning near-field optical microscope (SNOM) tip. The method is applied to a fiber Fabry-Perot microcavity. We demonstrate its use to determine the relative position of the antinodes at two different wavelengths. In addition, we use the SNOM tip as a point-like source allowing precise positioning of a microscope objective with respect to the cavity mode.


title: A note on nearly Platonic graphs with connectivity one

abstract: A k-regular planar graph G is nearly Platonic when all faces but one are of the same degree while the remaining face is of a different degree. We show that no such graphs with connectivity one can exist. This complements a recent result by Keith, Froncek, and Kreher on non-existence of 2-connected nearly Platonic graphs.


title: Fast Computation of Sep$_\lambda$ via Interpolation-based Globality Certificates

abstract: Given two square matrices $A$ and $B$, we propose a new approach for computing the smallest value $\varepsilon \geq 0$ such that $A+E$ and $A+F$ share an eigenvalue, where $\|E\|=\|F\|=\varepsilon$. In 2006, Gu and Overton proposed the first algorithm for computing this quantity, called $\mathrm{sep}_\lambda(A,B)$ ("sep-lambda"), using ideas inspired from an earlier algorithm of Gu for computing the distance to uncontrollability. However, the algorithm of Gu and Overton is extremely expensive, which limits it to the tiniest of problems, and until now, no other algorithms have been known. Our new algorithm can be orders of magnitude faster and can solve problems where $A$ and $B$ are of moderate size. Moreover, our method consists of many "embarrassingly parallel" computations, and so it can be further accelerated on multi-core hardware. Finally, we also propose the first algorithm to compute an earlier version of sep-lambda where $\|E\| + \|F\|=\varepsilon$.


title: Causes of the Diurnal Variation observed in $\gamma$-ray Spectrum using NaI (Tl) Detector

abstract: Present paper studies the $\gamma$-ray spectrum observed using NaI(Tl) scintillation detector over a year. The paper reports the presence of a distinct diurnal pattern in the total number of $\gamma$-ray counts detected by the NaI detector. The counts start decreasing after sunrise and show gradual recovery after sunset. The amplitude of this variation is quite significant ( > 10%) at the observation site Tirunelveli, South India. However, further investigation based on different energy ranges reveals that the mentioned diurnal pattern is actually present only in the energies related to the terrestrial background radioactivity. The study demonstrates that the pattern is associated with the radioactivity of isotopes of radon ($^{222}$Rn and $^{220}$Rn) and their daughter radionuclides. The amplitude of the diurnal variation is found to have seasonal dependence, with the lowest amplitude during June-October (~ 2%), and highest in April-May months (~ 14%). The dependence of the amplitude of the diurnal variation on various atmospheric parameters, namely temperature, relative humidity, and vertical wind is examined. The observed diurnal pattern is attributed to the generation of the inversion layer. The distribution of concentration of radon and thoron progenies in an air column changes with the vertical mixing and atmospheric boundary layer (ABL), resulting in the diurnal variation of $\gamma$-rays.


title: Varieties corresponding to classes of complemented posets

abstract: As algebraic semantics of the logic of quantum mechanics there are usually used orthomodular posets, i.e. bounded posets with a complementation which is an antitone involution and where the join of orthogonal elements exists and the orthomodular law is satisfied. When we omit the condition that the complementation is an antitone involution, then we obtain skew-orthomodular posets. To each such poset we can assign a bounded lambda-lattice in a non-unique way. Bounded lambda-lattices are lattice-like algebras whose operations are not necessarily associative. We prove that any of the following properties for bounded posets with a unary operation can be characterized by certain identities of an arbitrary assigned lambda-lattice: complementarity, orthogonality, almost skew-orthomodularity and skew-orthomodularity. It is shown that these identities are independent. Finally, we show that the variety of skew-orthomodular lambda-lattices is congruence permutable as well as congruence regular.


title: Causality-based tests to detect the influence of confounders on mobile health diagnostic applications: a comparison with restricted permutations

abstract: Machine learning practice is often impacted by confounders. Confounding can be particularly severe in remote digital health studies where the participants self-select to enter the study. While many different confounding adjustment approaches have been proposed in the literature, most of these methods rely on modeling assumptions, and it is unclear how robust they are to violations of these assumptions. This realization has recently motivated the development of restricted permutation methods to quantify the influence of observed confounders on the predictive performance of a machine learning models and evaluate if confounding adjustment methods are working as expected. In this paper we show, nonetheless, that restricted permutations can generate biased estimates of the contribution of the confounders to the predictive performance of a learner, and we propose an alternative approach to tackle this problem. By viewing a classification task from a causality perspective, we are able to leverage conditional independence tests between predictions and test set labels and confounders in order to detect confounding on the predictive performance of a classifier. We illustrate the application of our causality-based approach to data collected from mHealth study in Parkinson's disease.


title: Unsupervised Medical Image Segmentation with Adversarial Networks: From Edge Diagrams to Segmentation Maps

abstract: We develop and approach to unsupervised semantic medical image segmentation that extends previous work with generative adversarial networks. We use existing edge detection methods to construct simple edge diagrams, train a generative model to convert them into synthetic medical images, and construct a dataset of synthetic images with known segmentations using variations on extracted edge diagrams. This synthetic dataset is then used to train a supervised image segmentation model. We test our approach on a clinical dataset of kidney ultrasound images and the benchmark ISIC 2018 skin lesion dataset. We show that our unsupervised approach is more accurate than previous unsupervised methods, and performs reasonably compared to supervised image segmentation models. All code and trained models are available at https://github.com/kiretd/Unsupervised-MIseg.


title: A 2-Site for Continuous 2-Group Actions

abstract: Elmendorf's Theorem states that the category of continuous actions of a topological group is a Grothendieck topos in the sense that it is equivalent to a category of sheaves on a site. This paper offers a 2-dimensional generalization by showing that a certain 2-category of continuous actions of a topological 2-group is 2-equivalent to a 2-category of 2-sheaves on a suitable 2-site.


title: Incentivized Exploration for Multi-Armed Bandits under Reward Drift

abstract: We study incentivized exploration for the multi-armed bandit (MAB) problem where the players receive compensation for exploring arms other than the greedy choice and may provide biased feedback on reward. We seek to understand the impact of this drifted reward feedback by analyzing the performance of three instantiations of the incentivized MAB algorithm: UCB, $\varepsilon$-Greedy, and Thompson Sampling. Our results show that they all achieve $\mathcal{O}(\log T)$ regret and compensation under the drifted reward, and are therefore effective in incentivizing exploration. Numerical examples are provided to complement the theoretical analysis.


title: LRP2020: The Opportunity of Young Nearby Associations with the Advent of the Gaia Mission

abstract: This white paper proposes leveraging high-quality Gaia data available to the worldwide scientific community and complement it with support from Canadian-related facilities to place Canada as a leader in the fields of stellar associations and exoplanet science, and to train Canadian highly qualified personnel through graduate and post-graduate research grants. Gaia has sparked a new era in the study of stellar kinematics by measuring precise distances and proper motions for 1.3 billion stars. These data have already generated more than 1700 scientific papers and are guaranteed to remain the source of many more papers for the upcoming decades. More than 900 new age-calibrated young low-mass stars have already been discovered as a direct consequence of the second Gaia data release. Some of these may already be host stars to known exoplanet systems or may become so with the progress of the TESS mission that is expected to discover 10,000 nearby transiting exoplanets in the upcoming decade. This places Canada in a strategic position to leverage Gaia data because it has access to several high-resolution spectrometers on 1-4 m class telescopes (e.g. The ESPaDOnS, SPIRou and NIRPS), that would allow to quickly characterize this large number of low-mass stars and their exoplanet systems. This white paper describes the opportunity in such scientific projects that could place Canada as a leader in the fields of stellar associations and exoplanets.


title: PRESTvO: PRivacy Enabled Smartphone-based access To vehicle On-board units

abstract: Smartphones are quickly moving toward complementing or even replacing traditional car keys. We advocate a role-based access control policy mixed with attributes that facilitates access to various functionalities of vehicular on-board units from smartphones. We use a rights-based access control policy for in-vehicle functionalities similar to the case of a file allocation table of a contemporary OS, in which read, write or execute operations can be performed over various vehicle functions. Further, to assure the appropriate security, we develop a protocol suite using identity-based cryptography and we rely on group signatures that preserve the anonymity of group members for assuring privacy and traceability. To prove the feasibility of our approach, we develop a proof-of-concept implementation with modern smartphones, aftermarket Android head-units and test computational feasibility on a real-world in-vehicle controller. Our implementation relies on state-of-the-art cryptography, including traditional building blocks and more modern pairing-friendly curves, that facilitate the adoption of group signatures and identity-based cryptography in automotive-based scenarios.


title: Deconstructing Stellar Consensus (Extended Version)

abstract: Some of the recent blockchain proposals, such as Stellar and Ripple, allow for open membership while using quorum-like structures typical for classical Byzantine consensus with closed membership. This is achieved by constructing quorums in a decentralised way: each participant independently chooses whom to trust, and quorums arise from these individual decisions. Unfortunately, the consensus protocols underlying such blockchains are poorly understood, and their correctness has not been rigorously investigated. In this paper we rigorously prove correct the Stellar Consensus Protocol (SCP), with our proof giving insights into the protocol structure and its use of lower-level abstractions. To this end, we first propose an abstract version of SCP that uses as a black box Stellar's federated voting primitive (analogous to reliable Byzantine broadcast), previously investigated by Garc\'ia-P\'erez and Gotsman. The abstract consensus protocol highlights a modular structure in Stellar and can be proved correct by reusing the previous results on federated voting. However, it is unsuited for realistic implementations, since its processes maintain infinite state. We thus establish a refinement between the abstract protocol and the concrete SCP that uses only finite state, thereby carrying over the result about the correctness of former to the latter. Our results help establish the theoretical foundations of decentralised blockchains like Stellar and gain confidence in their correctness.


title: HyPar-Flow: Exploiting MPI and Keras for Scalable Hybrid-Parallel DNN Training using TensorFlow

abstract: To reduce training time of large-scale DNNs, scientists have started to explore parallelization strategies like data-parallelism, model-parallelism, and hybrid-parallelism. While data-parallelism has been extensively studied and developed, several problems exist in realizing model-parallelism and hybrid-parallelism efficiently. Four major problems we focus on are: 1) defining a notion of a distributed model across processes, 2) implementing forward/back-propagation across process boundaries that requires explicit communication, 3) obtaining parallel speedup on an inherently sequential task, and 4) achieving scalability without losing out on a model's accuracy. To address these problems, we create HyPar-Flow --- a model-size/-type agnostic, scalable, practical, and user-transparent system for hybrid-parallel training by exploiting MPI, Keras, and TensorFlow. HyPar-Flow provides a single API that can be used to perform data, model, and hybrid parallel training of any Keras model at scale. We create an internal distributed representation of the user-provided Keras model, utilize TF's Eager execution features for distributed forward/back-propagation across processes, exploit pipelining to improve performance and leverage efficient MPI primitives for scalable communication. Between model partitions, we use send and recv to exchange layer-data/partial-errors while allreduce is used to accumulate/average gradients across model replicas. Beyond the design and implementation of HyPar-Flow, we also provide comprehensive correctness and performance results on three state-of-the-art HPC systems including TACC Frontera (#5 on Top500.org). For ResNet-1001, an ultra-deep model, HyPar-Flow provides: 1) Up to 1.6x speedup over Horovod-based data-parallel training, 2) 110x speedup over single-node on 128 Stampede2 nodes, and 3) 481x speedup over single-node on 512 Frontera nodes.


title: On the Strength of Connectivity of Inhomogeneous Random K-out Graphs

abstract: Random graphs are an important tool for modelling and analyzing the underlying properties of complex real-world networks. In this paper, we study a class of random graphs known as the inhomogeneous random K-out graphs which were recently introduced to analyze heterogeneous sensor networks secured by the pairwise scheme. In this model, first, each of the $n$ nodes is classified as type-1 (respectively, type-2) with probability $0<\mu<1$ (respectively, $1-\mu)$ independently from each other. Next, each type-1 (respectively, type-2) node draws 1 arc towards a node (respectively, $K_n$ arcs towards $K_n$ distinct nodes) selected uniformly at random, and then the orientation of the arcs is ignored. From the literature on homogeneous K-out graphs wherein all nodes select $K_n$ neighbors (i.e., $\mu=0$), it is known that when $K_n \geq2$, the graph is $K_n$-connected asymptotically almost surely (a.a.s.) as $n$ gets large. In the inhomogeneous case (i.e., $\mu>0$), it was recently established that achieving even 1-connectivity a.a.s. requires $K_n=\omega(1)$. Here, we provide a comprehensive set of results to complement these existing results. First, we establish a sharp zero-one law for $k$-connectivity, showing that for the network to be $k$-connected a.a.s., we need to set $K_n = \frac{1}{1-\mu}(\log n +(k-2)\log\log n + \omega(1))$ for all $k=2, 3, \ldots$. Despite such large scaling of $K_n$ being required for $k$-connectivity, we show that the trivial condition of $K_n \geq 2$ for all $n$ is sufficient to ensure that inhomogeneous K-out graph has a connected component of size $n-O(1)$ whp.


title: Identifying predictive biomarkers of CIMAvaxEGF success in advanced Lung Cancer Patients

abstract: Objectives: To identify predictive biomarkers of CIMAvaxEGF success in the treatment of Non-Small Cell Lung Cancer Patients. Methods: Data from a clinical trial evaluating the effect on survival time of CIMAvax-EGF versus best supportive care were analyzed retrospectively following the causal inference approach. Pre-treatment potential predictive biomarkers included basal serum EGF concentration, peripheral blood parameters and immunosenescence biomarkers (The proportion of CD8 + CD28- T cells, CD4+ and CD8+ T cells, CD4 CD8 ratio and CD19+ B cells. The 33 patients with complete information were included. The predictive causal information (PCI) was calculated for all possible models. The model with a minimum number of predictors, but with high prediction accuracy (PCI>0.7) was selected. Good, rare and poor responder patients were identified using the predictive probability of treatment success. Results: The mean of PCI increased from 0.486, when only one predictor is considered, to 0.98 using the multivariate approach with all predictors. The model considering the proportion of CD4+ T cell, basal EGF concentration, NLR, Monocytes, and Neutrophils as predictors were selected (PCI>0.74). Patients predicted as good responders according to the pre-treatment biomarkers values treated with CIMAvax-EGF had a significant higher observed survival compared with the control group (p=0.03). No difference was observed for bad responders. Conclusions: Peripheral blood parameters and immunosenescence biomarkers together with basal EGF concentration in serum resulted in good predictors of the CIMAvax-EGF success in advanced NSCLC. The study illustrates the application of a new methodology, based on causal inference, to evaluate multivariate pre-treatment predictors


title: Hitting probabilities for L\'{e}vy processes on the real line

abstract: We prove sharp two-sided estimates on the tail probability of the first hitting time of bounded interval as well as its asymptotic behaviour for general non-symmetric processes which satisfy an integral condition \[ \int_0^{\infty} \frac{d\xi}{1+\operatorname{Re} \psi(\xi)}<\infty. \] To this end, we first prove and then apply the global scale invariant Harnack inequality. Results are obtained under certain conditions on the characteristic exponent. We provide a wide class of L\'{e}vy processs which satisfy these assumptions.


title: TESS Reveals HD 118203 b to be a Transiting Planet

abstract: The exoplanet HD 118203 b, orbiting a bright (V = 8.05) host star, was discovered using the radial velocity method by da Silva et al. (2006), but was not previously known to transit. TESS photometry has revealed that this planet transits its host star. Five planetary transits were observed by TESS, allowing us to measure the radius of the planet to be $1.133 \pm 0.031 R_J$, and to calculate the planet mass to be $2.173 \pm 0.078 M_J$. The host star is slightly evolved with an effective temperature of $T_{\rm eff} = 5692 \pm 83$ K and a surface gravity of ${\rm log}(g) = 3.891 \pm 0.019$. With an orbital period of $6.134980 \pm 0.000038$ days and an eccentricity of $0.316 \pm 0.021$, the planet occupies a transitional regime between circularized hot Jupiters and more dynamically active planets at longer orbital periods. The host star is among the ten brightest known to have transiting giant planets, providing opportunities for both planetary atmospheric and asteroseismic studies.


title: Learning from Data-Rich Problems: A Case Study on Genetic Variant Calling

abstract: Next Generation Sequencing can sample the whole genome (WGS) or the 1-2% of the genome that codes for proteins called the whole exome (WES). Machine learning approaches to variant calling achieve high accuracy in WGS data, but the reduced number of training examples causes training with WES data alone to achieve lower accuracy. We propose and compare three different data augmentation strategies for improving performance on WES data: 1) joint training with WES and WGS data, 2) warmstarting the WES model from a WGS model, and 3) joint training with the sequencing type specified. All three approaches show improved accuracy over a model trained using just WES data, suggesting the ability of models to generalize insights from the greater WGS data while retaining performance on the specialized WES problem. These data augmentation approaches may apply to other problem areas in genomics, where several specialized models would each see only a subset of the genome.


title: Cohomology with local coefficients and knotted manifolds

abstract: We show how the classical notions of cohomology with local coefficients, CW-complex, covering space, homeomorphism equivalence, simple homotopy equivalence, tubular neighbourhood, and spinning can be encoded on a computer and used to calculate ambient isotopy invariants of continuous embeddings $N\hookrightarrow M$ of one topological manifold into another. More specifically, we describe an algorithm for computing the homology $H_n(X,A)$ and cohomology $H^n(X,A)$ of a finite connected CW-complex X with local coefficients in a $\mathbb Z\pi_1X$-module $A$ when $A$ is finitely generated over $\mathbb Z$. It can be used, in particular, to compute the integral cohomology $H^n(\widetilde X_H,\mathbb Z)$ and induced homomorphism $H^n(X,\mathbb Z) \rightarrow H^n(\widetilde X_H,\mathbb Z)$ for the covering map $p\colon \widetilde X_H \rightarrow X$ associated to a finite index subgroup $H < \pi_1X$, as well as the corresponding homology homomorphism. We illustrate an open-source implementation of the algorithm by using it to show that: (i) the degree $2$ homology group $H_2(\widetilde X_H,\mathbb Z)$ distinguishes between the homotopy types of the complements $X\subset \mathbb R^4$ of the spun Hopf link and Satoh's tube map of the welded Hopf link (these two complements having isomorphic fundamental groups and integral homology); (ii) the degree $1$ homology homomorphism $H_1(p^{-1}(B),\mathbb Z) \rightarrow H_1(\widetilde X_H,\mathbb Z)$ distinguishes between the homeomorphism types of the complements $X\subset \mathbb R^3$ of the granny knot and the reef knot, where $B\subset X$ is the knot boundary (these two complements again having isomorphic fundamental groups and integral homology).


title: Improving Robustness of Task Oriented Dialog Systems

abstract: Task oriented language understanding in dialog systems is often modeled using intents (task of a query) and slots (parameters for that task). Intent detection and slot tagging are, in turn, modeled using sentence classification and word tagging techniques respectively. Similar to adversarial attack problems with computer vision models discussed in existing literature, these intent-slot tagging models are often over-sensitive to small variations in input -- predicting different and often incorrect labels when small changes are made to a query, thus reducing their accuracy and reliability. However, evaluating a model's robustness to these changes is harder for language since words are discrete and an automated change (e.g. adding `noise') to a query sometimes changes the meaning and thus labels of a query. In this paper, we first describe how to create an adversarial test set to measure the robustness of these models. Furthermore, we introduce and adapt adversarial training methods as well as data augmentation using back-translation to mitigate these issues. Our experiments show that both techniques improve the robustness of the system substantially and can be combined to yield the best results.


title: A LASSO-Inspired Approach for Localizing Power System Infeasibility

abstract: A method is proposed to identify and localize the cause of network collapse with augmented power flow analysis for a grid model with insufficient resources. Owing to heavy network loading, insufficient generation, component failures and unavoidable disturbances, power grid models can sometimes fail to converge to a feasible solution for a steady-state power flow study. For applications such as system expansion planning, it is desirable to locate the system buses that are contributing to network infeasibilities to facilitate corrective actions. This paper proposes a novel LASSO-inspired regularization of the power flow matrix that enforces sparsity to localize and quantify infeasibilities in the network. One of the examples demonstrates how the proposed method is capable of localizing a source of blackout to a single dominant bus on an 80k+ bus eastern interconnection model.


title: A Circuit-Theoretic Approach to State Estimation

abstract: Traditional state estimation (SE) methods that are based on nonlinear minimization of the sum of localized measurement error functionals are known to suffer from non-convergence and large residual errors. In this paper we propose an equivalent circuit formulation (ECF)-based SE approach that inherently considers the complete network topology and associated physical constraints. We analyze the mathematical differences between the two approaches and show that our approach produces a linear state-estimator that is mathematically a quadratic programming (QP) problem with closed-form solution. Furthermore, this formulation imposes additional topology-based constraints that provably shrink the feasible region and promote convergence to a more physically meaningful solution. From a probabilistic viewpoint, we show that our method applies prior knowledge into the estimate, thus converging to a more physics-based estimate than the traditional observation-driven maximum likelihood estimator (MLE). Importantly, incorporation of the entire system topology and underlying physics, while being linear, makes ECF-based SE advantageous for large-scale systems.


title: Large-scale optical interferometry in general spacetimes

abstract: We introduce a convenient formalism to evaluate the frequency-shift affecting a light signal propagating on a general curved background. Our formulation, which is based on the laws of geometric optics in a general relativistic setting, allows to obtain a transparent generalization of the Doppler frequency-shift without requiring to perform Local Lorentz transformations. It is easily applicable to stationary spacetimes, and in particular to the near-Earth experiments where geometry is described in the parametrized post-Newtonian approximation. We apply our recipe to evaluate the phase-shift arising in large-scale optical interferometric experiments, as the optical version of the Colella-Overhauser-Werner experiment.


title: Variants of epigroups and primary conjugacy

abstract: In a semigroup $S$ with fixed $c\in S$, one can construct a new semigroup $(S,\cdot_c)$ called a \emph{variant} by defining $x\cdot_c y:=xcy$. Elements $a,b\in S$ are \emph{primarily conjugate} if there exist $x,y\in S^1$ such that $a=xy, b=yx$. This coincides with the usual conjugacy in groups, but is not transitive in general semigroups. Ara\'{u}jo \emph{et al.} proved that transitivity holds in a variety $\mathcal{W}$ of epigroups containing all completely regular semigroups and their variants, and asked if transitivity holds for all variants of semigroups in $\mathcal{W}$. We answer this affirmatively as part of a study of varieties and variants of epigroups.


title: Cloudlet capture by Transitional Disk and FU Orionis stars

abstract: After its formation, a young star spends some time traversing the molecular cloud complex in which it was born. It is therefore not unlikely that, well after the initial cloud collapse event which produced the star, it will encounter one or more low mass cloud fragments, which we call "cloudlets" to distinguish them from full-fledged molecular clouds. Some of this cloudlet material may accrete onto the star+disk system, while other material may fly by in a hyperbolic orbit. In contrast to the original cloud collapse event, this process will be a "cloudlet flyby" and/or "cloudlet capture" event: A Bondi-Hoyle-Lyttleton type accretion event, driven by the relative velocity between the star and the cloudlet. As we will show in this paper, if the cloudlet is small enough and has an impact parameter similar or less than $GM_{*}/v^2_\infty$ (with $v_\infty$ being the approach velocity), such a flyby and/or capture event would lead to arc-shaped or tail-shaped reflection nebulosity near the star. Those shapes of reflection nebulosity can be seen around several transitional disks and FU Orionis stars. Although the masses in the those arcs appears to be much less than the disk masses in these sources, we speculate that higher-mass cloudlet capture events may also happen occasionally. If so, they may lead to the tilting of the outer disk, because the newly infalling matter will have an angular momentum orientation entirely unrelated to that of the disk. This may be one possible explanation for the highly warped/tilted inner/outer disk geometries found in several transitional disks. We also speculate that such events, if massive enough, may lead to FU Orionis outbursts.


title: Coordination Group Formation for OnLine Coordinated Routing Mechanisms

abstract: This study considers that the collective route choices of travelers en route represent a resolution of their competition on network routes. Well understanding this competition and coordinating their route choices help mitigate urban traffic congestion. Even though existing studies have developed such mechanisms (e.g., the CRM [1]), we still lack the quantitative method to evaluate the coordination penitential and identify proper coordination groups (CG) to implement the CRM. Thus, they hit prohibitive computing difficulty when implemented with many opt-in travelers. Motived by this view, this study develops mathematical approaches to quantify the coordination potential between two and among multiple travelers. Next, we develop the adaptive centroid-based clustering algorithm (ACCA), which splits travelers en route in a local network into CGs, each with proper size and strong coordination potential. Moreover, the ACCA is statistically secured to stop at a local optimal clustering solution, which balances the inner-cluster and inter-cluster coordination potential. It can be implemented by parallel computation to accelerate its computing efficiency. Furthermore, we propose a clustering based coordinated routing mechanism (CB-CRM), which implements a CRM on each individual CG. The numerical experiments built upon both Sioux Falls and Hardee city networks show that the ACCA works efficiently to form proper coordination groups so that as compared to the CRM, the CB-CRM significantly improves computation efficiency with minor system performance loss in a large network. This merit becomes more apparent under high penetration and congested traffic condition. Last, the experiments validate the good features of the ACCA as well as the value of implementing parallel computation.


title: Modeling The Temporally Constrained Preemptions of Transient Cloud VMs

abstract: Transient cloud servers such as Amazon Spot instances, Google Preemptible VMs, and Azure Low-priority batch VMs, can reduce cloud computing costs by as much as $10\times$, but can be unilaterally preempted by the cloud provider. Understanding preemption characteristics (such as frequency) is a key first step in minimizing the effect of preemptions on application performance, availability, and cost. However, little is understood about temporally constrained preemptions---wherein preemptions must occur in a given time window. We study temporally constrained preemptions by conducting a large scale empirical study of Google's Preemptible VMs (that have a maximum lifetime of 24 hours), develop a new preemption probability model, new model-driven resource management policies, and implement them in a batch computing service for scientific computing workloads. Our statistical and experimental analysis indicates that temporally constrained preemptions are not uniformly distributed, but are time-dependent and have a bathtub shape. We find that existing memoryless models and policies are not suitable for temporally constrained preemptions. We develop a new probability model for bathtub preemptions, and analyze it through the lens of reliability theory. To highlight the effectiveness of our model, we develop optimized policies for job scheduling and checkpointing. Compared to existing techniques, our model-based policies can reduce the probability of job failure by more than $2\times$. We also implement our policies as part of a batch computing service for scientific computing applications, which reduces cost by $5\times$ compared to conventional cloud deployments and keeps performance overheads under $3\%$.


title: All It Takes is 20 Questions!: A Knowledge Graph Based Approach

abstract: 20 Questions (20Q) is a two-player game. One player is the answerer, and the other is a questioner. The answerer chooses an entity from a specified domain and does not reveal this to the other player. The questioner can ask at most 20 questions to the answerer to guess the entity. The answerer can reply to the questions asked by saying yes/no/maybe. In this paper, we propose a novel approach based on the knowledge graph for designing the 20Q game on Bollywood movies. The system assumes the role of the questioner and asks questions to predict the movie thought by the answerer. It uses a probabilistic learning model for template-based question generation and answers prediction. A dataset of interrelated entities is represented as a weighted knowledge graph, which updates as the game progresses by asking questions. An evolutionary approach helps the model to gain a better understanding of user choices and predicts the answer in fewer questions over time. Experimental results show that our model was able to predict the correct movie in less than 10 questions for more than half of the times the game was played. This kind of model can be used to design applications that can detect diseases by asking questions based on symptoms, improving recommendation systems, etc.


title: Lightning Impulse Current Tests on some Electroconductive Fabrics

abstract: On the search of lightweight lightning protection materials that can be used as part of lightning protection systems, we investigate some types of electroconductive fabrics by applying several lightning impulse currents in laboratory. Samples of four commercially available electroconductive textiles were analyzed, two rip-stop, a plain-weave, a nonwoven, and additionally a carbon-impregnated polymeric film. Under laboratory conditions, each sample was subjected to several lengthwise subsequent lightning-like currents of 8/20 us standard waveform, recording both voltage and current signals. Optical and scanning electron microscope observations were performed after tests, revealing some patterns or morphological changes on the fabric surface. Despite these changes, the investigated conductive textiles withstand the several lightning impulse currents applied. Results suggest that some conductive fabrics could be used in personal mobile shelters, to protect human beings against the earth potential rise caused by a close lightning discharge.


title: A Multireference Quantum Krylov Algorithm for Strongly Correlated Electrons

abstract: We introduce a multireference selected quantum Krylov (MRSQK) algorithm suitable for quantum simulation of many-body problems. MRSQK is a low-cost alternative to the quantum phase estimation algorithm that generates a target state as a linear combination of non-orthogonal Krylov basis states. This basis is constructed from a set of reference states via real-time evolution avoiding the numerical optimization of parameters. An efficient algorithm for the evaluation of the off-diagonal matrix elements of the overlap and Hamiltonian matrices is discussed and a selection procedure is introduced to identify a basis of orthogonal references that ameliorates the linear dependency problem. Preliminary benchmarks on linear H$_6$, H$_8$, and BeH$_2$ indicate that MRSQK can predict the energy of these systems accurately using very compact Krylov bases.


title: A Reproducibility Study of "IP Spoofing Detection in Inter-Domain Traffic"

abstract: IP spoofing enables reflection and amplification attacks, which cause major threats to the current Internet infrastructure. Detecting IP packets with incorrect source addresses would help to improve the situation. This is easy at the attacker's network, but very challenging at Internet eXchange Points (IXPs) or in transit networks. In this reproducibility study, we revisit the paper \textit{Detection, Classification, and Analysis of Inter-Domain Traffic with Spoofed Source IP Addresses} published at ACM IMC 2017. Using data from a different IXP and from a different time, we were not able to reproduce the results. Unfortunately, our further analysis reveals structural problems of the state of the art methodology, which are not easy to overcome.


title: Classifying Relevant Social Media Posts During Disasters Using Ensemble of Domain-agnostic and Domain-specific Word Embeddings

abstract: The use of social media as a means of communication has significantly increased over recent years. There is a plethora of information flow over the different topics of discussion, which is widespread across different domains. The ease of information sharing has increased noisy data being induced along with the relevant data stream. Finding such relevant data is important, especially when we are dealing with a time-critical domain like disasters. It is also more important to filter the relevant data in a real-time setting to timely process and leverage the information for decision support. However, the short text and sometimes ungrammatical nature of social media data challenge the extraction of contextual information cues, which could help differentiate relevant vs. non-relevant information. This paper presents a novel method to classify relevant social media posts during disaster events by ensembling the features of both domain-specific word embeddings as well as more generic domain-agnostic word embeddings. Therefore, we develop and evaluate a hybrid feature engineering framework for integrating diverse semantic representations using a combination of word embeddings to efficiently classify a relevant social media post. The application of the proposed classification framework could help in filtering public posts at large scale, given the growing usage of social media posts in recent years.


title: Negative sampling in semi-supervised learning

abstract: We introduce Negative Sampling in Semi-Supervised Learning (NS3L), a simple, fast, easy to tune algorithm for semi-supervised learning (SSL). NS3L is motivated by the success of negative sampling/contrastive estimation. We demonstrate that adding the NS3L loss to state-of-the-art SSL algorithms, such as the Virtual Adversarial Training (VAT), significantly improves upon vanilla VAT and its variant, VAT with Entropy Minimization. By adding the NS3L loss to MixMatch, the current state-of-the-art approach on semi-supervised tasks, we observe significant improvements over vanilla MixMatch. We conduct extensive experiments on the CIFAR10, CIFAR100, SVHN and STL10 benchmark datasets.


title: Nonconvex Stochastic Nested Optimization via Stochastic ADMM

abstract: We consider the stochastic nested composition optimization problem where the objective is a composition of two expected-value functions. We proposed the stochastic ADMM to solve this complicated objective. In order to find an $\epsilon$ stationary point where the expected norm of the subgradient of corresponding augmented Lagrangian is smaller than $\epsilon$, the total sample complexity of our method is $\mathcal{O}(\epsilon^{-3})$ for the online case and $\mathcal{O} \Bigl((2N_1 + N_2) + (2N_1 + N_2)^{1/2}\epsilon^{-2}\Bigr)$ for the finite sum case. The computational complexity is consistent with proximal version proposed in \cite{zhang2019multi}, but our algorithm can solve more general problem when the proximal mapping of the penalty is not easy to compute.


title: Design and Implementation of a Three-Link Brachiation Robot with Optimal Control Based Trajectory Tracking Controller

abstract: This paper reports the design and implementation of a three-link brachiation robot. The robot is able to travel along horizontal monkey bars using continuous arm swings. We build a full order dynamics model for the robot and formulate each cycle of robot swing motion as an optimal control problem. The iterative Linear Quadratic Regulator (iLQR) algorithm is used to find the optimal control strategy during one swing. We select suitable robot design parameters by comparing the cost of robot motion generated by the iLQR algorithm for different robot designs. In particular, using this approach we show the importance of having a body link and low inertia arms for efficient brachiation. Further, we propose a trajectory tracking controller that combines a cascaded PID controller and an input-output linearization controller to enable the robot to track desired trajectory precisely and reject external disturbance during brachiation. Experiments on the simulated robot and the real robot demonstrate that the robot can robustly swing between monkey bars with same or different spacing of handholds.


title: New bounds on the spectral radius of graphs based on the moment problem

abstract: Let $\mathcal{G}$ be an undirected graph with adjacency matrix $A$ and spectral radius $\rho$. Let $w_k, \phi_k$ and $\phi_k^{(i)}$ be, respectively, the number walks of length $k$, closed walks of length $k$ and closed walks starting and ending at vertex $i$ after $k$ steps. In this paper, we propose a measure-theoretic framework which allows us to relate walks in a graph with its spectral properties. In particular, we show that $w_k, \phi_k$ and $\phi_k^{(i)}$ can be interpreted as the moments of three different measures, all of them supported on the spectrum of $A$. Building on this interpretation, we leverage results from the classical moment problem to formulate a hierarchy of new lower and upper bounds on $\rho$, as well as provide alternative proofs to several well-known bounds in the literature.


title: The Power of Two Choices for Random Walks

abstract: We apply the power-of-two-choices paradigm to a random walk on a graph: rather than moving to a uniform random neighbour at each step, a controller is allowed to choose from two independent uniform random neighbours. We prove that this allows the controller to significantly accelerate the hitting and cover times in several natural graph classes. In particular, we show that the cover time becomes linear in the number $n$ of vertices on discrete tori and bounded degree trees, of order $\mathcal{O}(n \log \log n)$ on bounded degree expanders, and of order $\mathcal{O}(n (\log \log n)^2)$ on the Erd\H{o}s-R\'{e}nyi random graph in a certain sparsely connected regime. We also consider the algorithmic question of computing an optimal strategy, and prove a dichotomy in efficiency between computing strategies for hitting and cover times.


title: Incentive Compatible Active Learning

abstract: We consider active learning under incentive compatibility constraints. The main application of our results is to economic experiments, in which a learner seeks to infer the parameters of a subject's preferences: for example their attitudes towards risk, or their beliefs over uncertain events. By cleverly adapting the experimental design, one can save on the time spent by subjects in the laboratory, or maximize the information obtained from each subject in a given laboratory session; but the resulting adaptive design raises complications due to incentive compatibility. A subject in the lab may answer questions strategically, and not truthfully, so as to steer subsequent questions in a profitable direction. We analyze two standard economic problems: inference of preferences over risk from multiple price lists, and belief elicitation in experiments on choice over uncertainty. In the first setting, we tune a simple and fast learning algorithm to retain certain incentive compatibility properties. In the second setting, we provide an incentive compatible learning algorithm based on scoring rules with query complexity that differs from obvious methods of achieving fast learning rates only by subpolynomial factors. Thus, for these areas of application, incentive compatibility may be achieved without paying a large sample complexity price.


title: Orbifolds from Modular Orbits

abstract: Given a two-dimensional conformal field theory with a global symmetry, we propose a method to implement an orbifold construction by taking orbits of the modular group. For the case of cyclic symmetries we find that this approach always seems to be consistent, even in asymmetric orbifold cases where the usual construction does not yield a modular invariant theory; our approach keeps modular invariance manifest but may give a result that is equivalent to the original theory. For the case that the symmetry is a subgroup of a continuous flavor symmetry, we can give explicit constructions of the spectrum, with twisted sectors corresponding to a non-standard group projection on an enlarged twisted sector Hilbert space.


title: Renormalization group flow of coupled tensorial group field theories: Towards the Ising model on random lattices

abstract: We introduce a new family of tensorial field theories by coupling different fields in a non-trivial way, with a view towards the investigation of the coupling between matter and gravity in the quantum regime. As a first step, we consider the simple case with two tensors of the same rank coupled together, with Dirac like kinetic kernel. We focus especially on rank-$3$ tensors, which lead to a power counting just-renormalizable model, and interpret Feynman graphs as Ising configurations on random lattices. We investigate the renormalization group flow for this model, using two different and complementary tools for approximations, namely, the effective vertex expansion method and finite-dimensional truncations for the flowing action. Due to the complicated structure of the resulting flow equations, we divided the work into two parts. In this first part we only investigate the fundamental aspects on the construction of the model and the different ways to get tractable renormalization group equations, while their numerical analysis will be addressed in a companion paper.


title: Star Formation Rate in late-type galaxies: I- The H(alpha) and FUV integrated values

abstract: The Star Formation Rate (SFR) was determined from a sample of 36 late-type galaxies (14 dS and 22 Sm), from the H(alpha) and Far Ultraviolet (FUV) flux. We found that the SFR(FUV) was in most cases higher than the SFR(H(alpha)).We also obtained that the SFR is larger for Sm galaxies and smaller for barred galaxies, for any morphological type, in both diagnostic methods. In addition, a study of the spatial distribution of star formation within these galaxies was made, concluding that there is not a preferential place for the star formation. Finally, we studied the role of the Diffuse Ionized Gas in the SFR value, finding that the H(alpha) flux contribution coming from this gas should not be included in the SFR determination unless it is verified that leaking photons are the only source of the neutral gas ionization.


title: Influence of structural discontinuities present in the core of red-giant stars on the observed mixed-mode pattern and characterization of their properties

abstract: The space-borne missions CoRoT and Kepler have provided seismic data of unprecedented quality. Among the observed stars, red giants show a complex oscillation pattern exhibiting pressure modes as well as mixed modes. The latter carry information on the radiative region properties of these stars. The very high precision of Kepler data provide enough accuracy to decipher the complex structure of the mixed-mode pattern and, therefore, deduce precise information on the structure of the stellar core. In this work, we studied the precise influence of the core structural discontinuities on the mixed-mode pattern. These phenomena have indeed an influence on the gravity waves that propagate inside the core of the star and, therefore, on the mixed-mode pattern. In this study, we aim to investigate the impact of core structural discontinuities on the observed mixed-mode pattern and measure their properties for several red giants. We identified several objects showing evidence of deviations in their mixed-mode frequency pattern that are characteristic of core structural discontinuities. We fitted these deviations and showed that they likely correspond to the influence of the inner convective core of these stars.


title: State Estimation for Legged Robots Using Contact-Centric Leg Odometry

abstract: Our goal is to send legged robots into challenging, unstructured terrains that wheeled systems cannot traverse. Moreover, precise estimation of the robot's position and orientation in rough terrain is especially difficult. To address this problem, we introduce a new state estimation algorithm which we term Contact-Centric Leg Odometry (COCLO). This new estimator uses a Square Root Unscented Kalman Filter (SR-UKF) to fuse multiple proprioceptive sensors available on a legged robot. In contrast to IMU-centric filtering approaches, COCLO formulates prediction and measurement models according to the contact status of legs. Additionally, COCLO has an indirect measurement model using joint velocities to estimate the robot body velocity. In rough terrain, when IMUs suffer from large amounts of noise, COCLO's contact-centric approach outperforms previous IMU-centric methods. To demonstrate improved state estimation accuracy, we compare COCLO with Visual Inertial Navigation System (VINS), a state-of-the-art visual inertial odometry, in three different environments: flat ground, ramps, and stairs. COCLO achieves better estimation precision than VINS in all three environments and is robust to unstable motion. Finally, we also show that COCLO and a modified VINS can work in tandem to improve each other's performance.


title: Numerical and experimental realization of analytical SLAM

abstract: Analytical approach to SLAM problem was introduced in the recent years. In our work we investigate the method numerically with the motivation of using the algorithm in a real hardware experiments. We perform a robustness test of the algorithm and apply it to the robotic hardware in two different setups. In one we try to recover a map of the environment using bearing angle measurements and radial distance measurements. The another setup utilizes only bearing angle information.


title: Spectra of Ni V and Fe V in the Vacuum Ultraviolet

abstract: This work presents 97 remeasured Fe V wavelengths (1200 \r{A} to 1600 \r{A}) and 123 remeasured Ni V wavelengths (1200 \r{A} to 1400 \r{A}) with uncertainties of approximately 2m\r{A}. An additional 67 remeasured Fe V wavelengths and 72 remeasured Ni V wavelengths with uncertainties greater than 2m\r{A} are also reported. A systematic calibration error is also identified in the previous Ni V wavelengths and is corrected in this work. Furthermore, a new energy level optimization of Ni V is presented that includes level values as well as Ritz wavelengths. This work improves upon the available data used for observations of quadruply ionized nickel (Ni V) in white dwarf stars. This compilation is specifically targeted towards observations of the G191-B2B white dwarf spectrum that has been used to test for variations in the fine structure constant, $\alpha$, in the presence of strong gravitational fields. The laboratory wavelengths for these ions were thought to be the cause of inconsistent conclusions regarding the variation limit of $\alpha$ as observed through the white dwarf spectrum. These inconsistencies can now be addressed with the improved laboratory data presented here.


title: WASP-52b. The effect of starspot correction on atmospheric retrievals

abstract: We perform atmospheric retrievals on the full optical to infrared ($0.3-5 \, \mu \mathrm{m}$) transmission spectrum of the inflated hot Jupiter WASP-52b by combining HST/STIS, WFC3 IR, and Spitzer/IRAC observations. As WASP-52 is an active star which shows both out-of-transit photometric variability and starspot crossings during transits, we account for the contribution of non-occulted active regions in the retrieval. We recover a $0.1-10\times$ solar atmospheric composition, in agreement with core accretion predictions for giant planets, and a weak contribution of aerosols. We also obtain a $<3000$ K temperature for the starspots, a measure which is likely affected by the models used to fit instrumental effects in the transits, and a 5% starspot fractional coverage, compatible with expectations for the host star's spectral type. Such constraints on the planetary atmosphere and on the activity of its host star will inform future JWST GTO observations of this target.


title: Turbulence forecasting via Neural ODE

abstract: Fluid turbulence is characterized by strong coupling across a broad range of scales. Furthermore, besides the usual local cascades, such coupling may extend to interactions that are non-local in scale-space. As such the computational demands associated with explicitly resolving the full set of scales and their interactions, as in the Direct Numerical Simulation (DNS) of the Navier-Stokes equations, in most problems of practical interest are so high that reduced modeling of scales and interactions is required before further progress can be made. While popular reduced models are typically based on phenomenological modeling of relevant turbulent processes, recent advances in machine learning techniques have energized efforts to further improve the accuracy of such reduced models. In contrast to such efforts that seek to improve an existing turbulence model, we propose a machine learning(ML) methodology that captures, de novo, underlying turbulence phenomenology without a pre-specified model form. To illustrate the approach, we consider transient modeling of the dissipation of turbulent kinetic energy, a fundamental turbulent process that is central to a wide range of turbulence models using a Neural ODE approach. After presenting details of the methodology, we show that this approach outperforms state-of-the-art approaches.


title: 92c/MFlops/s, Ultra-Large-Scale Neural-Network Training on a PIII Cluster

abstract: Artificial neural networks with millions of adjustable parameters and a similar number of training examples are a potential solution for difficult, large-scale pattern recognition problems in areas such as speech and face recognition, classification of large volumes of web data, and finance. The bottleneck is that neural network training involves iterative gradient descent and is extremely computationally intensive. In this paper we present a technique for distributed training of Ultra Large Scale Neural Networks (ULSNN) on Bunyip, a Linux-based cluster of 196 Pentium III processors. To illustrate ULSNN training we describe an experiment in which a neural network with 1.73 million adjustable parameters was trained to recognize machine-printed Japanese characters from a database containing 9 million training patterns. The training runs with a average performance of 163.3 GFlops/s (single precision). With a machine cost of \$150,913, this yields a price/performance ratio of 92.4c/MFlops/s (single precision). For comparison purposes, training using double precision and the ATLAS DGEMM produces a sustained performance of 70 MFlops/s or \$2.16 / MFlop/s (double precision).


title: A Proof of Principle: Multi-Modality Radiotherapy Optimization

abstract: Radiotherapy is used to treat cancer patients by damaging DNA of tumor cells using ionizing radiation. Photons are the most widely used radiation type for therapy, having been put into use soon after the first discovery of X-rays in 1895. However, there are emerging interests and developments of other radiation modalities such as protons and carbon ions, owing to their unique biological and physical characteristics that distinguish these modalities from photons. Current attempts to determine an optimal radiation modality or an optimal combination of multiple modalities are empirical and in the early stage of development. In this paper, we propose a mathematical framework to optimize full radiation dose distributions and fractionation schedules of multiple radiation modalities, aiming to maximize the damage to the tumor while limiting the damage to the normal tissue to the corresponding tolerance level. This formulation gives rise to a non-convex, mixed integer program and we propose a bilevel optimization algorithm, to efficiently solve it. The upper level problem is to optimize the fractionation schedule using the dose distribution optimized in the lower level. We demonstrate the feasibility of our novel framework and algorithms in a simple 2-dimensional phantom with two different radiation modalities, where clinical intuition can be easily drawn. The results of our numerical simulations agree with the clinical intuition, validating our approach and showing the promise of the framework for further clinical investigation.


title: Precursors to plastic failure in a numerical simulation of CuZr metallic glass

abstract: We deform, in pure shear, a thin sample of Cu$_{50}$Zr$_{50}$ metallic glass using a molecular dynamics simulation up to, and including, failure. The experiment is repeated ten times in order to have average values and standard deviations. Although failure occurs at the same value of the externally imposed strain for the ten samples, there is significant sample-to-sample variation in the specific microscopic material behavior. Failure can occur along one, two, or three planes, located at the boundaries of previously formed shear bands. These shear bands form shortly before failure. However, well before their formation and at external strains where plastic deformation just begins to be significant, non-affine displacement organizes itself along localized bands. The shear bands subsequently form at the edges of these non-affine-displacement-bands, and present an alternating rotation-quadrupole structure, as found previously by \c{S}opu et al. (2017) in the case of a notched sample loaded in tension. The thickness of shear bands is roughly determined by the available plastic energy. The onset of shear banding is accompanied by a sharp increase in the rate of change of the rotation angle localization, the strain localization, and the non-affine square displacement.


title: CHEETAH: An Ultra-Fast, Approximation-Free, and Privacy-Preserved Neural Network Framework based on Joint Obscure Linear and Nonlinear Computations

abstract: Machine Learning as a Service (MLaaS) is enabling a wide range of smart applications on end devices. However, such convenience comes with a cost of privacy because users have to upload their private data to the cloud. This research aims to provide effective and efficient MLaaS such that the cloud server learns nothing about user data and the users cannot infer the proprietary model parameters owned by the server. This work makes the following contributions. First, it unveils the fundamental performance bottleneck of existing schemes due to the heavy permutations in computing linear transformation and the use of communication intensive Garbled Circuits for nonlinear transformation. Second, it introduces an ultra-fast secure MLaaS framework, CHEETAH, which features a carefully crafted secret sharing scheme that runs significantly faster than existing schemes without accuracy loss. Third, CHEETAH is evaluated on the benchmark of well-known, practical deep networks such as AlexNet and VGG-16 on the MNIST and ImageNet datasets. The results demonstrate more than 100x speedup over the fastest GAZELLE (Usenix Security'18), 2000x speedup over MiniONN (ACM CCS'17) and five orders of magnitude speedup over CryptoNets (ICML'16). This significant speedup enables a wide range of practical applications based on privacy-preserved deep neural networks.


title: Pose estimation and bin picking for deformable products

abstract: Robotic systems in manufacturing applications commonly assume known object geometry and appearance. This simplifies the task for the 3D perception algorithms and allows the manipulation to be more deterministic. However, those approaches are not easily transferable to the agricultural and food domains due to the variability and deformability of natural food. We demonstrate an approach applied to poultry products that allows picking up a whole chicken from an unordered bin using a suction cup gripper, estimating its pose using a Deep Learning approach, and placing it in a canonical orientation where it can be further processed. Our robotic system was experimentally evaluated and is able to generalize to object variations and achieves high accuracy on bin picking and pose estimation tasks in a real-world environment.


title: TCT: A Cross-supervised Learning Method for Multimodal Sequence Representation

abstract: Multimodalities provide promising performance than unimodality in most tasks. However, learning the semantic of the representations from multimodalities efficiently is extremely challenging. To tackle this, we propose the Transformer based Cross-modal Translator (TCT) to learn unimodal sequence representations by translating from other related multimodal sequences on a supervised learning method. Combined TCT with Multimodal Transformer Network (MTN), we evaluate MTN-TCT on the video-grounded dialogue which uses multimodality. The proposed method reports new state-of-the-art performance on video-grounded dialogue which indicates representations learned by TCT are more semantics compared to directly use unimodality.