title: UQ-CHI: An Uncertainty Quantification-Based Contemporaneous Health Index for Degenerative Disease Monitoring

abstract: Developing knowledge-driven contemporaneous health index (CHI) that can precisely reflect the underlying patient across the course of the condition's progression holds a unique value, like facilitating a range of clinical decision-making opportunities. This is particularly important for monitoring degenerative condition such as Alzheimer's disease (AD), where the condition of the patient will decay over time. Detecting early symptoms and progression sign, and continuous severity evaluation, are all essential for disease management. While a few methods have been developed in the literature, uncertainty quantification of those health index models has been largely neglected. To ensure the continuity of the care, we should be more explicit about the level of confidence in model outputs. Ideally, decision-makers should be provided with recommendations that are robust in the face of substantial uncertainty about future outcomes. In this paper, we aim at filling this gap by developing an uncertainty quantification based contemporaneous longitudinal index, named UQ-CHI, with a particular focus on continuous patient monitoring of degenerative conditions. Our method is to combine convex optimization and Bayesian learning using the maximum entropy learning (MEL) framework, integrating uncertainty on labels as well. Our methodology also provides closed-form solutions in some important decision making tasks, e.g., such as predicting the label of a new sample. Numerical studies demonstrate the effectiveness of the propose UQ-CHI method in prediction accuracy, monitoring efficacy, and unique advantages if uncertainty quantification is enabled practice.


title: Tubes of finite Chen-type

abstract: In this paper, we consider surfaces in the 3-dimensional Euclidean space E3 which are of finite III-type, that is, they are of finite type, in the sense of B.-Y. Chen, corresponding to the third fundamental form. We present an important family of surfaces, namely, tubes in E3 .We show that tubes are of infinite III-type.


title: Ruled surfaces of finite type with respect to the second fundamental form

abstract: In this article, we consider surfaces in the 3-dimensional Euclidean space E3 without parabolic points which are of finite II-type, that is, they are of finite type, in the sense of B.-Y. Chen, corresponding to the second fundamental form. We study an important family of surfaces, namely, ruled surfaces in E3. We show that ruled surfaces are of infinite II-type.


title: A new stability test for linear neutral differential equations

abstract: We obtain new explicit exponential stability conditions for the linear scalar neutral equation with two bounded delays $ \dot{x}(t)-a(t)\dot{x}(g(t))+b(t)x(h(t))=0, $ where $ 0\leq a(t)\leq A_0<1$, $0<b_0\leq b(t)\leq B$, using the Bohl-Perron theorem and a transformation of the neutral equation into a differential equation with an infinite number of delays. The results are applied to the neutral logistic equation.


title: Adapting free-space fast multipole method for layered media Green's function: algorithm and analysis

abstract: In this paper, we present a numerical algorithm for the accurate and efficient computation of the convolution of the frequency domain layered media Green's function with a given density function. Instead of compressing the convolution matrix directly as in the classical fast multipole method, fast direct solvers, and fast H-matrix algorithms, the new algorithm considers a translated form of the original matrix so that many existing building blocks from the highly optimized free-space fast multipole method can be easily adapted to the Sommerfeld integral representations of the layered media Green's function. An asymptotic analysis is performed on the Sommerfeld integrals for large orders to provide an estimate of the decay rate in the new "multipole" and "local" expansions. In order to avoid the highly oscillatory integrand in the original Sommerfeld integral representations when the source and target are close to each other, or when they are both close to the interface in the scattered field, mathematically equivalent alternative direction integral representations are introduced. The convergence of the multipole and local expansions and formulas and quadrature rules for the original and alternative direction integral representations are numerically validated.


title: WebProt\'eg\'e: A Cloud-Based Ontology Editor

abstract: We present WebProt\'eg\'e, a tool to develop ontologies represented in the Web Ontology Language (OWL). WebProt\'eg\'e is a cloud-based application that allows users to collaboratively edit OWL ontologies, and it is available for use at https://webprotege.stanford.edu. WebProt\'ege\'e currently hosts more than 68,000 OWL ontology projects and has over 50,000 user accounts. In this paper, we detail the main new features of the latest version of WebProt\'eg\'e.


title: On stability of linear neutral differential equations in the Hale form

abstract: We present new explicit exponential stability conditions for the linear scalar neutral equation with two variable coefficients and delays $$ (x(t)-a(t)x(g(t)))'=-b(t)x(h(t)), $$ where $|a(t)|<1$, $b(t)\geq 0$, $h(t)\leq t$, $g(t)\leq t$, in the case when the delays $t-h(t)$, $t-g(t)$ are bounded, as well as an asymptotic stability condition, if the delays can be unbounded.


title: Cross helicity sign reversals in the dissipative scales of magnetohydrodynamic turbulence

abstract: We perform direct numerical simulations of magnetohydrodynamic (MHD) turbulence with kinetic energy and cross helicity injections at large scales. We observe that the cross helicity changes sign as we go from large and intermediate scales to small scales. In addition, the magnetic reconnections are strongest at the regions where the cross helicity changes sign and becomes smallest in magnitude. Thus, our simulations provide an important window to explore the regions of magnetic reconnections in nonlinear MHD.


title: Geometry free materials enabled by transformation optics for enhancing the intensity of electromagnetic waves in an arbitrary domain

abstract: The concentration of electromagnetic waves (EM) is of utmost importance in many engineering applications such as solar-cells. According to the transformation optics (TO) methodology, a feasible approach for obtaining arbitrary shape concentrators is proposed. In contrary to the previous works, the obtained materials are homogeneous and independent of the concentrator shape. That is regardless of the input geometry always one constant material, which is optic-null medium (ONM), is used and the performance of device will not alter. This competency will ease the design process and circumvent the sophisticated calculations of the necessitating materials. To authenticate the concept, several numerical full-wave simulations were performed for different shape concentrators. It was observed that numerical results exhibit strong agreement with the theoretical investigations and corroborate the generality and effectiveness of the proposed designing method.


title: Heat and water vapor transfer in the wake of a falling ice sphere and its implication for secondary ice formation in clouds

abstract: We perform direct numerical simulations of the settling of an ice sphere in an ambient fluid accounting for heat and mass transfer with the aim of studying in a meteorological context the case of falling graupel in humid air. The study is motivated by the fact that falling graupels in clouds are heated by the latent heat released during the accretion of liquid water droplets. They may therefore be considerably warmer than their surrounding and evaporate water vapor, which mixes with the surrounding air in the wake of the graupel, thereby creating transient zones of supersaturation there. The problem of a falling graupel is modeled as that of a heated sphere falling in a quiescent ambient fluid under the action of gravity. The coupling between the temperature and velocity fields is accounted for by the Boussinesq approximation. This problem can be parameterized by four parameters: the particle/fluid density ratio, the Galileo number, the Prandtl number, and the Richardson number. A separate scalar transport equation accounts for the vapor transport. Typical cloud conditions involve small temperature differences between the sphere and the surrounding, yielding relatively small Richardson numbers for both heat and mass transport. We give a special emphasis to the Galileo numbers 150, 170, 200 and 300 in order to analyze the specificities of each settling regime. The questions addressed in this study are mainly methodological and concern the influence of the settling regime and the mobility of the sphere on the structure of the scalar fields, the possible influence of modest Richardson numbers on the structure of the wake, and the possible application of this simulation framework to the investigation of the saturation in the wake of a falling graupel.


title: Patterned Supersolids in Dipolar Bose Systems

abstract: We study by means of first principle Quantum Monte Carlo simulations the ground state phase diagram of a system of dipolar bosons with aligned dipole moments, and with the inclusion of a two-body repulsive potential of varying range. The system is shown to display a supersolid phase in a relatively broad region of the phase diagram, featuring different crystalline patterns depending on the density and on the range of the repulsive part of the interaction (scattering length). The supersolid phase is sandwiched between a classical crystal of parallel filaments and a homogeneous superfluid phase. We show that a "roton" minimum appears in the elementary excitation spectrum of the superfluid as the system approaches crystallization.


title: Assessing the quantumness of the annealing dynamics via Leggett Gargs inequalities: a weak measurement approach

abstract: Adiabatic quantum computation (AQC) is a promising counterpart of universal quantum computation, based on the key concept of quantum annealing (QA). QA is claimed to be at the basis of commercial quantum computers and benefits from the fact that the detrimental role of decoherence and dephasing seems to have poor impact on the annealing towards the ground state. While many papers show interesting optimization results with a sizable number of qubits, a clear evidence of a full quantum coherent behavior during the whole annealing procedure is still lacking. In this paper we show that quantum non-demolition (weak) measurements of Leggett Garg inequalities can be used to efficiently assess the quantumness of the QA procedure. Numerical simulations based on a weak coupling Lindblad approach are compared with classical Langevin simulations to support our statements.


title: A Cold Atom Radio-Frequency Magnetometer

abstract: We propose and demonstrate a radio-frequency atomic magnetometer with sub-Doppler laser cooled rubidium-87. With a simple and compact design, our system demonstrates a sensitivity of $330~pT/\sqrt{Hz}$ in an unshielded environment, thus matching or surpassing previously reported cold atoms designs. By merging the multiple uses and robustness of radio-frequency atomic magnetometers with the detailed control of laser cooling, our cold atom radio-frequency magnetometer has the potential to move applications of atomic magnetometry to high spatial resolution. Direct impact in metrology for applied sciences, materials characterization, and nanotechnology can be anticipated.


title: On the Number of Discrete Chains

abstract: We study a generalization of Erd\H os's unit distances problem to chains of $k$ distances. Given $\mathcal P,$ a set of $n$ points, and a sequence of distances $(\delta_1,\ldots,\delta_k)$, we study the maximum possible number of tuples of distinct points $(p_1,\ldots,p_{k+1})\in \mathcal P^{k+1}$ satisfying $|p_j p_{j+1}|=\delta_j$ for every $1\leq j \leq k$. We study the problem in $\mathbb R^2$ and in $\mathbb R^3$, and derive upper and lower bounds for this family of problems.


