title: Comet 66P/du Toit: not a near Earth main belt comet

abstract: Main belt comets (MBCs) are a peculiar class of volatile-containing objects with comet-like morphology and asteroid-like orbits. However, MBCs are challenging targets to study remotely due to their small sizes and the relatively large distance they are from the Sun and the Earth. Recently, a number of weakly active short-period comets have been identified that might originate in the asteroid main belt. Among all of the known candidates, comet 66P/du Toit has been suggested to have one of the highest probabilities of coming from the main belt. We obtained medium and high-resolution spectra of 66P from 300-2500 nm with the X-shooter/VLT and the UVES/VLT instruments in July 2018. We also obtained a series of narrow-band images of 66P to monitor the gas and dust activity between May and July 2018 with TRAPPIST-South. In addition, we applied a dust model to characterize the dust coma of 66P and performed dynamical simulations to study the orbital evolution of 66P. We derive the OPR of ammonia (NH$_3$) in 66P to be 1.08$\pm$0.06, which corresponds to a nuclear spin temperature of $\sim$34 K. We computed the production rates of OH, NH, CN, C$_3,$ and C$_2$ radicals and measured the dust proxy, Af$\rho$. The dust analysis reveals that the coma can be best-fit with an anisotropic model and the peak dust production rate is about 55 kg s$^{-1}$ at the perihelion distance of 1.29 au. Dynamical simulations show that 66P is moderately asteroidal with the capture time, t$_{cap} \sim 10^4$ yr. Our observations demonstrate that the measured physical properties of 66P are consistent with other typical short-period comets and differ significantly from other MBCs. Therefore, 66P is unlikely to have a main belt origin.


title: Recent Advances in End-to-End Spoken Language Understanding

abstract: This work investigates spoken language understanding (SLU) systems in the scenario when the semantic information is extracted directly from the speech signal by means of a single end-to-end neural network model. Two SLU tasks are considered: named entity recognition (NER) and semantic slot filling (SF). For these tasks, in order to improve the model performance, we explore various techniques including speaker adaptation, a modification of the connectionist temporal classification (CTC) training criterion, and sequential pretraining.


title: Chow quotients of Grassmannians by diagonal subtori

abstract: The literature on maximal torus orbits in the Grassmannian is vast; in this paper we initiate a program to extend this to diagonal subtori. Our main focus is generalizing portions of Kapranov's seminal work on Chow quotient compactifications of these orbit spaces. This leads naturally to discrete polymatroids, generalizing the matroidal framework underlying Kapranov's results. By generalizing the Gelfand-MacPherson isomorphism, these Chow quotients are seen to compactify spaces of arrangements of parameterized linear subspaces, and a generalized Gale duality holds here. A special case is birational to the Chen-Gibney-Krashen moduli space of pointed trees of projective spaces, and we show that the question of whether this birational map is an isomorphism is a specific instance of a much more general question that hasn't previously appeared in the literature, namely, whether the geometric Borel transfer principle in non-reductive GIT extends to an isomorphism of Chow quotients.


title: Symplectic Recurrent Neural Networks

abstract: We propose Symplectic Recurrent Neural Networks (SRNNs) as learning algorithms that capture the dynamics of physical systems from observed trajectories. An SRNN models the Hamiltonian function of the system by a neural network and furthermore leverages symplectic integration, multiple-step training and initial state optimization to address the challenging numerical issues associated with Hamiltonian systems. We show that SRNNs succeed reliably on complex and noisy Hamiltonian systems. We also show how to augment the SRNN integration scheme in order to handle stiff dynamical systems such as bouncing billiards.


title: Angles of Random Simplices and Face Numbers of Random Polytopes

abstract: Pick $d+1$ points uniformly at random on the unit sphere in $\mathbb R^d$. What is the expected value of the angle sum of the simplex spanned by these points? Choose $n$ points uniformly at random in the $d$-dimensional ball. What is the expected number of faces of their convex hull? We answer these and some related questions of stochastic geometry. To this end, we compute expected internal angles of random simplices whose vertices are independent random points sampled from one of the following $d$-dimensional distributions: (i) the beta distribution with the density proportional to $(1-\|x\|^2)^{\beta}$, where $x$ is belongs to the unit ball in $\mathbb R^d$; (ii) the beta' distribution with the density proportional to $(1+\|x\|^2)^{-\beta}$, where $x\in\mathbb{R}^{d}$. These results imply explicit formulae for the expected face numbers of the following random polytopes: (a) the typical Poisson-Voronoi cell; (b) the zero cell of the Poisson hyperplane tessellation; (c) beta and beta' polytopes defined as convex hulls of i.i.d. samples from the corresponding distributions.


title: Phase diagram and thermal Hall conductivity of spin-liquid Kekul\'{e}-Kitaev model

abstract: In this work we study the phase diagram of Kekul\'{e}-Kitaev model. The model is defined on a honeycomb lattice with bond dependent anisotropic exchange interactions making it exactly solvable in terms of Majorana representation of spins in close analogy to the Kitaev model. However, the energy spectrum of Majorana fermions has a multi-band structure characterized by Chern numbers 0, $\pm$1, and $\pm2$. We obtained the phase diagram of the model in the plane of exchange couplings and in the presence of a magnetic field and found chiral topological and trivial spin-liquid ground states. In the absence of magnetic field most part of the phase diagram is a trivial gapped phase continuously connected to an Abelian phase, while in the presence of the magnetic field a topological phase arises. Furthermore, motivated by recent thermal measurements on the spin-liquid candidate $\alpha$-RuCl$_{3}$, we calculated the thermal Hall conductivity at different regimes of parameters and temperatures and found the latter is quantized over a wide range of temperatures.


title: Futures-Based Resource Trading and Fair Pricing in Real-Time IoT Networks

abstract: Technology advancements on sensing, communications, and computing directly accelerate the recent development of Internet of Things (IoT), leading to rich and diverse applications in industry and business processes. For enabling IoT applications, smart devices and networks rely on the capabilities for information collection, processing and tight collaboration with their neighbors. However, the rapid proliferation of IoT devices and real-time IoT applications result in dramatically increased demand for communication and computing resources, posing a significant challenge for the underlying resource constrained wireless networks. Due to their resource constraint, it is extremely difficult for such devices and networks to meet stringent requirements of many real-time IoT applications. To overcome these challenges, resource sharing among IoT devices and networks has been studied as an effective solution by exploiting unused resources.


title: Real space Berry curvature of itinerant electron systems with spin-orbit interaction

abstract: By considering an extended double-exchange model with spin-orbit coupling (SOC), we derive a general form of the Berry phase $\gamma$ that electrons pick up when moving around a closed loop. This form generalizes the well-known result valid for SU(2) invariant systems, $\gamma=\Omega/2$, where $\Omega$ is the solid angle subtended by the local magnetic moments enclosed by the loop. The general form of $\gamma$ demonstrates that collinear and coplanar magnetic textures can also induce a Berry phase different from 0 or $\pi$, smoothly connecting the result for SU(2) invariant systems with the well-known result of Karplus and Luttinger for collinear ferromagnets with finite SOC. By taking the continuum limit of the theory, we also derive the corresponding generalized form of the real space Berry curvature. The new expression is a generalization of the scalar spin chirality, which is presented in an explicitly covariant form. We finally show how these simple concepts can be used to understand the origin of the spontaneous topological Hall effect that has been recently reported in collinear and coplanar antiferromagnetic phases of correlated materials.


title: MMD-Bayes: Robust Bayesian Estimation via Maximum Mean Discrepancy

abstract: In some misspecified settings, the posterior distribution in Bayesian statistics may lead to inconsistent estimates. To fix this issue, it has been suggested to replace the likelihood by a pseudo-likelihood, that is the exponential of a loss function enjoying suitable robustness properties. In this paper, we build a pseudo-likelihood based on the Maximum Mean Discrepancy, defined via an embedding of probability distributions into a reproducing kernel Hilbert space. We show that this MMD-Bayes posterior is consistent and robust to model misspecification. As the posterior obtained in this way might be intractable, we also prove that reasonable variational approximations of this posterior enjoy the same properties. We provide details on a stochastic gradient algorithm to compute these variational approximations. Numerical simulations indeed suggest that our estimator is more robust to misspecification than the ones based on the likelihood.


title: Learning Everywhere: A Taxonomy for the Integration of Machine Learning and Simulations

abstract: We present a taxonomy of research on Machine Learning (ML) applied to enhance simulations together with a catalog of some activities. We cover eight patterns for the link of ML to the simulations or systems plus three algorithmic areas: particle dynamics, agent-based models and partial differential equations. The patterns are further divided into three action areas: Improving simulation with Configurations and Integration of Data, Learn Structure, Theory and Model for Simulation, and Learn to make Surrogates.


title: Rotation surfaces of constant Gaussian curvature as Riemannian approximation scheme in sub-Riemannian Heisenberg space $\mathbb H^1$

abstract: We verify if Gausssian curvature of surfaces and normal curvature of curves in surfaces introduced by Diniz-Veloso arXiv:1210.7110 and by Balogh-Tyson-Vecchi arXiv:1604.00180 to prove Gauss-Bonnet theorems in Heisenberg space $\mathbb H^1$ are equal. The authors in arXiv:1604.00180 utilize a limit of Gaussian and normal curvatures defined in the Riemannian approximations scheme $(\mathbb R^3,g_L)$ introduced by Gromov to study sub-Riemannian spaces. They show that these limits exist (unlike the limit of Riemannian surface area form or length form), and they obtain Gauss-Bonnet theorem in $\mathbb H^1$ as limit of Gauss-Bonnet theorems in $(\mathbb R^3,g_L)$ when $L$ goes to infinity. The approach in arXiv:1210.7110 uses an analogous of Gauss map defined in the unitary cylinder and they show that the curvature defined by the Gauss map corresponds to the curvature of covariant derivative which is null on left invariant vector fields of $\mathbb H^1$. The proof of Gauss-Bonnet theorem follows as the classical one. Curvatures do not coincide. For comparison sake we apply the same formalism of arXiv:1210.7110 to get the curvatures of arXiv:1604.00180 . With the obtained formulas, it is possible to prove the Gauss-Bonnet theorem in arXiv:1604.00180 as a straightforward application of Stokes theorem. To exemplify the Gaussian curvature of arXiv:1604.00180 we calculate the rotation surfaces of constant curvature in $\mathbb H^1$, which are only of three types.


