title: The Escape Fraction of Ionizing Photons During the Epoch of Reionization: observability with the Square Kilometre Array

abstract: One of the most important parameters in characterizing the Epoch of Reionization, the escape fraction of ionizing photons, $f_\mathrm{esc}$, remains unconstrained both observationally and theoretically. With recent work highlighting the impact of galaxy-scale feedback on the instantaneous value of $f_\mathrm{esc}$, it is important to develop a model in which reionization is self-consistently coupled to galaxy evolution. In this work, we present such a model and explore how physically motivated functional forms of $f_\mathrm{esc}$ affect the evolution of ionized hydrogen within the intergalactic medium. Using the $21$cm power spectrum evolution, we investigate the likelihood of observationally distinguishing between a constant $f_\mathrm{esc}$ and other models that depend upon different forms of galaxy feedback. We find that changing the underlying connection between $f_\mathrm{esc}$ and galaxy feedback drastically alters the large-scale $21$cm power. The upcoming Square Kilometre Array Low Frequency instrument possesses the sensitivity to differentiate between our models at a fixed optical depth, requiring only $200$ hours of integration time focused on redshifts $z = 7.5-8.5$. Generalizing these results to account for a varying optical depth will require multiple $800$ hour observations spanning redshifts $z = 7-10$. This presents an exciting opportunity to observationally constrain one of the most elusive parameters during the Epoch of Reionization.


title: Efficient conversion of chemical energy into mechanical work by Hsp70 chaperones

abstract: Hsp70 molecular chaperones are abundant ATP-dependent nanomachines that actively reshape non-native, misfolded proteins and assist a wide variety of essential cellular processes. Here we combine complementary computational/theoretical approaches to elucidate the structural and thermodynamic details of the chaperone-induced expansion of a substrate protein, with a particular emphasis on the critical role played by ATP hydrolysis. We first determine the conformational free-energy cost of the substrate expansion due to the binding of multiple chaperones using coarse-grained molecular simulations. We then exploit this result to implement a non-equilibrium rate model which estimates the degree of expansion as a function of the free energy provided by ATP hydrolysis. Our results are in quantitative agreement with recent single-molecule FRET experiments and highlight the stark non-equilibrium nature of the process, showing that Hsp70s are optimized to convert effectively chemical energy into mechanical work close to physiological conditions.


title: Low-latency job scheduling with preemption for the development of deep learning

abstract: One significant challenge in the job scheduling of computing clusters for the development of deep learning algorithms is the efficient scheduling of trial-and-error (TE) job, the type of job in which the users seek to conduct small-scale experiments while monitoring their processes. Unfortunately, the existing job schedulers to date do not feature well-balanced scheduling for the mixture of TE jobs and best-effort (BE) jobs, or they can handle the mixture in limited situations at most. To fill in this niche, we propose an algorithm that can significantly reduce the latency of TE jobs in versatile situations without greatly elongating the slowdown of the BE jobs. Our algorithm efficiently schedules both TE and BE jobs by selectively preempting the BE jobs that can be, when the time comes, resumed without much delay. In our simulation study with synthetic and real workloads, we were able to reduce the 95th percentile of the slowdown rates for the TE jobs in the standard FIFO strategy by 96.6%, while compromising the median of the BE slowdown rates by only 18.0% and the 95th percentile by only 23.9%.


title: On the linear independence constraint qualification in disjunctive programming

abstract: Mathematical programs with disjunctive constraints (MPDCs for short) cover several different problem classes from nonlinear optimization including complementarity-, vanishing-, cardinality-, and switching-constrained optimization problems. In this paper, we introduce an abstract but reasonable version of the prominent linear independence constraint qualification which applies to MPDCs. Afterwards, we derive first- and second-order optimality conditions for MPDCs under validity of this constraint qualification based on so-called strongly stationary points. Finally, we apply our findings to some popular classes of disjunctive programs and compare the obtained results to those ones available in the literature. Particularly, new second-order optimality conditions for mathematical programs with switching constraints are by-products of our approach.


title: Restructuring Conversations using Discourse Relations for Zero-shot Abstractive Dialogue Summarization

abstract: Dialogue summarization is a challenging problem due to the informal and unstructured nature of conversational data. Recent advances in abstractive summarization have been focused on data-hungry neural models and adapting these models to a new domain requires the availability of domain-specific manually annotated corpus created by linguistic experts. We propose a zero-shot abstractive dialogue summarization method that uses discourse relations to provide structure to conversations, and then uses an out-of-the-box document summarization model to create final summaries. Experiments on the AMI and ICSI meeting corpus, with document summarization models like PGN and BART, shows that our method improves the ROGUE score by up to 3 points, and even performs competitively against other state-of-the-art methods.


title: TDHF Theory and Its Extensions for the Multinucleon Transfer Reaction: a Mini Review

abstract: Time-dependent Hartree-Fock (TDHF) theory has been a powerful tool in describing a variety of complex nuclear dynamics microscopically without empirical parameters. In this contribution, recent advances in nuclear dynamics studies with TDHF and its extensions are briefly reviewed, along the line with the study of multinucleon transfer (MNT) reactions. The latter lies at the core of this Research Topic, whose application for production of extremely neutron-rich nuclei has been extensively discussed in recent years. Having in mind the ongoing theoretical developments, it is envisaged how microscopic theories may contribute to the future MNT study.


title: Towards Commercializing Vanadium Dioxide Films: Investigation of the Impact of Different Interface on the Deterioration Process for Largely Extended Service Life

abstract: Long term stability is the most pressing issue that impedes commercialization of Vanadium Dioxide (VO2) based functional films, which show a gradual loss of relative phase transition performance, especially in humid conditions when serving as smart windows. Here, we investigated the impact of different interface on the deterioration process of VO2 films and proposed a novel encapsulation structure for largely extended service life. Hydrophobic and stable hafnium dioxide (HfO2) layers have been incorporated with VO2 films for encapsulated surfaces and cross sections. With modified thickness and structure of HfO2 layers, the degradation process of VO2 can be effectively suppressed. The proposed films can retain stable phase transition performances under high relative humidity (90%) and temperature (60 Celsius) over 100 days, which is equal to about 16 years in the real environment. Improving the stability of VO2 materials is a necessary step towards commercializing production of high performance films for long term use.


title: Echo State Networks: analysis, training and predictive control

abstract: The goal of this paper is to investigate the theoretical properties, the training algorithm, and the predictive control applications of Echo State Networks (ESNs), a particular kind of Recurrent Neural Networks. First, a condition guaranteeing incremetal global asymptotic stability is devised. Then, a modified training algorithm allowing for dimensionality reduction of ESNs is presented. Eventually, a model predictive controller is designed to solve the tracking problem, relying on ESNs as the model of the system. Numerical results concerning the predictive control of a nonlinear process for pH neutralization confirm the effectiveness of the proposed algorithms for the identification, dimensionality reduction, and the control design for ESNs.


title: The nature of Els\"{a}sser variables in compressible MHD

abstract: The Els\"{a}sser variables are often used in studies of plasma turbulence, in helping differentiate between MHD waves propagating parallel or anti-parallel to the main magnetic field. While for pure Alfv\'en waves in a homogeneous plasma the method is strictly valid, we show that compressible, magnetoacoustic waves are in general described by both Els\"{a}sser variables. Furthermore, in a compressible and inhomogeneous plasma, the pure MHD waves (Alfv\'en, fast and slow) are no longer normal modes, but waves become linearly coupled or display mixed properties of Alfv\'en and magnetoacoustic nature. These waves are necessarily described by both Els\"{a}sser variables and therefore the Els\"{a}sser formalism cannot be used to strictly separate parallel and anti-parallel propagating waves. Nevertheless, even in an inhomogeneous plasma, for a highly Alfv\'enic wave the Els\"{a}sser variable corresponding to the propagation direction appears still dominating. We suggest that for Alfv\'enic waves, the relative amplitude of Els\"{a}sser variables depends on the local degree of inhomogeneity and other plasma and wave properties. This finding has implications for turbulence studies in inhomogeneous and compressible plasmas, such as the solar corona and solar wind.


title: A Primer on Weyl Semimetals: Down the Discovery of Topological Phases

abstract: Recently discovered Weyl semimetals (WSM) have found special place in topological condensed matter studies for they represent first example of massless Weyl fermions found in condensed matter systems. A WSM shows gapless bulk energy spectra with Dirac-like point degeneracies, famously called Weyl nodes, which carry with themselves well defined chiralities and topologically protected chiral charges. One finds the Berry curvature of the Bloch bands to become singular, like in a magnetic monopole, at these Weyl nodes. Moreover, these systems feature topological surface states in the form of open Fermi arcs. In this review, we undergo a concise journey from graphene based Dirac physics to Weyl semimetals: the underlying Hamiltonians, their basic features and their unique response to external electric and magnetic fields in order to provide a basic walk-through of how the Weyl physics unfolded with time starting from the discovery of Graphene.


title: Index and nullity of proper biharmonic maps in spheres

abstract: In recent years, the study of the bienergy functional has attracted the attention of a large community of researchers, but there are not many examples where the second variation of this functional has been thoroughly studied. We shall focus on this problem and, in particular, we shall compute the exact index and nullity of some known examples of proper biharmonic maps. Moreover, we shall analyse a case where the domain is not compact. More precisely, we shall prove that a large family of proper biharmonic maps $\varphi:\mathbb{R} \to \mathbb{S}^2$ is strictly stable with respect to compactly supported variations. In general, the computations involved in this type of problems are very long. For this reason, we shall also define and apply to specific examples a suitable notion of index and nullity with respect to equivariant variations.


title: A General Framework for Prediction in Time Series Models