title: Computational prediction of lattice thermal conductivity -- a comparison of molecular dynamics and Boltzmann transport approaches

abstract: The predictive modeling of lattice thermal conductivity is of fundamental importance for the understanding and design of materials for a wide range of applications. Two major approaches, namely molecular dynamics (MD) simulations and calculations solving approximately the Boltzmann transport equation (BTE), have been developed to compute the lattice thermal conductivity. We present a detailed direct comparison of these two approaches, using as prototypical cases MgO and PbTe. The comparison, carried out using empirical potentials, takes into account the effects of fourth order phonon scattering, temperature-dependent phonon frequencies (phonon renormalization), and investigates the effects of quantum vs. classical statistics. We clarify that equipartition, as opposed to Maxwell Boltzmann, govern the statistics of phonons in MD simulations. We find that lattice thermal conductivity values from MD and BTE show an apparent, satisfactory agreement; however such an agreement is the result of error cancellations. We also show that the primary effect of statistics on thermal conductivity is via the scattering rate dependence on phonon populations.


title: Latent Translation: Crossing Modalities by Bridging Generative Models

abstract: End-to-end optimization has achieved state-of-the-art performance on many specific problems, but there is no straight-forward way to combine pretrained models for new problems. Here, we explore improving modularity by learning a post-hoc interface between two existing models to solve a new task. Specifically, we take inspiration from neural machine translation, and cast the challenging problem of cross-modal domain transfer as unsupervised translation between the latent spaces of pretrained deep generative models. By abstracting away the data representation, we demonstrate that it is possible to transfer across different modalities (e.g., image-to-audio) and even different types of generative models (e.g., VAE-to-GAN). We compare to state-of-the-art techniques and find that a straight-forward variational autoencoder is able to best bridge the two generative models through learning a shared latent space. We can further impose supervised alignment of attributes in both domains with a classifier in the shared latent space. Through qualitative and quantitative evaluations, we demonstrate that locality and semantic alignment are preserved through the transfer process, as indicated by high transfer accuracies and smooth interpolations within a class. Finally, we show this modular structure speeds up training of new interface models by several orders of magnitude by decoupling it from expensive retraining of base generative models.


title: Stochastic control stabilizing unstable or chaotic maps

abstract: The paper considers a stabilizing stochastic control which can be applied to a variety of unstable and even chaotic maps. Compared to previous methods introducing control by noise, we relax assumptions on the class of maps, as well as consider a wider range of parameters for the same maps. This approach allows to stabilize unstable and chaotic maps by noise. The interplay between the map properties and the allowed neighbourhood where a solution can start to be stabilized is explored: as instability of the original map increases, the interval of allowed initial conditions narrows. A directed stochastic control aiming at getting to the target neighbourhood almost sure is combined with a controlling noise. Simulations illustrate that for a variety of problems, an appropriate bounded noise can stabilize an unstable positive equilibrium, without a limitation on the initial value.


title: One-Dimensional Flow Spectra and Cumulative Energy from Two Pipe Facilities

abstract: Experiments have been conducted to assess the sizes and energy fractions of structure in fully developed turbulent pipe flow regime in two pipe facilities, ColaPipe at BTU Cottbus-Senftenberg, and CICLoPE at University of Bologna, for shear Reynolds number in the range $2.5\cdot{10^3}\le{\mathrm{Re_{\tau}}}\le{{3.7\cdot{10^4}}}$, utilizing a single hot-wire probe. Considerations are given to the spectra of the streamwise velocity fluctuations, and to large scale motions and their energy contents from the pipe near-wall to centerline. The analysis of the velocity fluctuations revealed a Reynolds-number dependent inner peak at a fixed wall normal location, however, an outer peak seems not to appear that might be attributed either to low Reynolds number effect or not high enough spatial resolution of hot-wire probe, motivating further study utilizing nanoscale probes. Sizes of the large scale, and very large scale structures were estimated to have wavelengths of 3$R$, and 20$R$ at high Reynolds number, srespectively. The fractional energy contents in wavelengths associated with the large scale motions at various wall normal locations showed maximum contribution to the turbulent kinetic energy near the outer limit of the logarithmic layer.


title: The Structure of Integral Parabolic Subgroups of Orthogonal Groups

abstract: We determine the detailed structure of parabolic subgroups of orthogonal groups over $\mathbb{Z}$, and deduce the precise form of canonical boundary components in toroidal compactifications of orthogonal Shimura varieties.


title: Quantifying Perceptual Distortion of Adversarial Examples

abstract: Recent work has shown that additive threat models, which only permit the addition of bounded noise to the pixels of an image, are insufficient for fully capturing the space of imperceivable adversarial examples. For example, small rotations and spatial transformations can fool classifiers, remain imperceivable to humans, but have large additive distance from the original images. In this work, we leverage quantitative perceptual metrics like LPIPS and SSIM to define a novel threat model for adversarial attacks. To demonstrate the value of quantifying the perceptual distortion of adversarial examples, we present and employ a unifying framework fusing different attack styles. We first prove that our framework results in images that are unattainable by attack styles in isolation. We then perform adversarial training using attacks generated by our framework to demonstrate that networks are only robust to classes of adversarial perturbations they have been trained against, and combination attacks are stronger than any of their individual components. Finally, we experimentally demonstrate that our combined attacks retain the same perceptual distortion but induce far higher misclassification rates when compared against individual attacks.


title: Local Computation Algorithms for Spanners

abstract: A graph spanner is a fundamental graph structure that faithfully preserves the pairwise distances in the input graph up to a small multiplicative stretch. The common objective in the computation of spanners is to achieve the best-known existential size-stretch trade-off efficiently. Classical models and algorithmic analysis of graph spanners essentially assume that the algorithm can read the input graph, construct the desired spanner, and write the answer to the output tape. However, when considering massive graphs containing millions or even billions of nodes not only the input graph, but also the output spanner might be too large for a single processor to store. To tackle this challenge, we initiate the study of local computation algorithms (LCAs) for graph spanners in general graphs, where the algorithm should locally decide whether a given edge $(u,v) \in E$ belongs to the output spanner. Such LCAs give the user the `illusion' that a specific sparse spanner for the graph is maintained, without ever fully computing it. We present the following results: -For general $n$-vertex graphs and $r \in \{2,3\}$, there exists an LCA for $(2r-1)$-spanners with $\widetilde{O}(n^{1+1/r})$ edges and sublinear probe complexity of $\widetilde{O}(n^{1-1/2r})$. These size/stretch tradeoffs are best possible (up to polylogarithmic factors). -For every $k \geq 1$ and $n$-vertex graph with maximum degree $\Delta$, there exists an LCA for $O(k^2)$ spanners with $\widetilde{O}(n^{1+1/k})$ edges, probe complexity of $\widetilde{O}(\Delta^4 n^{2/3})$, and random seed of size $\mathrm{polylog}(n)$. This improves upon, and extends the work of [Lenzen-Levi, 2018]. We also complement our results by providing a polynomial lower bound on the probe complexity of LCAs for graph spanners that holds even for the simpler task of computing a sparse connected subgraph with $o(m)$ edges.


title: Convolutional Analysis Operator Learning: Dependence on Training Data

abstract: Convolutional analysis operator learning (CAOL) enables the unsupervised training of (hierarchical) convolutional sparsifying operators or autoencoders from large datasets. One can use many training images for CAOL, but a precise understanding of the impact of doing so has remained an open question. This paper presents a series of results that lend insight into the impact of dataset size on the filter update in CAOL. The first result is a general deterministic bound on errors in the estimated filters, and is followed by a bound on the expected errors as the number of training samples increases. The second result provides a high probability analogue. The bounds depend on properties of the training data, and we investigate their empirical values with real data. Taken together, these results provide evidence for the potential benefit of using more training data in CAOL.


title: Forecasting Angular Cross Correlations Between Diffuse X-ray Emission and the Thermal Sunyaev-Zel'dovich Effect

abstract: X-ray emission and the thermal Sunyaev-Zel'dovich distortion to the Cosmic Microwave Background are two important handles on the gas content of the Universe. The cross-correlation between these effects eliminates noise bias and reduces observational systematic effects. Using analytic models for the cluster profile, we develop a halo model formalism to study this cross-correlation and apply it to forecast the signal-to-noise of upcoming measurements from eROSITA and the Simons Observatory. In the soft X-ray band (0.5--2 keV), we forecast a signal-to-noise of 174 for the cross-power spectrum. Over a wide range of the scales, the X-rays will be signal-dominated, and so sample variance is important. In particular, non-Gaussian (4-point) contributions to the errors highlight the utility of masking massive clusters. Masking clusters down to $10^{14}\ M_\odot$ increases the signal-to-noise of the cross-spectrum to 201. We perform a Fisher Analysis on the fitting coefficients of the Battaglia et al. gas profiles and on cosmological parameters. We find that the cross spectrum is most sensitive to the overall scale of the profiles of pressure and electron density, as well as cosmological parameters $\sigma_8$ and $H_0$, but that the large number of parameters form a degenerate set, which makes extracting the information more challenging. Our modeling framework is flexible, and in the future, we can easily extend it to forecast the spatial cross-correlations of surveys of X-ray lines available to high-energy-resolution microcalorimetry, to studies of the Warm-Hot Intergalactic Medium, and other effects.


title: Photovoltaic Effect of Atomtronics Induced by Artificial Gauge Field

abstract: We investigate photovoltaic effect of atomtronics induced by artificial gauge field in four optical potentials. Effective magnetic flux gives rise to polarization of atom occupation probability which creates current of atomtronics. The relation between atomic current and magnetic flux behaves like the current-phase property in Josephson junction. The photovoltaic cell is well defined by the atomic opened system which have effective voltage and two different poles that correspond to two internal states of atomtronics. The atom flow is controllable by changing the direction of incident light and other system parameters. Detection of the atomic current intensity is available through light emission optical spectrum in experiments.


title: Torsion of elastic solids with sparse voids parallel to the twist axis

abstract: With the purpose of investigating a linear elastic solid containing a dilute distribution of cylindrical and prismatic holes parallel to the torsion axis, the full-field solution for an infinite elastic plane containing a single void and subject to torsion is derived. The obtained solution is exploited to derive the analytic expressions for the Stress Concentration Factor (SCF) related to the presence of an elliptical hole, for the Stress Intensity Factor (SIF) for hypocycloidal-shaped hole and star-shaped cracks, and for the Notch Stress Intensity Factor (NSIF) for star-shaped polygons. Special sets of the void location are obtained for which peculiar mechanical behaviours are displayed, such as the stress annihilation at some points along the boundary of elliptical voids and the stress singularity removal at the cusps/points of hypocycloidal shaped/isotoxal star-shaped polygonal voids. By means of finite element simulations it is finally shown that the presented closed-form expressions for the stress intensification provide reliable predictions even for finite domain realizations and, in particular, the infinite-plane solution remains highly accurate when the size of smooth and non-smooth external boundary is greater than twice and five times the void dimension, respectively. Under these geometrical conditions, the derived analytical expressions represent a valid 'guide tool' in mechanical design.


title: An IDEA: An Ingestion Framework for Data Enrichment in AsterixDB

abstract: Big Data today is being generated at an unprecedented rate from various sources such as sensors, applications, and devices, and it often needs to be enriched based on other reference information to support complex analytical queries. Depending on the use case, the enrichment operations can be compiled code, declarative queries, or machine learning models with different complexities. For enrichments that will be frequently used in the future, it can be advantageous to push their computation into the ingestion pipeline so that they can be stored (and queried) together with the data. In some cases, the referenced information may change over time, so the ingestion pipeline should be able to adapt to such changes to guarantee the currency and/or correctness of the enrichment results. In this paper, we present a new data ingestion framework that supports data ingestion at scale, enrichments requiring complex operations, and adaptiveness to reference data changes. We explain how this framework has been built on top of Apache AsterixDB and investigate its performance at scale under various workloads.


title: The computational power of parsing expression grammars

abstract: We study the computational power of parsing expression grammars (PEGs). We begin by constructing PEGs with unexpected behaviour, and surprising new examples of languages with PEGs, including the language of palindromes whose length is a power of two, and a binary-counting language. We then propose a new computational model, the scaffolding automaton, and prove that it exactly characterises the computational power of parsing expression grammars (PEGs). Using this characterisation we show that: (*) PEGs have unexpected power and semantics. We present several PEGs with surprising behaviour, and languages which, unexpectedly, have PEGs, including a PEG for the language of palindromes whose length is a power of two. (*) PEGs are computationally `universal', in the following sense: take any computable function $f:\{0,1\}^\ast\to \{0,1\}^\ast$; then there exists a computable function $g: \{0,1\}^\ast \to \mathbb{N}$ such that $\{ f(x) \#^{g(x)} x \mid x \in \{0,1\}^\ast \}$ has a PEG. (*) There can be no pumping lemma for PEGs. There is no total computable function $A$ with the following property: for every well-formed PEG $G$, there exists $n_0$ such that for every string $x \in \mathcal{L}(G)$ of size $|x| \ge n_0$, the output $y = A(G, x)$ is in $\mathcal{L}(G)$ and has $|y| > |x|$. (*) PEGs are strongly non real-time for Turing machines. There exists a language with a PEG, such that neither it nor its reverse can be recognised by any multi-tape online Turing machine which is allowed to do only $o(n/\log n)$ steps after reading each input symbol.


