title: An error indicator-based adaptive reduced order model for nonlinear structural mechanics -- application to high-pressure turbine blades

abstract: The industrial application motivating this work is the fatigue computation of aircraft engines' high-pressure turbine blades. The material model involves nonlinear elastoviscoplastic behavior laws, for which the parameters depend on the temperature. For this application, the temperature loading is not accurately known and can reach values relatively close to the creep temperature: important nonlinear effects occur and the solution strongly depends on the used thermal loading. We consider a nonlinear reduced order model able to compute, in the exploitation phase, the behavior of the blade for a new temperature field loading. The sensitivity of the solution to the temperature makes {the classical unenriched proper orthogonal decomposition method} fail. In this work, we propose a new error indicator, quantifying the error made by the reduced order model in computational complexity independent of the size of the high-fidelity reference model. In our framework, when the {error indicator} becomes larger than a given tolerance, the reduced order model is updated using one time step solution of the high-fidelity reference model. The approach is illustrated on a series of academic test cases and applied on a setting of industrial complexity involving 5 million degrees of freedom, where the whole procedure is computed in parallel with distributed memory.


title: Search for heavy neutrinos at CERN SPS

abstract: The phenomenology in the neutrino sector requires physics beyond the Standard Model. One possibility is the existence of new massive leptonic states which could be probed at the high intensity machines. The present results on heavy neutral leptons from the study of kaon decays in flight with the NA48/2 and NA62 experiments are presented and the future prospects for such searches at CERN SPS are discussed.


title: k-Spectra of weakly-c-Balanced Words

abstract: A word $u$ is a scattered factor of $w$ if $u$ can be obtained from $w$ by deleting some of its letters. That is, there exist the (potentially empty) words $u_1,u_2,..., u_n$, and $v_0,v_1,..,v_n$ such that $u = u_1u_2...u_n$ and $w = v_0u_1v_1u_2v_2...u_nv_n$. We consider the set of length-$k$ scattered factors of a given word w, called here $k$-spectrum and denoted $\ScatFact_k(w)$. We prove a series of properties of the sets $\ScatFact_k(w)$ for binary strictly balanced and, respectively, $c$-balanced words $w$, i.e., words over a two-letter alphabet where the number of occurrences of each letter is the same, or, respectively, one letter has $c$-more occurrences than the other. In particular, we consider the question which cardinalities $n= |\ScatFact_k(w)|$ are obtainable, for a positive integer $k$, when $w$ is either a strictly balanced binary word of length $2k$, or a $c$-balanced binary word of length $2k-c$. We also consider the problem of reconstructing words from their $k$-spectra.


title: Fugacity versus chemical potential in nonadditive generalizations of the ideal Fermi-gas

abstract: We compare two approaches to the generalization of the ordinary Fermi-statistics based on the nonadditive Tsallis $q$-exponential used in the Gibbs factor instead of the conventional exponential function. Both numerical and analytical calculations are made for the chemical potential, fugacity, energy, and the specific heat of the ideal gas obeying such generalized types of statistics. In the approach based on the Gibbs factor containing the chemical potential, high temperature behavior of the specific heat significantly deviates from the expected classical limit, while at low temperatures it resembles that of the ordinary ideal Fermi-gas. On the contrary, when the fugacity enters as a multiplier at the Gibbs factor, the high-temperature limit reproduces the classical ideal gas correctly. At low temperatures, however, some interesting results are observed, corresponding to non-zero specific heat at the absolute zero temperature or a finite (non-zero) minimal temperature. These results, though exotic from the first glance, might be applicable in effective modeling of physical phenomena in various domains.


title: The Sequential Test for Chaos

abstract: This paper reveals a novel numerical method, the sequential test, which approves chaos through sequences of numbers observations. The method alights alongside the Lyapunov exponent and bifurcation diagram test. Explicitly elucidation of the method application for both continuous and discrete systemswas given affiliated with the corresponding algorithms. The theoretical results are exemplified on systems satisfying different types of definitions of chaos or numerical methods. The results are supplemented with illustrative graphics.


title: High temperature convergence of the KMS boundary conditions: The Bose-Hubbard model on a finite graph

abstract: The Kubo-Martin-Schwinger condition is a widely studied fundamental property in quantum statistical mechanics which characterises the thermal equilibrium states of quantum systems. In the seventies, G. Gallavotti and E. Verboven, proposed an analogue to the KMS condition for classical mechanical systems and highlighted its relationship with the Kirkwood-Salzburg equations and with the Gibbs equilibrium measures. In the present article, we prove that in a certain limiting regime of high temperature the classical KMS condition can be derived from the quantum condition in the simple case of the Bose-Hubbard dynamical system on a finite graph. The main ingredients of the proof are Golden-Thompson inequality, Bogoliubov inequality and semiclassical analysis.


title: Structure of our Galactic Bulge from CN Measurements

abstract: The double red clumps (DRCs) are now dominantly believed to be the strong observational line of evidence of the so-called X-shaped Galactic bar structures. Recently, Lee et al. reported a subtle mean \ds\ difference between the DRCs and suggested a dichotomic picture that can be seen in globular clusters: the faint red clump is the first generation, while the bright red clump corresponds to the second generation (SG). They argued that the magnitude difference between the DRCs is due to different stellar populations, and is not due to the geometric difference between the DRCs. Our reanalysis shows that their data do not appear to support the idea of the multiple population-induced DRCs in our Galactic bulge. We perform fully empirical Monte Carlo simulations and find that the shape of the \ds\ distributions is the most stringent evidence to pursue. Our results strongly suggest that the CN distributions toward the Galactic bulge are qualitatively consistent with the X-shaped Galactic bulge with a minor fraction of the SG of about 2 -- 3\%.


title: The expectation of cosmic ray proton and helium energy spectrum below 4 PeV measured by LHAASO

abstract: Large High Altitude Air Shower Observatory(LHAASO) is a composite cosmic ray observatory consisting of three detector arrays: kilometer square array (KM2A) which includes the electromagnetic detector array and muon detector array, water Cherenkov detector array (WCDA) and wide field of view Cherenkov telescope array (WFCTA). One of the main scientific objectives of LHAASO is to precisely measure the cosmic rays energy spectrum of individual components from 1014 eV to 1018 eV. The hybrid observation will be employed by LHAASO experiment, in which the lateral and longitudinal distributions of the extensive air shower can be observed simultaneously. Thus many kinds of parameters can be used for primary nuclei identification. In this paper, high purity cosmic ray simulation samples of light nuclei component are obtained through Multi-Variable Analysis. The apertures of 1/4 LHAASO array for pure proton and mixed proton and helium (H&He) samples are 900 m2Sr and 1800 m2Sr respectively. A prospect of proton and H&He spectra from 100 TeV to 4 PeV is discussed.

categories: astro-ph.HE astro-ph.IM hep-ex physics.ins-det


title: OpenTapioca: Lightweight Entity Linking for Wikidata

abstract: We propose a simple Named Entity Linking system that can be trained from Wikidata only. This demonstrates the strengths and weaknesses of this data source for this task and provides an easily reproducible baseline to compare other systems against. Our model is lightweight to train, to run and to keep synchronous with Wikidata in real time.


title: Regularity for the fully nonlinear parabolic thin obstacle problem

abstract: We prove $C^{1, \alpha}$ regularity (in the parabolic sense) for the viscosity solution of a boundary obstacle problem with a fully nonlinear parabolic equation in the interior. Following the method which was first introduced for the harmonic case by L. Caffarelli in 1979, we extend the results of I. Athanasopoulos (1982) who studied the linear parabolic case and the results of E. Milakis and L. Silvestre (2008) who treated the fully nonlinear elliptic case.


title: Preservation of normality by transducers

abstract: We consider input-deterministic finite state transducers with infinite inputs and infinite outputs, and we consider the property of Borel normality on infinite words. When these transducers are given by a strongly connected set of states, and when the input is a Borel normal sequence, the output is an infinite word such that every word has a frequency given by a weighted automaton over the rationals. We prove that there is an algorithm that decides in cubic time whether an input-deterministic transducer preserves normality.


title: The Seventh Answer Set Programming Competition: Design and Results

abstract: Answer Set Programming (ASP) is a prominent knowledge representation language with roots in logic programming and non-monotonic reasoning. Biennial ASP competitions are organized in order to furnish challenging benchmark collections and assess the advancement of the state of the art in ASP solving. In this paper, we report on the design and results of the Seventh ASP Competition, jointly organized by the University of Calabria (Italy), the University of Genova (Italy), and the University of Potsdam (Germany), in affiliation with the 14th International Conference on Logic Programming and Non-Monotonic Reasoning (LPNMR 2017). (Under consideration for acceptance in TPLP).


title: Data Augmentation Using GANs

abstract: In this paper we propose the use of Generative Adversarial Networks (GAN) to generate artificial training data for machine learning tasks. The generation of artificial training data can be extremely useful in situations such as imbalanced data sets, performing a role similar to SMOTE or ADASYN. It is also useful when the data contains sensitive information, and it is desirable to avoid using the original data set as much as possible (example: medical data). We test our proposal on benchmark data sets using different network architectures, and show that a Decision Tree (DT) classifier trained using the training data generated by the GAN reached the same, (and surprisingly sometimes better), accuracy and recall than a DT trained on the original data set.


title: Numerical Analysis of Unsteady Implicitly Constituted Incompressible Fluids: Three-Field Formulation

abstract: In the classical theory of fluid mechanics a linear relationship between the shear stress and the symmetric velocity gradient tensor is often assumed. Even when a nonlinear relationship is assumed, it is typically formulated in terms of an explicit relation. Implicit constitutive models provide a theoretical framework that generalises this, allowing for general implicit constitutive relations. Since it is generally not possible to solve explicitly for the shear stress in the constitutive relation, a natural approach is to include the shear stress as a fundamental unknown in the formulation of the problem. In this work we present a mixed formulation with this feature, discuss its solvability and approximation using mixed finite element methods, and explore the convergence of the numerical approximations to a weak solution of the model.