abstract: In this paper we propose a general framework to analyze prediction in time series models and show how a wide class of popular time series models satisfies this framework. We postulate a set of high-level assumptions, and formally verify these assumptions for the aforementioned time series models. Our framework coincides with that of Beutner et al. (2019, arXiv:1710.00643) who establish the validity of conditional confidence intervals for predictions made in this framework. The current paper therefore complements the results in Beutner et al. (2019, arXiv:1710.00643) by providing practically relevant applications of their theory.


title: Topological phase transitions in glassy quantum matter

abstract: Amorphous systems have rapidly gained promise as novel platforms for topological matter. In this work we establish a scaling theory of amorphous topological phase transitions driven by the density of lattice points in two dimensions. By carrying out a finite-size scaling analysis of topological invariants averaged over discrete and continuum random geometries, we discover unique critical properties of Chern and $\mathbb{Z}_2$ glass transitions. Even for short-range hopping models the Chern glass phase may persist down to the fundamental lower bound given by the classical percolation threshold. While the topological indices accurately satisfy the postulated one-parameter scaling, they do not generally flow to the closest integer value in the thermodynamic limit. Furthermore, the value of the critical exponent describing the diverging localization length varies continuously along the phase boundary and is not fixed by the symmetry class of the Hamiltonian. We conclude that the critical behaviour of amorphous topological systems exhibit characteristic features not observed in disordered systems, motivating a wealth of new research directions.


title: Calculation of eigenvalues by Greens-functions and the Lippmann-Schwinger equation

abstract: We calculate eigenvalues of one-dimensional quantum-systems by the exact numerical solution of the Lippmann-Schwinger equation, analogous to the scattering problem. To illustrate our method, we treat elementary problems: the harmonic and quartic oscillator, a symmetric and a skew double-well potential, and potentials with finite and infinite depth. Furthermore, we show how our method can be used for eigenvalue-engineering.


title: Modularity-in-Design of Dynamical Network Systems: Retrofit Control Approach

abstract: In this paper, we develop a modular design method of decentralized controllers for linear dynamical network systems, where multiple subcontroller designers aim at individually regulating their local control performance with accessibility only to their respective subsystem models. First, we derive a constrained version of the Youla parameterization that characterizes all retrofit controllers for a single subcontroller, defined as an add-on type subcontroller that manages a subsystem. The resultant feedback system is kept robustly stable for any variation in the neighboring subsystems, other than the subsystem of interest, provided that the original system is stable prior to implementing the retrofit control. Subsequently, we find out a unique internal structure of the retrofit controllers, assuming that the interaction input signal from the neighboring subsystems is measurable. Furthermore, we show that the simultaneous implementation of multiple retrofit controllers, designed by individual subcontroller designers, can improve the upper bound of the overall control performance. Finally, the practical significance of the method is demonstrated via an illustrative example of frequency regulation using the IEEE 68-bus power system model.


title: EasyLabel: A Semi-Automatic Pixel-wise Object Annotation Tool for Creating Robotic RGB-D Datasets

abstract: Developing robot perception systems for recognizing objects in the real-world requires computer vision algorithms to be carefully scrutinized with respect to the expected operating domain. This demands large quantities of ground truth data to rigorously evaluate the performance of algorithms. This paper presents the EasyLabel tool for easily acquiring high quality ground truth annotation of objects at the pixel-level in densely cluttered scenes. In a semi-automatic process, complex scenes are incrementally built and EasyLabel exploits depth change to extract precise object masks at each step. We use this tool to generate the Object Cluttered Indoor Dataset (OCID) that captures diverse settings of objects, background, context, sensor to scene distance, viewpoint angle and lighting conditions. OCID is used to perform a systematic comparison of existing object segmentation methods. The baseline comparison supports the need for pixel- and object-wise annotation to progress robot vision towards realistic applications. This insight reveals the usefulness of EasyLabel and OCID to better understand the challenges that robots face in the real-world. Copyright 20XX IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.


title: Sharp gap edges in dense planetary rings: An axisymmetric diffusion model