title: Gravitational waves from compact dark matter objects in the solar system

abstract: Dark matter could be composed of compact dark objects (CDOs). A close binary of CDOs orbiting in the interior of solar system bodies can be a loud source of gravitational waves (GWs) for the LIGO and VIRGO detectors. We perform the first search ever for this type of signal and rule out close binaries, with separations of order 300 m, orbiting near the center of the Sun with GW frequencies (twice the orbital frequency) between 50 and 550 Hz and CDO masses above $\approx 10^{-9} M_\odot$. This mass limit is eight orders of magnitude lower than the mass probed in a LIGO search at extra galactic distances.


title: An Online Decision-Theoretic Pipeline for Responder Dispatch

abstract: The problem of dispatching emergency responders to service traffic accidents, fire, distress calls and crimes plagues urban areas across the globe. While such problems have been extensively looked at, most approaches are offline. Such methodologies fail to capture the dynamically changing environments under which critical emergency response occurs, and therefore, fail to be implemented in practice. Any holistic approach towards creating a pipeline for effective emergency response must also look at other challenges that it subsumes - predicting when and where incidents happen and understanding the changing environmental dynamics. We describe a system that collectively deals with all these problems in an online manner, meaning that the models get updated with streaming data sources. We highlight why such an approach is crucial to the effectiveness of emergency response, and present an algorithmic framework that can compute promising actions for a given decision-theoretic model for responder dispatch. We argue that carefully crafted heuristic measures can balance the trade-off between computational time and the quality of solutions achieved and highlight why such an approach is more scalable and tractable than traditional approaches. We also present an online mechanism for incident prediction, as well as an approach based on recurrent neural networks for learning and predicting environmental features that affect responder dispatch. We compare our methodology with prior state-of-the-art and existing dispatch strategies in the field, which show that our approach results in a reduction in response time with a drastic reduction in computational time.


title: In situ analysis of catalyst composition during gold catalyzed GaAs nanowire growth

abstract: Semiconductor nanowires offer the opportunity to incorporate novel structures and functionality into electronic and optoelectronic devices. A clear understanding of the nanowire growth mechanism is essential for well-controlled growth of structures with desired properties, but the understanding is currently limited by a lack of empirical measurements of important parameters during growth, such as catalyst particle composition. However, this is difficult to accurately determine by investigating post-growth. We report direct measurement of the catalyst composition of individual gold seeded GaAs nanowires inside an electron microscope as they grow. The Ga content in the catalyst during growth increased with both temperature and Ga precursor flux. A direct comparison of the calculated phase diagrams of the Au-Ga-As ternary system to the measured catalyst composition not only lets us estimate the As content in the catalyst but also indicates the relevance of phase diagrams to understanding nanowire growth.


title: End-to-End Jet Classification of Quarks and Gluons with the CMS Open Data

abstract: We describe the construction of end-to-end jet image classifiers based on simulated low-level detector data to discriminate quark- vs. gluon-initiated jets with high-fidelity simulated CMS Open Data. We highlight the importance of precise spatial information and demonstrate competitive performance to existing state-of-the-art jet classifiers. We further generalize the end-to-end approach to event-level classification of quark vs. gluon di-jet QCD events. We compare the fully end-to-end approach to using hand-engineered features and demonstrate that the end-to-end algorithm is robust against the effects of underlying event and pile-up.


title: Accelerated Steady-State Torque Computation for Induction Machines using Parallel-In-Time Algorithms

abstract: This paper focuses on efficient steady-state computations of induction machines. In particular, the periodic Parareal algorithm with initial-value coarse problem (PP-IC) is considered for acceleration of classical time-stepping simulations via non-intrusive parallelization in time domain, i.e., existing implementations can be reused. Superiority of this parallel-in-time method is in its direct applicability to time-periodic problems, compared to, e.g, the standard Parareal method, which only solves an initial-value problem, starting from a prescribed initial value. PP-IC is exploited here to obtain the steady state of several operating points of an induction motor, developed by Robert Bosch GmbH. Numerical experiments show that acceleration up to several dozens of times can be obtained, depending on availability of parallel processing units. Comparison of PP-IC with existing time-periodic explicit error correction method highlights better robustness and efficiency of the considered time-parallel approach.


title: Thresholding normally distributed data creates complex networks

abstract: Network data sets are often constructed by some kind of thresholding procedure. The resulting networks frequently possess properties such as heavy-tailed degree distributions, clustering, large connected components and short average shortest path lengths. These properties are considered typical of complex networks and appear in many contexts, prompting consideration of their universality. Here we introduce a simple model for correlated relational data and study the network ensemble obtained by thresholding it. We find that some, but not all, of the properties associated with complex networks can be seen after thresholding the correlated data, even though the underlying data are not "complex". In particular, we observe heavy-tailed degree distributions, a large numbers of triangles, and short path lengths, while we do not observe non-vanishing clustering or community structure.


title: Rational points in the moduli space of genus two

abstract: We build a database of genus 2 curves defined over $\mathbb Q$ which contains all curves with minimal absolute height $h \leq 5$, all curves with moduli height $\mathfrak h \leq 20$, and all curves with extra automorphisms in standard form $y^2=f(x^2)$ defined over $\mathbb Q$ with height $h \leq 101$. For each isomorphism class in the database, an equation over its minimal field of definition is provided, the automorphism group of the curve, Clebsch and Igusa invariants. The distribution of rational points in the moduli space $\mathcal M_2$ for which the field of moduli is a field of definition is discussed and some open problems are presented.


title: On Singularities of Flat Affine Systems With $n$ States and $n-1$ Controls

abstract: We study the set of intrinsic singularities of flat affine systems with $n-1$ controls and $n$ states using the notion of Lie-B\"acklund atlas, previously introduced by the authors. For this purpose, we prove two easily computable sufficient conditions to construct flat outputs as a set of independent first integrals of distributions of vector fields, the first one in a generic case, namely in a neighborhood of a point where the $n-1$ control vector fields are independent, and the second one at a degenerate point where $p-1$ control vector fields are dependent of the $n-p$ others, with $p>1$. We show that the set of intrinsic singularities includes the set of points where the system does not satisfy the strong accessibility rank condition and is included in the set where the distribution of vector fields, introduced in the generic case, is singular. We conclude this analysis by three examples of apparent singularites of flat systems in generic and non generic degenerate cases.