title: Interference-Precancelled Pilot Design for LMMSE Channel Estimation of GFDM

abstract: Generalized frequency division multiplexing (GFDM) is a promising candidate waveform for next-generation wireless communication systems. However, GFDM channel estimation is still challenging due to the inherent interference. In this paper, we formulate a pilot design framework with linear minimum mean square error (LMMSE) channel estimation for GFDM, and propose a novel pilot design to achieve interference precancellation during pilot generation with the fixed transmit sample values at selected frequency bins. Numerical results demonstrate that the proposed method reduces the channel estimation mean square error and the symbol error rate (SER) in high signal-to-noise ratio (SNR) regions, compared with the conventional methods.


title: ISTHMUS: Secure, Scalable, Real-time and Robust Machine Learning Platform for Healthcare

abstract: In recent times, machine learning (ML) and artificial intelligence (AI) based systems have evolved and scaled across different industries such as finance, retail, insurance, energy utilities, etc. Among other things, they have been used to predict patterns of customer behavior, to generate pricing models, and to predict the return on investments. But the successes in deploying machine learning models at scale in those industries have not translated into the healthcare setting. There are multiple reasons why integrating ML models into healthcare has not been widely successful, but from a technical perspective, general-purpose commercial machine learning platforms are not a good fit for healthcare due to complexities in handling data quality issues, mandates to demonstrate clinical relevance, and a lack of ability to monitor performance in a highly regulated environment with stringent security and privacy needs. In this paper, we describe Isthmus, a turnkey, cloud-based platform which addresses the challenges above and reduces time to market for operationalizing ML/AI in healthcare. Towards the end, we describe three case studies which shed light on Isthmus capabilities. These include (1) supporting an end-to-end lifecycle of a model which predicts trauma survivability at hospital trauma centers, (2) bringing in and harmonizing data from disparate sources to create a community data platform for inferring population as well as patient level insights for Social Determinants of Health (SDoH), and (3) ingesting live-streaming data from various IoT sensors to build models, which can leverage real-time and longitudinal information to make advanced time-sensitive predictions.


title: Modeling Single Picker Routing Problems in Classical and Modern Warehouses

abstract: The standard single picker routing problem (SPRP) seeks the cost-minimal tour to collect a set of given articles in a rectangular single-block warehouse with parallel picking aisles and a dedicated storage policy, i.e, each SKU is only available from one storage location in the warehouse. We present a compact formulation that forgoes classical subtour elimination constraints by directly exploiting two of the properties of an optimal picking tour used in the dynamic programming algorithm of Ratliff and Rosenthal (1983). We extend the formulation to three important settings prevalent in modern e-commerce warehouses: scattered storage, decoupling of picker and cart, and multiple end depots. In numerical studies, our formulation outperforms existing standard SPRP formulations from the literature and proves able to solve large instances within short runtimes. Realistically sized instances of the three problem extensions can also be solved with low computational effort. We find that decoupling of picker and cart can lead to substantial cost savings depending on the speed and capacity of the picker when traveling alone, whereas additional end depots have rather limited benefits in a single-block warehouse.


title: Parallel Machine Scheduling to Minimize Energy Consumption

abstract: Given n jobs with release dates, deadlines and processing times we consider the problem of scheduling them on m parallel machines so as to minimize the total energy consumed. Machines can enter a sleep state and they consume no energy in this state. Each machine requires Q units of energy to awaken from the sleep state and in its active state the machine can process jobs and consumes a unit of energy per unit time. We allow for preemption and migration of jobs and provide the first constant approximation algorithm for this problem.


title: Initial Assessment of Monocrystalline Silicon Solar Cells as Large-Area Sensors for Precise Flux Calibration

abstract: As the precision frontier of large-area survey astrophysics advances towards the one millimagnitude level, flux calibration of astronomical instrumentation remains an ongoing challenge. We describe initial testing of silicon solar cells as large-aperture precise calibration photodiodes. We present measurements of dark current, linearity, frequency response, spatial response uniformity, and noise characteristics of the Sunpower C60 solar cells, an interdigitated back-contact 125mm x 125mm monocrystalline solar cell. We find that these devices hold considerable promise as large-area flux calibration sensors and warrant further characterization.


title: Quantum Information Approaches to Quantum Gravity

abstract: In this thesis we apply techniques from quantum information theory to study quantum gravity within the framework of the anti-de Sitter / conformal field theory correspondence (AdS/CFT). Through AdS/CFT, progress has been made in understanding the structure of entanglement in quantum field theories, and in how gravitational physics can emerge from these structures. However, this understanding is far from complete and will require the development of new tools to quantify correlations in CFT. This thesis presents refinements of a duality between operator product expansion (OPE) blocks in the CFT, giving the contribution of a conformal family to the OPE, and geodesic integrated fields in AdS which are diffeomorphism invariant quantities. This duality was originally discovered in the maximally symmetric setting of pure AdS dual to the CFT ground state. In less symmetric states the duality must be modified. Working with excited states within AdS$_3$/CFT$_2$, this thesis shows how the OPE block decomposes into more fine-grained CFT observables that are dual to AdS fields integrated over non-minimal geodesics. Additionally, this thesis contains results on the dynamics of entanglement measures for general quantum systems. Results are presented for the family of quantum R\'enyi entropies and entanglement negativity. R\'enyi entropies are studied for general dynamics by imposing special initial conditions. Around pure, separable initial states, all R\'enyi entropies grow with the same timescale at leading, and next-to-leading order. Mathematical tools are developed for the differentiation of non-analytic matrix functions with respect to constrained arguments and are used to construct analytic expressions for derivatives of negativity. We establish bounds on the rate of change of state distinguishability and the rate of entanglement growth for closed systems. Note: Abstract shortened.


title: Wilf collapse in permutation classes

abstract: For a hereditary permutation class $\mathcal{C}$, we say that two permutations $\pi$ and $\sigma$ of $\mathcal{C}$ are Wilf-equivalent in $\mathcal{C}$, if $\mathcal{C}$ has the same number of permutations avoiding $\pi$ as those avoiding $\sigma$. We say that a permutation class $\mathcal{C}$ exhibits a Wilf collapse if the number of permutations of size $n$ in $\mathcal{C}$ is asymptotically larger than the number of Wilf-equivalence classes formed by these permutations. In this paper, we show that Wilf collapse is a surprisingly common phenomenon. Among other results, we show that Wilf collapse occurs in any permutation class with unbounded growth and finitely many sum-indecomposable permutations. Our proofs are based on encoding the elements of a permutation class $\mathcal{C}$ as words, and analyzing the structure of a random permutation in $\mathcal{C}$ using this representation.


title: On the Lagrangian Trajectories for the One-Dimensional Euler Alignment Model without Vacuum Velocity

abstract: A well-known result of Carrillo, Choi, Tadmor, and Tan states that the 1D Euler Alignment model with smooth interaction kernels possesses a 'critical threshold' criterion for the global existence or finite-time blowup of solutions, depending on the global nonnegativity (or lack thereof) of the quantity $e_0 = \partial_x u_0 + \phi*\rho_0$. In this note, we rewrite the 1D Euler Alignment model as a first-order system for the particle trajectories in terms of a certain primitive $\psi_0$ of $e_0$; using the resulting structure, we give a complete characterization of global-in-time existence versus finite-time blowup of regular solutions that does not require a velocity to be defined in the vacuum. We also prove certain upper and lower bounds on the separation of particle trajectories, valid for smooth and weakly singular kernels, and we use them to weaken the hypotheses of Tan sufficient for the global-in-time existence of a solution in the weakly singular case, when the order of the singularity lies in the range $s\in (0,\frac12)$.


title: An existence result for singular fractional Kirchhoff-Schr\"{o}dinger-Poisson system

abstract: In this paper, we study the existence of infinitely many weak solutions to a fractional Kirchhoff-Schr\"{o}dinger-Poisson system involving the weak singularity, i.e. when $0<\gamma<1$. Further, we obtain the existence of a solution with the strong singularity, i.e. when $\gamma>1$. We employ variational techniques to prove the existence and multiplicity results. Moreover, a $L^{\infty}$ estimate is obtained by using the Moser iteration method.


title: Density Functional Modeling and Total Scattering Analysis of the Atomic Structure of a Quaternary CaO-MgO-Al2O3-SiO2 (CMAS) Glass: Uncovering the Local Environment of Magnesium

abstract: Quaternary CaO-MgO-Al2O3-SiO2 glasses are important constituents of the Earth's lower crust and mantle, and they also have important industrial applications such as in metallurgical processes, concrete production and emerging low-CO2 cement technologies. In particular, these applications rely heavily on the composition-structure-reactivity relationships for CMAS glasses, which are not yet well established. In this study, we developed a robust method that combines force-field molecular dynamics (MD) simulations and density functional theory (DFT) calculations with X-ray/neutron scattering experiments to resolve the atomic structure of a CMAS glass. The final structural representation generated using this method is not only thermodynamically favorable (according to DFT calculations) but also agrees with experiments (including X-ray/neutron scattering data as well as literature data). Detailed analysis of the final structure (including partial pair distribution functions, coordination number, oxygen environment) enabled existing discrepancies in the literature to be reconciled and has revealed new structural information on the CMAS glass, specifically, (i) the unambiguous assignment of medium-range atomic ordering, (ii) the preferential role of Ca atoms as charge compensators and Mg atoms as network modifiers, (iii) the proximity of Mg atoms to free oxygen sites, and (iv) clustering of Mg atoms. Overall, this new structural information will enhance our mechanistic understanding on CMAS glass dissolution behavior, including dissolution-related mechanisms occurring during the formation of low-CO2 cements.


title: Representation-Optimal Multi-Robot Motion Planning using Conflict-Based Search