abstract: One of the most intriguing facets of Saturn's rings are the sharp edges of gaps in the rings where the surface density abruptly drops to zero. This is despite of the fact that the range over which a moon transfers angular momentum onto the ring material is much larger. Recent UVIS-scans of the edges of the Encke and Keeler gap show that this drop occurs over a range approximately equal to the rings' thickness. Borderies et al. (1982, 1989) show that this striking feature is likely related to the local reversal of the usually outward directed viscous transport of angular momentum in strongly perturbed regions. In this article we revise the Borderies et al. (1989) model using a granular flow model to define the shear and bulk viscosities and incorporate the angular momentum flux reversal effect into the axisymmetric diffusion model we developed for gaps in dense planetary rings (Gr\"atz et al. 2018). Finally, we apply our model to the Encke and Keeler division in order to estimate the shear and bulk viscosities in the vicinity of both gaps.


title: Analytic investigation of rotating holographic superconductors

abstract: In this paper we have investigated, in the probe limit, $s$-wave holographic superconductors in rotating $AdS_{3+1}$ spacetime using the matching method as well as the St{\"u}rm-Liouville eigenvalue approach. We have calculated the critical temperature using the matching technique in such a setting and our results are in agreement with previously reported results obtained using the St{\"u}rm-Liouville approach. We have then obtained the condensation operators using both analytical methods. The results obtained by both these techniques share the same features as found numerically. We observe that the rotation parameter of the black hole affects the critical temperature and the condensation operator in a non-trivial way.


title: Literature Survey on Finding Influential Communities in Large Scale Networks

abstract: Community or modular structure is considered to be a significant property of large scale real-world graphs such as social or information networks. Detecting influential clusters or communities in these graphs is a problem of considerable interest as it often accounts for the functionality of the system. We aim to provide a thorough exposition of the topic, including the main elements of the problem, a brief introduction of the existing research for both disjoint and overlapping community search, the idea of influential communities, its implications and the current state of the art and finally provide some insight on possible directions for future research.


title: A Generalization of Linear Positive Systems with Applications to Nonlinear Systems: Invariant Sets and the Poincar\'{e}-Bendixson Property

abstract: The dynamics of linear positive systems map the positive orthant to itself. In other words, it maps a set of vectors with zero sign variations to itself. This raises the following question: what linear systems map the set of vectors with $k$ sign variations to itself? We address this question using tools from the theory of cooperative dynamical systems and the theory of totally positive matrices. This yields a generalization of positive linear systems called $k$-positive linear systems, that reduces to positive systems for $k=1$. We describe applications of this new type of systems to the analysis of nonlinear dynamical systems. In particular, we show that such systems admit certain explicit invariant sets, and for the case $k=2$ establish the Poincar\'{e}-Bendixson property for any bounded trajectory.


title: Proximity exchange effects in MoSe$_2$ and WSe$_2$ heterostructures with CrI$_3$: twist angle, layer, and gate dependence

abstract: Proximity effects in two-dimensional (2D) van der Waals heterostructures offer controllable ways to tailor the electronic band structure of adjacent materials. Exchange proximity in particular is important for making materials magnetic without hosting magnetic ions. Such synthetic magnets could be used for studying magnetotransport in high-mobility 2D materials, or magneto-optics in highly absorptive nominally nonmagnetic semiconductors. Using first-principles calculations, we show that the proximity exchange in monolayer MoSe$_2$ and WSe$_2$ due to ferromagnetic monolayer CrI$_3$ can be tuned (even qualitatively) by twisting and gating. Remarkably, the proximity exchange remains the same when using antiferromagnetic CrI$_3$ bilayer, paving the way for optical and electrical detection of layered antiferromagnets. Interestingly, the proximity exchange is opposite to the exchange of the adjacent antiferromagnetic layer. Finally, we show that the exchange proximity is confined to the layer adjacent to CrI$_3$, and that adding a separating hBN barrier drastically reduces the proximity effect. We complement our it ab initio results with tight-binding modeling and solve the Bethe-Salpeter equation to provide experimentally verifiable optical signatures (in the exciton spectra) of the proximity exchange effects.


title: Minimum description length as an objective function for non-negative matrix factorization

abstract: Non-negative matrix factorization (NMF) is a dimensionality reduction technique which tends to produce a sparse representation of data. Commonly, the error between the actual and recreated matrices is used as an objective function, but this method may not produce the type of representation we desire as it allows for the complexity of the model to grow, constrained only by the size of the subspace and the non-negativity requirement. If additional constraints, such as sparsity, are imposed the question of parameter selection becomes critical. Instead of adding sparsity constraints in an ad-hoc manner we propose a novel objective function created by using the principle of minimum description length (MDL). Our formulation, MDL-NMF, automatically trades off between the complexity and accuracy of the model using a principled approach with little parameter selection or the need for domain expertise. We demonstrate our model works effectively on three heterogeneous data-sets and on a range of semi-synthetic data showing the broad applicability of our method.


title: An ancient and a primordial collisional family as the main sources of X-type asteroids of the inner Main Belt

abstract: The near-Earth asteroid population suggests the existence of an inner Main Belt source of asteroids that belongs to the spectroscopic X-complex and has moderate albedos. The identification of such a source has been lacking so far. We argue that the most probable source is one or more collisional asteroid families that escaped discovery up to now. We apply a novel method to search for asteroid families in the inner Main Belt population of asteroids belonging to the X-complex with moderate albedo. Instead of searching for asteroid clusters in orbital elements space, which could be severely dispersed when older than some billions of years, our method looks for correlations between the orbital semimajor axis and the inverse size of asteroids. This correlation is the signature of members of collisional families, which drifted from a common centre under the effect of the Yarkovsky thermal effect. We identify two previously unknown families in the inner Main Belt among the moderate-albedo X-complex asteroids. One of them, whose lowest numbered asteroid is (161) Author, is ~3 Gyrs-old, whereas the second one, whose lowest numbered object is (689) Zita, can be as old as the Solar System. Members of this latter family have orbital eccentricities and inclinations that spread them over the entire inner Main Belt, which is an indication that this family could be primordial, i.e. it formed before the giant planet orbital instability. The vast majority of moderate-albedo X-complex asteroids of the inner-Main Belt are genetically related, as they can be included into few asteroid families. Only nine X-complex asteroids with moderate albedo of the inner Main Belt cannot be included in asteroid families. We suggest that these bodies formed by direct accretion of the solids in the protoplanetary disk, and are thus surviving planetesimals.


title: CCD and DSLR maxima of RR Lyrae stars in 2016 and 2017

abstract: We present here maximum timings of RR Lyrae stars observed in 2016 and 2017 by VSOLJ member Kenji Hirosawa. The calculation of the time of maximum from the measurements provided by the observer have been made in the frame of the GEOS RR Lyr survey.


title: Riemannian optimization with a preconditioning scheme on the generalized Stiefel manifold

abstract: Optimization problems on the generalized Stiefel manifold (and products of it) are prevalent across science and engineering. For example, in computational science they arise in symmetric (generalized) eigenvalue problems, in nonlinear eigenvalue problems, and in electronic structures computations, to name a few problems. In statistics and machine learning, they arise, for example, in various dimensionality reduction techniques such as canonical correlation analysis. In deep learning, regularization and improved stability can be obtained by constraining some layers to have parameter matrices that belong to the Stiefel manifold. Solving problems on the generalized Stiefel manifold can be approached via the tools of Riemannian optimization. However, using the standard geometric components for the generalized Stiefel manifold has two possible shortcomings: computing some of the geometric components can be too expensive and convergence can be rather slow in certain cases. Both shortcomings can be addressed using a technique called Riemannian preconditioning, which amounts to using geometric components derived by a precoditioner that defines a Riemannian metric on the constraint manifold. In this paper we develop the geometric components required to perform Riemannian optimization on the generalized Stiefel manifold equipped with a non-standard metric, and illustrate theoretically and numerically the use of those components and the effect of Riemannian preconditioning for solving optimization problems on the generalized Stiefel manifold.


title: Search for heavy charged long-lived particles in the ATLAS detector in 36.1 fb$^{-1}$ of proton-proton collision data at $\sqrt{s} = 13$ TeV

abstract: A search for heavy charged long-lived particles is performed using a data sample of 36.1 fb$^{-1}$ of proton-proton collisions at $\sqrt{s} = 13$ TeV collected by the ATLAS experiment at the Large Hadron Collider. The search is based on observables related to ionization energy loss and time of flight, which are sensitive to the velocity of heavy charged particles traveling significantly slower than the speed of light. Multiple search strategies for a wide range of lifetimes, corresponding to path lengths of a few meters, are defined as model-independently as possible, by referencing several representative physics cases that yield long-lived particles within supersymmetric models, such as gluinos/squarks ($R$-hadrons), charginos and staus. No significant deviations from the expected Standard Model background are observed. Upper limits at 95% confidence level are provided on the production cross sections of long-lived $R$-hadrons as well as directly pair-produced staus and charginos. These results translate into lower limits on the masses of long-lived gluino, sbottom and stop $R$-hadrons, as well as staus and charginos of 2000 GeV, 1250 GeV, 1340 GeV, 430 GeV and 1090 GeV, respectively.


title: A Universal Algorithm for Variational Inequalities Adaptive to Smoothness and Noise

abstract: We consider variational inequalities coming from monotone operators, a setting that includes convex minimization and convex-concave saddle-point problems. We assume an access to potentially noisy unbiased values of the monotone operators and assess convergence through a compatible gap function which corresponds to the standard optimality criteria in the aforementioned subcases. We present a universal algorithm for these inequalities based on the Mirror-Prox algorithm. Concretely, our algorithm simultaneously achieves the optimal rates for the smooth/non-smooth, and noisy/noiseless settings. This is done without any prior knowledge of these properties, and in the general set-up of arbitrary norms and compatible Bregman divergences. For convex minimization and convex-concave saddle-point problems, this leads to new adaptive algorithms. Our method relies on a novel yet simple adaptive choice of the step-size, which can be seen as the appropriate extension of AdaGrad to handle constrained problems.


title: Diophantine Equations Involving the Euler Totient Function

abstract: We deal with various Diophantine equations involving the Euler totient function and various sequences of numbers, including factorials, powers, and Fibonacci sequences.


title: Exploiting locality in high-dimensional factorial hidden Markov models

abstract: We propose algorithms for approximate filtering and smoothing in high-dimensional Factorial hidden Markov models. The approximation involves discarding, in a principled way, likelihood factors according to a notion of locality in a factor graph associated with the emission distribution. This allows the exponential-in-dimension cost of exact filtering and smoothing to be avoided. We prove that the approximation accuracy, measured in a local total variation norm, is "dimension-free" in the sense that as the overall dimension of the model increases the error bounds we derive do not necessarily degrade. A key step in the analysis is to quantify the error introduced by localizing the likelihood function in a Bayes' rule update. The factorial structure of the likelihood function which we exploit arises naturally when data have known spatial or network structure. We demonstrate the new algorithms on synthetic examples and a London Underground passenger flow problem, where the factor graph is effectively given by the train network.


title: Dipolar Switching of Charge-Injection Barriers at Electrode/Semiconductor Interfaces as a Mechanism for Water-Induced Instabilities of Organic Devices

abstract: An electrode contact-related mechanism for the operational instability of organic electronic devices is proposed and confirmed via observation of a water-induced change in charge-injection barrier eights at the electrode/organic-semiconductor interfaces. Water molecules in air penetrate into the organic crystal via diffusion, and an external electric field orients the electric dipole of the water molecules at the electrode surfaces, leading to dipolar switching of the charge-injection barrier height. As a result of the switching, current-voltage curves of two-terminal Au-rubrene-Au devices change from symmetric to asymmetric, showing diode-like rectification and reversible switching of the diode polarity. The device shows the highest current switching ratio of 267 for the switching voltage of 3 V, corresponding to an electrode work function change of >144 meV. The mechanism proposed herein will be important especially for short-channel organic devices, which are indispensable for applications such as organic integrated circuits.


title: Rigidity theorems of Lagrangian submanifolds in the homogeneous nearly K\"ahler $\mathbb{S}^6(1)$

abstract: In this paper, we study Lagrangian submanifolds of the homogeneous nearly K\"ahler $6$-dimensional unit sphere $\mathbb{S}^6(1)$. As the main result, we derive a Simons' type integral inequality in terms of the second fundamental form for compact Lagrangian submanifolds of $\mathbb{S}^6(1)$. Moreover, we show that the equality sign occurs if and only if the Lagrangian submanifold is either the totally geodesic $\mathbb{S}^3(1)$ or the Dillen-Verstraelen-Vrancken's Berger sphere $S^3$ discribed in J Math Soc Japan, 42: 565-584, 1990.


title: Agent-Based Simulation Modelling for Reflecting on Consequences of Digital Mental Health

abstract: The premise of this working paper is based around agent-based simulation models and how to go about creating them from given incomplete information. Agent-based simulations are stochastic simulations that revolve around groups of agents that each have their own characteristics and can make decisions. Such simulations can be used to emulate real life situations and to create hypothetical situations without the need for real-world testing prior. Here we describe the development of an agent-based simulation model for studying future digital mental health scenarios. An incomplete conceptual model has been used as the basis for this development. To define differences in responses to stimuli we employed fuzzy decision making logic. The model has been implemented but not been used for structured experimentation yet. This is planned as our next step.


title: Three-body repulsive forces among identical bosons in one dimension

abstract: I consider non-relativistic bosons interacting via pairwise potentials with infinite scattering length and supporting no two-body bound states. To lowest order in effective field theory, these conditions lead to non-interacting bosons, since the coupling constant of the Lieb-Liniger model vanishes identically in this limit. Since any realistic pairwise interaction is not a mere delta function, the non-interacting picture is an idealisation indicating that the effect of interactions is weaker than in the case of off-resonant potentials. I show that the leading order correction to the ground state energy for more than two bosons is accurately described by the lowest order three-body force in effective field theory that arises due to the off-shell structure of the two-body interaction. For natural two-body interactions with a short-distance repulsive core and an attractive tail, the emergent three-body interaction is repulsive and, therefore, three bosons do not form any bound states. This situation is analogous to the two-dimensional repulsive Bose gas, when treated using the lowest-order contact interaction, where the scattering amplitude exhibits an unphysical Landau pole. The avoidance of this state in the three-boson problem proceeds in a way that parallels the two-dimensional case. These results pave the way for the experimental realisation of one-dimensional Bose gases with pure three-body interactions using ultracold atomic gases.


title: On the Convergence of Projected-Gradient Methods with Low-Rank Projections for Smooth Convex Minimization over Trace-Norm Balls and Related Problems

abstract: Smooth convex minimization over the unit trace-norm ball is an important optimization problem in machine learning, signal processing, statistics and other fields, that underlies many tasks in which one wishes to recover a low-rank matrix given certain measurements. While first-order methods for convex optimization enjoy optimal convergence rates, they require in worst-case to compute a full-rank SVD on each iteration, in order to compute the projection onto the trace-norm ball. These full-rank SVD computations however prohibit the application of such methods to large problems. A simple and natural heuristic to reduce the computational cost is to approximate the projection using only a low-rank SVD. This raises the question if, and under what conditions, this simple heuristic can indeed result in provable convergence to the optimal solution. In this paper we show that any optimal solution is a center of a Euclid. ball inside-which the projected-gradient mapping admits rank that is at most the multiplicity of the largest singular value of the gradient vector. Moreover, the radius of the ball scales with the spectral gap of this gradient vector. We show how this readily implies the local convergence (i.e., from a "warm-start" initialization) of standard first-order methods, using only low-rank SVD computations. We also quantify the effect of "over-parameterization", i.e., using SVD computations with higher rank, on the radius of this ball, showing it can increase dramatically with moderately larger rank. We extend our results also to the setting of optimization with trace-norm regularization and optimization over bounded-trace positive semidefinite matrices. Our theoretical investigation is supported by concrete empirical evidence that demonstrates the \textit{correct} convergence of first-order methods with low-rank projections on real-world datasets.


title: The Dynamics of Magnetism in Fe-Cr Alloys with Cr Clustering

abstract: The dynamics of magnetic moments in iron-chromium alloys with different levels of Cr clustering show unusual features resulting from the fact that even in a perfect body-centred cubic structure, magnetic moments experience geometric magnetic frustration resembling that of a spin glass. Due to the long range exchange coupling and configuration randomness, magnetic moments of Cr solutes remain non-collinear at all temperatures. To characterise magnetic properties of Fe-Cr alloys, we explore the temperature dependence of magnetisation, susceptibility, Curie temperature and spin-spin correlations with spatial resolution. The static and dynamic magnetic properties are correlated with the microstructure of Fe-Cr, where magnetisation and susceptibility are determined by the size of Cr precipitates at nominal Cr concentrations. The Curie temperature is always maximised when the solute concentration of Cr in the $\alpha$ phase is close to 5 to 6 at.\%, and the susceptibility of Fe atoms is always enhanced at the boundary between a precipitate and solid solution. Interaction between Cr and Fe stimulates magnetic disorder, lowering the effective Curie temperature. Dynamic simulation of evolution of magnetic correlations shows that the spin-spin relaxation time in Fe-Cr alloys is in the 20 to 40 ps range.


title: Basal Slip in Laves Phases: the Synchroshear Dislocation

abstract: Two different mechanisms have been reported in previous ab initio studies to describe basal slip in complex intermetallic Laves phases: synchroshear and undulating slip. To date, no clear answer has been given on which is the energetically favourable mechanism and whether either of them could effectively propagate as a dislocation. Using classical atomistic simulations supported by ab initio calculations, the present work removes the ambiguity and shows that the two mechanisms are, in fact, identical. Furthermore, we establish synchroshear as the mechanism for propagating dislocations within the basal plane in Laves phases.


title: Berry curvature unravelled by the Nernst effect in Mn$_3$Ge

abstract: The discovery of topological quantum materials represents a striking innovation in modern condensed matter physics with remarkable fundamental and technological implications. Their classification has been recently extended to topological Weyl semimetals, i.e., solid state systems which exhibit the elusive Weyl fermions as low-energy excitations. Here we show that the Nernst effect can be exploited as a sensitive probe for determining key parameters of the Weyl physics, applying it to the non-collinear antiferromagnet Mn$_3$Ge. This compound exhibits anomalous thermoelectric transport due to enhanced Berry curvature from Weyl points located extremely close to the Fermi level. We establish from our data a direct measure of the Berry curvature at the Fermi level and, using a minimal model of a Weyl semimetal, extract for the first time the Weyl point energy and their distance in momentum-space.


title: eMBB-URLLC Resource Slicing: A Risk-Sensitive Approach

abstract: Ultra Reliable Low Latency Communication (URLLC) is a 5G New Radio (NR) application that requires strict reliability and latency. URLLC traffic is usually scheduled on top of the ongoing enhanced Mobile Broadband (eMBB) transmissions (i.e., puncturing the current eMBB transmission) and cannot be queued due to its hard latency requirements. In this letter, we propose a risk-sensitive based formulation to allocate resources to the incoming URLLC traffic while minimizing the risk of the eMBB transmission (i.e., protecting the eMBB users with low data rate) and ensuring URLLC reliability. Specifically, the Conditional Value at Risk (CVaR) is introduced as a risk measure for eMBB transmission. Moreover, the reliability constraint of URLLC is formulated as a chance constraint and relaxed based on Markov's inequality. We decompose the formulated problem into two subproblems in order to transform it into a convex form and then alternatively solve them until convergence. Simulation results show that the proposed approach allocates resources to the incoming URLLC traffic efficiently while satisfying the reliability of both eMBB and URLLC.


title: Division of an angle into equal parts and construction of regular polygons by multi-fold origami

abstract: This article analyses geometric constructions by origami when up to $n$ simultaneous folds may be done at each step. It shows that any arbitrary angle can be $m$-sected if the largest prime factor of $m$ is $p\le n+2$. Also, the regular $m$-gon can be constructed if the largest prime factor of $\phi(m)$ is $q\le n+2$, where $\phi$ is Euler's totient function.


title: Energetics of magnetic transients in a solar active region plage

abstract: Densely packed coronal loops are rooted in plages in the vicinity of active regions on the Sun. The photospheric magnetic features underlying these plages are patches of mostly unidirectional magnetic field extending several arcsec on the solar surface. We aim to explore the transient nature of the magnetic field, its mixed-polarity characteristics, and the associated energetics in plages using high resolution observations and simulations. We used photospheric Fe I 6173 {\AA} spectropolarimetric observations of a decaying active region obtained from the SST. These data were inverted to retrieve the photospheric magnetic field underlying the plage as identified in the EUV emission maps obtained from the SDO/AIA. To obtain better insight into the evolution of plages, we performed 3D radiation magnetohydrodynamic simulations of magnetoconvection using the MURaM code. The observations show transient magnetic flux emergence and cancellation events within the plage on timescales of a few 100 s and on spatial scales comparable to granules. These transient events occur at the footpoints of active region plage loops. In one case the coronal response at the footpoints of these loops is clearly associated with the underlying transient. The simulations also reveal similar magnetic flux emergence and cancellation events that extend to even smaller spatial and temporal scales. Individual simulated transient events transfer an energy flux in excess of 1 MW m$^{-2}$ through the photosphere. We suggest that the magnetic transients could play an important role in the energetics of active region plage. Both in observations and simulations, the opposite-polarity magnetic field brought up by transient flux emergence cancels with the surrounding plage field. Magnetic reconnection associated with such transient events likely conduits magnetic energy to power the overlying chromosphere and coronal loops.


title: Long-range dispersion forces between molecules subject to ultra-short optical pulses from ab initio calculations

abstract: London-van der Waals dispersion forces are a fundamental component of condensed matter systems, biological processes, and self-assembly. In this letter we propose a method to calculate the C6 coefficients that characterize dispersion forces in the non-retarded regime for molecules in a coherent superposition of excited states. Several ultrafast femtosecond pump-probe schemes are investigated. We apply the method to LiH molecules and show that their London-van der Waals interaction can change dramatically after the interaction with the pump pulse. The pulse modulates the C6 coefficients, and the interplay between polarization, orientation of the molecules, and the dipole fields gives rise to a rich variety of combinations.


title: Numerical computation and new output bounds for time-limited balanced truncation of discrete-time systems

abstract: In this paper, balancing based model order reduction (MOR) for large-scale linear discrete-time time-invariant systems in prescribed finite time intervals is studied. The first main topic is the development of error bounds regarding the approximated output vector within the time limits. The influence of different components in the established bounds will be highlighted. After that, the second part of the article proposes strategies that enable an efficient numerical execution of time-limited balanced truncation for large-scale systems. Numerical experiments illustrate the performance of the proposed techniques.


title: Unconventional magnetism in the 4d$^{4}$ based ($S=1$) honeycomb system Ag$_{3}$LiRu$_{2}$O$_{6}$

abstract: We have investigated the thermodynamic and local magnetic properties of the Mott insulating system Ag$_{3}$LiRu$_{2}$O$_{6}$ containing Ru$^{4+}$ (4$d$$^{4}$) for novel magnetism. The material crystallizes in a monoclinic $C2/m$ structure with RuO$_{6}$ octahedra forming an edge-shared two-dimensional honeycomb lattice with limited stacking order along the $c$-direction. The large negative Curie-Weiss temperature ($\theta_{CW}$ = -57 K) suggests antiferromagnetic interactions among Ru$^{4+}$ ions though magnetic susceptibility and heat capacity show no indication of magnetic long-range order down to 1.8 K and 0.4 K, respectively. $^{7}$Li nuclear magnetic resonance (NMR) shift follows the bulk susceptibility between 120-300 K and levels off below 120 K. Together with a power-law behavior in the temperature dependent spin-lattice relaxation rate between 0.2 and 2 K, it suggest dynamic spin correlations with gapless excitations. Electronic structure calculations suggest an $S = 1$ description of the Ru-moments and the possible importance of further neighbour interactions as also bi-quadratic and ring-exchange terms in determining the magnetic properties. Analysis of our $\mu$SR data indicates spin freezing below 5 K but the spins remain on the borderline between static and dynamic magnetism even at 20 mK.


title: DVOLVER: Efficient Pareto-Optimal Neural Network Architecture Search

abstract: Automatic search of neural network architectures is a standing research topic. In addition to the fact that it presents a faster alternative to hand-designed architectures, it can improve their efficiency and for instance generate Convolutional Neural Networks (CNN) adapted for mobile devices. In this paper, we present a multi-objective neural architecture search method to find a family of CNN models with the best accuracy and computational resources tradeoffs, in a search space inspired by the state-of-the-art findings in neural search. Our work, called Dvolver, evolves a population of architectures and iteratively improves an approximation of the optimal Pareto front. Applying Dvolver on the model accuracy and on the number of floating points operations as objective functions, we are able to find, in only 2.5 days, a set of competitive mobile models on ImageNet. Amongst these models one architecture has the same Top-1 accuracy on ImageNet as NASNet-A mobile with 8% less floating point operations and another one has a Top-1 accuracy of 75.28% on ImageNet exceeding by 0.28% the best MobileNetV2 model for the same computational resources.


title: Sampling Free TDOA Localization in Millimeter Wave Networks

abstract: Time difference of arrival (TDOA) is a widely used technique for localizing a radio transmitter from the difference in signal arrival times at multiple receivers. For TDOA to work, the individual receivers must estimate the respective signal arrival times precisely, which requires sampling the signal at least double the rate of its highest frequency content, commonly known as the Nyquist rate. Such sampling is less practical for the millimeter wave band comprising of frequencies in 30-300 GHz range. In this paper, we propose an energy detection architecture for accurately estimating the time of arrival from a single picosecond Gaussian pulse, which enables TDOA localization without sampling at the receiver. We derive the closed form expression of the estimated time of arrival and validate it via simulation. We demonstrate that the proposed sampling-free TDOA can localize millimeter wave transmitters as accurately as the conventional sampling-based TDOA.


title: On distributive join-semilattices

abstract: Motivated by Gentzen disjunction elimination rule in his Natural Deduction calculus and reading inequalities with meet in a natural way, we conceive a notion of distributivity for join-semilattices. We prove that it is equivalent to a notion present in the literature. In the way, we prove that those notions are linearly ordered. We finally consider the notion of distributivity in join-semilattices with arrow, that is, the algebraic structure corresponding to the disjunction-conditional fragment of intuitionistic logic.


title: Probing restarting activity in hard X-ray selected giant radio galaxies

abstract: With their sizes larger than 0.7 Mpc, Giant Radio Galaxies (GRGs) are the largest individual objects in the Universe. To date, the reason why they reach such enormous extensions is still unclear. One of the proposed scenarios suggests that they are the result of multiple episodes of jet activity. Cross-correlating the INTEGRAL+Swift AGN population with radio catalogues (NVSS, FIRST, SUMSS), we found that 22% of the sources are GRG (a factor four higher than those selected from radio catalogues). Remarkably, 80% of the sample shows signs of restarting radio activity. The X-ray properties are consistent with this scenario, the sources being in a high-accretion, high-luminosity state with respect to the previous activity responsible for the radio lobes.


title: A desingularization of Kontsevich's compactification of twisted cubics in $V_5$

abstract: By definition, the del Pezzo $3$-fold $V_5$ is the intersection of $\mathrm{Gr}(2,5)$ with three hyperplanes in $\mathbb{P}^9$ under the Pl\"ucker embedding. Rational curves in $V_5$ have been studied in various contents of Fano geometry. In this paper, we propose an explicit birational relation of the Kontsevich and Simpson compactifications of twisted cubic curves in $V_5$. As a direct corollary, we obtain a desingularized model of Kontsevich compactification which induces the intersection cohomology group of Kontsevich's space.


title: Early Recognition of Sepsis with Gaussian Process Temporal Convolutional Networks and Dynamic Time Warping

abstract: Sepsis is a life-threatening host response to infection associated with high mortality, morbidity, and health costs. Its management is highly time-sensitive since each hour of delayed treatment increases mortality due to irreversible organ damage. Meanwhile, despite decades of clinical research, robust biomarkers for sepsis are missing. Therefore, detecting sepsis early by utilizing the affluence of high-resolution intensive care records has become a challenging machine learning problem. Recent advances in deep learning and data mining promise to deliver a powerful set of tools to efficiently address this task. This empirical study proposes two novel approaches for the early detection of sepsis: a deep learning model and a lazy learner based on time series distances. Our deep learning model employs a temporal convolutional network that is embedded in a Multi-task Gaussian Process Adapter framework, making it directly applicable to irregularly-spaced time series data. Our lazy learner, by contrast, is an ensemble approach that employs dynamic time warping. We frame the timely detection of sepsis as a supervised time series classification task. For this, we derive the most recent sepsis definition in an hourly resolution to provide the first fully accessible early sepsis detection environment. Seven hours before sepsis onset, our methods improve area under the precision--recall curve from 0.25 to 0.35/0.40 over the state of the art. This demonstrates that they are well-suited for detecting sepsis in the crucial earlier stages when management is most effective.


title: Explaining the decrease in ISM lithium at super-solar metallicities in the solar vicinity

abstract: We propose here that the lithium decrease at super-solar metallicities observed in high resolution spectroscopic surveys can be explained by the interplay of mixed populations, coming from the inner regions of the Milky Way disc. The lower lithium content of these stars is a consequence of inside-out disc formation, plus radial migration. In this framework, local stars with super-solar metallicities would have migrated to the solar vicinity and depleted their original lithium during their travel time. To arrive to such a result, we took advantage of the AMBRE catalog of lithium abundances combined with chemical evolution models which take into account the contribution to the lithium enrichment by different nucleosynthetic sources. A large proportion of migrated stars can explain the observed lower lithium abundance at super-solar metallicities. We stress that nowadays, there is no stellar model able to predict Li-depletion for such super-solar metallicity stars, and the Solar Li-depletion has to be assumed. In addition, it currently exists no solid quantitative estimate of the proportion of migrated stars in the Solar neighborhood and their travel time. Our results illustrate how important it is to properly include radial migration when comparing chemical evolution models to observations, and that in this case, the lithium decrease at larger metallicities does not necessarily imply that stellar yields have to be modified, contrary to previous claims in literature.


title: Nucleosynthetic Yields from Neutron Stars Accreting in Binary Common Envelopes

abstract: Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario for neutron star binaries, the system might pass through a phase where a neutron star spirals into the envelope of its giant star companion. These phases lead to mass accretion onto the neutron star. Accretion onto these common-envelope-phase neutron stars can eject matter that has undergone burning near to the neutron star surface. This paper presents nucleosynthetic yields of this ejected matter, using population synthesis models to study the importance of these nucleosynthetic yields in a galactic chemical evolution context. Depending on the extreme conditions in temperature and density found in the accreted material, both proton-rich and neutron-rich nucleosynthesis can be obtained, with efficient production of neutron rich isotopes of low Z material at the most extreme conditions, and proton rich isotopes, again at low Z, in lower density models. Final yields are found to be extremely sensitive to the physical modeling of the accretion phase. We show that neutron stars accreting in binary common envelopes might be a new relevant site for galactic chemical evolution, and therefore more comprehensive studies are needed to better constrain nucleosynthesis in these objects.


title: Design of materials properties and device performance in memristive systems

abstract: Future development of the modern nanoelectronics and its flagships internet of things and artificial intelligence as well as many related applications is largely associated with memristive elements. This technology offers a broad spectrum of functionalities, however, it follows predominantly a phenomenological approach and crucial challenge/limit for further development remains variability and lack of fundamental materials' design strategy. Here we demonstrate the vital importance of materials' purity for determining memristors' functionalities, showing that part per million foreign elements significantly change the performance. By appropriate choice of chemistry and amount of doping material we can selectively enhance desired operation mode. We highlight how dopant dependent structure and charge/potential distribution in the space charge layers and the cell capacitance determine the device kinetics and functions. We evidence for first time experimentally the relation between materials properties and switching/neuromorphic performance, thus providing rules and directions for a rational design of memristive devices.


title: Identification, Secrecy, Template, and Privacy-Leakage of Biometric Identification System Under Noisy Enrollment

abstract: In this study, we investigate fundamental trade-off among identification, secrecy, template, and privacy-leakage rates in biometric identification systems. Ignatenko and Willems (2015) studied this system assuming that the channel in the enrollment process of the system is noiseless and they did not consider the template rate. In the enrollment process, however, it is highly considered that noise occurs when bio-data is scanned. In this paper, we impose a noisy channel in the enrollment process and characterize the capacity region of the rate tuples. The capacity region is proved by a novel technique via two auxiliary random variables, which has never been seen in previous studies. As special cases, the obtained result shows that the characterization reduces to the one given by Ignatenko and Willems (2015) where the enrollment channel is noiseless and there is no constraint on the template rate, and it also coincides with the result derived by G\"unl\"u and Kramer (2018) where there is only one individual.


title: On the geometry of random polytopes

abstract: We present a simple proof to a fact recently established in [5]: let $\xi$ be a symmetric random variable that has variance $1$, let $\Gamma=(\xi_{ij})$ be an $N \times n$ random matrix whose entries are independent copies of $\xi$, and set $X_1,...,X_N$ to be the rows of $\Gamma$. Then under minimal assumptions on $\xi$ and as long as $N \geq c_1n$, $$ c_2 \bigl(B_\infty^n \cap \sqrt{\log(eN/n)} B_2^n \bigr) \subset {\rm absconv}(X_1,...,X_N) $$ with high probability.


title: What shapes a galaxy? - Unraveling the role of mass, environment and star formation in forming galactic structure

abstract: We investigate the dependence of galaxy structure on a variety of galactic and environmental parameters for ~500,000 galaxies at z<0.2, taken from the Sloan Digital Sky Survey data release 7 (SDSS-DR7). We utilise bulge-to-total stellar mass ratio, (B/T)_*, as the primary indicator of galactic structure, which circumvents issues of morphological dependence on waveband. We rank galaxy and environmental parameters in terms of how predictive they are of galaxy structure, using an artificial neural network approach. We find that distance from the star forming main sequence (Delta_SFR), followed by stellar mass (M_*), are the most closely connected parameters to (B/T)_*, and are significantly more predictive of galaxy structure than global star formation rate (SFR), or any environmental metric considered (for both central and satellite galaxies). Additionally, we make a detailed comparison to the Illustris hydrodynamical simulation and the LGalaxies semi-analytic model. In both simulations, we find a significant lack of bulge-dominated galaxies at a fixed stellar mass, compared to the SDSS. This result highlights a potentially serious problem in contemporary models of galaxy evolution.


title: Nonlinear electromagnetic response and Higgs mode excitation in BCS superconductors with impurities

abstract: We reveal that due to the presence of disorder oscillations of the order parameter amplitude called the Higgs mode can be effectively excited by the external electromagnetic radiation in usual BCS superconductors. This mechanism works for superconductors with both isotropic s-wave and anisotropic, such as d-wave, pairings. The non-linear response in the presence of impurities is captured by the quasiclassical formalism. We demonstrate that analytical solutions of the Eilenberger equation with impurity collision integral and external field drive coincide with the exact summation of ladder impurity diagrams. Using the developed formalism we show that resonant third-harmonic signal observed in recent experiments is naturally explained by the excitation of Higgs mode mediated by impurity scattering.


title: Against the Tyranny of Pure States in Quantum Theory

abstract: We argue that the notion of pure sate within Standard Quantum Mechanics is presently applied within the specialized literature in relation to two mutually inconsistent definitions. While the first (operational purity) provides a basis-dependent definition which makes reference to the certain prediction of measurement outcomes, the latter (trace-invariant purity) provides a purely abstract invariant definition which lacks operational content. In this work we derive a theorem which exposes the serious inconsistencies existent within these two incompatible definitions of purity.

categories: quant-ph math-ph math.MP physics.hist-ph


title: Expressive Power of Broadcast Consensus Protocols

abstract: Population protocols are a formal model of computation by identical, anonymous mobile agents interacting in pairs. Their computational power is rather limited: Angluin et al. have shown that they can only compute the predicates over $\mathbb{N}^k$ expressible in Presburger arithmetic. For this reason, several extensions of the model have been proposed, including the addition of devices called cover-time services, absence detectors, and clocks. All these extensions increase the expressive power to the class of predicates over $\mathbb{N}^k$ lying in the complexity class NL when the input is given in unary. However, these devices are difficult to implement, since they require that an agent atomically receives messages from all other agents in a population of unknown size; moreover, the agent must know that they have all been received. Inspired by the work of the verification community on Emerson and Namjoshi's broadcast protocols, we show that NL-power is also achieved by extending population protocols with reliable broadcasts, a simpler, standard communication primitive.


title: Rama: Controller Fault Tolerance in Software-Defined Networking Made Practical

abstract: In Software-Defined Networking (SDN), network applications use the logically centralized network view provided by the controller to remotely orchestrate the network switches. To avoid the controller being a single point of failure, traditional fault-tolerance techniques are employed to guarantee availability, a fundamental requirement in production environments. Unfortunately, these techniques fall short of ensuring correct network behaviour under controller failures. The problem of these techniques is that they deal with only part of the problem: guaranteeing that application and controller state remains consistent between replicas. However, in an SDN the switches maintain hard state that must also be handled consistently. Fault-tolerant SDN must therefore include switch state into the problem. A recently proposed fault-tolerant controller platform, Ravana, solves this problem by extending fault-tolerant SDN control with mechanisms that guarantee control messages to be processed transactionally and exactly once, at both the controllers and the switches. These guarantees are given even in the face of controller and switch crashes. The elegance of this solution comes at a cost. Ravana requires switches to be modified and OpenFlow to be extended with hitherto unforeseen additions to the protocol. In face of this challenge we propose Rama, a fault-tolerant SDN controller platform that offers the same strong guarantees as Ravana without requiring modifications to switches or to the OpenFlow protocol. Experiments with our prototype implementation show the additional overhead to be modest, making Rama the first fault-tolerant SDN solution that can be immediately deployable.


title: Interactively shaping robot behaviour with unlabeled human instructions

abstract: In this paper, we propose a framework that enables a human teacher to shape a robot behaviour by interactively providing it with unlabeled instructions. We ground the meaning of instruction signals in the task-learning process, and use them simultaneously for guiding the latter. We implement our framework as a modular architecture, named TICS (Task-Instruction-Contingency-Shaping) that combines different information sources: a predefined reward function, human evaluative feedback and unlabeled instructions. This approach provides a novel perspective for robotic task learning that lies between Reinforcement Learning and Supervised Learning paradigms. We evaluate our framework both in simulation and with a real robot. The experimental results demonstrate the effectiveness of our framework in accelerating the task-learning process and in reducing the number of required teaching signals.


title: Generation of very high-order high purity Gaussian modes via spatial light modulation

abstract: We experimentally demonstrate the conversion of a fundamental $\text{TEM}_{00}$ laser mode at 1064\,nm to higher order Hermite-Gaussian modes (HG) of arbitrary order via a commercially available liquid crystal Spatial Light Modulator (SLM). We particularly studied the $\text{HG}_{5,5}/\text{HG}_{10,10}/\text{HG}_{15,15}$ modes. A two-mirror plano-spherical cavity filters the higher-order modes spatially. We analyze the cleaned modes via a three-mirror diagnosis cavity and measure a mode purity of 96/93/78\% and a conversion efficiency of 6.6\%/3.7\%/1.7\% respectively. The generated high-purity Hermite-Gaussian modes can be employed for the mitigation of mirror thermal noise in optical cavities for both optical clocks and gravitational wave (GW) detectors. HG modes are then converted into high order LG modes which can be of particular interest in cold atom physics.


title: Optimized quasisymmetric stellarators are consistent with the Garren-Boozer construction

abstract: Most quasisymmetric stellarators to date have been designed by numerically optimizing the plasma boundary shape to minimize symmetry-breaking Fourier modes of the magnetic field strength $B$. At high aspect ratio, a faster approach is to directly construct the plasma shape from the equations of quasisymmetry near the magnetic axis derived by Garren and Boozer [Phys Fluids B 3, 2805 (1991)]. Here we show that the core shape and rotational transform of many optimization-based configurations can be accurately described by this direct-construction approach. This consistency supports use of the near-axis construction as an accurate analytical model for modern stellarator configurations.


title: On spatially irregular ordinary differential equations and a pathwise volatility modelling framework

abstract: This thesis develops a new framework for modelling price processes in finance, such as an equity price or foreign exchange rate. This can be related to the conventional Ito calculus-based framework through the time integral of a price's squared volatility, or `cumulative variance'. In the new framework, corresponding processes are strictly increasing, solve random ordinary differential equations (ODEs), and are composed with geometric Brownian motion. The new framework has no dependence on stochastic calculus, so processes can be studied on a pathwise basis using probability-free ODE techniques and functional analysis. The ODEs considered depend on continuous driving functions which are `spatially irregular', meaning they need not have any spatial regularity properties such as Holder continuity. They are however strictly increasing in time, thus temporally asymmetric. When sensible initial values are chosen, initial value problem (IVP) solutions are also strictly increasing, and the solution set of such IVPs is shown to contain all differentiable bijections on the non-negative reals. This enables the modelling of any non-negative volatility path which is not zero over intervals, via the time derivative of solutions. Despite this generality, new well-posedness results establish the uniqueness of solutions going forwards in time. Motivation to explore this framework comes from its connection with a time-changed Heston volatility model. The framework shows how Heston price processes can converge to a generalisation of the NIG Levy process, and reveals a deeper relationship between integrated CIR processes and the IG process. Within this framework, a `Riemann-Liouville-Heston' martingale model is defined which generalises these relationships to fractional counterparts. This model's implied volatilities are simulated, and exhibit features characteristic of leading volatility models.


title: Monte Carlo event generators for high energy particle physics event simulation

abstract: Monte Carlo event generators (MCEGs) are the indispensable workhorses of particle physics, bridging the gap between theoretical ideas and first-principles calculations on the one hand, and the complex detector signatures and data of the experimental community on the other hand. All collider physics experiments are dependent on simulated events by MCEG codes such as Herwig, Pythia, Sherpa, POWHEG, and MG5_aMC@NLO to design and tune their detectors and analysis strategies. The development of MCEGs is overwhelmingly driven by a vibrant community of academics at European Universities, who also train the next generations of particle phenomenologists. The new challenges posed by possible future collider-based experiments and the fact that the first analyses at Run II of the LHC are now frequently limited by theory uncertainties urge the community to invest into further theoretical and technical improvements of these essential tools. In this short contribution to the European Strategy Update, we briefly review the state of the art, and the further developments that will be needed to meet the challenges of the next generation.


title: The mass, location and heating of the dust in the Cassiopeia A supernova remnant

abstract: We model the thermal dust emission from dust grains heated by synchrotron radiation and by particle collisions, under conditions appropriate for four different shocked and unshocked gas components of the Cassiopeia A (Cas A) supernova remnant (SNR). By fitting the resulting spectral energy distributions (SEDs) to the observed SNR dust fluxes, we determine the required mass of dust in each component. We find the observed SED can be reproduced by $\sim 0.6 \, {\rm M_\odot}$ of silicate grains, the majority of which is in the unshocked ejecta and heated by the synchrotron radiation field. Warmer dust, located in the X-ray emitting reverse shock and blastwave regions, contribute to the shorter wavelength infrared emission but make only a small fraction of the total dust mass. Carbon grains can at most make up $\sim 25 \%$ of the total dust mass. Combined with estimates for the gas masses, we obtain dust-to-gas mass ratios for each component, which suggest that the condensation efficiency in the ejecta is high, and that dust in the shocked ejecta clumps is well protected from destruction by sputtering in the reverse shock.


title: Latitudinal differential rotation in the solar analogues 16 Cyg A and B

abstract: Asteroseismology has undergone a profound transformation as a scientific field following the CoRoT and Kepler space missions. The latter is now yielding the first measurements of latitudinal differential rotation obtained directly from oscillation frequencies. Differential rotation is a fundamental mechanism of the stellar dynamo effect. Our goal is to measure the amount of differential rotation in the solar analogues 16 Cyg A and B, which are the components of a binary system. These stars are the brightest observed by Kepler and have therefore been extensively observed, with exquisite precision on their oscillation frequencies. We modelled the acoustic power spectrum of 16 Cyg A and B using a model that takes into account the contribution of differential rotation to the rotational frequency splitting. The estimation was carried out in a Bayesian setting. We then inverted these results to obtain the rotation profile of both stars under the assumption of a solar-like functional form. We observe that the magnitude of latitudinal differential rotation has a strong chance of being solar-like for both stars, their rotation rates being higher at the equator than at the pole. The measured latitudinal differential rotation, defined as the difference of rotation rate between the equator and the pole, is $320\pm269$ nHz and $440^{+363}_{-383}$ nHz for 16 Cyg A and B, respectively, confirming that the rotation rates of these stars are almost solar-like. Their equatorial rotation rates are $535\pm75$ nHz and $565_{-129}^{+150}$ nHz. Our results are in good agreement with measurements obtained from spectropolarimetry, spectroscopy, and photometry. We present the first conclusive measurement of latitudinal differential rotation for solar analogues. Their rotational profiles are very close to those of the Sun. These results depend weakly on the uncertainties of the stellar parameters.


title: Stability results of properties related to the Bishop-Phelps-Bollob\'as property for operators

abstract: We prove that the class of Banach spaces $Y$ such that the pair $(\ell_1, Y)$ has the Bishop-Phelps-Bollob\'as property for operators is stable under finite products when the norm of the product is given by an absolute norm. We also provide examples showing that previous stability results obtained for that property are optimal.


title: Gravitational perturbations of non-singular black holes in conformal gravity

abstract: It is believed that in the near future, gravitational wave detections will become a promising tool not only to test gravity theories, but also to probe extremely curved spacetime regions in our universe, such as the surroundings of black holes. In this paper, we investigate the quasinormal modes (QNMs) of the axial gravitational perturbations of a class of non-singular black holes conformally related to the Schwarzschild black hole. These non-singular black holes can be regarded as the vacuum solution of a family of conformal gravity theories which are invariant under conformal transformations. After conformal symmetry is broken, these black holes produce observational signatures different from those of the Schwarzschild black hole, such as their QNM frequencies. We assume that the spacetime is described by the Einstein equation with the effective energy momentum tensor of an anisotropic fluid. The master equation describing the QNMs is derived, and the QNM frequencies are evaluated with the Wentzel-Kramers-Brillouin (WKB) method up to the 6th order. As expected, the QNM spectra of these non-singular black holes deviate from those of the Schwarzschild black hole, indicating the possibility of testing these black hole solutions with the help of future gravitational wave detections.


title: Practical scheme from QCD to phenomena via Dyson-Schwinger equations

abstract: We deliver a new scheme to compute the quark propagator and the quark-gluon interaction vertex through the coupled Dyson-Schwinger equations (DSEs) of QCD. We take the three-gluon vertex into account in our calculations, and implement the gluon propagator and the running coupling function fitted by the solutions of their respective DSEs. We obtain the momentum and current mass dependence of the quark propagator and the quark-gluon vertex, and the chiral quark condensate which agrees with previous results excellently. We also compute the quark-photon vertex within this scheme and give the anomalous chromo- and electro-magnetic moment of quark. The obtained results also agree with previous ones very well. These applications manifest that the new scheme is realistic and then practical for explaining the QCD-related phenomena.


title: Polynomial bounds on the Sobolev norms of the solutions of the nonlinear wave equation with time dependent potential

abstract: We consider the Cauchy problem for the nonlinear wave equation $u_{tt} - \Delta_x u +q(t, x) u + u^3 = 0$ with smooth potential $q(t, x) \geq 0$ having compact support with respect to $x$. The linear equation without the nonlinear term $u^3$ and potential periodic in $t$ may have solutions with exponentially increasing as $ t \to \infty$ norm $H^1({\mathbb R}^3_x)$. In [2] it was established that adding the nonlinear term $u^3$ the $H^1({\mathbb R}^3_x)$ norm of the solution is polynomially bounded for every choice of $q$. In this paper we show that $H^k({\mathbb R}^3_x)$ norm of this global solution is also polynomially bounded. To prove this we apply a different argument based on the analysis of a sequence $\{Y_k(n\tau_k)\}_{n = 0}^{\infty}$ with suitably defined energy norm $Y_k(t)$ and $0 < \tau_k <1.$


title: On two subclasses of Motzkin paths and their relation to ternary trees

abstract: Two subclasses of Motzkin paths, S-Motzkin and T-Motzkin paths, are introduced. We provide bijections between S-Motzkin paths and ternary trees, S-Motzkin paths and non-crossing trees, and T-Motzkin paths and ordered pairs of ternary trees. Symbolic equations for both paths, and thus generating functions for the paths, are provided. Using these, various parameters involving the two paths are analyzed.


title: Mobility-Aware Joint Service Placement and Routing in Space-Air-Ground Integrated Networks

abstract: People desire to be connected, no matter where they are. Recently, providing Internet access to on-board passengers has received a lot of attention from both industry and academia. However, in order to guarantee an acceptable Quality of Service (QoS) for the passenger services with low incurred cost, the path to route the services, as well as the datacenter (DC) to deploy the services should be carefully determined. This problem is challenging, due to different types of Air-to-Ground (A2G) connections, i.e., satellites and Direct-Air-to-Ground (DA2G) links. These A2G connection types differ in terms of cost, bandwidth, and latency. Furthermore, due to the flights' movements, it is important to consider adapting the service location accordingly. In this work, we formulate two Mixed Integer Linear Programs (MILPs) for the problem of Joint Service Placement and Routing (JSPR): i) Static (S-JSPR), and ii) Mobility-Aware~(MA-JSPR) in Space-Air-Ground Integrated Networks (SAGIN), with the objective of minimizing the total cost. We compare S-JSPR and MA-JSPR using comprehensive evaluations in a realistic European-based SAGIN. The obtained results show that the MA-JSPR model, by considering the future flight positions and using a service migration control, reduces the long-term total cost notably. Also, we show S-JSPR benefits from a low time-complexity and it achieves lower end-to-end delays comparing to MA-JSPR model.


title: From strong to weak temperature dependence of the two-photon entanglement resulting from the biexciton cascade inside a cavity

abstract: We investigate the degree of entanglement quantified by the concurrence of photon pairs that are simultaneously emitted in the biexciton-exciton cascade from a quantum dot in a cavity. Four dot-cavity configurations are compared that differ with respect to the detuning between the cavity modes and the quantum dot transitions, corresponding to different relative weights of direct two-photon and sequential single-photon processes. The dependence of the entanglement on the exciton fine-structure splitting $\delta$ is found to be significantly different for each of the four configurations. For a finite splitting and low temperatures, the highest entanglement is found when the cavity modes are in resonance with the two-photon transition between the biexciton and the ground state and, in addition, the biexciton has a finite binding energy of a few meV. However, this widely used configuration is rather strongly affected by phonons such that other dot-cavity configurations, that are commonly regarded as less suited for obtaining high degrees of entanglement, become more favorable already at temperatures on the order of 10 K and above. If the cavity modes are kept in resonance with one of the exciton-to-ground-state transitions and the biexciton binding energy is finite, the entanglement drastically drops for positive $\delta$ with rising temperatures when $T$ is below $\simeq$ 4 K, but is virtually independent of the temperature for higher $T$.


title: Hope4Genes: a Hopfield-like class prediction algorithm for transcriptomic data

abstract: After its introduction in 1982, the Hopfield model has been extensively applied for classification and pattern recognition. Recently, its great potential in gene expression patterns retrieval has also been shown. Following this line, we develop Hope4Genes a single-sample class prediction algorithm based on a Hopfield-like model. Differently from previous works, we here tested the performances of the algorithm for class prediction, a task of fundamental importance for precision medicine and therapeutic decision-making. Hope4Genes proved better performances than the state-of-art methodologies in the field independently of the size of the input dataset, its profiling platform, the number of classes and the typical class-imbalance present in biological data. Our results provide encoraging evidence that the Hopfield model, together with the use of its energy for the estimation of the false discoveries, is a particularly promising tool for precision medicine.


title: Some algebraic and geometric properties of hyper-K\"ahler fourfold symmetries

abstract: We complete the classification of order $5$ nonsymplectic automorphisms on hyper-K\"ahler fourfolds deformation equivalent to the Hilbert square of a K3 surface. We then compute the topological Lefschetz number of natural automorphisms of generalized Kummer fourfolds and we describe the geometry of their fix loci.


title: The Probabilistic Fault Tolerance of Neural Networks in the Continuous Limit

abstract: The loss of a few neurons in a brain rarely results in any visible loss of function. However, the insight into what "few" means in this context is unclear. How many random neuron failures will it take to lead to a visible loss of function? In this paper, we address the fundamental question of the impact of the crash of a random subset of neurons on the overall computation of a neural network and the error in the output it produces. We study fault tolerance of neural networks subject to small random neuron/weight crash failures in a probabilistic setting. We give provable guarantees on the robustness of the network to these crashes. Our main contribution is a bound on the error in the output of a network under small random Bernoulli crashes proved by using a Taylor expansion in the continuous limit, where close-by neurons at a layer are similar. The failure mode we adopt in our model is characteristic of neuromorphic hardware, a promising technology to speed up artificial neural networks, as well as of biological networks. We show that our theoretical bounds can be used to compare the fault tolerance of different architectures and to design a regularizer improving the fault tolerance of a given architecture. We design an algorithm achieving fault tolerance using a reasonable number of neurons. In addition to the theoretical proof, we also provide experimental validation of our results and suggest a connection to the generalization capacity problem.


title: Optimal Nonparametric Inference via Deep Neural Network

abstract: Deep neural network is a state-of-art method in modern science and technology. Much statistical literature have been devoted to understanding its performance in nonparametric estimation, whereas the results are suboptimal due to a redundant logarithmic sacrifice. In this paper, we show that such log-factors are not necessary. We derive upper bounds for the $L^2$ minimax risk in nonparametric estimation. Sufficient conditions on network architectures are provided such that the upper bounds become optimal (without log-sacrifice). Our proof relies on an explicitly constructed network estimator based on tensor product B-splines. We also derive asymptotic distributions for the constructed network and a relating hypothesis testing procedure. The testing procedure is further proven as minimax optimal under suitable network architectures.


title: Solvability in Gevrey classes of some linear functional equations

abstract: In this paper, we associate to each positive number k a new class of endomorphisms of the sheaf of germs of holomorphic functions on [-1,1] and prove the solvability in the Gevrey class G_k([-1,1]) of some linear functional equations related to endomorphisms.


title: Charged anisotropic strange stars in Finslerian geometry

abstract: We investigate a simplified model for the strange stars in the framework of Finslerian spacetime geometry, composed of charged fluid. It is considered that the fluid consisting of three flavor quarks including a small amount of non-interacting electrons to maintain the chemical equilibrium and assumed that the fluid is compressible by nature. To obtain the simplified form of charged strange star we considered constant flag curvature. Based on geometry, we have developed the field equations within the localized charge distribution. We considered that the strange quarks distributed within the stellar system are compiled with the MIT bag model type of equation of state (EOS) and the charge distribution within the system follows a power law. We represent the exterior spacetime by the Finslerian Ressiner-Nordstr{\"o}m space-time. The maximum anisotropic stress is obtained at the surface of the system. Whether the system is in equilibrium or not, has been examined with respect to the Tolman-Oppenheimer-Volkoff (TOV) equation, Herrera cracking concept, different energy conditions and adiabatic index. We obtain that the total charge is of the order of 10$^{20}$ C and the corresponding electric field is of around 10$^{22}$ V/m. The central density and central pressure vary inversely with the charge. Varying the free parameter (charge constant) of the model, we find the generalized mass-radius variation of strange stars and determine the maximum limited mass with the corresponding radius. Furthermore, we also considered the variation of mass and radius against central density respectively.


title: Topological pressure for conservative $C^1$-diffeomorphisms with no dominated splitting

abstract: We prove three formulas for computing topological pressure of $C^1$-generic conservative diffeomorphism and show the continuity of topological pressure with respect to these diffeomorphisms. We prove for these generic diffeomorphisms that there is no equilibrium states with positive measure theoretic entropy. In particular, for hyperbolic potentials, there is no equilibrium states. For $C^1$ generic conservative diffeomorphism on compact surfaces with no dominated splitting and $\phi_m(x):=-\frac{1}{m}\log \Vert D_x f^m\Vert, m \in \mathbb{N}$, we show that there exist equilibrium states with zero entropy and there exists a transition point $t_0$ for the family $\lbrace t \phi_m(x)\rbrace_{t\geq 0}$, such that there is no equilibrium states for $ t \in [0, t_0)$ and there is an equilibrium state for $t \in [t_0,+\infty)$.


title: Accuracy Evaluation of Overlapping and Multi-resolution Clustering Algorithms on Large Datasets

abstract: Performance of clustering algorithms is evaluated with the help of accuracy metrics. There is a great diversity of clustering algorithms, which are key components of many data analysis and exploration systems. However, there exist only few metrics for the accuracy measurement of overlapping and multi-resolution clustering algorithms on large datasets. In this paper, we first discuss existing metrics, how they satisfy a set of formal constraints, and how they can be applied to specific cases. Then, we propose several optimizations and extensions of these metrics. More specifically, we introduce a new indexing technique to reduce both the runtime and the memory complexity of the Mean F1 score evaluation. Our technique can be applied on large datasets and it is faster on a single CPU than state-of-the-art implementations running on high-performance servers. In addition, we propose several extensions of the discussed metrics to improve their effectiveness and satisfaction to formal constraints without affecting their efficiency. All the metrics discussed in this paper are implemented in C++ and are available for free as open-source packages that can be used either as stand-alone tools or as part of a benchmarking system to compare various clustering algorithms.


title: Demonstration of teleportation-based error correction in the IBM quantum computer

abstract: Quantum error correcting codes (QECC) are the key ingredients both for fault-tolerant quantum computation and quantum communication. Teleportation-based error correction (TEC) helps in detecting and correcting operational and erasure errors by performing X and Z measurements during teleportation. Here we demonstrate the TEC protocol for the detection and correction of a single bit-flip error by proposing a new quantum circuit. A single phase-flip error can also be detected and corrected using the above protocol. For the first time, we illustrate detection and correction of erasure error in the superconducting qubit-based IBM's 14-qubit quantum computer.


title: Bibliometrics for collaboration works

abstract: An important issue in bibliometrics is the weighing of co-authorship in the production of scientific collaborations, which are becoming the standard modality of research activity in many disciplines. The problem is especially relevant in the field of high-energy physics, where collaborations reach 3000 authors, but it can no longer be ignored also in other domains, like medicine or biology. We present theoretical and numerical arguments in favour of weighing the individual contributions as $1/N_{\rm aut}^\alpha$ where $N_{\rm aut}$ is the number of co-authors. When counting citations we suggest the exponent $\alpha\approx 1$, that corresponds to fractional counting. When counting the number of papers we suggest $\alpha \approx 1/3 - 1/2$, with the former (latter) value more appropriate for larger (smaller) collaborations. We expect and verify that the $h$ index scales as the square root of the average number of co-authors, and define a fractionalized $h$ index that does not scale with collaboration size.


title: The classical counterpart of the Aharonov-Bohm phase

abstract: The existence of the Aharonov-Bohm phase shows that the magnetic field may produce nonlocal effects in quantum mechanics. It is generally believed that such a nonlocal behavior of the magnetic field is not possible in classical physics and that this represents a clear difference between classical and quantum mechanics. Contrary to these beliefs, we argue that the classical counterpart of the Aharonov-Bohm phase, which is identified here with the electromagnetic angular momentum of the configuration formed by an electric charge moving around an infinitely-long solenoid enclosing a uniform magnetic flux, shows that the magnetic field may produce nonlocal effects in classical mechanics. We discuss this momentum in detail by putting special emphasis on its topological and nonlocal features. The momentum is topological because it depends on the number of windings the electric charge carries out around the solenoid and is nonlocal because the magnetic flux has no local consequences at any point on the charge trajectory. The topological feature allows us to introduce the concept of accumulated electromagnetic angular momentum and the nonlocal feature allows us to speak of a dynamical nonlocality attributable to the equations of classical electrodynamics.


title: Mapping the stellar age of the Milky Way bulge with the VVV. I. The method

abstract: Recent observational programmes are providing a global view of the Milky Way bulge that serves as template for detailed comparison with models and extragalactic bulges. A number of surveys (i.e. VVV, GIBS, GES, ARGOS, BRAVA, APOGEE) are producing comprehensive and detailed extinction, metallicity, kinematics and stellar density maps of the Galactic bulge with unprecedented accuracy. However, the still missing key ingredient is the distribution of stellar ages across the bulge. To overcome this limitation, we aim to age-date the stellar population in several bulge fields with the ultimate goal of deriving an age map of the Bulge. This paper presents the methodology and the first results obtained for a field along the Bulge minor axis, at $b=-6^\circ$. We use a new PSF-fitting photometry of the VISTA Variables in the V\'{i}a L\'{a}ctea (VVV) survey data to construct deep color-magnitude diagrams of the bulge stellar population down to $\sim$ 2 mag below the Main Sequence turnoff. We find the bulk of the bulge stellar population in the observed field along the minor axis to be at least older than $\sim$ 7.5 Gyr. In particular, when the metallicity distribution function spectroscopically derived by GIBS is used, the best fit to the data is obtained with a combination of synthetic populations with ages in between $\sim$ 7.5 Gyr and 11 Gyr. However, the fraction of stars younger than $\sim$ 10 Gyr strongly depends upon the number of Blue Straggler Stars present in the bulge. Simulations show that the observed color-magnitude diagram of the bulge in the field along the minor axis is incompatible with the presence of a conspicuous population of intermediate-age/young (i.e. $\lesssim 5$ Gyr) stars.


title: Intriguing & intuitive relation of diagonal Riemann tensor components to the corresponding 2-d Gaussian curvature for diagonal metrics of any dimensionality

abstract: Using Gauss's square-roots of the metric components, the diagonal Riemann tensor components for diagonal metrics are calculated. The result is a form which makes their source in the metric directly intuitive and displays an intriguing relation to Gaussian curvature. Several examples of calculation utilizing this formula are presented, and a speculative quesiton is raised about the possibility of an invariant characteristic to spaces amenable to orthogonal coordinates/diagonal metrics.


title: Learning to Solve Large-Scale Security-Constrained Unit Commitment Problems

abstract: Security-Constrained Unit Commitment (SCUC) is a fundamental problem in power systems and electricity markets. In practical settings, SCUC is repeatedly solved via Mixed-Integer Linear Programming, sometimes multiple times per day, with only minor changes in input data. In this work, we propose a number of machine learning (ML) techniques to effectively extract information from previously solved instances in order to significantly improve the computational performance of MIP solvers when solving similar instances in the future. Based on statistical data, we predict redundant constraints in the formulation, good initial feasible solutions and affine subspaces where the optimal solution is likely to lie, leading to significant reduction in problem size. Computational results on a diverse set of realistic and large-scale instances show that, using the proposed techniques, SCUC can be solved on average 4.3x faster with optimality guarantees, and 10.2x faster without optimality guarantees, but with no observed reduction in solution quality. Out-of-distribution experiments provides evidence that the method is somewhat robust against dataset shift.


title: Stochastic Enumeration with Importance Sampling

abstract: Many hard problems in the computational sciences are equivalent to counting the leaves of a decision tree, or, more generally, summing a cost function over the nodes. These problems include calculating the permanent of a matrix, finding the volume of a convex polyhedron, and counting the number of linear extensions of a partially ordered set. Many approximation algorithms exist to estimate such sums. One of the most recent is Stochastic Enumeration (SE), introduced in 2013 by Rubinstein. In 2015, Vaisman and Kroese provided a rigorous analysis of the variance of SE, and showed that SE can be extended to a fully polynomial randomized approximation scheme for certain cost functions on random trees. We present an algorithm that incorporates an importance function into SE, and provide theoretical analysis of its efficacy. We also present the results of numerical experiments to measure the variance of an application of the algorithm to the problem of counting linear extensions of a poset, and show that introducing importance sampling results in a significant reduction of variance as compared to the original version of SE.