title: Serre duality for Khovanov-Rozansky homology

abstract: We prove that the full twist is a Serre functor in the homotopy category of type A Soergel bimodules. As a consequence, we relate the top and bottom Hochschild degrees in Khovanov-Rozansky homology, categorifying a theorem of K\'alm\'an.


title: Constraints on Superconducting Cosmic Strings from the Global $21$-cm Signal before Reionization

abstract: Electromagnetic radiation from the cusp region of superconducting cosmic strings leads to a radio excess in the photon spectrum in the early universe and can produce a deep absorption feature in the global 21cm signal before the epoch of reionization. We study the constraints on the parameter space of superconducting strings which can be derived by demanding that the absorption feature is not larger in amplitude than what has recently been reported by the EDGES collaboration.


title: Capuchin: Causal Database Repair for Algorithmic Fairness

abstract: Fairness is increasingly recognized as a critical component of machine learning systems. However, it is the underlying data on which these systems are trained that often reflect discrimination, suggesting a database repair problem. Existing treatments of fairness rely on statistical correlations that can be fooled by statistical anomalies, such as Simpson's paradox. Proposals for causality-based definitions of fairness can correctly model some of these situations, but they require specification of the underlying causal models. In this paper, we formalize the situation as a database repair problem, proving sufficient conditions for fair classifiers in terms of admissible variables as opposed to a complete causal model. We show that these conditions correctly capture subtle fairness violations. We then use these conditions as the basis for database repair algorithms that provide provable fairness guarantees about classifiers trained on their training labels. We evaluate our algorithms on real data, demonstrating improvement over the state of the art on multiple fairness metrics proposed in the literature while retaining high utility.


title: Emission Measures and Emission-measure-weighted Temperatures of Shocked ISM and Ejecta in Supernova Remnants

abstract: A goal of supernova remnant (SNR) evolution models is to relate fundamental parameters of a supernova (SN) explosion and progenitor star to the current state of its SNR. The SNR hot plasma is characterized by its observed X-ray spectrum, which yields electron temperature, emission measure and abundances. Depending on their brightness, the properties of the plasmas heated by the SNR forward shock, reverse shock or both can be measured. The current work utilizes models which are spherically symmetric. One dimensional hydrodynamic simulations are carried out for SNR evolution prior to onset of radiative losses. From these, we derive dimensionless emission measures and emission-measure-weighted temperatures, and we present fitting formulae for these quantities as functions of scaled SNR time. These models allow one to infer SNR explosion energy, circumstellar medium density, age, ejecta mass and ejecta density profile from SNR observations. The new results are incorporated into the SNR modelling code SNRPy. The code is demonstrated with application to three historical SNRs: Kepler, Tycho and SN1006.


title: Bayes Optimal Early Stopping Policies for Black-Box Optimization

abstract: We derive an optimal policy for adaptively restarting a randomized algorithm, based on observed features of the run-so-far, so as to minimize the expected time required for the algorithm to successfully terminate. Given a suitable Bayesian prior, this result can be used to select the optimal black-box optimization algorithm from among a large family of algorithms that includes random search, Successive Halving, and Hyperband. On CIFAR-10 and ImageNet hyperparameter tuning problems, the proposed policies offer up to a factor of 13 improvement over random search in terms of expected time to reach a given target accuracy, and up to a factor of 3 improvement over a baseline adaptive policy that terminates a run whenever its accuracy is below-median.


title: From Mott to not: phenomenology of overdoped cuprates

abstract: Recently, we have argued that experimental data on superfluid density and terahertz conductivity of overdoped LSCO are compatible with a Landau Fermi liquid/Bardeen-Cooper-Schrieffer description of these samples, provided dopants are treated within "dirty $d$-wave theory" as weak scatterers. Here we test these ideas by comparing to specific heat and thermal conductivity data on LSCO, showing that the theory works extremely well across the overdoped region for similar disorder parameters. We then study the same properties in another overdoped cuprate, Tl-2201, thought to be quite "clean" since it exhibits quantum oscillations, low residual resistivities and small superconducting state Sommerfeld coefficients. Our results are consistent with the Tl-2201 system being $\approx 3$ times cleaner due in part to the dopant atoms' being located further from the CuO$_2$ plane. We conclude that cuprates can be described semiquantitatively in the overdoped regime by "dirty $d$-wave" theory, subject to significant Fermi liquid renormalizations, without introducing physics beyond the Landau-BCS paradigm.


title: Can a negative-mass cosmology explain dark matter and dark energy?