abstract: Multi-Agent Motion Planning (MAMP) is the problem of computing feasible paths for a set of agents given individual start and goal states. Given the hardness of MAMP, most of the research related to multi-agent systems has focused on multi-agent pathfinding (MAPF), which simplifies the problem by assuming a shared discrete representation of the space for all agents. The Conflict-Based Search algorithm (CBS) has proven a tractable means of generating optimal solutions in discrete spaces. However, neither CBS nor other discrete MAPF techniques can be directly applied to solve MAMP problems because of the assumption of the shared discrete representation of the agents' state space. In this work, we solve MAMP problems by adapting the techniques discovered in the MAPF scenario by CBS to the more general problem with heterogeneous agents in a continuous space. We demonstrate the scalability teams of up to 32 agents, demonstrate the ability to handle up to 8 high DOF manipulators, and plan for heterogeneous teams. In all scenarios, our approach plans significantly faster while providing higher quality solutions.


title: Random Phase Approximation Applied to Many-Body Noncovalent Systems

abstract: The random phase approximation (RPA) has received a considerable interest in the field of modeling systems where noncovalent interactions are important. Its advantages over widely used density functional theory (DFT) approximations are the exact treatment of exchange and the description of long-range correlation. In this work we address two open questions related to RPA. First, how accurately RPA describes nonadditive interactions encountered in many-body expansion of a binding energy. We consider three-body nonadditive energies in molecular and atomic clusters. Second, how does the accuracy of RPA depend on input provided by different DFT models, without resorting to selfconsistent RPA procedure which is currently impractical for calculations employing periodic boundary conditions. We find that RPA based on the SCAN0 and PBE0 models, i.e., hybrid DFT, achieves an overall accuracy between CCSD and MP3 on a dataset of molecular trimers of \v{R}ez\'{a}\v{c} et al. (J. Chem. Theory. Comput. 2015, 11, 3065) Finally, many-body expansion for molecular clusters and solids often leads to a large number of small contributions that need to be calculated with a high precision. We therefore present a cubic-scaling (or SCF-like) implementation of RPA in atomic basis set, which is designed for calculations with a high numerical precision.


title: GACNN: Training Deep Convolutional Neural Networks with Genetic Algorithm

abstract: Convolutional Neural Networks (CNNs) have gained a significant attraction in the recent years due to their increasing real-world applications. Their performance is highly dependent to the network structure and the selected optimization method for tuning the network parameters. In this paper, we propose novel yet efficient methods for training convolutional neural networks. The most of current state of the art learning method for CNNs are based on Gradient decent. In contrary to the traditional CNN training methods, we propose to optimize the CNNs using methods based on Genetic Algorithms (GAs). These methods are carried out using three individual GA schemes, Steady-State, Generational, and Elitism. We present new genetic operators for crossover, mutation and also an innovative encoding paradigm of CNNs to chromosomes aiming to reduce the resulting chromosome's size by a large factor. We compare the effectiveness and scalability of our encoding with the traditional encoding. Furthermore, the performance of individual GA schemes used for training the networks were compared with each other in means of convergence rate and overall accuracy. Finally, our new encoding alongside the superior GA-based training scheme is compared to Backpropagation training with Adam optimization.


title: Siamese Neural Networks for Wireless Positioning and Channel Charting

abstract: Neural networks have been proposed recently for positioning and channel charting of user equipments (UEs) in wireless systems. Both of these approaches process channel state information (CSI) that is acquired at a multi-antenna base-station in order to learn a function that maps CSI to location information. CSI-based positioning using deep neural networks requires a dataset that contains both CSI and associated location information. Channel charting (CC) only requires CSI information to extract relative position information. Since CC builds on dimensionality reduction, it can be implemented using autoencoders. In this paper, we propose a unified architecture based on Siamese networks that can be used for supervised UE positioning and unsupervised channel charting. In addition, our framework enables semisupervised positioning, where only a small set of location information is available during training. We use simulations to demonstrate that Siamese networks achieve similar or better performance than existing positioning and CC approaches with a single, unified neural network architecture.


title: A fast boundary integral method for high-order multiscale mesh generation

abstract: In this work we present an algorithm to construct an infinitely differentiable smooth surface from an input consisting of a (rectilinear) triangulation of a surface of arbitrary shape. The original surface can have non-trivial genus and multiscale features, and our algorithm has computational complexity which is linear in the number of input triangles. We use a smoothing kernel to define a function $\Phi$ whose level set defines the surface of interest. Charts are subsequently generated as maps from the original user-specified triangles to $\mathbb R^3$. The degree of smoothness is controlled locally by the kernel to be commensurate with the fineness of the input triangulation. The expression for~$\Phi$ can be transformed into a boundary integral, whose evaluation can be accelerated using a fast multipole method. We demonstrate the effectiveness and cost of the algorithm with polyhedral and quadratic skeleton surfaces obtained from CAD and meshing software.


title: On inverse spectral problems for Sturm-Liouville differential operators on closed sets

abstract: We study Sturm-Liouville operators on closed sets of a special structure, which are sometimes referred as time scales and often appear in modelling various real processes. Depending on the set structure, such operators unify both differential and difference operators. We obtain properties of their spectral characteristics and study inverse problems with respect to them. We prove that the spectral characteristics uniquely determine the operator.


title: Enhancing Electrical Properties and Seebeck Effects of WO$_3$ Thin Films Using Spray Pyrolysis: Insights into the Conductivity and Carrier Type

abstract: This study focuses on enhancing the electrical and thermoelectric properties of tungsten trioxide (WO$_3$) thin films using the spray pyrolysis technique. Three samples, namely A$_1$, A$_2$, and A$_3$, were prepared by depositing the films onto glass substrates at different deposition rates: R$_1=1$ mL/min, R$_2=3$ mL/min, and R$_3=7$ mL/min, respectively. Subsequently, the films underwent annealing at 500 $^\circ$C for 2 hours in air. A detailed investigation was conducted to analyze the structural, morphological, optical, and electrical properties of the $\mathrm{WO}_3$ thin films. The FE-SEM images revealed the formation of circular strings with varying diameters for all layers. The diameter of each string decreased as the deposition rate increased. XRD structural analysis indicated that, before annealing, sample R$_1$ exhibited a polycrystalline nature with mixed tetragonal-hexagonal phases, while the other two samples were amorphous. After annealing, sample A$_3$ also became polycrystalline. Moreover, the UV-Vis spectra analysis demonstrated a monotonic decrease in the band gap of the layers with increasing spray rate, potentially attributed to the attenuation of the quantum confinement effect. The Seebeck effect was used to authenticate the $n$-type conductivity of the WO$_3$ layers, as the Seebeck effect and conductivity exhibited a direct connection. Furthermore, the current-voltage variation was found to align with the oxygen vacancies sites. Remarkably, the sample A$_3$ with the highest deposition rate exhibited the lowest resistivity.


title: End-to-End Deep Convolutional Active Contours for Image Segmentation

abstract: The Active Contour Model (ACM) is a standard image analysis technique whose numerous variants have attracted an enormous amount of research attention across multiple fields. Incorrectly, however, the ACM's differential-equation-based formulation and prototypical dependence on user initialization have been regarded as being largely incompatible with the recently popular deep learning approaches to image segmentation. This paper introduces the first tight unification of these two paradigms. In particular, we devise Deep Convolutional Active Contours (DCAC), a truly end-to-end trainable image segmentation framework comprising a Convolutional Neural Network (CNN) and an ACM with learnable parameters. The ACM's Eulerian energy functional includes per-pixel parameter maps predicted by the backbone CNN, which also initializes the ACM. Importantly, both the CNN and ACM components are fully implemented in TensorFlow, and the entire DCAC architecture is end-to-end automatically differentiable and backpropagation trainable without user intervention. As a challenging test case, we tackle the problem of building instance segmentation in aerial images and evaluate DCAC on two publicly available datasets, Vaihingen and Bing Huts. Our reseults demonstrate that, for building segmentation, the DCAC establishes a new state-of-the-art performance by a wide margin.


title: Library network, a possible path to explainable neural networks

abstract: Deep neural networks (DNNs) may outperform human brains in complex tasks, but the lack of transparency in their decision-making processes makes us question whether we could fully trust DNNs with high stakes problems. As DNNs' operations rely on a massive number of both parallel and sequential linear/nonlinear computations, predicting their mistakes is nearly impossible. Also, a line of studies suggests that DNNs can be easily deceived by adversarial attacks, indicating that their decisions can easily be corrupted by unexpected factors. Such vulnerability must be overcome if we intend to take advantage of DNNs' efficiency in high stakes problems. Here, we propose an algorithm that can help us better understand DNNs' decision-making processes. Our empirical evaluations suggest that this algorithm can effectively trace DNNs' decision processes from one layer to another and detect adversarial attacks.


title: A Longitudinal Framework for Predicting Nonresponse in Panel Surveys

abstract: Nonresponse in panel studies can lead to a substantial loss in data quality due to its potential to introduce bias and distort survey estimates. Recent work investigates the usage of machine learning to predict nonresponse in advance, such that predicted nonresponse propensities can be used to inform the data collection process. However, predicting nonresponse in panel studies requires accounting for the longitudinal data structure in terms of model building, tuning, and evaluation. This study proposes a longitudinal framework for predicting nonresponse with machine learning and multiple panel waves and illustrates its application. With respect to model building, this approach utilizes information from multiple waves by introducing features that aggregate previous (non)response patterns. Concerning model tuning and evaluation, temporal cross-validation is employed by iterating through pairs of panel waves such that the training and test sets move in time. Implementing this approach with data from a German probability-based mixed-mode panel shows that aggregating information over multiple panel waves can be used to build prediction models with competitive and robust performance over all test waves.


title: Language-Agnostic Syllabification with Neural Sequence Labeling

abstract: The identification of syllables within phonetic sequences is known as syllabification. This task is thought to play an important role in natural language understanding, speech production, and the development of speech recognition systems. The concept of the syllable is cross-linguistic, though formal definitions are rarely agreed upon, even within a language. In response, data-driven syllabification methods have been developed to learn from syllabified examples. These methods often employ classical machine learning sequence labeling models. In recent years, recurrence-based neural networks have been shown to perform increasingly well for sequence labeling tasks such as named entity recognition (NER), part of speech (POS) tagging, and chunking. We present a novel approach to the syllabification problem which leverages modern neural network techniques. Our network is constructed with long short-term memory (LSTM) cells, a convolutional component, and a conditional random field (CRF) output layer. Existing syllabification approaches are rarely evaluated across multiple language families. To demonstrate cross-linguistic generalizability, we show that the network is competitive with state of the art systems in syllabifying English, Dutch, Italian, French, Manipuri, and Basque datasets.