title: From Static to Dynamic Anomaly Detection with Application to Power System Cyber Security

abstract: Developing advanced diagnosis tools to detect cyber attacks is the key to security of power systems. It has been shown that multivariate data injection attacks can bypass bad data detection schemes typically built on static behavior of the systems, which misleads operators to disruptive decisions. In this article, we depart from the existing static viewpoint to develop a diagnosis filter that captures the dynamics signatures of such a multivariate intrusion. To this end, we introduce a dynamic residual generator approach formulated as robust optimization programs in order to detect a class of disruptive multivariate attacks that potentially remain stealthy in view of a static bad data detector. We investigate two possible desired features: (i) a non-zero transient and (ii) a non-zero steady-state behavior of the residual generator in the presence of an attack. In case (i), the problem is reformulated as a finite, but possibly non-convex, optimization program. We further develop a linear programming relaxation that improves the scalability, and as such practicality, of the diagnosis filter design. In case (ii), it turns out that the resulting robust program admits an exact convex reformulation, yielding a Nash equilibrium between the attacker and the residual generator. This assertion has an interesting implication: the proposed approach is not conservative in the sense that the additional knowledge of the worst-case attack does not improve the diagnosis performance. To illustrate our theoretical results, we implement the proposed diagnosis filter to detect multivariate attacks on the system measurements deployed to generate the so-called Automatic Generation Control signals in a three-area IEEE 39-bus system.


title: Heterometallic Metal-Organic Frameworks of MOF-5 and UiO-66 Families: Insight from Computational Chemistry

abstract: We study the energetic stability and structural features of bimetallic metal-organic frameworks. Such heterometallic MOFs, which can result from partial substitutions between two types of cations, can have specific physical or chemical properties used for example in catalysis or gas adsorption. We work here to provide through computational chemistry a microscopic understanding of bimetallic MOFs and the distribution of cations within their structure. We develop a methodology based on a systematic study of possible cation distributions at all cation ratios by means of quantum chemistry calculations at the density functional theory level. We analyze the energies of the resulting bimetallic frameworks and correlate them with various disorder descriptors (functions of the bimetallic framework topology, regardless of exact atomic positions). We apply our methodology to two families of MOFs known for heterometallicity: MOF-5 (with divalent metal ions) and UiO-66 (with tetravalent metal ions). We observe that bimetallicity is overall more favorable for pairs of cations with sizes very close to each other, owing to a charge transfer mechanism inside secondary building units. For cation pairs with significant mutual size difference, metal mixing is globally less favorable, and the energy signifantly correlates with the coordination environment of linkers, determining their ability to adapt the mixing-induced strains. This effect is particularly strong in the UiO-66 family because of high cluster coordination number.


title: Magnetization density distribution of Sr$_2$IrO$_4$: Deviation from a local $j_\text{eff}=1/2$ picture

abstract: $5d$ iridium oxides are of huge interest due to the potential for new quantum states driven by strong spin-orbit coupling. The strontium iridate Sr$_2$IrO$_4$ is particularly in the spotlight because of the so-called $j_\text{eff}=1/2$ state consisting of a quantum superposition of the three local $t_{2g}$ orbitals with -- in its most simple version -- nearly equal population, which stabilizes an unconventional Mott insulating state. Here, we report an anisotropic and aspherical magnetization density distribution measured by polarized neutron diffraction in a magnetic field up to 5~T at 4~K, which strongly deviates from a local \jeffHalf picture even when distortion-induced deviations from the equal weights of the orbital populations are taken into account. Once reconstructed by the maximum entropy method and multipole expansion model refinement, the magnetization density shows cross-shaped positive four lobes along the crystallographic tetragonal axes with a large spatial extent, showing that the $xy$ orbital contribution is dominant. The analogy to the superconducting copper oxide systems might then be weaker than commonly thought.


title: Simple yet efficient real-time pose-based action recognition

abstract: Recognizing human actions is a core challenge for autonomous systems as they directly share the same space with humans. Systems must be able to recognize and assess human actions in real-time. In order to train corresponding data-driven algorithms, a significant amount of annotated training data is required. We demonstrated a pipeline to detect humans, estimate their pose, track them over time and recognize their actions in real-time with standard monocular camera sensors. For action recognition, we encode the human pose into a new data format called Encoded Human Pose Image (EHPI) that can then be classified using standard methods from the computer vision community. With this simple procedure we achieve competitive state-of-the-art performance in pose-based action detection and can ensure real-time performance. In addition, we show a use case in the context of autonomous driving to demonstrate how such a system can be trained to recognize human actions using simulation data.


title: Conserved currents and $\text{T}\bar{\text{T}}_s$ irrelevant deformations of 2D integrable field theories

abstract: It has been recently discovered that the $\text{T}\bar{\text{T}}$ deformation is closely-related to Jackiw-Teitelboim gravity. At classical level, the introduction of this perturbation induces an interaction between the stress-energy tensor and space-time and the deformed EoMs can be mapped, through a field-dependent change of coordinates, onto the corresponding undeformed ones. The effect of this perturbation on the quantum spectrum is non-perturbatively described by an inhomogeneous Burgers equation. In this paper, we point out that there exist infinite families of models where the geometry couples instead to generic combinations of local conserved currents labelled by the Lorentz spin. In spirit, these generalisations are similar to the $\text{J}\bar{\text{T}}$ model as the resulting theories and the corresponding scattering phase factors are not Lorentz invariant. The link with the $\text{J}\bar{\text{T}}$ model is discussed in detail. While the classical setup described here is very general, we shall use the sine-Gordon model and its CFT limit as explanatory quantum examples. Most of the final equations and considerations are, however, of broader validity or easily generalisable to more complicated systems.


title: An Efficient Algorithm for the Fast Delivery Problem

abstract: We study a problem where k autonomous mobile agents are initially located on distinct nodes of a weighted graph (with n nodes and m edges). Each autonomous mobile agent has a predefined velocity and is only allowed to move along the edges of the graph. We are interested in delivering a package, initially positioned in a source node s, to a destination node y. The delivery is achieved by the collective effort of the autonomous mobile agents, which can carry and exchange the package among them. The objective is to compute a delivery schedule that minimizes the delivery time of the package. In this paper, we propose an O(kn log n + km) time algorithm for this problem. This improves the previous state-of-the-art O(k^2 m + k n^2 + APSP) time algorithm for this problem, where APSP stands for the running-time of an algorithm for the All-Pairs Shortest Paths problem.


title: Search for ultralight scalar dark matter with NANOGrav pulsar timing arrays

abstract: An ultralight scalar field is a candidate for the dark matter. The ultralight scalar dark matter with mass around $10^{-23}\,{\rm eV}$ induces oscillations of the pulse arrival time in the sensitive frequency range of the pulsar timing arrays. We search for the ultralight scalar dark matter using the North American Nanohertz Observatory for Gravitational Waves 11-year Data Set. We give the 95% confidence upper limit for the signal induced by the ultralight scalar dark matter. In comparison with the published Bayesian upper limits on the amplitude of the ultralight scalar dark matter obtained by Bayesian analysis using the Parkes Pulsar Timing Array 12-year data set (Porayko et al. 2018), we find three times stronger upper limit in the frequency range from $10^{-8.34}$ to $10^{-8.19}\,{ \rm Hz}$ which corresponds to the mass range from $9.45\times10^{-24}$ to $1.34\times10^{-23}\,{\rm eV}$. In terms of the energy density of the dark matter, we find that the energy density near the Earth is less than $7\,{\rm GeV/cm^3}$ in the range from $10^{-8.55}$ to $10^{-8.01}\,{ \rm Hz}$ (from $5.83\times10^{-24}$ to $2.02\times10^{-23}\,{\rm eV}$). The strongest upper limit on the the energy density is given by $2\,{\rm GeV/cm^3}$ at a frequency $10^{-8.28}\,{ \rm Hz}$ (corresponding to a mass $1.09\times10^{-23}\,{\rm eV}$). We find that the signal of the ultralight scalar dark matter can be explained by the solar system ephemeris effect. Also, we reveal that the model of the solar system ephemeris effect prefers parameters which are contrary to the expectation that noise will be reduced on all pulsars.


title: The peculiar statistical mechanics of Optimal Learning Machines

abstract: Optimal Learning Machines (OLM) are systems that extract maximally informative representation of the environment they are in contact with, or of the data they are presented. It has recently been suggested that these systems are characterised by an exponential distribution of energy levels. In order to understand the peculiar properties of OLM within a broader framework, I consider an ensemble of optimisation problems over functions of many variables, part of which describe a sub-system and the rest account for its interaction with a random environment. The number of states of the sub-system with a given value of the objective function obeys a stretched exponential distribution, with exponent $\gamma$, and the interaction part is drawn at random from the same distribution, independently for each configuration of the whole system. Systems with $\gamma=1$ then correspond to OLM, and we find that they sit at the boundary between two regions with markedly different properties. For all $\gamma>0$ the system exhibits a freezing phase transition. The transition is discontinuous for $\gamma<1$ and it is continuous for $\gamma>1$. The region $\gamma>1$ corresponds to learnable energy landscapes and the behaviour of the sub-system becomes predictable as the size of the environment exceeds a critical threshold. For $\gamma<1$, instead, the energy landscape is unlearnable and the behaviour of the system becomes more and more unpredictable as the size of the environment increases. Sub-systems with $\gamma=1$ (OLM) feature a behaviour which is independent of the relative size of the environment. This is consistent with the expectation that efficient representations should be largely independent of the level of detail of the description of the environment.


title: Universality for critical KCM: infinite number of stable directions