abstract: A recent study by Farnes (2018) proposed an alternative cosmological model in which both dark matter and dark energy are replaced with a single fluid of negative mass. This paper presents a critical review of that model. A number of problems and discrepancies with observations are identified. For instance, the predicted shape and density of galactic dark matter halos are incorrect. Also, halos would need to be less massive than the baryonic component otherwise they would become gravitationally unstable. Perhaps the most challenging problem in this theory is the presence of a large-scale version of the `runaway' effect, which would result in all galaxies moving in random directions at nearly the speed of light. Other more general issues regarding negative mass in general relativity are discussed, such as the possibility of time-travel paradoxes.


title: L2 Observers for a Class of Nonlinear Systems with Unknown Inputs

abstract: We consider the problem of estimating the state and unknown input for a large class of nonlinear systems subject to unknown exogenous inputs. The exogenous inputs themselves are modeled as being generated by a nonlinear system subject to unknown inputs. The nonlinearities considered in this work are characterized by multiplier matrices that include many commonly encountered nonlinearities. We obtain a linear matrix inequality (LMI), that, if feasible, provides the gains for an observer which results in certified L2 performance of the error dynamics associated with the observer. We also present conditions which guarantee that the L2 norm of the error can be made arbitrarily small and investigate conditions for feasibility of the proposed LMIs.


title: Using an online sample to learn about an offline population

abstract: Online data sources offer tremendous promise to demography and other social sciences, but researchers worry that the group of people who are represented in online datasets can be different from the general population. We show that by sampling and anonymously interviewing people who are online, researchers can learn about both people who are online and people who are offline. Our approach is based on the insight that people everywhere are connected through in-person social networks, such as kin, friendship, and contact networks. We illustrate how this insight can be used to derive an estimator for tracking the *digital divide* in access to the internet, an increasingly important dimension of population inequality in the modern world. We conducted a large-scale empirical test of our approach, using an online sample to estimate internet adoption in five countries ($n \approx 15,000$). Our test embedded a randomized experiment whose results can help design future studies. Our approach could be adapted to many other settings, offering one way to overcome some of the major challenges facing demographers in the information age.


title: Manifold valued data analysis of samples of networks, with applications in corpus linguistics

abstract: Networks arise in many applications, such as in the analysis of text documents, social interactions and brain activity. We develop a general framework for extrinsic statistical analysis of samples of networks, motivated by networks representing text documents in corpus linguistics. We identify networks with their graph Laplacian matrices, for which we define metrics, embeddings, tangent spaces, and a projection from Euclidean space to the space of graph Laplacians. This framework provides a way of computing means, performing principal component analysis and regression, and carrying out hypothesis tests, such as for testing for equality of means between two samples of networks. We apply the methodology to the set of novels by Jane Austen and Charles Dickens.


title: How I Learned to Stop Worrying and Love Re-optimization

abstract: Cost-based query optimizers remain one of the most important components of database management systems for analytic workloads. Though modern optimizers select plans close to optimal performance in the common case, a small number of queries are an order of magnitude slower than they could be. In this paper we investigate why this is still the case, despite decades of improvements to cost models, plan enumeration, and cardinality estimation. We demonstrate why we believe that a re-optimization mechanism is likely the most cost-effective way to improve end-to-end query performance. We find that even a simple re-optimization scheme can improve the latency of many poorly performing queries. We demonstrate that re-optimization improves the end-to-end latency of the top 20 longest running queries in the Join Order Benchmark by 27%, realizing most of the benefit of perfect cardinality estimation.


title: Magnetic force theory combined with quasi-particle self-consistent GW method

abstract: We report a successful combination of magnetic force linear response theory with quasiparticle self-consistent GW method. The self-consistently determined wavefunctions and eigenvalues can just be used for the conventional magnetic force calculations. While its formulation is straightforward, this combination provides a way to investigate the effect of GW self-energy on the magnetic interactions which can hardly be quantified due to the limitation of current GW methodology in calculating the total energy difference in between different magnetic phases. In ferromagnetic $3d$ elements, GW self-energy slightly reduces the $d$ bandwidth and enhances the interactions while the same long-range feature is maintained. In antiferromagnetic transition-metal monoxides, QSGW significantly reduces the interaction strengths by enlarging the gap. Orbital-dependent magnetic force calculations show that the coupling between $e_g$ and the nominally-empty $4s$ orbital is noticeably large in MnO which is reminiscent of the discussion for cuprates regarding the role of Cu-$4s$ state. This combination of magnetic force theory with quasiparticle self-consistent GW can be a useful tool to study various magnetic materials.


title: X-ray spectroscopy of the candidate AGN in Henize 2-10 and NGC 4178: Likely supernova remnants

abstract: Black holes in dwarf/bulgeless galaxies play a crucial role in studying the co-evolution of galaxies and their central black holes. Identifying massive black holes in dwarf galaxies suggests that the growth of black holes could precede that of galaxies. However, some of the most intriguing candidate active galactic nuclei (AGN) in small galaxies have such low luminosities that the sample is vulnerable to contamination by other sources, such as supernova remnants. We re-analysed Chandra X-ray Observatory observations of candidate AGN in Henize 2-10 and NGC 4178, considering the potential signals of emission lines in the minimally-binned X-ray spectra. We find that hot plasma models, which are typical of supernova remnants, explain the observed spectra much better than simple power-law models, which are appropriate for AGN. We identify clear signals of X-ray lines in the faint X-ray source identified with the radio source in Henize 2-10 by Reines et al. 2016. Combining our work with the MUSE measurement of the ionization parameter in this region by Cresci et al. 2017 indicates that this radio and X-ray source is more likely a supernova remnant than an AGN. A similar analysis of the low-count X-ray spectrum of a candidate AGN in NGC 4178 shows that a hot plasma model is about seventeen times more probable than a simple power-law model. Our results indicate that investigation of X-ray spectra, even in a low-count regime, can be a crucial tool to identify thermally-dominated supernova remnants among AGN candidates.


title: Hybrid Block Successive Approximation for One-Sided Non-Convex Min-Max Problems: Algorithms and Applications

abstract: The min-max problem, also known as the saddle point problem, is a class of optimization problems which minimizes and maximizes two subsets of variables simultaneously. This class of problems can be used to formulate a wide range of signal processing and communication (SPCOM) problems. Despite its popularity, most existing theory for this class has been mainly developed for problems with certain special convex-concave structure. Therefore, it cannot be used to guide the algorithm design for many interesting problems in SPCOM, where various kinds of non-convexity arise. In this work, we consider a block-wise one-sided non-convex min-max problem, in which the minimization problem consists of multiple blocks and is non-convex, while the maximization problem is (strongly) concave. We propose a class of simple algorithms named Hybrid Block Successive Approximation (HiBSA), which alternatingly perform gradient descent-type steps for the minimization blocks and gradient ascent-type steps for the maximization problem. A key element in the proposed algorithm is the use of certain regularization and penalty sequences, which stabilize the algorithm and ensure convergence. We show that HiBSA converges to some properly defined first-order stationary solutions with quantifiable global rates. To validate the efficiency of the proposed algorithms, we conduct numerical tests on a number of problems, including the robust learning problem, the non-convex min-utility maximization problems, and certain wireless jamming problem arising in interfering channels.


title: Lingvo: a Modular and Scalable Framework for Sequence-to-Sequence Modeling

abstract: Lingvo is a Tensorflow framework offering a complete solution for collaborative deep learning research, with a particular focus towards sequence-to-sequence models. Lingvo models are composed of modular building blocks that are flexible and easily extensible, and experiment configurations are centralized and highly customizable. Distributed training and quantized inference are supported directly within the framework, and it contains existing implementations of a large number of utilities, helper functions, and the newest research ideas. Lingvo has been used in collaboration by dozens of researchers in more than 20 papers over the last two years. This document outlines the underlying design of Lingvo and serves as an introduction to the various pieces of the framework, while also offering examples of advanced features that showcase the capabilities of the framework.


title: On the propagation of regularity for solutions of the fractional Korteweg-de Vries equation

abstract: We consider the initial value problem (IVP) for the fractional Korteweg-de Vries equation (fKdV) \begin{equation}\label{abstracteq1} \left\{ \begin{array}{ll} \partial_{t}u-D_{x}^{\alpha}\partial_{x}u+u\partial_{x}u=0, & x,t\in\mathbb{R},\,0<\alpha<1, \\ u(x,0)=u_{0}(x).& \\ \end{array} \right. \end{equation} It has been shown that the solutions to certain dispersive equations satisfy the propagation of regularity phenomena. More precisely, it deals in determine whether regularity of the initial data on the right hand side of the real line is propagated to the left hand side by the flow solution. This property was found originally in solutions of Korteweg-de Vries (KdV) equation and it has been also verified in other dispersive equations as the Benjamin-Ono (BO) equation. Recently, it has been shown that the solutions of the dispersive generalized Benjamin-Ono (DGBO) equation, this is $\alpha\in (2,3)$ in \eqref{abstracteq1}; also satisfy the propagation of regularity phenomena. This is achieved by introducing a commutator decomposition to handle the dispersive part in the equation. Following the approach used in the DGBO, we prove that the solutions of the fKdV also satisfies the propagation of regularity phenomena. Consequently, this type of regularity travels with infinite speed to its left as time evolves.


title: Solving a Class of Non-Convex Min-Max Games Using Iterative First Order Methods

abstract: Recent applications that arise in machine learning have surged significant interest in solving min-max saddle point games. This problem has been extensively studied in the convex-concave regime for which a global equilibrium solution can be computed efficiently. In this paper, we study the problem in the non-convex regime and show that an \varepsilon--first order stationary point of the game can be computed when one of the player's objective can be optimized to global optimality efficiently. In particular, we first consider the case where the objective of one of the players satisfies the Polyak-{\L}ojasiewicz (PL) condition. For such a game, we show that a simple multi-step gradient descent-ascent algorithm finds an \varepsilon--first order stationary point of the problem in \widetilde{\mathcal{O}}(\varepsilon^{-2}) iterations. Then we show that our framework can also be applied to the case where the objective of the "max-player" is concave. In this case, we propose a multi-step gradient descent-ascent algorithm that finds an \varepsilon--first order stationary point of the game in \widetilde{\cal O}(\varepsilon^{-3.5}) iterations, which is the best known rate in the literature. We applied our algorithm to a fair classification problem of Fashion-MNIST dataset and observed that the proposed algorithm results in smoother training and better generalization.


title: Good Wild Harmonic Bundles and Good Filtered Higgs Bundles

abstract: We prove the Kobayashi-Hitchin correspondence between good wild harmonic bundles and polystable good filtered $\lambda$-flat bundles satisfying a vanishing condition. We also study the correspondence for good wild harmonic bundles with the homogeneity with respect to a group action, which is expected to provide another way to construct Frobenius manifolds.


title: The Power of Self-Reducibility: Selectivity, Information, and Approximation

abstract: This chapter provides a hands-on tutorial on the important technique known as self-reducibility. Through a series of "Challenge Problems" that are theorems that the reader will---after being given definitions and tools---try to prove, the tutorial will ask the reader not to read proofs that use self-reducibility, but rather to discover proofs that use self-reducibility. In particular, the chapter will seek to guide the reader to the discovery of proofs of four interesting theorems---whose focus areas range from selectivity to information to approximation---from the literature, whose proofs draw on self-reducibility. The chapter's goal is to allow interested readers to add self-reducibility to their collection of proof tools. The chapter simultaneously has a related but different goal, namely, to provide a "lesson plan" (and a coordinated set of slides is available online to support this use [Hem19]) for a lecture to a two-lecture series that can be given to undergraduate students---even those with no background other than basic discrete mathematics and an understanding of what polynomial-time computation is---to immerse them in hands-on proving, and by doing that, to serve as an invitation to them to take courses on Models of Computation or Complexity Theory.


title: Localized model reduction for parameterized problems

abstract: In this contribution we present a survey of concepts in localized model order reduction methods for parameterized partial differential equations. The key concept of localized model order reduction is to construct local reduced spaces that have only support on part of the domain and compute a global approximation by a suitable coupling of the local spaces. In detail, we show how optimal local approximation spaces can be constructed and approximated by random sampling. An overview of possible conforming and non-conforming couplings of the local spaces is provided and corresponding localized a posteriori error estimates are derived. We introduce concepts of local basis enrichment, which includes a discussion of adaptivity. Implementational aspects of localized model reduction methods are addressed. Finally, we illustrate the presented concepts for multiscale, linear elasticity and fluid-flow problems, providing several numerical experiments. This work has been accepted as a chapter in P. Benner, S. Grivet-Talocia, A. Quarteroni, G. Rozza, W.H.A. Schilders, L.M. Sileira. Handbook on Model Order Reduction. Walter De Gruyter GmbH, Berlin, 2019+.


title: Correlation in transport coefficients of hole-doped CuRhO$_2$ single crystals

abstract: To clarify the origin of the nontrivial thermoelectric properties observed in the delafossite oxide CuRhO$_2$ polycrystals, we have performed the systematic transport measurements on the single-crystalline CuRhO$_2$ samples. In the parent compound, we find a pronounced peak structure due to a phonon-drag effect in the temperature dependence of the Seebeck coefficient, which is also confirmed by the size effect experiments. In the Mg-substituted crystals, in contrast to the results of the polycrystals, both the resistivity and the Seebeck coefficient decrease with increasing Mg content $y$. In particular, the coefficient $A$ for the $T^2$ term of the resistivity and the $T$-linear coefficient for the Seebeck coefficient at low temperatures are well described within a simple relationship expected for metals, which is also applicable to the correlated materials with low carrier densities.


title: Very Large Scale Integration of Josephson-Junction-Based Superconductor Random Access Memories

abstract: Arrays of Vortex Transitional (VT) memory cells with functional density up to $1 Mbit/cm^2$ have been designed, fabricated, and successfully demonstrated. This progress is due to recent advances in design optimization and in superconductor electronics fabrication achieved at MIT Lincoln Laboratory. As a starting point, we developed a demo array of VT cells for the 100-${\mu}A/{\mu}m^2$ MIT LL fabrication process SFQ5ee with 8 niobium layers. The studied two-junction memory cell with a two-junction nondestructive readout occupied $168 {\mu}m^2$, resulting in an over $0.5 Mbit/cm^2$ functional density. Then, we reduced the cell area down to $99 {\mu}m^2$ (corresponding to over $0.9 Mbit/cm^2$ functional density) by utilizing self-shunted Josephson Junctions (JJs) with critical current density, $J_c$ of $600 {\mu}A/{\mu}m^2$ and eliminating shunt resistors. The fabricated high-$J_c$ memory cells were fully operational and possessed wide Read/Write current margins, quite close to the theoretically predicted values. We discuss approaches to further increasing the integration scale of superconductor memory and logic circuits: a) miniaturization of superconducting transformers by using soft magnetic materials; b) reduction of JJ area by using planar high-$J_c$ junctions similar to variable thickness bridges.


title: Bottom-up strategy for data retrieval and data entry over front-end application Software

abstract: Some people implement pattern and best practices without analyzing its efficiency on their projects. Consequently, our goal in this article is to convince software developers that it is worth to make an earnest effort to evaluate the use of best practices and software patterns. For such purpose, in this study we took a concrete case system for geographical locations inputs through user interfaces. Then, we performed a comparative study on a traditional method against our approach, named reverse logistic to retrieve results, by measuring the time that a user spends to perform actions when entering data into a system. Surprisingly, we had a decrease of 59% in the amount of time spent in comparison to the time spent on the traditional method. This result lays a foundation for feeding data from the typical final step and search based on string matching algorithms, speeding up the interaction between people and computer response


title: A Dictionary-Based Generalization of Robust PCA with Applications to Target Localization in Hyperspectral Imaging

abstract: We consider the decomposition of a data matrix assumed to be a superposition of a low-rank matrix and a component which is sparse in a known dictionary, using a convex demixing method. We consider two sparsity structures for the sparse factor of the dictionary sparse component, namely entry-wise and column-wise sparsity, and provide a unified analysis, encompassing both undercomplete and the overcomplete dictionary cases, to show that the constituent matrices can be successfully recovered under some relatively mild conditions on incoherence, sparsity, and rank. We leverage these results to localize targets of interest in a hyperspectral (HS) image based on their spectral signature(s) using the a priori known characteristic spectral responses of the target. We corroborate our theoretical results and analyze target localization performance of our approach via experimental evaluations and comparisons to related techniques.


title: Non-Iterative Characteristics Analysis for Ramp Loading with a Window

abstract: Ramp compression experiment are used to deduce the relation between compression and normal stress in a material, by measuring how a compression wave evolves as it propagates through different thicknesses of the sample material. The compression wave is generally measured by Doppler velocimetry from a surface that can be observed with optical or near-optical photons. For high-pressure ramp loading, the reflectivity of a free surface often decreases as it is accelerated by the ramp wave, and window materials transparent to the probing photons are used to keep the surface flatter and preserve its reflectivity. We previously described a method of analyzing ramp-wave data measured at the free surface which did not require numerical iteration. However, this method breaks down when the pressure at the surface changes and hence cannot be used for data taken with a finite-impedance window. We have now generalized this non-iterative analysis method to apply to measurements taken through a window. Free surfaces can be treated seamlessly, and the need for sampling at uniform intervals of velocity has been removed. These calculations require interpolation of partially-released states using the partially-constructed stress-compression relation, making them slower than the previous free-surface scheme, but they are still much more robust and fast than iterative analysis.


title: Interfacial closure of contacting surfaces

abstract: Understanding the contact between solid surfaces is a long standing problem which has a strong impact on the physics of many processes such as adhesion, friction, lubrication and wear. Experimentally, the investigation of solid/solid interfaces remains challenging today, due to the lack of experimental techniques able to provide sub-nanometer scale information on interfaces buried between millimeters of materials. Yet, a strong interest exists improving the modeling of contact mechanics of materials in order to adjust their interface properties (e.g. thermal transport, friction). We show here that the essential features of the residual gap between contacting surfaces can be measured using high energy X-ray synchrotron reflectivity. The presence of this nano-gap is general to the contact of solids. In some special case however, it can be removed when attractive forces take over repulsive contributions, depending on both height and wavelength of asperity distributions (roughness). A criterion for this instability is established in the standard case of van der Waals attractive forces and elastic asperity compression repulsive forces (Hertz model). This collapse instability is confirmed experimentally in the case of silicon direct bonding, using high-energy X-ray synchrotron reflectivity and adhesion energy measurements. The possibility to achieve fully closed interfaces at room temperature opens interesting perspectives to build stronger assemblies with smaller thermal budgets.


title: An\'alise t\'ermica\-fluidodin\^amica de um transformador de energia a seco atrav\'es da din\^amica dos fluidos computacional

abstract: The cooling system is a key factor in power transformer designs. The frequent need for optimized projects further increases the importance of development through of increasingly developed tools. This work aims to evaluate a transformer design, in which the refrigeration is performed through of exhaust fans and heat exchangers installed on the lower side of the transformer. The model is elaborated applying a symmetry, reducing the model to be analyzed. In order to evaluate the best configuration in terms of positioning, flow rate and quantity of understand the main thermal and fluid dynamics characteristics of the flow inside the transformer, the Computational Fluid Dynamics tool is used. The methodology employed is to solve the numerically, through the ANSYS CFX computational package, the flow and energy equations inside the transformer tank.


title: Simplify Power Flow Calculation Using Terminal Circuit and PMU Measurements

abstract: Power flow calculation methods have been developed in decades using power injections and Newton-Raphson method. The nonlinear characteristics of the power flow to the bus voltage require Jacobian matrix reformation and refactorization in each iteration. Power network is composed by resistors, reactors, and capacitors which is a linear circuit when investigating the node voltages with the node current injections. To take the advantage of the linearity, this paper proposed to use terminal circuit model and PMU voltage phase angle measurements to simplify power flow calculation. When updating current injections at the righthand side of power flow equations and using PMU voltage phase angle measurements representing PV buses voltage phase angle, the Jacobian matrix is constant during the iteration. The simplification reduces the computational efforts and improves the computation efficiency. The proposed method is tested on the IEEE 14-bus and IEEE 118-bus standard systems. The results are validated by traditional power flow solution and the computation efficiency is demonstrated.


title: The Laboratory of Material Analysis with Ion Beams, LAMFI-USP

abstract: LAMFI is a laboratory dedicated to the development and application of ion beam techniques for the analysis of bulk materials and thin films. Its main facilities comprise an 1.7MV Pelletron tandem accelerator and two analytical setups, one mainly for Rutherford Backscattering Spectrometry (RBS), and channeling, and another, for Particle Induced X-ray Emission (PIXE) analysis. Since its beginning in March 1992, over 10 institutions involving about 68 researchers, 42 graduate and 9 undergraduate students have been direct or indirectly users at LAMFI.


title: An Effective-Current Approach for Hall\'{e}n's Equation in Center-Fed Dipole Antennas with Finite Conductivity

abstract: We propose a remedy for the unphysical oscillations arising in the current distribution of carbon nanotube and imperfectly conducting antennas center-driven by a delta-function generator when the approximate kernel is used. We do so by formulating an effective current, which was studied in detail in a 2011 and a 2013 paper for a perfectly conducting linear cylindrical antenna of infinite length, with application to the finite-length antenna. We discuss our results in connection with the perfectly conducting antenna, providing perturbative corrections to the current distribution for a large conductance, as well as presenting a delta-sequence and the field of a Hertzian dipole for the effective current in the limit of vanishing conductance. To that end, we employ both analytical tools and numerical methods to compare with experimental results.


title: Kinetic models of BGK type and their numerical integration

abstract: This minicourse contains a description of recent results on the modelling of rarefied gases in weakly non equilibrium regimes, and the numerical methods used to approximate the resulting equations. Therefore this work focuses on BGK type approximations, rather than on full Boltzmann models. Within this framework, models for polyatomic gases and for mixtures will be considered. We will also address numerical issues characteristic of the difficulties one encounters when integrating kinetic equations. In particular, we will consider asymptotic preserving schemes, which are designed to approximate equilibrium solutions, without resolving the fast scales of the approach to equilibrium.


title: Asteroseismology of massive stars with the TESS mission: the runaway Beta Cep pulsator PHL 346 = HN Aqr

abstract: We report an analysis of the first known Beta Cep pulsator observed by the TESS mission, the runaway star PHL 346 = HN Aqr. The star, previously known as a singly-periodic pulsator, has at least 34 oscillation modes excited, 12 of those in the g-mode domain and 22 p modes. Analysis of archival data implies that the amplitude and frequency of the dominant mode and the stellar radial velocity were variable over time. A binary nature would be inconsistent with the inferred ejection velocity from the Galactic disc of 420 km/s, which is too large to be survivable by a runaway binary system. A kinematic analysis of the star results in an age constraint (23 +- 1 Myr) that can be imposed on asteroseismic modelling and that can be used to remove degeneracies in the modelling process. Our attempts to match the excitation of the observed frequency spectrum resulted in pulsation models that were too young. Hence, asteroseismic studies of runaway pulsators can become vital not only in tracing the evolutionary history of such objects, but to understand the interior structure of massive stars in general. TESS is now opening up these stars for detailed asteroseismic investigation.


title: Reduced-Rank Local Distance Metric Learning for k-NN Classification

abstract: We propose a new method for local distance metric learning based on sample similarity as side information. These local metrics, which utilize conical combinations of metric weight matrices, are learned from the pooled spatial characteristics of the data, as well as the similarity profiles between the pairs of samples, whose distances are measured. The main objective of our framework is to yield metrics, such that the resulting distances between similar samples are small and distances between dissimilar samples are above a certain threshold. For learning and inference purposes, we describe a transductive, as well as an inductive algorithm; the former approach naturally befits our framework, while the latter one is provided in the interest of faster learning. Experimental results on a collection of classification problems imply that the new methods may exhibit notable performance advantages over alternative metric learning approaches that have recently appeared in the literature.


title: The NIGENS General Sound Events Database

abstract: Computational auditory scene analysis is gaining interest in the last years. Trailing behind the more mature field of speech recognition, it is particularly general sound event detection that is attracting increasing attention. Crucial for training and testing reasonable models is having available enough suitable data -- until recently, general sound event databases were hardly found. We release and present a database with 714 wav files containing isolated high quality sound events of 14 different types, plus 303 `general' wav files of anything else but these 14 types. All sound events are strongly labeled with perceptual on- and offset times, paying attention to omitting in-between silences. The amount of isolated sound events, the quality of annotations, and the particular general sound class distinguish NIGENS from other databases.