title: An unbiased approach to low rank recovery

abstract: Low rank recovery problems have been a subject of intense study in recent years. While the rank function is useful for regularization it is difficult to optimize due to its non-convexity and discontinuity. The standard remedy for this is to exchange the rank function for the convex nuclear norm, which is known to favor low rank solutions under certain conditions. On the downside the nuclear norm exhibits a shrinking bias that can severely distort the solution in the presence of noise, which motivates the use of stronger non-convex alternatives. In this paper we study two such formulations. We characterize the critical points and give sufficient conditions for a low rank stationary point to be unique. Moreover, we derive conditions that ensure global optimality of the low rank stationary point and show that these hold under moderate noise levels.


title: Hunting white dwarfs at the center of planetary nebulae

abstract: We present Gemini-South observations of nine faint and extended planetary nebulae. Using direct images taken with the spectrograph GMOS, we built the $(u' - g')$ vs. $(g' - r')$ diagrams of the stars in the observed areas which allowed us, also considering their geometrical positions, to identify the probable central stars of the nebulae. Our stellar spectra of seven stars, also taken with GMOS, indicate that four (and probably two more) objects are white dwarfs of the DAO subtype. Moreover, the white dwarf status of the four stars is confirmed by the parameters $ T_{\mathrm{eff}}$ and $ \log g$ derived with the help of theoretical stellar spectra. Given this evidence, we propose that these hot stars are the central ionizing sources of the nebulae. With this work we hope to help improve the current scarce statistics on central white dwarfs in planetary nebulae.


title: Shaping Field Gradients for Nanolocalization

abstract: Deep sub-wavelength localization and displacement sensing of optical nanoantennas have emerged as extensively pursued objectives in nanometrology, where focused beams serve as high-precision optical rulers while the scattered light provides an optical readout. Here, we show that in these schemes using an optical excitation as a position gauge implies that the sensitivity to displacements of a nanoantenna depends on the spatial gradients of the excitation field. Specifically, we explore one of such novel localization schemes based on appearance of transversely spinning fields in strongly confined optical beams, resulting in far-field segmentation of left- and right-hand circular polarizations of the scattered light, an effect known as the giant spin-Hall effect of light. We construct vector beams with augmented transverse spin density gradient in the focal plane and experimentally confirm enhanced sensitivity of the far-field spin-segmentation to lateral displacements of an electric-dipole nanoantenna. We conclude that sculpturing of electromagnetic field gradients and intelligent design of scatterers pave the way towards future improvements in displacement sensitivity.


title: Not so Particular about Calibration: Smile Problem Resolved