abstract: Kinetically constrained models (KCM) are reversible interacting particle systems on $\mathbb{Z}^d$ with continuous-time constrained Glauber dynamics. They are a natural non-monotone stochastic version of the family of cellular automata with random initial state known as $\mathcal{U}$-bootstrap percolation. KCM have an interest in their own right, owing to their use for modelling the liquid-glass transition in condensed matter physics. In two dimensions there are three classes of models with qualitatively different scaling of the infection time of the origin as the density of infected sites vanishes. Here we study in full generality the class termed `critical'. Together with the companion paper by Martinelli and two of the authors we establish the universality classes of critical KCM and determine within each class the critical exponent of the infection time as well as of the spectral gap. In this work we prove that for critical models with an infinite number of stable directions this exponent is twice the one of their bootstrap percolation counterpart. This is due to the occurrence of `energy barriers', which determine the dominant behaviour for these KCM but which do not matter for the monotone bootstrap dynamics. Our result confirms the conjecture of Martinelli, Morris and the last author, who proved a matching upper bound.

categories: math.PR cond-mat.stat-mech math-ph math.MP


title: Salient Object Detection in the Deep Learning Era: An In-Depth Survey

abstract: As an essential problem in computer vision, salient object detection (SOD) has attracted an increasing amount of research attention over the years. Recent advances in SOD are predominantly led by deep learning-based solutions (named deep SOD). To enable in-depth understanding of deep SOD, in this paper, we provide a comprehensive survey covering various aspects, ranging from algorithm taxonomy to unsolved issues. In particular, we first review deep SOD algorithms from different perspectives, including network architecture, level of supervision, learning paradigm, and object-/instance-level detection. Following that, we summarize and analyze existing SOD datasets and evaluation metrics. Then, we benchmark a large group of representative SOD models, and provide detailed analyses of the comparison results. Moreover, we study the performance of SOD algorithms under different attribute settings, which has not been thoroughly explored previously, by constructing a novel SOD dataset with rich attribute annotations covering various salient object types, challenging factors, and scene categories. We further analyze, for the first time in the field, the robustness of SOD models to random input perturbations and adversarial attacks. We also look into the generalization and difficulty of existing SOD datasets. Finally, we discuss several open issues of SOD and outline future research directions.


title: Epitaxial Stabilisation of ${\bf \mathrm{Ge_{1-x}Sn_x}}$ Alloys

abstract: The thermodynamic stability of germanium tin $\mathrm{Ge_{1-x}Sn_x}$ alloys is investigated across the composition range $0 \le x \le 1$ by applying density functional theory (DFT) together with the cluster expansion formalism (CE). It is known that GeSn alloys are immiscible and that non-equilibrium growth techniques are required to produce metastable films and nanostructures. Insight into the driving forces behind component segregation is gained by investigating the equilibrium thermodynamics of GeSn systems. The alloy free energy of mixing is computed by combining enthalpies from CE with entropy terms for configurational and vibrational degrees of freedom. Volume deformations due to the large mismatch in ionic radii are readily found to be the key driving force for immiscibility at all temperatures of relevance. This leads to a study of epitaxial stabilisation by employing latticed matched substrates to favour growth of alloys with fractional compositions of $\mathrm{x=0}$, approximately $\mathrm{x=0.5}$ and $\mathrm{x=1}$. Reduction of the free energy of mixing due to epitaxial strain in thin films is quantified for each substrate leading to indicators for growth of kinetically stable films.


title: The Douglas-Rachford Algorithm for Convex and Nonconvex Feasibility Problems

abstract: The Douglas-Rachford method, a projection algorithm designed to solve continuous optimization problems, forms the basis of a useful heuristic for solving combinatorial optimization problems. In order to successfully use the method, it is necessary to formulate the problem at hand as a feasibility problem with constraint sets having efficiently computable nearest points. In this self-contained tutorial, we develop the convergence theory of projection algorithms within the framework of fixed point iterations, explain how to devise useful feasibility problem formulations, and demonstrate the application of the Douglas-Rachford method to said formulations. The paradigm is then illustrated on two concrete problems: a generalization of the "eight queens puzzle" known as the "(m,n)-queens problem", and the problem of constructing a probability distribution with prescribed moments.


title: Knowledge Distillation via Route Constrained Optimization

abstract: Distillation-based learning boosts the performance of the miniaturized neural network based on the hypothesis that the representation of a teacher model can be used as structured and relatively weak supervision, and thus would be easily learned by a miniaturized model. However, we find that the representation of a converged heavy model is still a strong constraint for training a small student model, which leads to a high lower bound of congruence loss. In this work, inspired by curriculum learning we consider the knowledge distillation from the perspective of curriculum learning by routing. Instead of supervising the student model with a converged teacher model, we supervised it with some anchor points selected from the route in parameter space that the teacher model passed by, as we called route constrained optimization (RCO). We experimentally demonstrate this simple operation greatly reduces the lower bound of congruence loss for knowledge distillation, hint and mimicking learning. On close-set classification tasks like CIFAR100 and ImageNet, RCO improves knowledge distillation by 2.14% and 1.5% respectively. For the sake of evaluating the generalization, we also test RCO on the open-set face recognition task MegaFace.


title: A Scalable Handwritten Text Recognition System

abstract: Many studies on (Offline) Handwritten Text Recognition (HTR) systems have focused on building state-of-the-art models for line recognition on small corpora. However, adding HTR capability to a large scale multilingual OCR system poses new challenges. This paper addresses three problems in building such systems: data, efficiency, and integration. Firstly, one of the biggest challenges is obtaining sufficient amounts of high quality training data. We address the problem by using online handwriting data collected for a large scale production online handwriting recognition system. We describe our image data generation pipeline and study how online data can be used to build HTR models. We show that the data improve the models significantly under the condition where only a small number of real images is available, which is usually the case for HTR models. It enables us to support a new script at substantially lower cost. Secondly, we propose a line recognition model based on neural networks without recurrent connections. The model achieves a comparable accuracy with LSTM-based models while allowing for better parallelism in training and inference. Finally, we present a simple way to integrate HTR models into an OCR system. These constitute a solution to bring HTR capability into a large scale OCR system.


title: Chromospheric polarimetry through multi-line observations of the 850 nm spectral region III: Chromospheric jets driven by twisted magnetic fields

abstract: We investigate the diagnostic potential of the spectral lines at 850 nm for understanding the magnetism of the lower atmosphere. For that purpose, we use a newly developed 3D simulation of a chromospheric jet to check the sensitivity of the spectral lines to this phenomenon as well as our ability to infer the atmospheric information through spectropolarimetric inversions of noisy synthetic data. We start comparing the benefits of inverting the entire spectrum at 850 nm versus only the Ca II 8542 A spectral line. We found a better match of the input atmosphere for the former case, mainly at lower heights. However, the results at higher layers were not accurate. After several tests, we determined that we need to weight more the chromospheric lines than the photospheric ones in the computation of the goodness of the fit. The new inversion configuration allows us to obtain better fits and consequently more accurate physical parameters. Therefore, to extract the most from multi-line inversions, a proper set of weights needs to be estimated. Besides that, we conclude again that the lines at 850 nm, or a similar arrangement with Ca II 8542 A plus Zeeman sensitive photospheric lines, poses the best observing configuration for examining the thermal and magnetic properties of the lower solar atmosphere.


title: A microscopic description for polarization in particle scatterings

abstract: We propose a microscopic description for the polarization from the first principle through the spin-orbit coupling in particle collisions. The model is different from previous ones based on local equilibrium assumptions for the spin degree of freedom. It is based on scatterings of particles as wave packets, an effective method to deal with particle scatterings at specified impact parameters. The polarization is then the consequence of particle collisions in a non-equilibrium state of spins. The spin-vorticity coupling naturally emerges from the spin-orbit one encoded in polarized scattering amplitudes of collisional integrals when one assumes local equilibrium in momentum but not in spin.


title: Multiplication and division of the orbital angular momentum of light with diffractive transformation optics

abstract: We present a method to efficiently multiply or divide the orbital angular momentum (OAM) of light beams using a sequence of two optical elements. The key-element is represented by an optical transformation mapping the azimuthal phase gradient of the input OAM beam onto a circular sector. By combining multiple circular-sector transformations into a single optical element, it is possible to perform the multiplication of the value of the input OAM state by splitting and mapping the phase onto complementary circular sectors. Conversely, by combining multiple inverse transformations, the division of the initial OAM value is achievable, by mapping distinct complementary circular sectors of the input beam into an equal number of circular phase gradients. The optical elements have been fabricated in the form of phase-only diffractive optics with high-resolution electron-beam lithography. Optical tests confirm the capability of the multiplier optics to perform integer multiplication of the input OAM, while the designed dividers are demonstrated to correctly split up the input beam into a complementary set of OAM beams. These elements can find applications for the multiplicative generation of higher-order OAM modes, optical information processing based on OAM-beams transmission, and optical routing/switching in telecom.


title: Efficient Blind Deblurring under High Noise Levels

abstract: The goal of blind image deblurring is to recover a sharp image from a motion blurred one without knowing the camera motion. Current state-of-the-art methods have a remarkably good performance on images with no noise or very low noise levels. However, the noiseless assumption is not realistic considering that low light conditions are the main reason for the presence of motion blur due to requiring longer exposure times. In fact, motion blur and high to moderate noise often appear together. Most works approach this problem by first estimating the blur kernel $k$ and then deconvolving the noisy blurred image. In this work, we first show that current state-of-the-art kernel estimation methods based on the $\ell_0$ gradient prior can be adapted to handle high noise levels while keeping their efficiency. Then, we show that a fast non-blind deconvolution method can be significantly improved by first denoising the blurry image. The proposed approach yields results that are equivalent to those obtained with much more computationally demanding methods.


title: Physics potential of the International Axion Observatory (IAXO)