title: Residual Spectrum of Power Bounded Operators

abstract: The residual spectrum of a power bounded operator lies in the open unit disk.


title: Mixed temperature-dependent order parameters in the extended Hubbard model

abstract: The extended Hubbard model can host s-wave, d-wave and p-wave superconducting phases depending on the values of the on-site and nearest-neighbour interactions. Upon detailed examination of the free energy functional of the gap in this model, we show that these symmetries are often dependent on temperature. The critical points of this functional are constrained by symmetry and allow us to formulate stringent conditions on the temperature profile of the gap function, applicable to other models as well. We discuss the finite temperature phase diagram of the extended Hubbard model, and point out the existence of symmetry transitions below $T_c$. Understanding the nature of these transitions is crucial to assessing the symmetry of unconventional superconductors.


title: Learned multi-stability in mechanical networks

abstract: We contrast the distinct frameworks of materials design and physical learning in creating elastic networks with desired stable states. In design, the desired states are specified in advance and material parameters can be optimized on a computer with this knowledge. In learning, the material physically experiences the desired stable states in sequence, changing the material so as to stabilize each additional state. We show that while designed states are stable in networks of linear Hookean springs, sequential learning requires specific non-linear elasticity. We find that such non-linearity stabilizes states in which strain is zero in some springs and large in others, thus playing the role of Bayesian priors used in sparse statistical regression. Our model shows how specific material properties allow continuous learning of new functions through deployment of the material itself.