abstract: We present a novel Monte Carlo based LSV calibration algorithm that applies to all stochastic volatility models, including the non-Markovian rough volatility family. Our framework overcomes the limitations of the particle method proposed by Guyon and Henry-Labord\`ere (2012) and theoretically guarantees a variance reduction without additional computational complexity. Specifically, we obtain a closed-form and exact calibration method that allows us to remove the dependency on both the kernel function and bandwidth parameter. This makes the algorithm more robust and less prone to errors or instabilities in a production environment. We test the efficiency of our algorithm on various hybrid (rough) local stochastic volatility models.


title: From Turbulence to Landscapes: Universality of Logarithmic Mean Profiles in Bounded Complex Systems

abstract: The logarithmic mean-velocity profile is a key experimental and theoretical result in wall-bounded turbulence. Similarly, here we show that the topographic surface emerging between parallel zero-elevation boundaries presents an intermediate region with a logarithmic mean-elevation profile. We use model simulations, which account for growth, erosion, and smoothing processes and give rise to complex topography with channel branching and fractal river networks, as well as data from a physical landscape-evolution experiment. Dimensional and self-similarity arguments are used to corroborate this finding. Our results suggest a universality of the logarithmic scaling in bounded complex systems out of equilibrium, of which landscape topography and turbulence are quintessential examples.


title: Threshold-Secure Coding with Shared Key

abstract: Cryptographic protocols are often implemented at upper layers of communication networks, while error-correcting codes are employed at the physical layer. In this paper, we consider utilizing readily-available physical layer functions, such as encoders and decoders, together with shared keys to provide a threshold-type security scheme. To this end, we first consider a scenario where the effect of the physical layer is omitted and all the channels between the involved parties are assumed to be noiseless. We introduce a model for threshold-secure coding, where the legitimate parties communicate using a shared key such that an eavesdropper does not get any information, in an information-theoretic sense, about the key as well as about any subset of the input symbols of size up to a certain threshold. Then, a framework is provided for constructing threshold-secure codes from linear block codes while characterizing the requirements to satisfy the reliability and security conditions. Moreover, we propose a threshold-secure coding scheme, based on Reed-Muller (RM) codes, that meets security and reliability conditions. It is shown that the encoder and the decoder of the scheme can be implemented efficiently with quasi-linear time complexity. In particular, a successive cancellation decoder is shown for the RM-based coding scheme. Then we extend the setup to the scenario where the channel between the legitimate parties is no longer noiseless. The reliability condition for noisy channels is then modified accordingly, and a method is described to construct codes attaining threshold security as well as desired reliability. Also, we propose a coding scheme based on RM codes for threshold security and robustness designed for binary erasure channels along with a unified successive cancellation decoder. The proposed threshold-secure coding schemes are flexible and can be adapted for different key lengths.


title: Information Transfer in Dynamical Systems and Optimal Placement of Actuators and Sensors for Control of Non-equilibrium Dynamics

abstract: In this paper we develop the concept of information transfer between the Borel-measurable sets for a dynamical system described by a measurable space and a non-singular transformation. The concept is based on how Shannon entropy is transferred between the measurable sets, as the dynamical system evolves. We show that the proposed definition of information transfer satisfies the usual notions of information transfer and causality, namely, zero transfer and transfer asymmetry. Furthermore, we show how the information transfer measure can be used to classify ergodicity and mixing. We also develop the computational methods for information transfer computation and apply the framework for optimal placements of actuators and sensors for control of non-equilibrium dynamics.


title: Rigid automorphisms of linking systems

abstract: A rigid automorphism of a linking system is an automorphism which restricts to the identity on the Sylow subgroup. A rigid inner automorphism is conjugation by an element in the center of the Sylow subgroup. At odd primes, it is known that each rigid automorphism of a centric linking system is inner. We prove that the group of rigid outer automorphisms of a linking system at the prime $2$ is elementary abelian, and that it splits over the subgroup of rigid inner automorphisms. In a second result, we show that if an automorphism of a finite group $G$ restricts to the identity on the centric linking system for $G$, then it is of $p'$-order modulo the group of inner automorphisms, provided $G$ has no nontrivial normal $p'$-subgroups. We present two applications of this last result, one to tame fusion systems.


title: Gradient Descent: The Ultimate Optimizer

abstract: Working with any gradient-based machine learning algorithm involves the tedious task of tuning the optimizer's hyperparameters, such as its step size. Recent work has shown how the step size can itself be optimized alongside the model parameters by manually deriving expressions for "hypergradients" ahead of time. We show how to automatically compute hypergradients with a simple and elegant modification to backpropagation. This allows us to easily apply the method to other optimizers and hyperparameters (e.g. momentum coefficients). We can even recursively apply the method to its own hyper-hyperparameters, and so on ad infinitum. As these towers of optimizers grow taller, they become less sensitive to the initial choice of hyperparameters. We present experiments validating this for MLPs, CNNs, and RNNs. Finally, we provide a simple PyTorch implementation of this algorithm (see people.csail.mit.edu/kach/gradient-descent-the-ultimate-optimizer).


title: The model companions of set theory

abstract: This is an introductory paper to a series of results linking generic absoluteness results for second and third order number theory to the model theoretic notion of model companionship. Specifically we develop here a general framework linking Woodin's generic absoluteness results for second order number theory and the theory of universally Baire sets to model companionship and show that (with the required care in details) a $\Pi_2$-property formalized in an appropriate language for second order number theory is forcible from some $T\supseteq\mathsf{ZFC}+$large cardinals if and only if it is consistent with the universal fragment of $T$ if and only if it is realized in the model companion of $T$. In particular we show that the first order theory of $H_{\omega_1}$ is the model companion of the first order theory of the universe of sets assuming the existence of class many Woodin cardinals, and working in a signature with predicates for $\Delta_0$-properties and for all universally Baire sets of reals. We will extend these results also to the theory of $H_{\aleph_2}$ in a follow up of this paper.


title: Generation of Optical Chirality Patterns with Plane Waves, Evanescent Waves and Surface Plasmon Waves

abstract: We systematically investigate the generation of optical chirality patterns by applying the superposition of two waves in three scenarios, namely plane waves in free space, evanescent waves of totally reflected light at dielectric interface and propagating surface plasmon waves on a metallic surface. In each scenario, the general analytical solution of the optical chirality pattern is derived for different polarization states and propagating directions of the two waves. The analytical solutions are verified by numerical simulations. Spatially structured optical chirality patterns can be generated in all scenarios if the incident polarization states and propagation directions are correctly chosen. Optical chirality enhancement can be obtained from the constructive interference of free-space circularly polarized light or enhanced evanescent waves of totally reflected light. Surface plasmon waves do not provide enhanced optical chirality unless the near-field intensity enhancement is sufficiently high. The structured optical chirality patterns may find applications in chirality sorting, chiral imaging and circular dichroism spectroscopy.

categories: physics.optics cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.other physics.chem-ph


title: Deep k-NN Defense against Clean-label Data Poisoning Attacks

abstract: Targeted clean-label data poisoning is a type of adversarial attack on machine learning systems in which an adversary injects a few correctly-labeled, minimally-perturbed samples into the training data, causing a model to misclassify a particular test sample during inference. Although defenses have been proposed for general poisoning attacks, no reliable defense for clean-label attacks has been demonstrated, despite the attacks' effectiveness and realistic applications. In this work, we propose a simple, yet highly-effective Deep k-NN defense against both feature collision and convex polytope clean-label attacks on the CIFAR-10 dataset. We demonstrate that our proposed strategy is able to detect over 99% of poisoned examples in both attacks and remove them without compromising model performance. Additionally, through ablation studies, we discover simple guidelines for selecting the value of k as well as for implementing the Deep k-NN defense on real-world datasets with class imbalance. Our proposed defense shows that current clean-label poisoning attack strategies can be annulled, and serves as a strong yet simple-to-implement baseline defense to test future clean-label poisoning attacks. Our code is available at https://github.com/neeharperi/DeepKNNDefense


title: A Simple and Effective Model for Answering Multi-span Questions

abstract: Models for reading comprehension (RC) commonly restrict their output space to the set of all single contiguous spans from the input, in order to alleviate the learning problem and avoid the need for a model that generates text explicitly. However, forcing an answer to be a single span can be restrictive, and some recent datasets also include multi-span questions, i.e., questions whose answer is a set of non-contiguous spans in the text. Naturally, models that return single spans cannot answer these questions. In this work, we propose a simple architecture for answering multi-span questions by casting the task as a sequence tagging problem, namely, predicting for each input token whether it should be part of the output or not. Our model substantially improves performance on span extraction questions from DROP and Quoref by 9.9 and 5.5 EM points respectively.


title: Towards Races in Linear Logic

abstract: Process calculi based in logic, such as $\pi$DILL and CP, provide a foundation for deadlock-free concurrent programming, but exclude non-determinism and races. HCP is a reformulation of CP which addresses a fundamental shortcoming: the fundamental operator for parallel composition from the $\pi$-calculus does not correspond to any rule of linear logic, and therefore not to any term construct in CP. We introduce non-deterministic HCP, which extends HCP with a novel account of non-determinism. Our approach draws on bounded linear logic to provide a strongly-typed account of standard process calculus expressions of non-determinism. We show that our extension is expressive enough to capture many uses of non-determinism in untyped calculi, such as non-deterministic choice, while preserving HCP's meta-theoretic properties, including deadlock freedom.


title: Lane Attention: Predicting Vehicles' Moving Trajectories by Learning Their Attention over Lanes

abstract: Accurately forecasting the future movements of surrounding vehicles is essential for safe and efficient operations of autonomous driving cars. This task is difficult because a vehicle's moving trajectory is greatly determined by its driver's intention, which is often hard to estimate. By leveraging attention mechanisms along with long short-term memory (LSTM) networks, this work learns the relation between a driver's intention and the vehicle's changing positions relative to road infrastructures, and uses it to guide the prediction. Different from other state-of-the-art solutions, our work treats the on-road lanes as non-Euclidean structures, unfolds the vehicle's moving history to form a spatio-temporal graph, and uses methods from Graph Neural Networks to solve the problem. Not only is our approach a pioneering attempt in using non-Euclidean methods to process static environmental features around a predicted object, our model also outperforms other state-of-the-art models in several metrics. The practicability and interpretability analysis of the model shows great potential for large-scale deployment in various autonomous driving systems in addition to our own.


title: Physical parameters of red supergiants in dwarf irregular galaxies in the Local Group

abstract: Increasing the statistics of evolved massive stars in the Local Group enables investigating their evolution at different metallicities. During the late stages of stellar evolution, the physics of some phenomena, such as episodic and systematic mass loss, are not well constrained. For example, the physical properties of red supergiants (RSGs) in different metallicity regimes remain poorly understood. Thus, we initiated a systematic study of RSGs in dwarf irregular galaxies (dIrrs) in the Local Group. The target selection is based on 3.6 $\mu$m and 4.5 $\mu$m photometry from archival Spitzer Space Telescope images of nearby galaxies. We selected 46 targets in the dIrrs IC 10, IC 1613, Sextans B, and the Wolf-Lundmark-Melotte (WLM) galaxy that we observed with the GTC-OSIRIS and VLT-FORS2 instruments. We used several photometric techniques together with a spectral energy distribution analysis to derive the luminosities and effective temperatures of known and newly discovered RSGs. We identified and spectroscopically confirmed 4 new RSGs, 5 previously known RSGs, and 5 massive asymptotic giant branch (AGB) stars. We added known objects from previous observations. In total, we present spectral classification and fundamental physical parameters of 25 late-type massive stars in the following dIrrs: Sextans A, Sextans B, IC 10, IC 1613, Pegasus, Phoenix, and WLM. This includes 17 RSGs and 8 AGB stars that have been identified here and previously. Based on our observational results and PARSEC evolutionary models, we draw the following conclusions: (i) a trend to higher minimum effective temperatures at lower metallicities and (ii) the maximum luminosity of RSGs appears to be constant at $log$($L/L$$_{\odot}$) $\approx$ $5.5$, independent of the metallicity of the host environment (up to $\mathrm{[Fe/H]}$ $\approx$ $-1$ dex).


title: The telescope conjecture at height 2 and the tmf resolution

abstract: Mahowald proved the height 1 telescope conjecture at the prime 2 as an application of his seminal work on bo-resolutions. In this paper we study the height 2 telescope conjecture at the prime 2 through the lens of tmf-resolutions. To this end we compute the structure of the tmf-resolution for a specifc type 2 complex Z. We find that, analogous to the height 1 case, the E1-page of the tmf-resolution possesses a decomposition into a v2-periodic summand, and an Eilenberg-MacLane summand which consists of bounded v2-torsion. However, unlike the height 1 case, the E2-page of the tmf-resolution exhibits unbounded v2-torsion. We compare this to the work of Mahowald-Ravenel-Shick, and discuss how the validity of the telescope conjecture is connected to the fate of this unbounded v2-torsion: either the unbounded v2-torsion kills itself off in the spectral sequence, and the telescope conjecture is true, or it persists to form v2-parabolas and the telescope conjecture is false. We also study how to use the tmf-resolution to effectively give low dimensional computations of the homotopy groups of Z. These computations allow us to prove a conjecture of the second author and Egger: the E(2)-local Adams-Novikov spectral sequence for Z collapses.


title: Scalable and modular material point method for large-scale simulations

abstract: In this paper, we describe a new scalable and modular material point method (MPM) code developed for solving large-scale problems in continuum mechanics. The MPM is a hybrid Eulerian-Lagrangian approach, which uses both moving material points and computational nodes on a background mesh. The MPM has been successfully applied to solve large-deformation problems such as landslides, failure of slopes, concrete flows, etc. Solving these large-deformation problems result in the material points actively moving through the mesh. Developing an efficient parallelisation scheme for the MPM code requires dynamic load-balancing techniques for both the material points and the background mesh. This paper describes the data structures and algorithms employed to improve the performance and portability of the MPM code. An object-oriented programming paradigm is adopted to modularise the MPM code. The Unified Modelling Language (UML) diagram of the MPM code structure is shown in Figure 1.


title: Explainable Clustering and Application to Wealth Management Compliance

abstract: Many applications from the financial industry successfully leverage clustering algorithms to reveal meaningful patterns among a vast amount of unstructured financial data. However, these algorithms suffer from a lack of interpretability that is required both at a business and regulatory level. In order to overcome this issue, we propose a novel two-steps method to explain clusters. A classifier is first trained to predict the clusters labels, then the Single Feature Introduction Test (SFIT) method is run on the model to identify the statistically significant features that characterise each cluster. We describe a real wealth management compliance use-case that highlights the necessity of such an interpretable clustering method. We illustrate the performance of our method using simulated data and through an experiment on financial ratios of U.S. companies.


title: A classification of finite simple amenable Z-stable C*-algebras, II, --C*-algebras with rational generalized tracial rank one

abstract: A classification theorem is obtained for a class of unital simple separable amenable Z-stable C*-algebras which exhausts all possible values of the Elliott invariant for unital stably finite simple separable amenable Z-stable C*-algebras. Moreover, it contains all unital simple separable amenable C*-algebras which satisfy the UCT and have finite rational tracial rank


title: Uniqueness of tangent cones to boundary points of two-dimensional almost-minimizing currents

abstract: We prove that tangent cones at singular boundary points of a two-dimensional current almost area minimizing are unique. Following the ideas exposed by White in [8], the result is achieved by combining a suitable epiperimetric inequality and an almost-monotonicity formula for the mass at boundary points.


title: Optimal Sketching for Kronecker Product Regression and Low Rank Approximation

abstract: We study the Kronecker product regression problem, in which the design matrix is a Kronecker product of two or more matrices. Given $A_i \in \mathbb{R}^{n_i \times d_i}$ for $i=1,2,\dots,q$ where $n_i \gg d_i$ for each $i$, and $b \in \mathbb{R}^{n_1 n_2 \cdots n_q}$, let $\mathcal{A} = A_1 \otimes A_2 \otimes \cdots \otimes A_q$. Then for $p \in [1,2]$, the goal is to find $x \in \mathbb{R}^{d_1 \cdots d_q}$ that approximately minimizes $\|\mathcal{A}x - b\|_p$. Recently, Diao, Song, Sun, and Woodruff (AISTATS, 2018) gave an algorithm which is faster than forming the Kronecker product $\mathcal{A}$ Specifically, for $p=2$ their running time is $O(\sum_{i=1}^q \text{nnz}(A_i) + \text{nnz}(b))$, where nnz$(A_i)$ is the number of non-zero entries in $A_i$. Note that nnz$(b)$ can be as large as $n_1 \cdots n_q$. For $p=1,$ $q=2$ and $n_1 = n_2$, they achieve a worse bound of $O(n_1^{3/2} \text{poly}(d_1d_2) + \text{nnz}(b))$. In this work, we provide significantly faster algorithms. For $p=2$, our running time is $O(\sum_{i=1}^q \text{nnz}(A_i) )$, which has no dependence on nnz$(b)$. For $p<2$, our running time is $O(\sum_{i=1}^q \text{nnz}(A_i) + \text{nnz}(b))$, which matches the prior best running time for $p=2$. We also consider the related all-pairs regression problem, where given $A \in \mathbb{R}^{n \times d}, b \in \mathbb{R}^n$, we want to solve $\min_{x} \|\bar{A}x - \bar{b}\|_p$, where $\bar{A} \in \mathbb{R}^{n^2 \times d}, \bar{b} \in \mathbb{R}^{n^2}$ consist of all pairwise differences of the rows of $A,b$. We give an $O(\text{nnz}(A))$ time algorithm for $p \in[1,2]$, improving the $\Omega(n^2)$ time needed to form $\bar{A}$. Finally, we initiate the study of Kronecker product low rank and low $t$-rank approximation. For input $\mathcal{A}$ as above, we give $O(\sum_{i=1}^q \text{nnz}(A_i))$ time algorithms, which is much faster than computing $\mathcal{A}$.


title: Steady state programming of controlled nonlinear systems via deep dynamic mode decomposition

abstract: This paper describes the optimal selection of a control policy to program the steady state of controlled nonlinear systems with hyperbolic fixed points. This work is motivated by the field of synthetic biology, in which saddle points are common (along with limit cycles), and the aim is to program cells to perform both digital and analog computation, though developing genetic digital computation has been the main focus. We frame the analog computing challenge of generating a steady state input-output function inside living cells. To program the steady state, a data-driven approach is taken wherein an approximation of the Koopman operator, identified via deep dynamic mode decomposition, is used to describe the dynamics of the system linearly. The new representation of the dynamics are then used to solve an optimization problem for the input which maximizes a direction in state space. Some added structure on the Koopman operator learning process for controlled systems is given for dynamics that are separable in the state and input. Finally, the methods are demonstrated on simulation examples of an incoherent feedforward loop and a combinatorial promoter system, two common network architectures seen in the field of synthetic biology.


title: Liouville-Riemann-Roch theorems on abelian coverings

abstract: The classical Riemann-Roch theorem has been extended by N. Nadirashvili and then M. Gromov and M. Shubin to computing indices of elliptic operators on compact (as well as non-compact) manifolds, when a divisor mandates a finite number of zeros and allows a finite number of poles of solutions. On the other hand, Liouville type theorems count the number of solutions that are allowed to have a "pole at infinity." Usually these theorems do not provide the exact dimensions of the spaces of such solutions (only finite-dimensionality, possibly with estimates or asymptotics of the dimension. An important case has been discovered by M. Avellaneda and F. H. Lin and advanced further by J. Moser and M. Struwe. It pertains periodic elliptic operators of divergent type, where, surprisingly, exact dimensions can be computed. This study has been extended by P. Li and Wang and brought to its natural limit for the case of periodic elliptic operators on co-compact abelian coverings by P. Kuchment and Pinchover. Comparison of the results and techniques of Nadirashvili and Gromov and Shubin with those of Kuchment and Pinchover shows significant similarities, as well as appearance of the same combinatorial expressions in the answers. Thus a natural idea was considered that possibly the results could be combined somehow in the case of co-compact abelian coverings, if the infinity is "added to the divisor." This work shows that such results indeed can be obtained, although they come out more intricate than a simple-minded expectation would suggest. Namely, the interaction between the finite divisor and the point at infinity is non-trivial.


title: FaSNet: Low-latency Adaptive Beamforming for Multi-microphone Audio Processing

abstract: Beamforming has been extensively investigated for multi-channel audio processing tasks. Recently, learning-based beamforming methods, sometimes called \textit{neural beamformers}, have achieved significant improvements in both signal quality (e.g. signal-to-noise ratio (SNR)) and speech recognition (e.g. word error rate (WER)). Such systems are generally non-causal and require a large context for robust estimation of inter-channel features, which is impractical in applications requiring low-latency responses. In this paper, we propose filter-and-sum network (FaSNet), a time-domain, filter-based beamforming approach suitable for low-latency scenarios. FaSNet has a two-stage system design that first learns frame-level time-domain adaptive beamforming filters for a selected reference channel, and then calculate the filters for all remaining channels. The filtered outputs at all channels are summed to generate the final output. Experiments show that despite its small model size, FaSNet is able to outperform several traditional oracle beamformers with respect to scale-invariant signal-to-noise ratio (SI-SNR) in reverberant speech enhancement and separation tasks. Moreover, when trained with a frequency-domain objective function on the CHiME-3 dataset, FaSNet achieves 14.3\% relative word error rate reduction (RWERR) compared with the baseline model. These results show the efficacy of FaSNet particularly in reverberant and noisy signal conditions.


title: Separation probabilities and analogues of a Zagier-Stanley formula

abstract: In this paper, we first obtain some analogues of a formula of Zagier (1995) and Stanley (2011). For instance, we prove that the number of pairs of $n$-cycles whose product has $k$ cycles and has $m$ given elements contained in distinct cycles (or separated) is given by $$ \frac{2 (n-1)! C_m(n+1,k)}{(n+m)(n+1-m)} $$ when $n-k$ is even, where $C_m(n,k)$ is the number of permutations of $n$ elements having $k$ cycles and separating $m$ given elements. As consequences, we obtain the formulas for certain separation probabilities due to Du and Stanley, answering a call of Stanley for simple combinatorial proofs. Furthermore, we obtain the expectation and variance of the number of fixed points in the product of two random $n$-cycles.


title: Modeling irregular boundaries using isoparametric elements in the Material Point Method

abstract: The Material Point Method (MPM) is a hybrid Eulerian-Lagrangian approach capable of simulating large deformation problems of history-dependent materials. While the MPM can represent complex and evolving material domains by using Lagrangian points, boundary conditions are often applied to the Eulerian nodes of the background mesh nodes. Hence, the use of a structured mesh may become prohibitively restrictive for modeling complex boundaries such as a landslide topography. We study the suitability of unstructured background mesh with isoparametric elements to model irregular boundaries in the MPM. An inverse mapping algorithm is used to transform the material points from the global coordinates to the local natural coordinates. Dirichlet velocity and frictional boundary conditions are applied in the local coordinate system at each boundary node. This approach of modeling complex boundary conditions is validated by modeling the dynamics of a gravity-driven rigid block sliding on an inclined plane. This method is later applied to a flume test of controlled debris flow on an inclined plane conducted by the United States Geological Survey (USGS).


title: Cosmology as a CFT$_1$

abstract: We show that the simplest FLRW cosmological system consisting in the homogeneous and isotropic massless Einstein-Scalar system enjoys a hidden conformal symmetry under the 1D conformal group $SL(2,\mathbb{R})$ acting as Mobius transformations in proper time. This invariance is made explicit through the mapping of FLRW cosmology onto conformal mechanics. On the one hand, we identify the corresponding conformal Noether charges, as combinations of the Hamiltonian scalar constraint, the extrinsic curvature and the 3D volume, which form a closed $sl(2,\mathbb{R})$ Lie algebra. On the other hand, this approach allows to write FLRW cosmology in terms of a AdS${}_{2}$ phase space and a Schwarzian action. Preserving this conformal structure at the quantum level fixes the ordering ambiguities in the Wheeler-de Witt quantization and allows to formulate FLRW quantum cosmology as a CFT${}_{1}$. We show that the CFT two-points correlator is realized as the overlap of the evolution in proper time of cosmological coherent wave-packets. In particular, the two-points function is built from a vacuum state which, although not conformally invariant, coincides with the cosmological vacuum annihilated by the scalar constraint. These results suggest new perspectives in classical and quantum cosmology, among which the possibility to apply the conformal bootstrap program to quantize cosmological backgrounds.


title: Distributed SGD Generalizes Well Under Asynchrony

abstract: The performance of fully synchronized distributed systems has faced a bottleneck due to the big data trend, under which asynchronous distributed systems are becoming a major popularity due to their powerful scalability. In this paper, we study the generalization performance of stochastic gradient descent (SGD) on a distributed asynchronous system. The system consists of multiple worker machines that compute stochastic gradients which are further sent to and aggregated on a common parameter server to update the variables, and the communication in the system suffers from possible delays. Under the algorithm stability framework, we prove that distributed asynchronous SGD generalizes well given enough data samples in the training optimization. In particular, our results suggest to reduce the learning rate as we allow more asynchrony in the distributed system. Such adaptive learning rate strategy improves the stability of the distributed algorithm and reduces the corresponding generalization error. Then, we confirm our theoretical findings via numerical experiments.


title: Learning from Observations Using a Single Video Demonstration and Human Feedback

abstract: In this paper, we present a method for learning from video demonstrations by using human feedback to construct a mapping between the standard representation of the agent and the visual representation of the demonstration. In this way, we leverage the advantages of both these representations, i.e., we learn the policy using standard state representations, but are able to specify the expected behavior using video demonstration. We train an autonomous agent using a single video demonstration and use human feedback (using numerical similarity rating) to map the standard representation to the visual representation with a neural network. We show the effectiveness of our method by teaching a hopper agent in the MuJoCo to perform a backflip using a single video demonstration generated in MuJoCo as well as from a real-world YouTube video of a person performing a backflip. Additionally, we show that our method can transfer to new tasks, such as hopping, with very little human feedback.


title: Stationary scalar hairy configurations supported by Neumann compact stars

abstract: We study stationary scalar field hairy configurations supported by asymptotically flat horizonless compact stars. At the star surface, we impose Neumann boundary conditions for the scalar field. With analytical methods, we obtain bounds on the frequency of scalar fields. For certain discrete frequency satisfying the bounds, we numerically get solutions of scalar hairy stars. We also disclose effects of model parameters on the discrete frequency of scalar fields.


title: Extended Hadamard expansions for the Airy functions

abstract: A new series expansion for the the Airy function is presented here that stems from the method of steepest descents and can be related to the Hadamard expansions as presented in prevous works cited in the manuscript, and which is convergent for all values of the complex variable. Hadamard expansions were introduced as an extension of the method of steepest descents and are defined in terms of a large number of non-systematic integration path subdivisions. Unlike them, the expansions in the present work originate in the splitting of the steepest descent in a number of segments that is not only finite but very small, and which are defined on the basis of the location of the branch points. One of the segments reaches to infinity and this gives rise to the presence of upper incomplete Gamma functions. This is one of the most important differences with the Hadamard series as defined in the aforementioned references, where all the incomplete Gamma functions are of the lower type. The theoretical interest of the new series expansion is twofold. First of all, it shows how to convert an asymptotic series into a convergent one with a finite splitting of the steepest descent path. Secondly, the inverse of the phase function that is part of the Laplace-type equation is Taylor-expanded around branch points to produce Puiseux series when necessary. In addition to this, the proposed analysis shows again how the Stokes phenomenon for the Airy function is related to the transition of the steepest descent paths at $\arg z = \pm 2 \pi/3$ from one to two. In regard to its computational application, these series expansions require a relatively small number of terms for each of them to reach a very high precision.


title: SteReFo: Efficient Image Refocusing with Stereo Vision

abstract: Whether to attract viewer attention to a particular object, give the impression of depth or simply reproduce human-like scene perception, shallow depth of field images are used extensively by professional and amateur photographers alike. To this end, high quality optical systems are used in DSLR cameras to focus on a specific depth plane while producing visually pleasing bokeh. We propose a physically motivated pipeline to mimic this effect from all-in-focus stereo images, typically retrieved by mobile cameras. It is capable to change the focal plane a posteriori at 76 FPS on KITTI images to enable real-time applications. As our portmanteau suggests, SteReFo interrelates stereo-based depth estimation and refocusing efficiently. In contrast to other approaches, our pipeline is simultaneously fully differentiable, physically motivated, and agnostic to scene content. It also enables computational video focus tracking for moving objects in addition to refocusing of static images. We evaluate our approach on the publicly available datasets SceneFlow, KITTI, CityScapes and quantify the quality of architectural changes.


title: REQ-YOLO: A Resource-Aware, Efficient Quantization Framework for Object Detection on FPGAs

abstract: Deep neural networks (DNNs), as the basis of object detection, will play a key role in the development of future autonomous systems with full autonomy. The autonomous systems have special requirements of real-time, energy-efficient implementations of DNNs on a power-constrained system. Two research thrusts are dedicated to performance and energy efficiency enhancement of the inference phase of DNNs. The first one is model compression techniques while the second is efficient hardware implementation. Recent works on extremely-low-bit CNNs such as the binary neural network (BNN) and XNOR-Net replace the traditional floating-point operations with binary bit operations which significantly reduces the memory bandwidth and storage requirement. However, it suffers from non-negligible accuracy loss and underutilized digital signal processing (DSP) blocks of FPGAs. To overcome these limitations, this paper proposes REQ-YOLO, a resource-aware, systematic weight quantization framework for object detection, considering both algorithm and hardware resource aspects in object detection. We adopt the block-circulant matrix method and propose a heterogeneous weight quantization using the Alternating Direction Method of Multipliers (ADMM), an effective optimization technique for general, non-convex optimization problems. To achieve real-time, highly-efficient implementations on FPGA, we present the detailed hardware implementation of block circulant matrices on CONV layers and develop an efficient processing element (PE) structure supporting the heterogeneous weight quantization, CONV dataflow and pipelining techniques, design optimization, and a template-based automatic synthesis framework to optimally exploit hardware resource. Experimental results show that our proposed REQ-YOLO framework can significantly compress the YOLO model while introducing very small accuracy degradation.


title: Ion exchange gels enhance organic electrochemical transistor performance in aqueous solution

abstract: Conjugated polymer-based organic electrochemical transistors (OECTs) are being studied for applications ranging from biochemical sensing to neural interfaces. While new conjugated polymers are being developed that can interface digital electronics with the aqueous chemistry of life, the vast majority of high-performance, high-mobility organic transistor materials developed over the past decades are extremely poor at taking up biologically-relevant ions. Here we incorporate an ion exchange gel into an OECT, demonstrating that this structure is capable of taking up biologically-relevant ions from solution and injecting larger, more hydrophobic ions into the underlying polymer semiconductor active layer in multiple hydrophobic conjugated polymers. Using poly[2,5-bis(3-tetradecylthiophen-2-yl) thieno[3,2-b]thiophene] (PBTTT) as a model semiconductor active layer and a blend of the ionic liquid 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (BMIM TFSI) and poly(vinylidene fluoride-co-hexafluoropropylene) (PVDF-HFP) as the ion exchange gel, we demonstrate more than a four order of magnitude improvement in OECT device transconductance and a one hundred-fold increase in ion injection kinetics. We demonstrate the ability of the ion exchange gel OECT to record biological signals by measuring the action potentials of a Venus flytrap plant. These results show the possibility of using interface engineering to open up a wider palette of organic semiconductor materials as OECTs that can be gated by aqueous solutions.


title: Fractional-Order Model Predictive Control for Neurophysiological Cyber-Physical Systems: A Case Study using Transcranial Magnetic Stimulation

abstract: Fractional-order dynamical systems are used to describe processes that exhibit temporal long-term memory and power-law dependence of trajectories. There has been evidence that complex neurophysiological signals like electroencephalogram (EEG) can be modeled by fractional-order systems. In this work, we propose a model-based approach for closed-loop Transcranial Magnetic Stimulation (TMS) to regulate brain activity through EEG data. More precisely, we propose a model predictive control (MPC) approach with an underlying fractional-order system (FOS) predictive model. Furthermore, MPC offers, by design, an additional layer of robustness to compensate for system-model mismatch, which the more traditional strategies lack. To establish the potential of our framework, we focus on epileptic seizure mitigation by computational simulation of our proposed strategy upon seizure-like events. We conclude by empirically analyzing the effectiveness of our method, and compare it with event-triggered open-loop strategies.


title: A new bound for smooth spline spaces

abstract: For a planar simplicial complex Delta contained in R^2, Schumaker proved that a lower bound on the dimension of the space C^r_k(Delta) of planar splines of smoothness r and polynomial degree at most k on Delta is given by a polynomial P_Delta(r,k), and Alfeld-Schumaker showed this polynomial gives the correct dimension when k >= 4r+1. Examples due to Morgan-Scott, Tohaneanu, and Yuan show that the equality dim C^r_k(Delta) = P_Delta(r,k) can fail when k = 2r or 2r+1. We prove that the equality dim C^r_k(Delta)= P_Delta(r,k) cannot hold in general for k <= (22r+7)/10.


title: Nested Distributed Gradient Methods with Stochastic Computation Errors

abstract: In this work, we consider the problem of a network of agents collectively minimizing a sum of convex functions. The agents in our setting can only access their local objective functions and exchange information with their immediate neighbors. Motivated by applications where computation is imperfect, including, but not limited to, empirical risk minimization (ERM) and online learning, we assume that only noisy estimates of the local gradients are available. To tackle this problem, we adapt a class of Nested Distributed Gradient methods (NEAR-DGD) to the stochastic gradient setting. These methods have minimal storage requirements, are communication aware and perform well in settings where gradient computation is costly, while communication is relatively inexpensive. We investigate the convergence properties of our method under standard assumptions for stochastic gradients, i.e. unbiasedness and bounded variance. Our analysis indicates that our method converges to a neighborhood of the optimal solution with a linear rate for local strongly convex functions and appropriate constant steplengths. We also show that distributed optimization with stochastic gradients achieves a noise reduction effect similar to mini-batching, which scales favorably with network size. Finally, we present numerical results to demonstrate the effectiveness of our method.


title: Spin-dependent refraction at the interface of lateral heterostructures of 2$H$-type transition-metal dichalcogenide monolayers

abstract: We study the refraction effect of electronic wave in hole-doped lateral heterojunctions of metallic and semiconducting transition-metal dichalcogenide monolayers. This effect is theoretically investigated in 2$H$-type MoSe$_2$-NbS$_2$ and WSe$_2$-NbS$_2$ junctions by combining the first-principles calculation and the lattice Green's function method. We show that the electronic waves change the direction of motion at the interface and collimate the velocity along two different directions depending on the spin. We find that the transmission probability increases with the charge density and that the direction of refracted electron beams is close to $\pm30^\circ$ with respect to the perpendicular axis to the interface. The metallic transition-metal dichalcogenide is essential for the refraction effect because of the strong trigonal-warping effect, the large Fermi surface, and the Zeeman-type spin-orbit coupling. The refraction effect enables to generate the spin-polarized electronic current by using a simple fabrication of transition-metal dichalcogenide monolayers.


title: On generalization of classical Hurwitz stability criteria for matrix polynomials

abstract: In this paper, we associate a class of Hurwitz matrix polynomials with Stieltjes positive definite matrix sequences. This connection leads to an extension of two classical criteria of Hurwitz stability for real polynomials to matrix polynomials: tests for Hurwitz stability via positive definiteness of block-Hankel matrices built from matricial Markov parameters and via matricial Stieltjes continued fractions. We obtain further conditions for Hurwitz stability in terms of block-Hankel minors and quasiminors, which may be viewed as a weak version of the total positivity criterion.


title: Using GANs for Sharing Networked Time Series Data: Challenges, Initial Promise, and Open Questions

abstract: Limited data access is a longstanding barrier to data-driven research and development in the networked systems community. In this work, we explore if and how generative adversarial networks (GANs) can be used to incentivize data sharing by enabling a generic framework for sharing synthetic datasets with minimal expert knowledge. As a specific target, our focus in this paper is on time series datasets with metadata (e.g., packet loss rate measurements with corresponding ISPs). We identify key challenges of existing GAN approaches for such workloads with respect to fidelity (e.g., long-term dependencies, complex multidimensional relationships, mode collapse) and privacy (i.e., existing guarantees are poorly understood and can sacrifice fidelity). To improve fidelity, we design a custom workflow called DoppelGANger (DG) and demonstrate that across diverse real-world datasets (e.g., bandwidth measurements, cluster requests, web sessions) and use cases (e.g., structural characterization, predictive modeling, algorithm comparison), DG achieves up to 43% better fidelity than baseline models. Although we do not resolve the privacy problem in this work, we identify fundamental challenges with both classical notions of privacy and recent advances to improve the privacy properties of GANs, and suggest a potential roadmap for addressing these challenges. By shedding light on the promise and challenges, we hope our work can rekindle the conversation on workflows for data sharing.


title: Towards modular and programmable architecture search

abstract: Neural architecture search methods are able to find high performance deep learning architectures with minimal effort from an expert. However, current systems focus on specific use-cases (e.g. convolutional image classifiers and recurrent language models), making them unsuitable for general use-cases that an expert might wish to write. Hyperparameter optimization systems are general-purpose but lack the constructs needed for easy application to architecture search. In this work, we propose a formal language for encoding search spaces over general computational graphs. The language constructs allow us to write modular, composable, and reusable search space encodings and to reason about search space design. We use our language to encode search spaces from the architecture search literature. The language allows us to decouple the implementations of the search space and the search algorithm, allowing us to expose search spaces to search algorithms through a consistent interface. Our experiments show the ease with which we can experiment with different combinations of search spaces and search algorithms without having to implement each combination from scratch. We release an implementation of our language with this paper.


title: G-0.02-0.07, the Compact HII Region Complex nearest to the Galactic Center with ALMA

abstract: We have observed the compact HII region complex nearest to the dynamical center of the Galaxy, G-0.02-0.07, using ALMA in the H42a recombination line, CS J=2-1, H13CO+ J=1-0, and SiO v=0, J=2-1 emission lines, and 86 GHz continuum emission. The HII regions HII-A to HII-C in the cluster are clearly resolved into a shell-like feature with a bright-half and a dark-half in the recombination line and continuum emission. The absorption features in the molecular emission lines show that HII-A, B and C are located on the near side of the 50 km/s Molecular Cloud (50MC) but HII-D is located on the far side. The electron temperatures and densities range Te=5150-5920 K and ne=950-2340 cm-3, respectively. The electron temperatures on the bright-half are slightly lower than those on the dark-half, while the electron densities on the bright-half are slightly higher than those on the dark-half. The HII regions are located on the molecular filaments in the 50MC. They have already broken through the filaments and are growing in the surrounding molecular gas. There are some shocked molecular gas components around the HII regions. From line width of the H42a recombination line, the expansion velocities from HII-A to HII-D are estimated to be Vexp=16.7, 11.6, 11.1, and 12.1 km/s, respectively. The expansion timescales from HII-A to HII-D are estimated to be Tage~1.4x0^4, 1.7x10^4, 2.0x10^4, and 0.7x10^4 years, respectively. The spectral types of the central stars from HII-A to HII-D are estimated to be O8V, O9.5V, O9V, and B0V, respectively. The positional relation among the HII regions, the SiO molecule enhancement area, and Class-I maser spots suggest that the shock wave caused by a cloud-cloud collision propagated along the line from HII-C to HII-A in the 50MC. The shock wave would trigger the massive star formation.


title: Embedding Dimension Phenomena in Intersection Complete Codes

abstract: Two tantalizing invariants of a combinatorial code $\mathcal C\subseteq 2^{[n]}$ are cdim$(\mathcal C)$ and odim$(\mathcal C)$, the smallest dimension in which $\mathcal C$ can be realized by convex closed or open sets, respectively. Cruz, Giusti, Itskov, and Kronholm showed that for intersection complete codes $\mathcal C$ with $m+1$ maximal codewords, odim$(\mathcal C)$ and cdim$(\mathcal C)$ are both bounded above by $\max\{2,m\}$. Results of Lienkaemper, Shiu, and Woodstock imply that odim and cdim may differ, even for intersection complete codes. We add to the literature on open and closed embedding dimensions of intersection complete codes with the following results: (*) If $\mathcal C$ is a simplicial complex, then cdim$(\mathcal{C}) = \mbox{odim}(\mathcal C)$, (*) If $\mathcal C$ is intersection complete, then cdim$(\mathcal C)\le \mbox{odim}(\mathcal C)$, (*) If $\mathcal C\subseteq 2^{[n]}$ is intersection complete with $n\ge 2$, then cdim$(\mathcal C) \le \min \{2d+1, n-1\}$, where $d$ is the dimension of the simplicial complex of $\mathcal C$, and (*) For each simplicial complex $\Delta\subseteq 2^{[n]}$ with $m\ge 2$ facets, the code $\mathcal S_\Delta := (\Delta \ast (n+1)) \cup \{[n]\}\subseteq 2^{[n+1]}$ is intersection complete, has $m+1$ maximal codewords, and satisfies odim$(\mathcal S_\Delta)=m$. In particular, for each $n\ge 3$ there exists an intersection complete code $\mathcal C\subseteq 2^{[n]}$ with odim$(\mathcal C) = \binom{n-1}{\lfloor (n-1)/2\rfloor}$. A key tool in our work is the study of sunflowers: arrangements of convex open sets in which the sets simultaneously meet in a central region, and nowhere else. We use Tverberg's theorem to study the structure of "$k$-flexible" sunflowers, and consequently obtain new lower bounds on odim$(\mathcal C)$ for intersection complete codes $\mathcal C$.


title: A perturbative method for resolving contact interactions in quantum mechanics

abstract: Long-range effective methods are ubiquitous in physics and in quantum theory, in particular. Furthermore, the reliability of such methods is higher when the nature of short-ranged interactions need not be modeled explicitly. This may be necessary for two reasons: (1) there are interactions that occur over a short range that cannot be accurately modeled with a potential function and/or (2) the entire Hamiltonian loses its reliability when applied at short distances. This work is an investigation of the utility and consequences of omitting a finite region of space from quantum mechanical analysis, accomplished by imposition of an artificial boundary behind which obscured short-ranged physical effects may operate. With this method, a free function of integration that depends on momentum is interpreted as a function encoding information needed to match a long-distance wavefunction to an appropriate state function on the other side of the boundary. Omitting part of the space from analysis implies that the strict unitarity requirement of quantum mechanics must be relaxed, since particles can actually propagate beyond the boundary. Strict orthogonality of eigenmodes and hermiticity of the Hamiltonian must also be relaxed in this method; however, all of these canonical relations are obeyed when averaged over sufficiently long times. What is achieved, therefore, appears to be an effective long-wavelength theory, at least for stationary systems. As examples, the quantum defect theory of the one-dimensional Coulomb interaction is recovered, as well as a new perspective of the inverse-square potential and the free particle, as well as the Wigner time delay associated with contact interactions. Potential applications of this method may include three-dimensional atomic systems and two-dimensional systems, such as graphene.

categories: quant-ph cond-mat.mes-hall hep-th physics.atom-ph


title: Multi-classifier prediction of knee osteoarthritis progression from incomplete imbalanced longitudinal data

abstract: Conventional inclusion criteria used in osteoarthritis clinical trials are not very effective in selecting patients who would benefit from a therapy being tested. Typically majority of selected patients show no or limited disease progression during a trial period. As a consequence, the effect of the tested treatment cannot be observed, and the efforts and resources invested in running the trial are not rewarded. This could be avoided, if selection criteria were more predictive of the future disease progression. In this article, we formulated the patient selection problem as a multi-class classification task, with classes based on clinically relevant measures of progression (over a time scale typical for clinical trials). Using data from two long-term knee osteoarthritis studies OAI and CHECK, we tested multiple algorithms and learning process configurations (including multi-classifier approaches, cost-sensitive learning, and feature selection), to identify the best performing machine learning models. We examined the behaviour of the best models, with respect to prediction errors and the impact of used features, to confirm their clinical relevance. We found that the model-based selection outperforms the conventional inclusion criteria, reducing by 20-25% the number of patients who show no progression. This result might lead to more efficient clinical trials.


title: LCAO-TDDFT-$k$-$\omega$: Spectroscopy in the Optical Limit

abstract: Understanding, optimizing, and controlling the optical absorption process, exciton gemination, and electron-hole separation and conduction in low dimensional systems is a fundamental problem in materials science. However, robust and efficient methods capable of modelling the optical absorbance of low dimensional macromolecular systems and providing physical insight into the processes involved have remained elusive. We employ a highly efficient linear combination of atomic orbitals (LCAOs) representation of the Kohn--Sham (KS) orbitals within time dependent density functional theory (TDDFT) in the reciprocal space ($k$) and frequency ($\omega$) domains, as implemented within our LCAO-TDDFT-$k$-$\omega$ code, and apply the derivative discontinuity correction of the exchange functional $\Delta_x$ to the KS eigenenergies. In so doing we are able to provide a semi-quantitative description of the optical absorption, conductivity, and polarizability spectra for prototypical 0D, 1D, 2D, and 3D systems within the optical limit ($\|\bf{q}\|\to0^+$) as compared to both available measurements and from solving the Bethe$-$Salpeter equation with quasiparticle $G_0 W_0$ eigenvalues ($G_0 W_0$-BSE). Specifically, we consider 0D fullerene (C$_{60}$), 1D metallic (10,0) and semiconducting (10,10) single-walled carbon nanotubes (SWCNTs), 2D graphene (GR) and phosphorene (PN), and 3D rutile (R-TiO$_2$) and anatase (A-TiO$_2$). For each system, we also employ the spatially resolved electron-hole density to provide direct physical insight into the nature of their optical excitations. These results demonstrate the reliability, applicability, efficiency, and robustness of our LCAO-TDDFT-$k$-$\omega$ code, and open the pathway to the computational design of macromolecular systems for optoelectronic, photovoltaic, and photocatalytic applications $in$ $silico$.


title: A Method for Geodesic Distance on Subdivision of Trees with Arbitrary Orders and Their Applications

abstract: Geodesic distance, sometimes called shortest path length, has proven useful in a great variety of applications, such as information retrieval on networks including treelike networked models. Here, our goal is to analytically determine the exact solutions to geodesic distances on two different families of growth trees which are recursively created upon an arbitrary tree $\mathcal{T}$ using two types of well-known operations, first-order subdivision and ($1,m$)-star-fractal operation. Different from commonly-used methods, for instance, spectral techniques, for addressing such a problem on growth trees using a single edge as seed in the literature, we propose a novel method for deriving closed-form solutions on the presented trees completely. Meanwhile, our technique is more general and convenient to implement compared to those previous methods mainly because there are not complicated calculations needed. In addition, the closed-form expression of mean first-passage time ($MFPT$) for random walk on each member in tree families is also readily obtained according to connection of our obtained results to effective resistance of corresponding electric networks. The results suggest that the two topological operations above are sharply different from each other due to $MFPT$ for random walks, and, however, have likely to show the similar performance, at least, on geodesic distance.