abstract: We review the physics potential of a next generation search for solar axions: the International Axion Observatory (IAXO). Endowed with a sensitivity to discover axion-like particles (ALPs) with a coupling to photons as small as $g_{a\gamma}\sim 10^{-12}$ GeV$^{-1}$, or to electrons $g_{ae}\sim$10$^{-13}$, IAXO has the potential to find the QCD axion in the 1 meV$\sim$1 eV mass range where it solves the strong CP problem, can account for the cold dark matter of the Universe and be responsible for the anomalous cooling observed in a number of stellar systems. At the same time, IAXO will have enough sensitivity to detect lower mass axions invoked to explain: 1) the origin of the anomalous "transparency" of the Universe to gamma-rays, 2) the observed soft X-ray excess from galaxy clusters or 3) some inflationary models. In addition, we review string theory axions with parameters accessible by IAXO and discuss their potential role in cosmology as Dark Matter and Dark Radiation as well as their connections to the above mentioned conundrums.


title: New insight on Young Stellar Objects accretion shocks -- a claim for NLTE opacities

abstract: Context. Accreted material onto CTTSs is expected to form a hot quasi-periodic plasma structure that radiates in X-rays. Simulations of this phenomenon only partly match with observations. They all rely on a static model for the chromosphere model and on the assumption that radiation and matter are decoupled. Aims. We explore the effects on the structure and on the dynamics of the accretion flow of both a shock-heated chromosphere and of the coupling between radiation and hydrodynamics. Methods. We simulate accretion columns falling onto a stellar chromosphere using the 1D ALE code AstroLabE. This code solves the hydrodynamics equations along with the two first momenta equations for radiation transfer, with the help of a dedicated opacity table for the coupling between matter and radiation. We derive the total electron and ions densities from collisional-radiative NLTE ionisation equilibrium. Results. The chromospheric acoustic heating has an impact on the duration of the cycle and on the structure of the heated slab. In addition, the coupling between radiation and hydrodynamics leads to a heating of the accretion flow and the chromosphere, inducing a possible unburial of the whole column. These two last conclusions are in agreement with the computed monochromatic intensity. Both effects (acoustic heating and radiation coupling) have an influence on the amplitude and temporal variations of the net X-ray luminosity, which varies between 30 and 94% of the incoming mechanical energy flux, depending on the model considered.


title: New results on pseudosquare avoidance

abstract: We start by considering binary words containing the minimum possible numbers of squares and antisquares (where an antisquare is a word of the form $x \overline{x}$), and we completely classify which possibilities can occur. We consider avoiding $x p(x)$, where $p$ is any permutation of the underlying alphabet, and $x t(x)$, where $t$ is any transformation of the underlying alphabet. Finally, we prove the existence of an infinite binary word simultaneously avoiding all occurrences of $x h(x)$ for every nonerasing morphism $h$ and all sufficiently large words $x$.


title: Weyl points in systems of multiple semiconductor-superconductor quantum dots

abstract: As an analogy to the Weyl point in k-space, we search for energy levels which close at a single point as a function of a three dimensional parameter space. Such points are topologically protected in the sense that any perturbation which acts on the two level subsystem can be corrected by tuning the control parameters. We find that parameter controlled Weyl points are ubiquitous in semiconductor-superconductor quantum dots and that they are deeply related to Majorana zero modes. In this paper, we present several semiconductor-superconductor quantum dot devices which host parameter controlled Weyl points. Further, we show how these points can be observed experimentally via conductance measurements.


title: Assessing the Sharpness of Satellite Images: Study of the PlanetScope Constellation

abstract: New micro-satellite constellations enable unprecedented systematic monitoring applications thanks to their wide coverage and short revisit capabilities. However, the large volumes of images that they produce have uneven qualities, creating the need for automatic quality assessment methods. In this work, we quantify the sharpness of images from the PlanetScope constellation by estimating the blur kernel from each image. Once the kernel has been estimated, it is possible to compute an absolute measure of sharpness which allows to discard low quality images and deconvolve blurry images before any further processing. The method is fully blind and automatic, and since it does not require the knowledge of any satellite specifications it can be ported to other constellations.


title: Laminar-turbulent coexistence in annular Couette flow

abstract: Annular Couette flow is the flow between two coaxial cylinders driven by the axial translation of the inner cylinder. It is investigated using direct numerical simulation in long domains, with an emphasis on the laminar-turbulent coexistence regime found for marginally low values of the Reynolds number. Three distinct flow regimes are demonstrated as the radius ratio $\eta$ is decreased from 0.8 to 0.5 and finally to 0.1. The high-$\eta$ regime features helically-shaped turbulent patches coexisting with laminar flow, as in planar shear flows. The moderate-$\eta$ regime does not feature any marked laminar-turbulent coexistence. In an effort to discard confinement effects, proper patterning is however recovered by artificially extending the azimuthal span beyond 2$\pi$. Eventually, the low-$\eta$ regime features localised turbulent structures different from the puffs commonly encountered in transitional pipe flow. In this new coexistence regime, turbulent fluctuations are surprisingly short-ranged. Implications are discussed in terms of phase transition and critical scaling.


title: Necessary and Sufficient Conditions on Measurements of Quantum Channels

abstract: Quantum supermaps are a higher-order generalization of quantum maps, taking quantum maps to quantum maps. It is known that any completely positive, trace non-increasing (CPTNI) map can be performed as part of a quantum measurement. By providing an explicit counterexample we show that, instead, not every quantum supermap sending a quantum channel to a CPTNI map can be realized in a measurement on quantum channels. We find that the supermaps that can be implemented in this way are exactly those transforming quantum channels into CPTNI maps even when tensored with the identity supermap. We link this result to the fact that the principle of causality fails in the theory of quantum supermaps.


title: PLOTS: Procedure Learning from Observations using Subtask Structure

abstract: In many cases an intelligent agent may want to learn how to mimic a single observed demonstrated trajectory. In this work we consider how to perform such procedural learning from observation, which could help to enable agents to better use the enormous set of video data on observation sequences. Our approach exploits the properties of this setting to incrementally build an open loop action plan that can yield the desired subsequence, and can be used in both Markov and partially observable Markov domains. In addition, procedures commonly involve repeated extended temporal action subsequences. Our method optimistically explores actions to leverage potential repeated structure in the procedure. In comparing to some state-of-the-art approaches we find that our explicit procedural learning from observation method is about 100 times faster than policy-gradient based approaches that learn a stochastic policy and is faster than model based approaches as well. We also find that performing optimistic action selection yields substantial speed ups when latent dynamical structure is present.


title: Transfer Entropy: where Shannon meets Turing

abstract: Transfer entropy is capable of capturing nonlinear source-destination relations between multi-variate time series. It is a measure of association between source data that are transformed into destination data via a set of linear transformations between their probability mass functions. The resulting tensor formalism is used to show that in specific cases, e.g., in the case the system consists of three stochastic processes, bivariate analysis suffices to distinguish true relations from false relations. This allows us to determine the causal structure as far as encoded in the probability mass functions of noisy data. The tensor formalism was also used to derive the Data Processing Inequality for transfer entropy.


title: Understanding electron-doped cuprate superconductors as hole superconductors

abstract: Since their experimental discovery in 1989, the electron-doped cuprate superconductors have presented both a major challenge and a major opportunity. The major challenge has been to determine whether these materials are fundamentally different from or essentially similar to their hole-doped counterparts; a major opportunity because answering this question would strongly constrain the possible explanations for what is the essential physics that leads to high temperature superconductivity in the cuprates, which is still not agreed upon. Here we argue that experimental results over the past 30 years on electron-doped cuprate materials have provided conclusive answers to these fundamental questions, by establishing that both in hole- and electron-doped cuprates, superconductivity originates in pairing of hole carriers in the same band. We discuss a model to describe this physics that is different from the generally accepted ones, and calculate physical observables that agree with experiment, in particular tunneling characteristics. We argue that our model is simpler, more natural and more compelling than other models. Unlike other models, ours was originally proposed before rather than after many key experiments were performed.


title: Location-Sector Analysis of International Profit Shifting on a Multilayer Ownership-Tax Network

abstract: Currently all countries including developing countries are expected to utilize their own tax revenues and carry out their own development for solving poverty in their countries. However, developing countries cannot earn tax revenues like developed countries partly because they do not have effective countermeasures against international tax avoidance. Our analysis focuses on treaty shopping among various ways to conduct international tax avoidance because tax revenues of developing countries have been heavily damaged through treaty shopping. To analyze the location and sector of conduit firms likely to be used for treaty shopping, we constructed a multilayer ownership-tax network and proposed multilayer centrality. Because multilayer centrality can consider not only the value owing in the ownership network but also the withholding tax rate, it is expected to grasp precisely the locations and sectors of conduit firms established for the purpose of treaty shopping. Our analysis shows that firms in the sectors of Finance & Insurance and Wholesale & Retail trade etc. are involved with treaty shopping. We suggest that developing countries make a clause focusing on these sectors in the tax treaties they conclude.


title: Emergence of magnetism in bulk amorphous palladium

abstract: Magnetism in palladium has been the subject of much work and speculation. Bulk crystalline palladium is paramagnetic with a high magnetic susceptibility. Palladium under pressure and palladium nanoclusters have generated interest to scrutinize its magnetic properties. Here we report another possibility: Palladium may become an itinerant ferromagnet in the amorphous bulk phase at atmospheric pressure. Atomic palladium is a d$^{10}$ element, whereas bulk crystalline Pd is a d$^{10-x}$(sp)$^{x}$ material; this, together with the possible presence of 'unsaturated bonds' in amorphous materials, may explain the remnant magnetism reported herein. This work presents and discusses magnetic effects in bulk amorphous palladium.


title: On a semismooth* Newton method for solving generalized equations

abstract: In the paper, a Newton-type method for the solution of generalized equations (GEs) is derived, where the linearization concerns both the single-valued and the multi-valued part of the considered GE. The method is based on the new notion of semismoothness${}^*$ which, together with a suitable regularity condition, ensure the local superlinear convergence. An implementable version of the new method is derived for a class of GEs, frequently arising in optimization and equilibrium models.