title: Parsing Gigabytes of JSON per Second

abstract: JavaScript Object Notation or JSON is a ubiquitous data exchange format on the Web. Ingesting JSON documents can become a performance bottleneck due to the sheer volume of data. We are thus motivated to make JSON parsing as fast as possible. Despite the maturity of the problem of JSON parsing, we show that substantial speedups are possible. We present the first standard-compliant JSON parser to process gigabytes of data per second on a single core, using commodity processors. We can use a quarter or fewer instructions than a state-of-the-art reference parser like RapidJSON. Unlike other validating parsers, our software (simdjson) makes extensive use of Single Instruction, Multiple Data (SIMD) instructions. To ensure reproducibility, simdjson is freely available as open-source software under a liberal license.


title: Multi-marginal Schrodinger bridges

abstract: We consider the problem to identify the most likely flow in phase space, of (inertial) particles under stochastic forcing, that is in agreement with spatial (marginal) distributions that are specified at a set of points in time. The question raised generalizes the classical Schrodinger Bridge Problem (SBP) which seeks to interpolate two specified end-point marginal distributions of overdamped particles driven by stochastic excitation. While we restrict our analysis to second-order dynamics for the particles, the data represents partial (i.e., only positional) information on the flow at {\em multiple} time-points. The solution sought, as in SBP, represents a probability law on the space of paths this closest to a uniform prior while consistent with the given marginals. We approach this problem as an optimal control problem to minimize an action integral a la Benamou-Brenier, and derive a time-symmetric formulation that includes a Fisher information term on the velocity field. We underscore the relation of our problem to recent measure-valued splines in Wasserstein space, which is akin to that between SBP and Optimal Mass Transport (OMT). The connection between the two provides a Sinkhorn-like approach to computing measure-valued splines. We envision that interpolation between measures as sought herein will have a wide range of applications in signal/images processing as well as in data science in cases where data have a temporal dimension.


title: Environment mediated multipartite and multidimensional entanglement

abstract: Quantum entanglement is usually considered a fragile quantity and decoherence through coupling to an external environment, such as a thermal reservoir, can quickly destroy the entanglement resource. This doesn't have to be the case and the environment can be engineered to assist in the formation of entanglement. We investigate a system of qubits and higher dimensional spins interacting only through their mutual coupling to a reservoir. We explore the entanglement of multipartite and multidimensional system as mediated by the bath and show that at low temperatures and intermediate coupling strengths multipartite entanglement may form between qubits and between higher spins, i.e., qudits. We characterise the multipartite entanglement using an entanglement witness based upon the structure factor and demonstrate its validity versus the directly calculated entanglement of formation, suggesting possible experiments for its measure.


title: Comparison of Deep Neural Networks and Deep Hierarchical Models for Spatio-Temporal Data

abstract: Spatio-temporal data are ubiquitous in the agricultural, ecological, and environmental sciences, and their study is important for understanding and predicting a wide variety of processes. One of the difficulties with modeling spatial processes that change in time is the complexity of the dependence structures that must describe how such a process varies, and the presence of high-dimensional complex data sets and large prediction domains. It is particularly challenging to specify parameterizations for nonlinear dynamic spatio-temporal models (DSTMs) that are simultaneously useful scientifically and efficient computationally. Statisticians have developed deep hierarchical models that can accommodate process complexity as well as the uncertainties in the predictions and inference. However, these models can be expensive and are typically application specific. On the other hand, the machine learning community has developed alternative "deep learning" approaches for nonlinear spatio-temporal modeling. These models are flexible yet are typically not implemented in a probabilistic framework. The two paradigms have many things in common and suggest hybrid approaches that can benefit from elements of each framework. This overview paper presents a brief introduction to the deep hierarchical DSTM (DH-DSTM) framework, and deep models in machine learning, culminating with the deep neural DSTM (DN-DSTM). Recent approaches that combine elements from DH-DSTMs and echo state network DN-DSTMs are presented as illustrations.


title: Buried double CuO chains in YBa$_2$Cu$_4$O$_8$ uncovered by nano-ARPES

abstract: The electron dynamics in the CuO chains has been elusive in Y-Ba-Cu-O cuprate systems by means of standard angle-resolved photoemission spectroscopy (ARPES); cleaved sample exhibits areas terminated by both CuO-chain or BaO layers, and the size of a typical beam results in ARPES signals that are superposed from both terminations. Here, we employ spatially-resolved ARPES with submicrometric beam (nano-ARPES) to reveal the surface-termination-dependent electronic structures of the double CuO chains in YBa$_2$Cu$_4$O$_8$. We present the first observation of sharp metallic dispersions and Fermi surfaces of the double CuO chains buried underneath the CuO$_2$-plane block on the BaO terminated surface. While the observed Fermi surfaces of the CuO chains are highly one-dimensional, the electrons in the CuO-chains do not undergo significant electron correlations and no signature of a Tomonaga-Luttinger liquid nor a marginal Fermi liquid is found. Our works represent an important experimental step toward understanding of the charge dynamics and provides a starting basis for modelling the high-$T_c$ superconductivity in YBCO cuprate systems.


title: Singularities in Reissner-Nordstr\"om black holes

abstract: We study black holes produced by the collapse of a spherically symmetric charged scalar field in asymptotically flat space. We employ a late time expansion and show decaying fluxes of radiation through the event horizon imply the black hole must contain a null singularity on the Cauchy horizon and a central spacelike singularity.


title: A pattern recognition algorithm for quantum annealers

abstract: The reconstruction of charged particles will be a key computing challenge for the high-luminosity Large Hadron Collider (HL-LHC) where increased data rates lead to large increases in running time for current pattern recognition algorithms. An alternative approach explored here expresses pattern recognition as a Quadratic Unconstrained Binary Optimization (QUBO) using software and quantum annealing. At track densities comparable with current LHC conditions, our approach achieves physics performance competitive with state-of-the-art pattern recognition algorithms. More research will be needed to achieve comparable performance in HL-LHC conditions, as increasing track density decreases the purity of the QUBO track segment classifier.


title: The Spiral Structure of the Milky Way Galaxy from Radio Observations

abstract: The study of the Milky Way structure is vital importance for a better understanding of the universe. This study digs deeper into the structure of Milky Way Galaxys with emphasis to its spiral arms. Multiple data sources were sought including the NASA websites. As found out, the Milky Way is a disk shaped galaxy composed of four spiral arms, two large and the rest are smaller. The spiral arms are regions of actively forming new stars dominated by young stars, dust and gas. Dust and gas are primary materials for the formation of new stars. This dust is also responsible for the reddish appearance of the stars since it absorbs more blue light than red. Besides the disk, there are numerous older stars, which appear white in color according to the picture by COBE. Additionally, the Milky Way contains a dark halo, extensive but void of luminous stars. The extensive halo was discovered due to the gravitational pull it exerts on existing visible matter.


title: A systematic study of ULIRGs using near-infrared absorption bands reveals a strong UV environment in their star-forming regions

abstract: We present a systematic study of the 3.0 um H2O ice and the 3.4 um aliphatic carbon absorption features toward 48 local ultraluminous infrared galaxies (ULIRGs) using spectra obtained by the AKARI Infrared Camera to investigate the UV environment in their star-forming regions. All the ULIRGs in our sample exhibit a ratio of optical depth of H2O ice to silicate dust (tau3.0/tau9.7) that is lower than that in the Taurus dark cloud. This implies that ULIRGs cannot be described as an ensemble of low-mass star-forming regions and that a significant amount of high-mass star-forming regions contribute to star-forming clouds in local ULIRGs. The results also show that the ratios of optical depth of aliphatic carbon to silicate dust, tau3.4/tau9.7, exhibit diverse values. We investigate two effects that can affect this ratio: the geometric temperature gradient (which increases the ratio) and the intense UV environment (which decreases it). The geometric temperature gradient is typically considered as a sign of active galactic nuclei (AGN). ULIRGs with AGN signs (optical classification, NIR color, and a PAH emission strength of 3.3 um) indeed tend to exhibit a large tau3.4/tau9.7 ratio. However, we find that the presence of buried AGN is not the only cause of the geometric temperature gradient, because the enhancement of the ratio is also evident in pure starburst-like ULIRGs without these AGN signs. Regarding the intense UV environment in star-forming regions, the correlation between the aliphatic carbon ratio and the ratio of the [C II] 158 um line luminosity to the far-infrared luminosity (L[CII]/LFIR), which represents the UV environment in photodissociation regions, implies that the intense UV environment causes the decrease of the aliphatic carbon ratio. We find that an intense UV environment (G/nH > 3) in star-forming regions is needed for the aliphatic carbon ratio to be suppressed.


title: On How Users Edit Computer-Generated Visual Stories

abstract: A significant body of research in Artificial Intelligence (AI) has focused on generating stories automatically, either based on prior story plots or input images. However, literature has little to say about how users would receive and use these stories. Given the quality of stories generated by modern AI algorithms, users will nearly inevitably have to edit these stories before putting them to real use. In this paper, we present the first analysis of how human users edit machine-generated stories. We obtained 962 short stories generated by one of the state-of-the-art visual storytelling models. For each story, we recruited five crowd workers from Amazon Mechanical Turk to edit it. Our analysis of these edits shows that, on average, users (i) slightly shortened machine-generated stories, (ii) increased lexical diversity in these stories, and (iii) often replaced nouns and their determiners/articles with pronouns. Our study provides a better understanding on how users receive and edit machine-generated stories,informing future researchers to create more usable and helpful story generation systems.


title: Comparison between continuous- and discrete-mode coherent feedback for the Jaynes-Cummings model

abstract: Using the example of the Jaynes-Cummings model, we present a comparison between time-delayed coherent feedback mediated by reservoirs with continuous and discrete mode structures and work out their qualitative differences. In contrast to the discrete-mode case, the continuous-mode case results in the well-known single-delay dynamics which can, e.g., stabilize Rabi oscillations. The discrete-mode case, however, shows population trapping, not present in the continuous-mode model. Given these differences, we discuss the cavity output spectra and show how these characteristic properties are spectrally identifiable. This work demonstrates the fundamental difference between the continuous-mode case, which represents a truly dissipative mechanism, and the discrete-mode case that is in principle based on a coherent excitation exchange process.


title: E-LSTM-D: A Deep Learning Framework for Dynamic Network Link Prediction

abstract: Predicting the potential relations between nodes in networks, known as link prediction, has long been a challenge in network science. However, most studies just focused on link prediction of static network, while real-world networks always evolve over time with the occurrence and vanishing of nodes and links. Dynamic network link prediction thus has been attracting more and more attention since it can better capture the evolution nature of networks, but still most algorithms fail to achieve satisfied prediction accuracy. Motivated by the excellent performance of Long Short-Term Memory (LSTM) in processing time series, in this paper, we propose a novel Encoder-LSTM-Decoder (E-LSTM-D) deep learning model to predict dynamic links end to end. It could handle long term prediction problems, and suits the networks of different scales with fine-tuned structure. To the best of our knowledge, it is the first time that LSTM, together with an encoder-decoder architecture, is applied to link prediction in dynamic networks. This new model is able to automatically learn structural and temporal features in a unified framework, which can predict the links that never appear in the network before. The extensive experiments show that our E-LSTM-D model significantly outperforms newly proposed dynamic network link prediction methods and obtain the state-of-the-art results.


title: A weighted endpoint weak-type estimate for multilinear Calder\'on-Zygmund operators

abstract: Two proofs of a weighted weak-type $\left(1,\ldots,1;\frac{1}{m}\right)$ estimate for multilinear Calder\'on-Zygmund operators are given. The ideas are motivated by different proofs of the classical weak-type $(1,1)$ estimate for Calder\'on-Zygmund operators. One proof uses the Calder\'on-Zygmund decomposition, and the other proof is motivated by ideas of Nazarov, Treil, and Volberg.


title: Ample line bundles, global generation and $K_0$ on quasi-projective derived schemes

abstract: The purpose of this note is to extend some classical results on quasi-projective schemes to the setting of derived algebraic geometry. Namely, we want to show that any vector bundle on a derived scheme admitting an ample line bundle can be twisted to be globally generated. Moreover, we provide a presentation of $K^0(X)$ as the Grothendieck group of vector bundles modulo exact sequences on any quasi-projective derived scheme $X$.


title: On the Mean Radiative Efficiency of Accreting Massive Black Holes in AGNs And QSOs

abstract: Radiative efficiency is an important physical parameter that describes the fraction of accretion material converted to radiative energy for accretion onto massive black holes (MBHs). With the simplest So{\l}tan argument, the radiative efficiency of MBHs can be estimated by matching the mass density of MBHs in the local universe to the accreted mass density by MBHs during AGN/QSO phases. In this paper, we estimate the local MBH mass density through a combination of various determinations of the correlations between the masses of MBHs and the properties of MBH host galaxies, with the distribution functions of those galaxy properties. We also estimate the total energy density radiated by AGNs and QSOs by using various AGN/QSO X-ray luminosity functions in the literature. We then obtain several hundred estimates of the mean radiative efficiency of AGNs/QSOs. Under the assumption that those estimates are independent of each other and free of systematic effects, we apply the median statistics as described by Gott et al.\cite{got01} and find the mean radiative efficiency of AGNs/QSOs is $\epsilon=0.105^{+0.006}_{-0.008}$, which is consistent with the canonical value $\sim0.1$. Considering that about $20\%$ Compton-thick objects may be missed from current available X-ray surveys, the true mean radiative efficiency may be actually $\sim 0.12$.


title: Roton-induced Bose polaron in the presence of synthetic spin-orbit coupling

abstract: We theoretically investigate the quasiparticle (polaron) properties of an impurity immersing in a Bose-Einstein condensate with equal Rashba and Dresselhaus spin-orbit coupling at zero temperature. In the presence of spin-orbit coupling, all bosons can condense into a single plane-wave state with finite momentum, and the corresponding excitation spectrum shows an intriguing roton minimum. We find that the polaron properties are strongly modified by this roton minimum, where the ground state of attractive polaron acquires a nonzero momentum and anisotropic effective mass. Across the resonance of the interaction between impurity and atoms, the polaron evolves into a tight-binding dimer. We show that the evolution is not smooth when the roton structure of the condensate becomes apparent, and a first-order phase transition from a phonon-induced polaron to a roton-induced polaron is observed at a critical interaction strength.


title: The Absolute Orders on the Coxeter Groups $A_n$ and $B_n$ are Sperner

abstract: Over 50 years ago, Rota posted the following celebrated `Research Problem': prove or disprove that the partial order of partitions on an $n$-set (i.e., the refinement order) is Sperner. A counterexample was eventually discovered by Canfield in 1978. However, Harper and Kim recently proved that a closely related order --- i.e., the refinement order on the symmetric group --- is not only Sperner, but strong Sperner. Equivalently, the well-known absolute order on the symmetric group is strong Sperner. In this paper, we extend these results by giving a concise, elegant proof that the absolute orders on the Coxeter groups $A_n$ and $B_n$ are strong Sperner.


title: Polarization amplitude near quantum critical points

abstract: We discuss the polarization amplitude of quantum spin systems in one dimension. In particular, we closely investigate it in gapless phases of those systems based on the two-dimensional conformal field theory. The polarization amplitude is defined as the ground-state average of a twist operator which induces a large gauge transformation attaching the unit amount of the U(1) flux to the system. We show that the polarization amplitude under the periodic boundary condition is sensitive to perturbations around the fixed point of the renormalization-group flow rather than the fixed point itself even when the perturbation is irrelevant. This dependence is encoded into the scaling law with respect to the system size. In this paper, we show how and why the scaling law of the polarization amplitude encodes the information of the renormalization-group flow. In addition, we show that the polarization amplitude under the antiperiodic boundary condition is determined fully by the fixed point in contrast to that under the periodic one and that it visualizes clearly the nontriviality of spin systems in the sense of the Lieb-Schultz-Mattis theorem.


title: On the Sensitivity of Adversarial Robustness to Input Data Distributions

abstract: Neural networks are vulnerable to small adversarial perturbations. Existing literature largely focused on understanding and mitigating the vulnerability of learned models. In this paper, we demonstrate an intriguing phenomenon about the most popular robust training method in the literature, adversarial training: Adversarial robustness, unlike clean accuracy, is sensitive to the input data distribution. Even a semantics-preserving transformations on the input data distribution can cause a significantly different robustness for the adversarial trained model that is both trained and evaluated on the new distribution. Our discovery of such sensitivity on data distribution is based on a study which disentangles the behaviors of clean accuracy and robust accuracy of the Bayes classifier. Empirical investigations further confirm our finding. We construct semantically-identical variants for MNIST and CIFAR10 respectively, and show that standardly trained models achieve comparable clean accuracies on them, but adversarially trained models achieve significantly different robustness accuracies. This counter-intuitive phenomenon indicates that input data distribution alone can affect the adversarial robustness of trained neural networks, not necessarily the tasks themselves. Lastly, we discuss the practical implications on evaluating adversarial robustness, and make initial attempts to understand this complex phenomenon.


title: Provably size-guaranteed mesh generation with superconvergence

abstract: The properties and applications of superconvergence on size-guaranteed Delaunay triangulation generated by bubble placement method (BPM), are studied in this paper. First, we derive a mesh condition that the difference between the actual side length and the desired length $h$ is as small as ${\cal O}(h^{1+{\alpha}})$ $({\alpha}>0)$. Second, the superconvergence estimations are analyzed on linear and quadratic finite element for elliptic boundary value problem based on the above mesh condition. In particular, the mesh condition is suitable for many known superconvergence estimations of different equations. Numerical tests are provided to verify the theoretical findings and to exhibit the superconvergence property on BPM-based grids.


title: Topological insulator and particle pumping in a one-dimensional shaken optical lattice

abstract: We propose a simple method to simulate and detect topological insulators with cold atoms trapped in a one-dimensional bichromatic optical lattice subjected to a time-periodic modulation. The tight-binding form of this shaken system is equivalent to the periodically driven Aubry-Andre model. We demonstrate that this model can be mapped into a two-dimensional Chern insulator model, whose energy spectrum hosts a topological phase within an experimentally accessible parameter regime. By tuning the laser phase adiabatically, such one-dimensional system constitutes a natural platform to realize topological particle pumping. We show that the Chern number characterizing the topological features of this system can be measured by detecting the density shift after one cycle of pumping.