title: Magnetization in the zig-zag layered Ising model and orthogonal polynomials

abstract: We discuss the magnetization $M_m$ in the $m$-th column of the zig-zag layered 2D Ising model on a half-plane using Kadanoff-Ceva fermions and orthogonal polynomials techniques. Our main result gives an explicit representation of $M_m$ via $m\times m$ Hankel determinants constructed from the spectral measure of a certain Jacobi matrix which encodes the interaction parameters between the columns. We also illustrate our approach by giving short proofs of the classical Kaufman-Onsager-Yang and McCoy-Wu theorems in the homogeneous setup and expressing $M_m$ as a Toeplitz+Hankel determinant for the homogeneous sub-critical model in presence of a boundary magnetic field.


title: Realistic Hair Simulation Using Image Blending

abstract: In this presented work, we propose a realistic hair simulator using image blending for dermoscopic images. This hair simulator can be used for benchmarking and validation of the hair removal methods and in data augmentation for improving computer aided diagnostic tools. We adopted one of the popular implementation of image blending to superimpose realistic hair masks to hair lesion. This method was able to produce realistic hair masks according to a predefined mask for hair. Thus, the produced hair images and masks can be used as ground truth for hair segmentation and removal methods by inpainting hair according to a pre-defined hair masks on hairfree areas. Also, we achieved a realism score equals to 1.65 in comparison to 1.59 for the state-of-the-art hair simulator.


title: Axi-symmetrization near point vortex solutions for the 2D Euler equation

abstract: We prove a definitive theorem on the asymptotic stability of point vortex solutions to the full Euler equation in 2 dimensions. More precisely, we show that a small, Gevrey smooth, and compactly supported perturbation of a point vortex leads to a global solution of the Euler equation in 2D, which converges weakly as $t\to\infty$ to a radial profile with respect to the vortex. The position of the point vortex, which is time dependent, stabilizes rapidly and becomes the center of the final, radial profile. The mechanism that leads to stabilization is mixing and inviscid damping.


title: Critically Examining the "Neural Hype": Weak Baselines and the Additivity of Effectiveness Gains from Neural Ranking Models

abstract: Is neural IR mostly hype? In a recent SIGIR Forum article, Lin expressed skepticism that neural ranking models were actually improving ad hoc retrieval effectiveness in limited data scenarios. He provided anecdotal evidence that authors of neural IR papers demonstrate "wins" by comparing against weak baselines. This paper provides a rigorous evaluation of those claims in two ways: First, we conducted a meta-analysis of papers that have reported experimental results on the TREC Robust04 test collection. We do not find evidence of an upward trend in effectiveness over time. In fact, the best reported results are from a decade ago and no recent neural approach comes close. Second, we applied five recent neural models to rerank the strong baselines that Lin used to make his arguments. A significant improvement was observed for one of the models, demonstrating additivity in gains. While there appears to be merit to neural IR approaches, at least some of the gains reported in the literature appear illusory.


title: Video Object Segmentation and Tracking: A Survey

abstract: Object segmentation and object tracking are fundamental research area in the computer vision community. These two topics are diffcult to handle some common challenges, such as occlusion, deformation, motion blur, and scale variation. The former contains heterogeneous object, interacting object, edge ambiguity, and shape complexity. And the latter suffers from difficulties in handling fast motion, out-of-view, and real-time processing. Combining the two problems of video object segmentation and tracking (VOST) can overcome their respective difficulties and improve their performance. VOST can be widely applied to many practical applications such as video summarization, high definition video compression, human computer interaction, and autonomous vehicles. This article aims to provide a comprehensive review of the state-of-the-art tracking methods, and classify these methods into different categories, and identify new trends. First, we provide a hierarchical categorization existing approaches, including unsupervised VOS, semi-supervised VOS, interactive VOS, weakly supervised VOS, and segmentation-based tracking methods. Second, we provide a detailed discussion and overview of the technical characteristics of the different methods. Third, we summarize the characteristics of the related video dataset, and provide a variety of evaluation metrics. Finally, we point out a set of interesting future works and draw our own conclusions.


title: On the first frequency of reinforced partially hinged plates

abstract: We consider a partially hinged rectangular plate and its normal modes. The dynamical properties of the plate are influenced by the spectrum of the associated eingenvalue problem. In order to improve the stability of the plate, it seems reasonable to place a certain amount of stiffening material in appropriate regions. If we look at the partial differential equation appearing in the model, this corresponds to insert a suitable weight coefficient inside the equation. A possible way to locate such regions is to study the eigenvalue problem associated to the aforementioned weighted equation. In the present paper we focus our attention essentially on the first eigenvalue and on its minimization in terms of the weight. We prove the existence of minimizing weights inside special classes and we try to describe them together with the corresponding eigenfunctions.


title: The methods of thermal field theory for degenerate quantum plasmas in astrophysical compact objects

abstract: In the study of degenerate plasmas contained within compact astrophysical objects, both special relativity and general relativity play important roles. After reviewing the existing treatment in the literature, here we employ the methods of relativistic thermal quantum field theory to compute the equation of states of degenerate matter for compact astrophysical objects such as the white dwarfs and the neutron stars. In particular, we compute the equation of states that include leading order corrections due to the finite temperature, the fine-structure constant as well as the effect of gravitational time dilation. We show that the fine-structure constant correction remains well-defined even in the non-relativistic regime in contrast to the existing treatment in the literature.


title: Terrestrial deuterium-to-hydrogen ratio in water in hyperactive comets

abstract: The D/H ratio in cometary water has been shown to vary between 1 and 3 times the Earth's oceans value, in both Oort cloud comets and Jupiter-family comets originating from the Kuiper belt. We present new sensitive spectroscopic observations of water isotopologues in the Jupiter-family comet 46P/Wirtanen carried out using the GREAT spectrometer aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The derived D/H ratio of $(1.61 \pm 0.65) \times 10^{-4}$ is the same as in the Earth's oceans. Although the statistics are limited, we show that interesting trends are already becoming apparent in the existing data. A clear anti-correlation is seen between the D/H ratio and the active fraction, defined as the ratio of the active surface area to the total nucleus surface. Comets with an active fraction above 0.5 typically have D/H ratios in water consistent with the terrestrial value. These hyperactive comets, such as 46P/Wirtanen, require an additional source of water vapor in their coma, explained by the presence of subliming icy grains expelled from the nucleus. The observed correlation may suggest that hyperactive comets belong to a population of ice-rich objects that formed just outside the snow line, or in the outermost regions of the solar nebula, from water thermally reprocessed in the inner disk that was transported outward during the early disk evolution. The observed anti-correlation between the active fraction and the nucleus size seems to argue against the first interpretation, as planetesimals near the snow line are expected to undergo rapid growth. Alternatively, isotopic properties of water outgassed from the nucleus and icy grains may be different due to fractionation effects at sublimation. In this case, all comets may share the same Earth-like D/H ratio in water, with profound implications for the early solar system and the origin of Earth's oceans.


title: The degeneration of critical point in Z(3)spin system. A proposal for QCD confinement-deconfinement phase transition in the color space

abstract: The problem of the phase transition of a Z(3) spin system is a complex issue. A numerical simulation in the framework of the mean field theory using the Metropolis algorithm reveals: (a) the existence of second order phase transition with a degenerated critical point which could be considered as a state of resonance (b) hysteresis phenomena which are accompanied with tricritical crossover around the end point of first order transition. Based on these results we propose a scenario about confinement-deconfinement phase transition in SU(3) symmetry of color QCD which, as it is known, has the same center with Z(3) spin system symmetry.


title: Monitoring dynamic collagen reorganization during skin stretching with fast polarization-resolved second harmonic generation imaging

abstract: The mechanical properties of biological tissues are strongly correlated to the specific distribution of their collagen fibers. Monitoring the dynamic reorganization of the collagen network during mechanical stretching is however a technical challenge because it requires mapping orientation of collagen fibers in a thick and deforming sample. In this work, a fast polarization-resolved SHG microscope is implemented to map collagen orientation during mechanical assays. This system is based on line-to-line switching of polarization using an electro-optical modulator and works in epidetection geometry. After proper calibration, it successfully highlights the collagen dynamic alignment along the traction direction in ex vivo murine skin dermis. This microstructure reorganization is quantified by the entropy of the collagen orientation distribution as a function of the stretch ratio. It exhibits a linear behavior, whose slope is measured with a good accuracy. This approach can be generalized to probe a variety of dynamic processes in thick tissues.


title: A strong order $3/4$ method for SDEs with discontinuous drift coefficient

abstract: In this paper we study strong approximation of the solution of a scalar stochastic differential equation (SDE) at the final time in the case when the drift coefficient may have discontinuities in space. Recently it has been shown in [M\"uller-Gronbach, T., and Yaroslavtseva, L., On the performance of the Euler-Maruyama scheme for SDEs with discontinuous drift coefficient, arXiv:1809.08423 (2018)] that for scalar SDEs with a piecewise Lipschitz drift coefficient and a Lipschitz diffusion coefficient that is non-zero at the discontinuity points of the drift coefficient the classical Euler-Maruyama scheme achieves an $L_p$-error rate of at least $1/2$ for all $p\in [1,\infty)$. Up to now this was the best $L_p$-error rate available in the literature for equations of that type. In the present paper we construct a method based on finitely many evaluations of the driving Brownian motion that even achieves an $L_p$-error rate of at least $3/4$ for all $p\in [1,\infty)$ under additional piecewise smoothness assumptions on the coefficients. To obtain this result we prove in particular that a quasi-Milstein scheme achieves an $L_p$-error rate of at least $3/4$ in the case of coefficients that are both Lipschitz continuous and piecewise differentiable with Lipschitz continuous derivatives, which is of interest in itself.


title: Credential Masquerading and OpenSSL Spy: Exploring ROS 2 using DDS security

abstract: The trend toward autonomous robot deployments is on an upward growth curve. These robots are undertaking new tasks and are being integrated into society. Examples of this trend are autonomous vehicles, humanoids, and eldercare. The movement from factory floors to streets and homes has also increased the number of vulnerabilities that adversaries can utilize. To improve security, Robot Operating System (ROS) 2 has standardized on using Data Distributed Services (DDS) as the messaging layer, which supports a security standard for protecting messages between parties with access control enforcement. DDS security is dependent on the OpenSSL and a security configuration file that specifies sensitive data location. DSS Security assumes that the underlining Operating System (OS) is secure and that the dependencies are consistent, but ongoing integrity checks are not performed. This paper looks at two vulnerabilities that we exploit using an OpenSSL spy process and a security property file manipulation. An overview of each exploit is provided with an evaluation of mitigation technologies that may be employed in client computers, servers, and other areas. Since, ROS 2 and DDS run in user space, these processes are prone to vulnerabilities. We provide recommendations about mitigation technology, as currently autonomous platforms are being deployed without safe-guards for on or off-line threats. The Trust Platform Module (TPM) is new to robotic systems, but the standard usage model does not provide risk mitigation above the OS layer for the types of attacks we discuss.


title: Symmetry Breaking Creates Electro-Momentum Coupling in Piezoelectric Metamaterials

abstract: Metamaterials posses microstructure designed to acquire properties not found in nature. An epitome in acoustics and solid mechanics is Willis coupling, which refers to the particle velocity-stress coupling, and of great significance since it controls mechanical waves. We here reveal new couplings, analogous to Willis coupling, when considering patterns of materials that mechanically interact with magnetic, electric or thermal fields. To this end, we develop a rigorous homogenization method for the effective properties of such responsive metamaterials. As an example, we apply the scheme to piezoelectric materials, and unveil coupling of the velocity and electric fields. Hence, Willis-like couplings in responsive metamaterials open new avenues for active wave control by modulation of external stimuli.


title: Theoretical study of phase transitions in Sb2S3, Bi2S3 and Sb2Se3 under compression

abstract: We report a theoretical study of Sb2S3, Sb2Se3 and Bi2S3 sesquichalcogenides at hydrostatic pressures up to 60 GPa. We explore the possibility that the R-3m, C2/m, C2/c and Im-3m phases observed in sesquichalcogenides with heavier cations, viz. Bi2Se3, Bi2Te3 and Sb2Te3, could also be formed in the former compounds as suggested by recent experiments. Our calculations show that the C2/m and C2/c phases are energetically unstable in any of the three compounds and over the entire range of pressures examined. In contrast, the disordered bcc-like Im-3m phase is energetically stable at high pressures, but only the Sb2Se3 phase does not show dynamical instabilities below 60 GPa. Our calculations show that the Pnma structure is the most energetically stable Sb2S3 and Bi2S3 phase at ambient pressure whereas, surprisingly, the R-3m phase of Sb2Se3 is predicted to have the lowest enthalpy energy at 0 GPa, in contradiction to experimental evidence. Our calculations show that both Pnma and R-3m phases are dynamical and mechanically stable at room pressure. These results suggest that the formation of the R-3m phase in Sb2Se3 could be feasible at close to ambient conditions. To aid its identification, we provide a theoretical crystal structure (lattice and atomic parameters) and complete infrared and Raman spectra for this phase.


title: Long time behavior of the NLS-Szeg{\"o} equation

abstract: We are interested in the influence of filtering the positive Fourier modes to the integrable non linear Schr{\"o}dinger equation. Equivalently, we want to study the effect of dispersion added to the cubic Szeg{\"o} equation, leading to the NLS-Szeg{\"o} equation on the circle $\mathbb{S}^1$\begin{equation*}i \partial\_t u + \epsilon^{\alpha}\partial\_x^2 u = \Pi(|u|^2 u),\qquad 0<\epsilon<1, \qquad \alpha\geq 0.\end{equation*}There are two sets of results in this paper. The first result concerns the long time Sobolev estimates for small data. The second set of results concerns the orbital stability of plane wave solutions. Some instability results are also obtained, leading to the wave turbulence phenomenon.


title: Magnonic Analogue Black/White Hole Horizon in Superfluid $^3$He-B: experiment

abstract: We provide experimental details of the first experiment made in zero temperature limit ($\sim$ 600\,$\mu$K) studying the magnonic black/white hole horizon analogue using absolutely pure physical system based on the spin superfluidity in superfluid $^3$He-B. We show that spin precession waves propagating on the background of the spin super-currents in a channel between two Bose-Einstein condensates of magnons in form of homogeneously precessing domains mimic the properties of the black/white horizon. Once the white hole horizon is formed and blocks the propagation of the spin-precession waves between two domains, we observed an amplification effect, i.e. when the energy of the spin precession waves reflected from the horizon is higher than the energy of the excited spin precession waves before horizon was formed. Moreover, the estimated temperature of the spontaneous Hawking radiation in this model system is about four orders of magnitude lower than the system's background temperature what makes it a promising tool to study the effect of spontaneous Hawking radiation.


title: Undecidability of future timeline-based planning over dense temporal domains

abstract: Planning is one of the most studied problems in computer science. In this paper, we consider the timeline-based approach, where the domain is modeled by a set of independent, but interacting, components, identified by a set of state variables, whose behavior over time (timelines) is governed by a set of temporal constraints (synchronization rules). Timeline-based planning in the dense-time setting has been recently shown to be undecidable in the general case, and undecidability relies on the high expressiveness of the trigger synchronization rules. In this paper, we strengthen the previous negative result by showing that undecidability already holds under the future semantics of the trigger rules which limits the comparison to temporal contexts in the future with respect to the trigger.


title: Radial basis function network using Lambert-Tsallis Wq function

abstract: The present work brings two applications of the Lambert-Tsallis Wq function in radial basis function networks (RBFN). Initially, a RBFN is used to discriminate between entangled and disentangled bipartite of qubit states. The kernel used is based on the Lambert-Tsallis Wq function for q = 2 and the quantum relative disentropy is used as distance measure between quantum states. Following, a RBFN with the same kernel is used to estimate the probability density function of a set of data samples.


title: Super-resolution of near-colliding point sources

abstract: We consider the problem of stable recovery of sparse signals of the form $$F(x)=\sum_{j=1}^d a_j\delta(x-x_j),\quad x_j\in\mathbb{R},\;a_j\in\mathbb{C}, $$ from their spectral measurements, known in a bandwidth $\Omega$ with absolute error not exceeding $\epsilon>0$. We consider the case when at most $p\le d$ nodes $\{x_j\}$ of $F$ form a cluster whose extent is smaller than the Rayleigh limit ${1\over\Omega}$, while the rest of the nodes are well separated. Provided that $\epsilon \lessapprox SRF^{-2p+1}$, where $SRF=(\Omega\Delta)^{-1}$ and $\Delta$ is the minimal separation between the nodes, we show that the minimax error rate for reconstruction of the cluster nodes is of order ${1\over\Omega}SRF^{2p-1}\epsilon$, while for recovering the corresponding amplitudes $\{a_j\}$ the rate is of the order $SRF^{2p-1}\epsilon$. Moreover, the corresponding minimax rates for the recovery of the non-clustered nodes and amplitudes are ${\epsilon\over\Omega}$ and $\epsilon$, respectively. These results suggest that stable super-resolution is possible in much more general situations than previously thought. Our numerical experiments show that the well-known Matrix Pencil method achieves the above accuracy bounds.


title: Continual Learning for Sentence Representations Using Conceptors

abstract: Distributed representations of sentences have become ubiquitous in natural language processing tasks. In this paper, we consider a continual learning scenario for sentence representations: Given a sequence of corpora, we aim to optimize the sentence encoder with respect to the new corpus while maintaining its accuracy on the old corpora. To address this problem, we propose to initialize sentence encoders with the help of corpus-independent features, and then sequentially update sentence encoders using Boolean operations of conceptor matrices to learn corpus-dependent features. We evaluate our approach on semantic textual similarity tasks and show that our proposed sentence encoder can continually learn features from new corpora while retaining its competence on previously encountered corpora.


title: The digest2 NEO Classification Code

abstract: We describe the digest2 software package, a fast, short-arc orbit classifier for small Solar System bodies. The digest2 algorithm has been serving the community for more than 13 years. The code provides a score, D2, which represents a pseudo-probability that a tracklet belongs to a given Solar System orbit type. digest2 is primarily used as a classifier for Near-Earth Object (NEO) candidates, to identify those to be prioritized for follow-up observation. We describe the historical development of digest2 and demonstrate its use on real and synthetic data. We find that digest2 can accurately and precisely distinguish NEOs from non- NEOs. At the time of detection, 14% of NEO tracklets and 98.5% of non-NEOs tracklets have D2 below the critical value of D2 = 65. 94% of our simulated NEOs achieved the maximum D2 = 100 and 99.6% of NEOs achieved $D2 \ge 65$ at least once during the simulated 10-year timeframe. We demonstrate that D2 varies as a function of time, rate of motion, magnitude and sky-plane location, and show that NEOs tend to have lower D2 at low Solar elongations close to the ecliptic. We use our findings to recommend future development directions for the digest2 code.


title: Semi-Lagrangian Exponential Integration with application to the rotating shallow water equations

abstract: In this paper we propose a novel way to integrate time-evolving partial differential equations that contain nonlinear advection and stiff linear operators, combining exponential integration techniques and semi-Lagrangian methods. The general formulation is built from the solution of an integration factor problem with respect to the problem written with a material derivative, so that the exponential integration scheme naturally incorporates the nonlinear advection. Semi-Lagrangian techniques are used to treat the dependence of the exponential integrator on the flow trajectories. The formulation is general, as many exponential integration techniques could be combined with different semi-Lagrangian methods. This formulation allows an accurate solution of the linear stiff operator, a property inherited by the exponential integration technique. It also provides a sufficiently accurate representation of the nonlinear advection, even with large time-step sizes, a property inherited by the semi-Lagrangian method. Aiming for application in weather and climate modeling, we discuss possible combinations of well established exponential integration techniques and state-of-the-art semi-Lagrangian methods used operationally in the application. We show experiments for the planar rotating shallow water equations. When compared to traditional exponential integration techniques, the experiments reveal that the coupling with semi-Lagrangian allows stabler integration with larger time-step sizes. From the application perspective, which already uses semi-Lagrangian methods, the exponential treatment could improve the solution of wave-dispersion when compared to semi-implicit schemes.


title: D\'ecompositions \`a la Steinberg sur une cat\'egorie additive

abstract: We give a description of simple functors taking finitely generated values, from a small additive category to the category of vector spaces over a field. This result is analogous to Steinberg's tensor product theorems in group representation theory. Our results rest on the notion of polynomial functor introduced by Eilenberg and Mac Lane. We give applications to representations of general linear groups or to finiteness properties of functor categories.


title: Learning Manipulation Skills Via Hierarchical Spatial Attention

abstract: Learning generalizable skills in robotic manipulation has long been challenging due to real-world sized observation and action spaces. One method for addressing this problem is attention focus -- the robot learns where to attend its sensors and irrelevant details are ignored. However, these methods have largely not caught on due to the difficulty of learning a good attention policy and the added partial observability induced by a narrowed window of focus. This article addresses the first issue by constraining gazes to a spatial hierarchy. For the second issue, we identify a case where the partial observability induced by attention does not prevent Q-learning from finding an optimal policy. We conclude with real-robot experiments on challenging pick-place tasks demonstrating the applicability of the approach.


title: Square-root nuclear norm penalized estimator for panel data models with approximately low-rank unobserved heterogeneity

abstract: This paper considers a nuclear norm penalized estimator for panel data models with interactive effects. The low-rank interactive effects can be an approximate model and the rank of the best approximation unknown and grow with sample size. The estimator is solution of a well-structured convex optimization problem and can be solved in polynomial-time. We derive rates of convergence, study the low-rank properties of the estimator, estimation of the rank and of annihilator matrices when the number of time periods grows with the sample size. Two-stage estimators can be asymptotically normal. None of the procedures require knowledge of the variance of the errors.


title: Cantor-Bernstein implies Excluded Middle

abstract: We prove in constructive logic that the statement of the Cantor-Bernstein theorem implies excluded middle. This establishes that the Cantor-Bernstein theorem can only be proven assuming the full power of classical logic. The key ingredient is a theorem of Mart\'in Escard\'o stating that quantification over a particular subset of the Cantor space $2^{\mathbb{N}}$, the so-called one-point compactification of $\mathbb{N}$, preserves decidable predicates.


title: Champagne Bubbles : Isolation and Characterization of amphiphilic macromolecules responsible for the stability of the collar at the Champagne / air interface

abstract: The effervescence and the stability of the ring of fine bubbles crowning the surface of a champagne glass, the "collar", constitute one of the hallmarks of Champagne. Defects in the stability of this collar are not well understood and account for a significant proportion of bottle return. This study aims to better understand the surface properties of champagne wine such that the foaming properties can be controlled more effectively. Early studies on Champagne foaming properties using the "Mosalux" measurement of the foam level formed by air flow in champagne through fritted glass pointed to a link between protein concentration and foam level [5] but no satisfactory correlation between protein content and foam stability was established. Later measurements were conducted with either ultra-filtrates or ultra-concentrates [4]. The authors demonstrated that macromolecule concentration was an essential parameter in the foam stability. The stability of bubbles is usually ascribed to the presence of an adsorption layer formed at the gas/liquid interface and its properties [3]. Thus surface properties of Champagne were analysed by ellipsometry and tensiometry. Measurements conducted on base wine, on ultra-filtered base wines and degassed champagne samples showed the presence of an adsorption layer formed at the air/champagne wine interface [6] and that adsorption layer being composed of macromolecules in a 104 to 105 molecular range [7]. Previous studies on champagne wine macromolecules had shown wine macromolecules to be mostly proteins and polysaccharides [8] with very little insight as to the chemical constitution. The present study describes the isolation and characterization of these macromolecules and their link with the adsorption layer.


title: Smoothing of the slowly extracted coasting beam from a synchrotron

abstract: Slow extraction of beam from synchrotrons or storage rings as required by many fixed target experiments is performed by controlled excitation and feeding of a structural lattice resonance. Due to the sensitive nature of this resonant extraction process, the temporal structure of the extracted beam is modulated by the minuscule current fluctuations present on the quadrupole magnet power supplies. Such a modulation lead to pile-ups in detectors and a significant reduction in accumulated event statistics. This contribution proposes and experimentally demonstrates that by an introduction of further modulation on quadrupole currents with a specific amplitude and frequency, the inherent power supply fluctuations are mitigated leading to a smoothening of the beam temporal structure. The slow extraction beam dynamics associated with this method are explained along with the operational results.

categories: physics.acc-ph physics.app-ph physics.ins-det physics.med-ph


title: Improved algorithms for left factorial residues

abstract: We present improved algorithms for computing the left factorial residues $!p=0!+1!+...+(p-1)! \!\mod p$. We use these algorithms for the calculation of the residues $!p\!\mod p$, for all primes $p$ up to $2^{40}$. Our results confirm that Kurepa's left factorial conjecture is still an open problem, as they show that there are no odd primes $p<2^{40}$ such that $p$ divides $!p$. Additionally, we confirm that there are no socialist primes $p$ with $5<p<2^{40}$.


title: Microwave to optical conversion with atoms on a superconducting chip

abstract: We describe a scheme to coherently convert a microwave photon of a superconducting co-planar waveguide resonator to an optical photon emitted into a well-defined temporal and spatial mode. The conversion is realized by a cold atomic ensemble trapped above the surface of the superconducting atom chip, near the antinode of the microwave cavity. The microwave photon couples to a strong Rydberg transition of the atoms that are also driven by a pair of laser fields with appropriate frequencies and wavevectors for an efficient wave-mixing process. With only few thousand atoms in an ensemble of moderate density, the microwave photon can be completely converted into an optical photon emitted with high probability into the phase matched direction and, e.g., fed into a fiber waveguide. This scheme operates in a free-space configuration, without requiring strong coupling of the atoms to a resonant optical cavity.


title: Semi-implicit methods for the dynamics of elastic sheets

abstract: Recent applications (e.g. active gels and self-assembly of elastic sheets) motivate the need to efficiently simulate the dynamics of thin elastic sheets. We present semi-implicit time stepping algorithms to improve the time step constraints that arise in explicit methods while avoiding much of the complexity of fully-implicit approaches. For a triangular lattice discretization with stretching and bending springs, our semi-implicit approach involves discrete Laplacian and biharmonic operators, and is stable for all time steps in the case of overdamped dynamics. For a more general finite-difference formulation that can allow for general elastic constants, we use the analogous approach on a square grid, and find that the largest stable time step is two to three orders of magnitude greater than for an explicit scheme. For a model problem with a radial traveling wave form of the reference metric, we find transitions from quasi-periodic to chaotic dynamics as the sheet thickness is reduced, wave amplitude is increased, and damping constant is reduced.


title: Spreading dynamics of reactive surfactants driven by Marangoni convection

abstract: We consider the spreading dynamics of some insoluble surface-active species along an aqueous interface. The model includes both diffusion, Marangoni convection and first-order reaction kinetics. An exact solution of the nonlinear transport equations is derived in the regime of large Schmidt number, where viscous effects are dominant. We demonstrate that the variance of the surfactant distribution increases linearly with time, providing an unambiguous definition for the enhanced diffusion coefficient observed in the experiments. The model thus presents new insight regarding the actuation of camphor grains at the water-air interface.


title: Re-enterant efficiency of phototaxis in Chlamydomonas reinhardtii cells

abstract: Phototaxis is one of the most fundamental stimulus-response behaviors in biology wherein motile micro-organisms sense light gradients to swim towards the light source. Apart from single cell survival and growth, it plays a major role at the global scale of aquatic ecosystem and bio-reactors. We study photoaxis of single celled algae Chalmydomonas reinhardtii as a function of cell number density and light stimulus using high spatio-temporal video microscopy. Surprisingly, the phototactic efficiency has a minimum at a well-defined number density, for a given light gradient, above which the phototaxis behaviour of collection of cells can even exceed the performance obtainable from single isolated cells. We show that the origin of enhancement of performance above the critical concentration lies in the slowing down of the cells which enables them to sense light more effectively. We also show that this steady state phenomenology is well captured by a modelling the phototactic response as a density dependent torque acting on an active Brownian particle.


title: Visualizing the decision-making process in deep neural decision forest

abstract: Deep neural decision forest (NDF) achieved remarkable performance on various vision tasks via combining decision tree and deep representation learning. In this work, we first trace the decision-making process of this model and visualize saliency maps to understand which portion of the input influence it more for both classification and regression problems. We then apply NDF on a multi-task coordinate regression problem and demonstrate the distribution of routing probabilities, which is vital for interpreting NDF yet not shown for regression problems. The pre-trained model and code for visualization will be available at https://github.com/Nicholasli1995/VisualizingNDF


title: Absence of inter-layer tunnel coupling of $K$-valley electrons in bilayer MoS$_{2}$

abstract: In Bernal stacked bilayer graphene interlayer coupling significantly affects the electronic bandstructure compared to monolayer graphene. Here we present magnetotransport experiments on high-quality $n$-doped bilayer MoS$_{2}$. By measuring the evolution of the Landau levels as a function of electron density and applied magnetic field we are able to investigate the occupation of conduction band states, the interlayer coupling in pristine bilayer MoS$_{2}$, and how these effects are governed by electron-electron interactions. We find that the two layers of the bilayer MoS$_{2}$ behave as two independent electronic systems where a two-fold Landau level's degeneracy is observed for each MoS$_{2}$ layer. At the onset of the population of the bottom MoS$_{2}$ layer we observe a large negative compressibility caused by the exchange interaction. These observations, enabled by the high electronic quality of our samples, demonstrate weak interlayer tunnel coupling but strong interlayer electrostatic coupling in pristine bilayer MoS$_{2}$. The conclusions from the experiments may be relevant also to other transition metal dichalcogenide materials.


title: Linear Bounded Composition of Tree-Walking Tree Transducers: Linear Size Increase and Complexity

abstract: Compositions of tree-walking tree transducers form a hierarchy with respect to the number of transducers in the composition. As main technical result it is proved that any such composition can be realized as a linear bounded composition, which means that the sizes of the intermediate results can be chosen to be at most linear in the size of the output tree. This has consequences for the expressiveness and complexity of the translations in the hierarchy. First, if the computed translation is a function of linear size increase, i.e., the size of the output tree is at most linear in the size of the input tree, then it can be realized by just one, deterministic, tree-walking tree transducer. For compositions of deterministic transducers it is decidable whether or not the translation is of linear size increase. Second, every composition of deterministic transducers can be computed in deterministic linear time on a RAM and in deterministic linear space on a Turing machine, measured in the sum of the sizes of the input and output tree. Similarly, every composition of nondeterministic transducers can be computed in simultaneous polynomial time and linear space on a nondeterministic Turing machine. Their output tree languages are deterministic context-sensitive, i.e., can be recognized in deterministic linear space on a Turing machine. The membership problem for compositions of nondeterministic translations is nondeterministic polynomial time and deterministic linear space. The membership problem for the composition of a nondeterministic and a deterministic tree-walking tree translation (for a nondeterministic IO macro tree translation) is log-space reducible to a context-free language, whereas the membership problem for the composition of a deterministic and a nondeterministic tree-walking tree translation (for a nondeterministic OI macro tree translation) is possibly NP-complete.


title: Optimal Recovery of Precision Matrix for Mahalanobis Distance from High Dimensional Noisy Observations in Manifold Learning

abstract: Motivated by establishing theoretical foundations for various manifold learning algorithms, we study the problem of Mahalanobis distance (MD), and the associated precision matrix, estimation from high-dimensional noisy data. By relying on recent transformative results in covariance matrix estimation, we demonstrate the sensitivity of \MD~and the associated precision matrix to measurement noise, determining the exact asymptotic signal-to-noise ratio at which MD fails, and quantifying its performance otherwise. In addition, for an appropriate loss function, we propose an asymptotically optimal shrinker, which is shown to be beneficial over the classical implementation of the MD, both analytically and in simulations. The result is extended to the manifold setup, where the nonlinear interaction between curvature and high-dimensional noise is taken care of. The developed solution is applied to study a multiscale reduction problem in the dynamical system analysis.


title: Plasma Wakefield Accelerator Research 2019 - 2040: A community-driven UK roadmap compiled by the Plasma Wakefield Accelerator Steering Committee (PWASC)

abstract: The acceleration gradients generated in a laser- or beam-driven plasma wakefield accelerator are typically three orders of magnitude greater than those produced by a conventional accelerator, and hence plasma accelerators can open a route to a new generation of very compact machines. In addition, plasma-based accelerators can generate beams with unique properties, such as tens of kiloamp peak currents, attosecond bunch duration, ultrahigh brightness and intrinsic particle beam-laser pulse synchronization. In this roadmap we review the status of plasma accelerator research in the UK. We outline potential applications, describe the research and development required to enable those applications, and discuss synergies with related areas of research. We also set-out the resources required to realise these ambitions and provide a timeline for advances in the key areas.

categories: physics.acc-ph physics.app-ph physics.optics physics.plasm-ph


title: Crystal plasticity analysis of deformation anisotropy of lamellar TiAl alloy: 3D microstructure-based modelling and in-situ micro-compression

abstract: Detailed microstructure characterisation and in-situ micropillar compression were coupled with crystal plasticity-based finite element modelling (CP-FEM) to study the micro-mechanisms of plastic anisotropy in lamellar TiAl alloys. The consideration of microstructure in both simulation and in-situ experiments enables in-depth understanding of micro-mechanisms responsible for the highly anisotropic deformation response of TiAl on the intra-lamella and inter-lamella scales. This study focuses on two specific configurations of {\gamma}/{\alpha}_2 lamellar microstructure with the {\gamma}/{\alpha}_2 interfaces being aligned {25}^o and {55}^o to the loading direction. Microstructure-based CP-FEM shows that longituginal slip of super and ordinary dislocations are most responsible for the plastic anisotropy in the {25}^o micropillar while the anisotropy of the {55}^o micropillar is due to longitudinal superdislocations and longitudinal twins. In addition, transversal superdislocations were more active, making the deformation in the {25}^o micropillar less localised than that in the {55}^o micropillar. Moreover, the CP-FEM model successfully predicted substantial build-up of internal stresses at {\gamma}/{\alpha}_2 interfaces, which is believed to be detrimental to the ductility in TiAl. However, as evidenced by the model, the detrimental internal stresses can be significantly relieved by the activation of transverse deformation twinning, suggesting that the ductility of TiAl can be improved by promoting transverse twins.


title: Quandle Cocycle Quivers

abstract: We incorporate quandle cocycle information into the quandle coloring quivers we defined in arXiv:1807.10465 to define weighted directed graph-valued invariants of oriented links we call \textit{quandle cocycle quivers}. This construction turns the quandle cocycle invariant into a small category, yielding a categorification of the quandle cocycleinvariant. From these graphs we define several new link invariants including a 2-variable polynomial which specializes to the usual quandle cocycle invariant. Examples and computations are provided.


title: Mean Force Kinetic Theory: a Convergent Kinetic Theory for Weakly and Strongly Coupled Plasmas

abstract: A new closure of the BBGKY hierarchy is developed, which results in a convergent kinetic equation that provides a rigorous extension of plasma kinetic theory into the regime of strong Coulomb coupling. The approach is based on a single expansion parameter which enforces that the exact equilibrium limit is maintained at all orders. Because the expansion parameter does not explicitly depend on the range or the strength of the interaction potential, the resulting kinetic theory does not suffer from the typical divergences at short and long length scales encountered when applying the standard kinetic equations to Coulomb interactions. The approach demonstrates that particles effectively interact via the potential of mean force and that the range of this force determines the size of the collision volume. When applied to a plasma, the collision operator is shown to be related to the effective potential theory [Baalrud and Daligault, Phys. Rev. Lett 110, 235001 (2013)]. The relationship between this and previous kinetic theories is discussed.


title: A New Class of Scaling Matrices for Scaled Trust Region Algorithms

abstract: A new class of affine scaling matrices for the interior point Newton-type methods is considered to solve the nonlinear systems with simple bounds. We review the essential properties of a scaling matrix and consider several well-known scaling matrices proposed in the literature. We define a new scaling matrix that is the convex combination of these matrices. The proposed scaling matrix inherits those interesting properties of the individual matrices and satisfies additional desired requirements. The numerical experiments demonstrate the superiority of the new scaling matrix in solving several important test problems.


title: Mechanisms to Splay-Bend Nematic Phases

abstract: While twist-bend nematic phases have been extensively studied, the experimental observation of two dimensional, oscillating splay-bend phases is recent. We consider two theoretical models that have been used to explain the formation of twist-bend phases -- flexoelectricity and bond orientational order -- as mechanisms to induce splay-bend phases. Flexoelectricity is a viable mechanism, and splay and bend flexoelectric couplings can lead to splay-bend phases with different modulations. We show that while bond orientational order circumvents the need for higher order terms in the free energy, the important role of nematic symmetry and phase chirality rules it out as a basic mechanism.


title: Analyzing the benefits of communication channels between deep learning models

abstract: As artificial intelligence systems spread to more diverse and larger tasks in many domains, the machine learning algorithms, and in particular the deep learning models and the databases required to train them are getting bigger themselves. Some algorithms do allow for some scaling of large computations by leveraging data parallelism. However, they often require a large amount of data to be exchanged in order to ensure the shared knowledge throughout the compute nodes is accurate. In this work, the effect of different levels of communications between deep learning models is studied, in particular how it affects performance. The first approach studied looks at decentralizing the numerous computations that are done in parallel in training procedures such as synchronous and asynchronous stochastic gradient descent. In this setting, a simplified communication that consists of exchanging low bandwidth outputs between compute nodes can be beneficial. In the following chapter, the communication protocol is slightly modified to further include training instructions. Indeed, this is studied in a simplified setup where a pre-trained model, analogous to a teacher, can customize a randomly initialized model's training procedure to accelerate learning. Finally, a communication channel where two deep learning models can exchange a purposefully crafted language is explored while allowing for different ways of optimizing that language.


title: Risk Convergence of Centered Kernel Ridge Regression with Large Dimensional Data

abstract: This paper carries out a large dimensional analysis of a variation of kernel ridge regression that we call \emph{centered kernel ridge regression} (CKRR), also known in the literature as kernel ridge regression with offset. This modified technique is obtained by accounting for the bias in the regression problem resulting in the old kernel ridge regression but with \emph{centered} kernels. The analysis is carried out under the assumption that the data is drawn from a Gaussian distribution and heavily relies on tools from random matrix theory (RMT). Under the regime in which the data dimension and the training size grow infinitely large with fixed ratio and under some mild assumptions controlling the data statistics, we show that both the empirical and the prediction risks converge to a deterministic quantities that describe in closed form fashion the performance of CKRR in terms of the data statistics and dimensions. Inspired by this theoretical result, we subsequently build a consistent estimator of the prediction risk based on the training data which allows to optimally tune the design parameters. A key insight of the proposed analysis is the fact that asymptotically a large class of kernels achieve the same minimum prediction risk. This insight is validated with both synthetic and real data.