title: On a $q-$analog of a singularly perturbed problem of irregular type with two complex time variables

abstract: Analytic solutions and their formal asymptotic expansions for a family of the singularly perturbed $q-$difference-differential equations in the complex domain are constructed. They stand for a $q-$analog of the singularly perturbed partial differential equations considered in our recent work [A. Lastra, S. Malek, Boundary layer expansions for initial value problems with two complex time variables, submitted 2019]. In the present work, we construct outer and inner analytic solutions of the main equation, each of them showing asymptotic expansions of essentially different nature with respect to the perturbation parameter. The appearance of the $-1$-branch of Lambert $W$ function will be crucial in this respect.


title: A homogenized limit for the 2D Euler equations in a perforated domain

abstract: We study the motion of an ideal incompressible fluid in a perforated domain. The porous medium is composed of inclusions of size $a$ separated by distances $\tilde d$ and the fluid fills the exterior. We analyse the asymptotic behavior of the fluid when $(a,\tilde d) \to (0,0)$. If the inclusions are distributed on the unit square, this issue is studied recently when $\frac{\tilde d}a$ tends to zero or infinity, leaving aside the critical case where the volume fraction of the porous medium is below its possible maximal value but non-zero. In this paper, we provide the first result in this regime. In contrast with former results, we obtain an Euler type equation where a homogenized term appears in the elliptic problem relating the velocity and the vorticity. Our analysis is based on the so-called method of reflections whose convergence provides novel estimates on the solutions to the div-curl problem which is involved in the 2D-Euler equations.


title: Linear MIM-Width of Trees

abstract: We provide an $O(n \log n)$ algorithm computing the linear maximum induced matching width of a tree and an optimal layout.


title: Rapid Node Cardinality Estimation in Heterogeneous Machine-to-Machine Networks

abstract: Machine-to-Machine (M2M) networks are an emerging technology with applications in various fields, including smart grids, healthcare, vehicular telematics and smart cities. Heterogeneous M2M networks contain different types of nodes, e.g., nodes that send emergency, periodic, and normal type data. An important problem is to rapidly estimate the number of active nodes of each node type in every time frame in such a network. In this paper, we design two schemes for estimating the active node cardinalities of each node type in a heterogeneous M2M network with $T$ types of nodes, where $T \ge 2$ is an arbitrary integer. Our schemes consist of two phases-- in phase 1, coarse estimates are computed, and in phase 2, these estimates are used to compute the final estimates to the required accuracy. We analytically derive a condition for one of our schemes that can be used to decide as to which of two possible approaches should be used in phase 2 to minimize its execution time. The expected number of time slots required to execute and the expected energy consumption of each active node under one of our schemes are analysed. Using simulations, we show that our proposed schemes require significantly fewer time slots to execute compared to estimation schemes designed for a heterogeneous M2M network in prior work, and also, compared to separately executing a well-known estimation protocol designed for a homogeneous network in prior work $T$ times to estimate the cardinalities of the $T$ node types, even though all these schemes obtain estimates with the same accuracy.


title: A Bridge Anchored on Both Sides: Formal Deduction in Introductory CS, and Code Proofs in Discrete Math

abstract: There is a sharp disconnect between the programming and mathematical portions of the standard undergraduate computer science curriculum, leading to student misunderstanding about how the two are related. We propose connecting the subjects early in the curriculum---specifically, in CS1 and the introductory discrete mathematics course---by using formal reasoning about programs as a bridge between them. This article reports on Haverford and Grinnell College's experience in constructing the end points of this bridge between CS1 and discrete mathematics. Haverford's long-standing "3-2-1" curriculum introduces code reasoning in conjunction with introductory programming concepts, and Grinnell's discrete mathematics introduces code reasoning as a motivation for logic and formal deduction. Both courses present code reasoning in a style based on symbolic code execution techniques from the programming language community, but tuned to address the particulars of each course. These courses rely primarily on traditional means of proof authoring with pen-and-paper. This is unsatisfactory for students who receive no feedback until grading on their work and instructors who must shoulder the burden of interpreting students' proofs and giving useful feedback. To this end, we also describe the current state of Orca, an in-development proof assistant for undergraduate education that we are developing to address these issues in our courses. Finally, in teaching our courses, we have discovered a number of educational research questions about the effectiveness of code reasoning in bridging the gap between programming and mathematics, and the ability of tools like \orca to support this pedagogy. We pose these research questions as next steps to formalize our initial experiences in our courses with the hope of eventually generalizing our approaches for wider adoption.


title: The What-If Tool: Interactive Probing of Machine Learning Models

abstract: A key challenge in developing and deploying Machine Learning (ML) systems is understanding their performance across a wide range of inputs. To address this challenge, we created the What-If Tool, an open-source application that allows practitioners to probe, visualize, and analyze ML systems, with minimal coding. The What-If Tool lets practitioners test performance in hypothetical situations, analyze the importance of different data features, and visualize model behavior across multiple models and subsets of input data. It also lets practitioners measure systems according to multiple ML fairness metrics. We describe the design of the tool, and report on real-life usage at different organizations.


title: Emergency Management Systems and Algorithms: a Comprehensive Survey

abstract: Owing to the increasing frequency and destruction of natural and manmade disasters to modern highly-populated societies, emergency management, which provides solutions to prevent or address disasters, have drawn considerable research over the last few decades and become a multidisciplinary area. Because of its open and inclusive nature, new technologies always tend to influence, change or even revolutionise this research area. Hence, it is imperative to consolidate the state-of-the-art studies and knowledge to meet the research needs and identify the future research directions. The paper presents a comprehensive and systemic review of the existing research in the field of emergency management from both the system design aspect and algorithm engineering aspect. We begin with the history and evolution of the emergency management research. Then the two main research topics of this area, "emergency navigation" and "emergency search and rescue planning", are introduced and discussed. Finally, we suggest the emerging challenges and opportunities from system optimisation, evacuee behaviour modelling and optimisation, computing patterns, data analysis, energy and cyber security aspects.


title: Geometric inequivalence of metric and Palatini formulations of General Relativity

abstract: Projective invariance is a symmetry of the Palatini version of General Relativity which is not present in the metric formulation. The fact that the Riemann tensor changes nontrivially under projective transformations implies that, unlike in the usual metric approach, in the Palatini formulation this tensor is subject to a gauge freedom, which allows some ambiguities even in its scalar contractions. In this sense, we show that for the Schwarzschild solution there exists a projective gauge in which the (affine) Kretschmann scalar, $K\equiv {R^\alpha}_{\beta\mu\nu}{R_\alpha}^{\beta\mu\nu}$, can be set to vanish everywhere. This puts forward that the divergence of curvature scalars may, in some cases, be avoided by a gauge transformation of the connection.


title: Characterization of Overlap in Observational Studies

abstract: Overlap between treatment groups is required for non-parametric estimation of causal effects. If a subgroup of subjects always receives the same intervention, we cannot estimate the effect of intervention changes on that subgroup without further assumptions. When overlap does not hold globally, characterizing local regions of overlap can inform the relevance of causal conclusions for new subjects, and can help guide additional data collection. To have impact, these descriptions must be interpretable for downstream users who are not machine learning experts, such as policy makers. We formalize overlap estimation as a problem of finding minimum volume sets subject to coverage constraints and reduce this problem to binary classification with Boolean rule classifiers. We then generalize this method to estimate overlap in off-policy policy evaluation. In several real-world applications, we demonstrate that these rules have comparable accuracy to black-box estimators and provide intuitive and informative explanations that can inform policy making.


title: To save the Environmental degradation: Please pay for ecosystem services!

abstract: Project planners should consider the monetary value of ecosystem services in the regional service category and add it to the cost calculation of the project as environmental cost. Managers should monitor and regularly maintain the ecological environment within the project area to maintain the benefits of environmental costs.With regard to the problem of predicting how the model changes over time, we use LSTM to realize dynamic evaluation. The analysis that has been made shows that the monetary value of urban ecological services has increased. The environment is still being damaged. Through neural network memory and learning, the problem can be found, and the model can be modified in time.


title: Wandering domains arising from Lavaurs maps with Siegel disks

abstract: The classification of Fatou components for rational functions was concluded with Sullivan's proof of the No Wandering Domains Theorem in 1985. In 2016 it was shown, in joint work of the first and last author with Buff, Dujardin and Raissy, that wandering domains do exist in higher dimensions. In fact, wandering domains arise even for a seemingly simple class of maps: polynomial skew products. While the construction gives an infinite dimensional class of examples, and has been extended to polynomial automorphisms of $\mathbb C^4$ by Hahn and the last author, the currently known wandering domains are essentially unique. Our goal in this paper is to construct a second example, arising from similar techniques, but with distinctly different dynamical behavior. Instead of wandering domains arising from a Lavaurs map with an attracting fixed point, we construct a domain arising from a Lavaurs map with a fixed point of Siegel type. Siegel disks are not robust under perturbations, as opposed to attracting fixed points. We prove a necessary and sufficient condition for the existence of a so-called trapping domain for non-autonomous dynamical systems given by sequences of maps converging parabolically towards a Siegel type limit map. Guaranteeing that this condition is satisfied in our current construction requires a reconsideration of the proof of the original wandering domain, as more precise estimates on the rate of convergence towards the Lavaurs map are required. By adapting ideas introduced recently by Bedford, Smillie and Ueda, and by proving the existence of parabolic curves, with control on their domains of definition, we prove that the convergence rate is parabolic.


title: Spatial Reuse in IEEE 802.11ax WLANs

abstract: Dealing with massively crowded scenarios is one of the most ambitious goals of next-generation wireless networks. With this goal in mind, the IEEE 802.11ax amendment includes, among other techniques, the Spatial Reuse (SR) operation. The SR operation encompasses a set of unprecedented techniques {that are expected to significantly boost Wireless Local Area Networks (WLANs) performance in dense environments}. In particular, the main objective of the SR operation is to maximize the utilization of the medium by increasing the number of parallel transmissions. Nevertheless, due to the novelty of the operation, its performance gains remain largely unknown. In this paper, we first provide a gentle tutorial of the SR operation included in the IEEE 802.11ax. Then, we analytically model SR and delve into the new kinds of MAC-level interactions among network devices. Finally, we provide a simulation-driven analysis to showcase the potential of SR in various deployments, comprising different network densities and traffic loads. Our results show that the SR operation can significantly improve the medium utilization, especially in scenarios under high interference conditions. Moreover, our results demonstrate the non-intrusive design characteristic of SR, which allows enhancing the number of simultaneous transmissions with a low impact on the environment. We conclude the paper by giving some thoughts on the main challenges and limitations of the IEEE 802.11ax SR operation, including research gaps and future directions.


title: Energy stable schemes for gradient flows based on novel auxiliary variable with energy bounded above

abstract: In this paper, we consider a novel auxiliary variable method to obtain energy stable schemes for gradient flows. The auxiliary variable based on energy bounded above does not limited to the hypothetical conditions adopted in previous approaches. We proved the unconditional energy stability for all the semi-discrete schemes carefully and rigorously. The novelty of the proposed schemes is that the computed values for the functional in square root are guaranteed to be positive. This method, termed novel auxiliary energy variable (NAEV) method does not consider any bounded below restrictions any longer. However, these restrictions are necessary in invariant energy quadratization (IEQ) and scalar auxiliary variable (SAV) approaches which are very popular methods recently. This property of guaranteed positivity is not available in previous approaches. A comparative study of classical SAV and NAEV approaches is considered to show the accuracy and efficiency. Finally, we present various 2D numerical simulations to demonstrate the stability and accuracy.


title: On the Tate and Standard Conjectures over Finite Fields

abstract: For an abelian variety over a finite field, Clozel (1999) showed that l-homological equivalence coincides with numerical equivalence for infinitely many l, and the author (1999) gave a criterion for the Tate conjecture to follow from Tate's theorem on divisors. We generalize both statements to motives, and apply them to other varieties including K3 surfaces.


title: From rotating fluid masses and Ziegler's paradox to Pontryagin- and Krein spaces and bifurcation theory

abstract: Three classical systems, the Kelvin gyrostat, the Maclaurin spheroids, and the Ziegler pendulum have directly inspired development of the theory of Pontryagin and Krein spaces with indefinite metric and singularity theory as independent mathematical topics, not to mention stability theory and nonlinear dynamics. From industrial applications in shipbuilding, turbomachinery, and artillery to fundamental problems of astrophysics, such as asteroseismology and gravitational radiation --- that is the range of phenomena involving the Krein collision of eigenvalues, dissipation-induced instabilities, and spectral and geometric singularities on the neutral stability surfaces, such as the famous Whitney's umbrella.

categories: math.DS math-ph math.MP nlin.CD physics.class-ph


title: Feasibility of Cosmic Microwave Background Observations using Radiometers based on Whispering Gallery Mode Resonators

abstract: The fundamentals of the whispering gallery mode (WGM) resonators are well established in the literature, with several successful proof-of-concept experiments. One remarkable benefit of this technology is the room-temperature operation. This characteristic could be used to build a new generation of radiometers that do not need to be cooled down to cryogenic temperatures to reach high sensitivities. In this article, a study of the viability of technological transfer is undertaken, beginning with a brief review of the theoretical background that will be applied and leading to a proposal for a novel spectro-polarimeter design. Simulations for a radiometer based on WGM resonance are analyzed and compared with state-of-the-art coherent receivers. The results are then discussed, finding that although promising, WGM technology needs R\&D in several directions in order to be competitive, whose are suggested by the authors with the idea of inspiring the work of the researchers in the field towards a new direction or approach.


title: Large, tunable valley splitting and single-spin relaxation mechanisms in a Si/Si$_x$Ge$_{1-x}$ quantum dot

abstract: Valley splitting is a key figure of silicon-based spin qubits. Quantum dots in Si/SiGe heterostructures reportedly suffer from a relatively low valley splitting, limiting the operation temperature and the scalability of such qubit devices. Here, we demonstrate a robust and large valley splitting exceeding 200 $\mu$eV in a gate-defined single quantum dot, hosted in molecular-beam epitaxy-grown $^{28}$Si/SiGe. The valley splitting is monotonically and reproducibly tunable up to 15 % by gate voltages, originating from a 6 nm lateral displacement of the quantum dot. We observe static spin relaxation times $T_1>1$ s at low magnetic fields in our device containing an integrated nanomagnet. At higher magnetic fields, $T_1$ is limited by the valley hotspot and by phonon noise coupling to intrinsic and artificial spin-orbit coupling, including phonon bottlenecking.


title: Adaptive inference for a semiparametric generalized autoregressive conditional heteroskedasticity model

abstract: This paper considers a semiparametric generalized autoregressive conditional heteroskedasticity (S-GARCH) model. For this model, we first estimate the time-varying long run component for unconditional variance by the kernel estimator, and then estimate the non-time-varying parameters in GARCH-type short run component by the quasi maximum likelihood estimator (QMLE). We show that the QMLE is asymptotically normal with the parametric convergence rate. Next, we construct a Lagrange multiplier test for linear parameter constraint and a portmanteau test for model checking, and obtain their asymptotic null distributions. Our entire statistical inference procedure works for the non-stationary data with two important features: first, our QMLE and two tests are adaptive to the unknown form of the long run component; second, our QMLE and two tests share the same efficiency and testing power as those in variance targeting method when the S-GARCH model is stationary.


title: Multiple membership multilevel models

abstract: Multiple membership multilevel models are an extension of standard multilevel models for non-hierarchical data that have multiple membership structures. Traditional multilevel models involve hierarchical data structures whereby lower-level units such as students are nested within higher-level units such as schools and where these higher-level units may in turn be nested within further groupings or clusters such as school districts, regions, and countries. With hierarchical data structures, there is an exact nesting of each lower-level unit in one and only one higher-level unit. For example, each student attends one school, each school is located within one school district, and so on. However, social reality is more complicated than this, and so social and behavioural data often do not follow pure or strict hierarchies. Two types of non-hierarchical data structures which often appear in practice are cross-classified and multiple membership structures. In this article, we describe multiple membership data structures and multiple membership models which can be used to analyse them.


title: Answer Extraction for Why Arabic Questions Answering Systems: EWAQ

abstract: With the increasing amount of web information, questions answering systems becomes very important to allow users to access to direct answers for their requests. This paper presents an Arabic Questions Answering Systems based on entailment metrics. The type of questions which this paper focuses on is why questions. There are many reasons lead us to develop this system: generally, the lack of Arabic Questions Answering Systems and scarcity Arabic Questions Answering Systems which focus on why questions. The goal of the proposed system in this research is to extract answers from re-ranked retrieved passages which are retrieved by search engines. This system extracts the answer only to why questions. This system is called by EWAQ: Entailment based Why Arabic Questions Answering. Each answer is scored with entailment metrics and ranked according to their scores in order to determine the most possible correct answer. EWAQ is compared with search engines: yahoo, google and ask.com, the well-established web-based Questions Answering systems, using manual test set. In EWAQ experiments, it is showed that the accuracy is increased by implementing the textual entailment in re-raking the retrieved relevant passages by search engines and deciding the correct answer. The obtained results show that using entailment based similarity can help significantly to tackle the why Answer Extraction module in Arabic language.


title: Nonnegative Matrix Factorization with Local Similarity Learning

abstract: Existing nonnegative matrix factorization methods focus on learning global structure of the data to construct basis and coefficient matrices, which ignores the local structure that commonly exists among data. In this paper, we propose a new type of nonnegative matrix factorization method, which learns local similarity and clustering in a mutually enhancing way. The learned new representation is more representative in that it better reveals inherent geometric property of the data. Nonlinear expansion is given and efficient multiplicative updates are developed with theoretical convergence guarantees. Extensive experimental results have confirmed the effectiveness of the proposed model.


title: Charged pion condensation and duality in dense and hot chirally and isospin asymmetric quark matter in the framework of NJL$_2$ model

abstract: In this paper we investigate in the large-$N_c$ limit ($N_c$ is the number of colored quarks) the phase structure of a massless (1+1)-dimensional quark model with four-quark interaction and in the presence of baryon ($\mu_B$), isospin ($\mu_I$) and chiral isospin ($\mu_{I5}$) chemical potentials as well as at nonzero temperature. It is established that chiral isospin chemical potential leads to the generation of charged pion condensation (PC) in dense (nonzero baryon density) and chiral asymmetric quark matter for a wide range of isospin densities. It is shown that there exists a duality correspondence between the chiral symmetry breaking and the charged PC phenomena at any values of temperature even for very hot quark gluon plasma. Moreover, it is shown that charged PC phase with nonzero baryon density can be induced in the model at comparatively high temperatures. This opens up new possible physical systems, where it can be of importance, such as heavy ion collisions, just born neutron stars (proto-neutron stars), supernovas as well as neutron star mergers.


title: Interpretable Segmentation of Medical Free-Text Records Based on Word Embeddings

abstract: Is it true that patients with similar conditions get similar diagnoses? In this paper we show NLP methods and a unique corpus of documents to validate this claim. We (1) introduce a method for representation of medical visits based on free-text descriptions recorded by doctors, (2) introduce a new method for clustering of patients' visits and (3) present an~application of the proposed method on a corpus of 100,000 visits. With the proposed method we obtained stable and separated segments of visits which were positively validated against final medical diagnoses. We show how the presented algorithm may be used to aid doctors during their practice.


title: Non-homogeneous extensions of Cantor minimal systems

abstract: Floyd gave an example of a minimal dynamical system which was an extension of an odometer and the fibres of the associated factor map were either singletons or intervals. Gjerde and Johansen showed that the odometer could be replaced by any Cantor minimal system. Here, we show further that the intervals can be generalized to cubes of arbitrary dimension and to attractors of certain iterated function systems. We discuss applications.


title: EU H2020 Gauss project. Geo-Fencing Software System

abstract: The Geofencing system is the key to operate the Unmanned Aerial Vehicle (UAV) within the safe and appropriate zone to avoid public concerns and other privacy issues. The system is designed to keep the UAV away from geofenced obstacles using the onboard GNSS and IMU location. The Geofencing system is part of the H2020 GAUSS project and facilities other subsystems, for instance, to support the command and control link, which is the security measure to secure the UAV from hijacking and signal spoofing. The regulatory authorities expressed the concern of having UAVs flying in the no-fly zone and causing troubles from offending private privacy to hazards at airport airspace. Hence the geofence system shall provide guidance message, which enables the UAV to evacuate from no-fly-zone, based on real-time updated location. This thesis aims to first illustrate the generation of geofence and then apply the geofence system on UAV operation. This application enables UAV to fly in the designated area without human intervention. The project is built with JAVA using GIS-enabled Database Management System and Open Soured Map data powered by OpenStreetMap and OS map. This method has been tested by simulations which had results of high accuracy.


title: GP-VAE: Deep Probabilistic Time Series Imputation

abstract: Multivariate time series with missing values are common in areas such as healthcare and finance, and have grown in number and complexity over the years. This raises the question whether deep learning methodologies can outperform classical data imputation methods in this domain. However, naive applications of deep learning fall short in giving reliable confidence estimates and lack interpretability. We propose a new deep sequential latent variable model for dimensionality reduction and data imputation. Our modeling assumption is simple and interpretable: the high dimensional time series has a lower-dimensional representation which evolves smoothly in time according to a Gaussian process. The non-linear dimensionality reduction in the presence of missing data is achieved using a VAE approach with a novel structured variational approximation. We demonstrate that our approach outperforms several classical and deep learning-based data imputation methods on high-dimensional data from the domains of computer vision and healthcare, while additionally improving the smoothness of the imputations and providing interpretable uncertainty estimates.


title: Private key encryption and recovery in blockchain

abstract: The disruptive technology of blockchain can deliver secure solutions without the need for a central authority. In blockchain protocols, assets that belong to a participant are controlled through the private key of an asymmetric key pair that is owned by the participant. Although, this lets blockchain network participants to have sovereignty on their assets, it comes with the responsibility of managing their own keys. Currently, there exists two major bottlenecks in managing keys; $a)$ users don't have an efficient and secure way to store their keys, $b)$ no efficient recovery mechanism exists in case the keys are lost. In this study, we propose secure methods to efficiently store and recover keys. For the first, we introduce an efficient encryption mechanism to securely encrypt and decrypt the private key using the owner's biometric signature. For the later, we introduce an efficient recovery mechanism using biometrics and secret sharing scheme. By applying the proposed key encryption and recovery mechanism, asset owners are able to securely store their keys on their devices and recover the keys in case they are lost.


title: Long-lived nuclear coherences inside helium nanodroplets

abstract: Much of our knowledge about dynamics and functionality of molecular systems has been achieved with femtosecond time-resolved spectroscopy. Despite extensive technical developments over the past decades, some classes of systems have eluded dynamical studies so far. Here, we demonstrate that superfluid helium nanodroplets, acting as thermal bath of 0.4 K temperature to stabilize weakly bound or reactive systems, are well suited for time-resolved studies of single molecules solvated in the droplet interior. By observing vibrational wavepacket motion of indium dimers (In$_2$) for over 50 ps, we demonstrate that the perturbation imposed by this quantum liquid can be lower by a factor of 10-100 compared to any other solvent, which uniquely allows to study processes depending on long nuclear coherence in a dissipative environment. Furthermore, tailor-made microsolvation environments inside droplets will enable to investigate the solvent influence on intramolecular dynamics in a wide tuning range from molecular isolation to strong molecule-solvent coupling.


title: Well-posedness of Stochastic Port-Hamiltonian Systems on Infinite-dimensional Spaces

abstract: Stochastic port-Hamiltonian systems on infinite-dimensional spaces governed by It\^o stochastic differential equations (SDEs) are introduced and some properties of this new class of systems are studied. They are an extension of stochastic port-Hamiltonian systems defined on a finite-dimensional state space. The concept of well-posedness in the sense of Weiss and Salamon is generalized to the stochastic context. Under this extended definition, stochastic port-Hamiltonian systems are shown to be well-posed. The theory is illustrated on an example of a vibrating string subject to a Hilbert space-valued Gaussian white noise process.


title: Quantum wave-particle superposition in a delayed-choice experiment

abstract: Wave-particle duality epitomizes the counterintuitive character of quantum physics. A striking illustration is the quantum delay-choice experiment, which is based on Wheeler's classic delayed-choice gedanken experiment, but with the addition of a quantum-controlled device enabling wave-to-particle transitions. Here we realize a quantum delayed-choice experiment in which we control the wave and the particle states of photons and in particular the phase between them, thus directly establishing the created quantum nature of the wave-particle. We generate three-photon entangled states and inject one photon into a Mach--Zehnder interferometer embedded in a 187-m-long two-photon Hong-Ou-Mandel interferometer. The third photon is sent 141m away from the interferometers and remotely prepares a two-photon quantum gate according to independent active choices under Einstein locality conditions. We have realized transitions between wave and particle states in both classical and quantum scenarios, and therefore tests of the complementarity principle that go fundamentally beyond earlier implementations.


title: Learning in Competitive Network with Haeusslers Equation adapted using FIREFLY algorithm

abstract: Many of the competitive neural network consists of spatially arranged neurons. The weigh matrix that connects cells represents local excitation and long-range inhibition. They are known as soft-winner-take-all networks and shown to exhibit desirable information-processing. The local excitatory connections are many times predefined hand-wired based depending on spatial arrangement which is chosen using the previous knowledge of data. Here we present learning in recurrent network through Haeusslers equation and modified wiring scheme based on biologically based Firefly algorithm. Following results show learning in such network from input patterns without hand-wiring with fixed topology.


title: A Note on String Excitations and the Higuchi Bound

abstract: In this brief note we consider the interaction between high spin excitations in string theory along the Regge trajectory and the Higuchi bound in de Sitter space. There is always a point along the Regge trajectory where the Higuchi bound is violated. However, this point precisely coincides with a string whose length is of order the de Sitter Hubble scale. String theory therefore manifests no immediate inconsistency as long as the string scale $M_s$ is above the Hubble scale $H$. However, an implication is that the Regge trajectory must be significantly modified at some ultraviolet scale. Insisting that this modification should occur no earlier than the Planck scale would lead to a bound on the string scale of $M_s > \sqrt{H M_p}$.


title: Optimality of impulse control problem in refracted L\'evy model with Parisian ruin and transaction costs

abstract: In this paper we investigate an optimal dividend problem with transaction costs, where the surplus process is modelled by a refracted L\'evy process and the ruin time is considered with Parisian delay. Presence of the transaction costs implies that one need to consider the impulse control problem as a control strategy in such model. An impulse policy $(c_1,c_2)$, which is to reduce the reserves to some fixed level $c_1$ whenever they are above another level $c_2$ is an important strategy for the impulse control problem. Therefore, we give sufficient conditions under which the above described impulse policy is optimal. Further, we give the new analytical formulas for the Parisian refracted $q$-scale functions in the case of the linear Brownian motion and the Cr\'amer-Lundberg process with exponential claims. Using these formulas we show that for these models there exists a unique $(c_1, c_2)$ policy which is optimal for the impulse control problem. Numerical examples are also provided.


title: Approximately Stable Matchings with General Constraints

abstract: This paper focuses on two-sided matching where one side (a hospital or firm) is matched to the other side (a doctor or worker) so as to maximize a cardinal objective under general feasibility constraints. In a standard model, even though multiple doctors can be matched to a single hospital, a hospital has a responsive preference and a maximum quota. However, in practical applications, a hospital has some complicated cardinal preference and constraints. With such preferences (e.g., submodular) and constraints (e.g., knapsack or matroid intersection), stable matchings may fail to exist. This paper first determines the complexity of checking and computing stable matchings based on preference class and constraint class. Second, we establish a framework to analyze this problem on packing problems and the framework enables us to access the wealth of online packing algorithms so that we construct approximately stable algorithms as a variant of generalized deferred acceptance algorithm. We further provide some inapproximability results.


title: Which Algorithmic Choices Matter at Which Batch Sizes? Insights From a Noisy Quadratic Model

abstract: Increasing the batch size is a popular way to speed up neural network training, but beyond some critical batch size, larger batch sizes yield diminishing returns. In this work, we study how the critical batch size changes based on properties of the optimization algorithm, including acceleration and preconditioning, through two different lenses: large scale experiments, and analysis of a simple noisy quadratic model (NQM). We experimentally demonstrate that optimization algorithms that employ preconditioning, specifically Adam and K-FAC, result in much larger critical batch sizes than stochastic gradient descent with momentum. We also demonstrate that the NQM captures many of the essential features of real neural network training, despite being drastically simpler to work with. The NQM predicts our results with preconditioned optimizers, previous results with accelerated gradient descent, and other results around optimal learning rates and large batch training, making it a useful tool to generate testable predictions about neural network optimization.


title: Global Cardinality Constraints Make Approximating Some Max-2-CSPs Harder

abstract: Assuming the Unique Games Conjecture, we show that existing approximation algorithms for some Boolean Max-2-CSPs with cardinality constraints are optimal. In particular, we prove that Max-Cut with cardinality constraints is UG-hard to approximate within \approx 0.858, and that Max-2-Sat with cardinality constraints is UG-hard to approximate within \approx 0.929. In both cases, the previous best hardness results were the same as the hardness of the corresponding unconstrained Max-2-CSP (\approx 0.878 for Max-Cut, and \approx 0.940 for Max-2-Sat). The hardness obtained for Max-2-Sat applies to monotone Max-2-Sat instances, meaning that we also obtain tight inapproximability for the Max-k-Vertex-Cover problem.


title: Mu2e calorimeter readout system

abstract: The Mu2e electromagnetic calorimeter is made of two disks of un-doped parallelepiped CsI crystals readout by SiPM. There are 674 crystals in one disk and each crystal is readout by an array of two SiPM. The readout electronics is composed of two types of modules: 1) the front-end module hosts the shaping amplifier and the high voltage linear regulator; since one front-end module is interfaced to one SiPM, a total of 2696 modules are needed for the entire calorimeter; 2) a waveform digitizer provides a further level of amplification and digitizes the SiPM signal at the sampling frequency of $200\ \text{M}\text{Hz}$ with 12-bits ADC resolution; since one board digitizes the data received from 20 SiPMs, a total of 136 boards are needed. The readout system operational conditions are hostile: ionization dose of $20\ \text{krads}$, neutron flux of $10^{12}\ \mathrm{n}(1\ \text{MeVeq})/\text{cm}^2$, magnetic field of $1\ \text{T}$ and in vacuum level of $10^{-4}\ \text{Torr}$. A description of the readout system and qualification tests is reported.


title: Some analytic results on interpolating sesqui-harmonic maps

abstract: In this article we study various analytic aspects of interpolating sesqui-harmonic maps between Riemannian manifolds where we mostly focus on the case of a spherical target. The latter are critical points of an energy functional that interpolates between the functionals for harmonic and biharmonic maps. In the case of a spherical target we will derive a conservation law and use it to show the smoothness of weak solutions. Moreover, we will obtain several classification results for interpolating sesqui-harmonic maps.


title: AMS-100: The Next Generation Magnetic Spectrometer in Space -- An International Science Platform for Physics and Astrophysics at Lagrange Point 2

abstract: The next generation magnetic spectrometer in space, AMS-100, is designed to have a geometrical acceptance of $100\,\mathrm{m}^2\,\mathrm{sr}$ and to be operated for at least ten years at the Sun-Earth Lagrange Point 2. Compared to existing experiments, it will improve the sensitivity for the observation of new phenomena in cosmic rays, and in particular in cosmic antimatter, by at least a factor of 1000. The magnet design is based on high temperature superconductor tapes, which allow the construction of a thin solenoid with a homogeneous magnetic field of 1 Tesla inside. The inner volume is instrumented with a silicon tracker reaching a maximum detectable rigidity of 100 TV and a calorimeter system that is 70 radiation lengths deep, equivalent to four nuclear interaction lengths, which extends the energy reach for cosmic-ray nuclei up to the PeV scale, i.e. beyond the cosmic-ray knee. Covering most of the sky continuously, AMS-100 will detect high-energy gamma rays in the calorimeter system and by pair conversion in the thin solenoid, reconstructed with excellent angular resolution in the silicon tracker.


title: Symmetries and reductions on the noncommutative Kadomtsev-Petviashvili and Gelfand-Dickey hierarchies

abstract: In this paper, we construct the additional flows of the noncommutative Kadomtsev-Petviashvili(KP) hierarchy and the additional symmetry flows constitute an infinite dimensional Lie algebra $W_{1+\infty}$. In addition, the generating function of the additional symmetries can also be proved to have a nice form in terms of wave functions and this generating symmetry is used to construct the noncommutative KP hierarchy with self-consistent sources and the constrained noncommutative KP hierarchy. The above results will be further generalized to the noncommutative Gelfand-Dickey hierarchies which contains many interesting noncommutative integrable systems such as the noncommutative KdV hierarchy and noncommutative Boussinesq hierarchy. Meanwhile, we construct two new noncommutative systems including odd noncommutative C type Gelfand-Dickey and even noncommutative C type Gelfand-Dickey hierarchies. Also using the symmetry, we can construct a new noncommutative Gelfand-Dickey hierarchy with self-consistent sources. Basing on the natural differential Lax operator of the noncommutative Gelfand-Dickey hierarchy, the string equations of the noncommutative Gelfand-Dickey hierarchy are also derived.


title: Using data-reduction techniques to analyse biomolecular trajectories

abstract: This chapter discusses the way in which dimensionality reduction algorithms such as diffusion maps and sketch-map can be used to analyze molecular dynamics trajectories. The first part discusses how these various algorithms function, as well as practical issues such as landmark selection and how these algorithms can be used when the data to be analyzed, comes from enhanced sampling trajectories. In the later parts, a comparison between the results obtained by applying various algorithms to two sets of sample data is performed and discussed. This section is then followed by a summary of how one algorithm, in particular, sketch-map, has been applied to a range of problems. The chapter concludes with a discussion on the directions that we believe this field is currently moving.


title: A Surface Radio Array for the Enhancement of IceTop and its Science Prospects

abstract: Radio detection of air showers in the current era has progressed immensely to effectively extract the properties of these air showers. Primary cosmic rays with energies of hundreds of PeV have been successfully measured with the method of radio detection. There are also attempts to observe high-energy neutrinos with this technique. Current radio experiments measuring cosmic-ray air showers mostly operate in the frequency range of 30-80 MHz. An optimization of the frequency band of operation can be done for maximizing the signal-to-noise ratio that can be achieved by an array of radio antennas at the South Pole, operated along with IceTop. Such an array can improve the reconstruction of air showers performed with IceTop. The prospect of using such an optimized radio array for measuring gamma rays of PeV energies from the Galactic Center is discussed.


title: Enumerating path diagrams in connection with $q$-tangent and $q$-secant numbers

abstract: We enumerate height-restricted path diagrams associated with $q$-tangent and $q$-secant numbers by considering convergents of continued fractions, leading to expressions involving basic hypergeometric functions. Our work generalises some results by M. Josuat-Verg\'es for unrestricted path diagrams [European Journal of Combinatorics 31 (2010) 1892].


title: Thermopower of the electron-doped manganese pnictide LaMnAsO

abstract: Upon chemical substitution of oxygen with fluor, LaMnAsO has been electron-doped in experiments, resulting in samples with remarkably high Seebeck coefficients of around -300$\mu VK^{-1}$ at room temperature and 3% doping. Within the framework of density functional theory plus dynamical mean-field theory we are not only able to reproduce these experimental observations, but can also provide a thorough investigation of the underlying mechanisms. By considering electronic correlations in the half-filled Mn-3d shells, we trace the high Seebeck coefficient back to an asymmetry in the spectral function, which is due to emergence of incoherent spectral weight under doping and a strong renormalisation of the unoccupied states. This is only possible in correlated systems and cannot be explained by DFT-based band structure calculations.


title: First-principles study of UO$_2$ lattice thermal-conductivity: A simple description

abstract: Modeling the high-$T$ paramagnetic state of bulk UO$_2$ by a non-spin-polarized calculation and neglecting the Hubbard-U correction for the $f$ electrons in U atoms, the lattice thermal conductivity of bulk UO$_2$ is investigated by the exact solution of the Boltzmann transport equation for the steady-state phonon distribution function. The results show that TA branches corresponding to U-atoms vibrations have the largest lifetimes and therefore have dominant role in thermal conductivity, while the optical branches corresponding mainly to O-atoms vibrations have the shortest lifetimes. Using this simple model, our results for the thermal conductivity show a very good agreement with the experiments. The calculations are repeated for bulk UO$_2$ with different U-235 concentrations of 3\%, 5\%, 7\%, and 20\%, and the results show a small decrease of thermal conductivity which arise from scattering of phonons by impurities.


title: A method for computing the Perron root for primitive matrices

abstract: Following the Perron-Frobenius theorem, the spectral radius of a primitive matrix is a simple eigenvalue. It is shown that for a primitive matrix $A$, there is a positive rank one matrix $X$ such that $B = A \circ X$, where $\circ$ denotes the Hadamard product of matrices, and such that the row (column) sums of matrix $B$ are the same and equal to the Perron root. An iterative algorithm is presented to obtain matrix $B$ without an explicit knowledge of $X$. The convergence rate of this algorithm is similar to that of the power method but it uses less computational load. A byproduct of the proposed algorithm is a new method for calculating the first eigenvector.


title: Multi-beam-energy operation for the continuous-wave x-ray free electron laser

abstract: The parallel operation of multiple undulator lines with a wide spectral range is an important way to increase the efficiency of x-ray free electron laser (XFEL) facilities, especially for machines with high-repetition-rate. In this paper, a delay system based on four double bend achromats is proposed to delay electron beams, thereby changing the arrival time of those delayed electron beams in the accelerating structure behind the system. Combined with kickers, the delay system can be used to generate bunch-to-bunch energy changed electron beams in a continuous wave XFEL facility. Start-to-end simulations based on the Shanghai high-repetition-rate XFEL and extreme light facility parameters are performed to demonstrate that the delay system can flexibly control electron beam energy from 1.48 to 8.74 GeV at the end of the linac.


title: S-folds and (non-)supersymmetric Janus solutions

abstract: The S-fold description of Janus-type solutions of type IIB supergravity is investigated. This is done by first studying a $\,\textrm{U}(1) \times \textrm{U}(1)\,$ invariant sector of the four-dimensional dyonically-gauged $\,{[\,\textrm{SO}(1,1) \times \textrm{SO}(6)\,] \ltimes \mathbb{R}^{12}}\,$ maximal supergravity that arises upon reduction of type IIB supergravity on $\,\mathbb{R} \, \times \, \textrm{S}^{5}\,$. Two AdS$_{4}$ solutions preserving $\,\textrm{SU}(3)\,$ and $\,\textrm{SO}(6)\,$ gauge symmetry together with $\,\mathcal{N}=1\,$ and $\,\mathcal{N}=0\,$ supersymmetry are found within this sector. Fetching techniques from the E$_{7(7)}$ exceptional field theory, these solutions are uplifted to ten-dimensional S-folds of type IIB Janus-type solutions of the form $\,\textrm{AdS}_{4} \times \mathbb{R} \times \textrm{M}_{5}\,$. The solutions presented here are natural candidates for the holographic duals of three-dimensional $\,\mathcal{N}=1\,$ and $\,\mathcal{N}=0\,$ interface super-Yang--Mills theories with $\,\textrm{SU}(3)\,$ and $\,\textrm{SU}(4)\,$ internal symmetry.


title: Lectures on Faster-than-Light Travel and Time Travel

abstract: These lecture notes were prepared for a 25-hour course for advanced undergraduate students participating in Perimeter Institute's Undergraduate Summer Program. The lectures cover some of what is currently known about the possibility of superluminal travel and time travel within the context of established science, that is, general relativity and quantum field theory. Previous knowledge of general relativity at the level of a standard undergraduate-level introductory course is recommended, but all the relevant material is included for completion and reference. No previous knowledge of quantum field theory, or anything else beyond the standard undergraduate curriculum, is required. Advanced topics in relativity, such as causal structures, the Raychaudhuri equation, and the energy conditions are presented in detail. Once the required background is covered, concepts related to faster-than-light travel and time travel are discussed. After introducing tachyons in special relativity as a warm-up, exotic spacetime geometries in general relativity such as warp drives and wormholes are discussed and analyzed, including their limitations. Time travel paradoxes are also discussed in detail, including some of their proposed resolutions.


title: Spot the Difference: Accuracy of Numerical Simulations via the Human Visual System

abstract: Comparative evaluation lies at the heart of science, and determining the accuracy of a computational method is crucial for evaluating its potential as well as for guiding future efforts. However, metrics that are typically used have inherent shortcomings when faced with the under-resolved solutions of real-world simulation problems. We show how to leverage crowd-sourced user studies in order to address the fundamental problems of widely used classical evaluation metrics. We demonstrate that such user studies, which inherently rely on the human visual system, yield a very robust metric and consistent answers for complex phenomena without any requirements for proficiency regarding the physics at hand. This holds even for cases away from convergence where traditional metrics often end up inconclusive results. More specifically, we evaluate results of different essentially non-oscillatory (ENO) schemes in different fluid flow settings. Our methodology represents a novel and practical approach for scientific evaluations that can give answers for previously unsolved problems.


title: Universality Classes of Stabilizer Code Hamiltonians

abstract: Stabilizer code quantum Hamiltonians have been introduced with the intention of physically realizing a quantum memory because of their resilience to decoherence. In order to analyze their finite temperature thermodynamics, we show how to generically solve their partition function using duality techniques. By unveiling each model's universality class and effective dimension, insights may be gained on their finite temperature dynamics and robustness. Our technique is demonstrated in particular on the 4D Toric Code and Haah's Code -- we find that the former falls into the 4D Ising universality class, whereas Haah's Code exhibits dimensional reduction and falls into the 1D Ising universality class.


title: Resource theory of entanglement for bipartite quantum channels

abstract: The traditional perspective in quantum resource theories concerns how to use free operations to convert one resourceful quantum state to another one. For example, a fundamental and well known question in entanglement theory is to determine the distillable entanglement of a bipartite state, which is equal to the maximum rate at which fresh Bell states can be distilled from many copies of a given bipartite state by employing local operations and classical communication for free. It is the aim of this paper to take this kind of question to the next level, with the main question being: What is the best way of using free channels to convert one resourceful quantum channel to another? Here we focus on the the resource theory of entanglement for bipartite channels and establish several fundamental tasks and results regarding it. In particular, we establish bounds on several pertinent information processing tasks in channel entanglement theory, and we define several entanglement measures for bipartite channels, including the logarithmic negativity and the $\kappa$-entanglement. We also show that the max-Rains information of [B\"auml et al., Physical Review Letters, 121, 250504 (2018)] has a divergence interpretation, which is helpful for simplifying the results of this earlier work.


title: 24 rational curves on K3 surfaces

abstract: Given d in IN, we prove that all smooth K3 surfaces (over any field of characteristic p other than 2,3) of degree greater than 84d^2 contain at most 24 rational curves of degree at most d. In the exceptional characteristics, the same bounds hold for non-unirational K3 surfaces, and we develop analogous results in the unirational case. For d at least 3, we also construct K3 surfaces of any degree greater than 4d(d+1) with 24 rational curves of degree exactly d, thus attaining the above bounds.


title: Ricci-Based Gravity theories and their impact on Maxwell and nonlinear electromagnetic models

abstract: We extend the correspondence between metric-affine Ricci-Based Gravity theories and General Relativity (GR) to the case in which the matter sector is represented by linear and nonlinear electromagnetic fields. This complements previous studies focused on fluids and scalar fields. We establish the general algorithm that relates the matter fields in the GR and RBG frames and consider some applications. In particular, we find that the so-called Eddington-inspired Born-Infeld gravity theory coupled to Maxwell electromagnetism is in direct correspondence with GR coupled to Born-Infeld electromagnetism. We comment on the potential phenomenological implications of this relation.


title: Evolution of ONeMg Core in Super-AGB Stars towards Electron-Capture Supernovae: Effects of Updated Electron-Capture Rate

abstract: Stars with $\sim 8-10~{M}_{\odot}$ evolve to form a strongly degenerate ONeMg core. When the core mass becomes close to the Chandrasekhar mass, the core undergoes electron captures on $^{24}$Mg and $^{20}$Ne, which induce the electron-capture supernova (ECSN). In order to clarify whether the ECSN leads to a collapse or thermonuclear explosion, we calculate the evolution of an $8.4~M_\odot$ star from the main sequence until the oxygen ignition in the ONeMg core. We apply the latest electron-capture rate on $^{20}$Ne including the second forbidden transition, and investigate how the location of the oxygen ignition (center or off-center) and the $Y_e$ distribution depend on the input physics and the treatment of the semiconvection and convection. The central density when the oxygen deflagration is initiated, $\rho_{\rm c,def}$, can be significantly higher than that of the oxygen ignition thanks to the convection, and we estimate $\log_{10}(\rho_{\rm c,def}/\mathrm{g~cm^{-3}})>10.10$. We perform two-dimensional simulations of the flame propagation to examine how the final fate of the ONeMg core depends on the $Y_e$ distribution and $\rho_{\rm c,def}$. We find that the deflagration starting from $\log_{10}(\rho_{\rm c,def}/\mathrm{g~cm^{-3}})>10.01 (< 10.01)$ leads to a collapse (thermonuclear explosion). Since our estimate of $\rho_{\rm c,def}$ exceeds this critical value, the ONeMg core is likely to collapse, although further studies of the convection and semiconvection before the deflagration are important.


title: Predictively Consistent Prior Effective Sample Sizes

abstract: Determining the sample size of an experiment can be challenging, even more so when incorporating external information via a prior distribution. Such information is increasingly used to reduce the size of the control group in randomized clinical trials. Knowing the amount of prior information, expressed as an equivalent prior effective sample size (ESS), clearly facilitates trial designs. Various methods to obtain a prior's ESS have been proposed recently. They have been justified by the fact that they give the standard ESS for one-parameter exponential families. However, despite being based on similar information-based metrics, they may lead to surprisingly different ESS for non-conjugate settings, which complicates many designs with prior information. We show that current methods fail a basic predictive consistency criterion, which requires the expected posterior-predictive ESS for a sample of size $N$ to be the sum of the prior ESS and $N$. The expected local-information-ratio ESS is introduced and shown to be predictively consistent. It corrects the ESS of current methods, as shown for normally distributed data with a heavy-tailed Student-t prior and exponential data with a generalized Gamma prior. Finally, two applications are discussed: the prior ESS for the control group derived from historical data, and the posterior ESS for hierarchical subgroup analyses.


title: Bandwidth Is Not Enough: "Hidden" Outlier Noise and Its Mitigation

abstract: In addition to ever-present thermal noise, communication and sensor systems can contain significant amounts of interference with outlier (e.g. impulsive) characteristics. Such outlier interference (including that caused by nonlinear signal distortions, e.g. clipping) can be efficiently mitigated in real-time using intermittently nonlinear filters. Depending on the interference nature and composition, improvements in the quality of the signal of interest achieved by such filtering will vary from "no harm" to substantial. In this tutorial, we explain in detail why the underlying outlier nature of interference often remains obscured, discussing the many challenges and misconceptions associated with state-of-art analog and/or digital nonlinear mitigation techniques, especially when addressing complex practical interference scenarios. We then focus on the methodology and tools for real-time outlier noise mitigation, demonstrating how the "excess band" observation of outlier noise enables its efficient in-band mitigation. We introduce the basic real-time nonlinear components that are used for outlier noise filtering and provide examples of their implementation. We further describe complementary nonlinear filtering arrangements for wide- and narrow-band outlier noise reduction, providing illustrations of their performance and the effect on channel capacity. Finally, we outline "effectively analog" digital implementations of these filtering structures, discuss their broader applications, and comment on the ongoing development of the platform for their demonstration and testing. To emphasize the effectiveness and versatility of this approach, in our examples we use particularly challenging waveforms that severely obscure low-amplitude outlier noise, such as broadband chirp signals (e.g. used in radar, sonar, and spread-spectrum communications) and "bursty," high crest factor signals (e.g. OFDM).


title: Magnetic dipole ordering in resonant dielectric metasurfaces

abstract: Artificial magnetism at optical frequencies can be realized in metamaterials composed of periodic arrays of subwavelength elements, also called "meta-atoms". Optically-induced magnetic moments can be arranged in both unstaggered structures, naturally associated with ferromagnetic (FM) order, or staggered structures, linked correspondingly to antiferromagnetic (AFM) order. Here we demonstrate that such magnetic dipole orders of the lattices of meta-atoms can appear in low-symmetry Mie-resonant metasurfaces where each asymmetric dielectric (non-magnetic) meta-atom supports a localized trapped mode. We reveal that these all-dielectric resonant metasurfaces possess not only strong optical magnetic response but also they demonstrate a significant polarization rotation of the propagating electromagnetic waves at both FM and AFM resonances. We confirm these findings experimentally by measuring directly the spectral characteristics of different modes excited in all-dielectric metasurfaces, and mapping near-field patterns of the electromagnetic fields at the microwave frequencies.


title: Geometric Matrix Midranges

abstract: We define geometric matrix midranges for positive definite Hermitian matrices and study the midrange problem from a number of perspectives. Special attention is given to the midrange of two positive definite matrices before considering the extension of the problem to $N > 2$ matrices. We compare matrix midrange statistics with the scalar and vector midrange problem and note the special significance of the matrix problem from a computational standpoint. We also study various aspects of geometric matrix midrange statistics from the viewpoint of linear algebra, differential geometry and convex optimization.


title: Properties of Nanocrystalline Silicon Probed by Optomechanics

abstract: Nanocrystalline materials exhibit properties that can differ substantially from those of their single crystal counterparts. As such, they provide ways to enhance and optimise their functionality for devices and applications. Here we report on the optical, mechanical and thermal properties of nanocrystalline silicon probed by means of optomechanical nanobeams to extract information of the dynamics of optical absorption, mechanical losses, heat generation and dissipation. The optomechanical nanobeams are fabricated using nanocrystalline films prepared by annealing amorphous silicon layers at different temperatures. The resulting crystallite sizes and the stress in the films can be controlled by the annealing temperature and time and, consequently, the properties of the films can be tuned relatively freely, as demonstrated here by means of electron microscopy and Raman scattering. We show that the nanocrystallite size and the volume fraction of the grain boundaries play a key role in the dissipation rates through non-linear optical and thermal processes. Promising optical (13000) and mechanical (1700) quality factors were found in the optomechanical cavity realised in the nanocrystalline Si resulting from annealing at 950 C. The enhanced absorption and recombination rates via the intra-gap states and the reduced thermal conductivity boost the potential to exploit these non-linear effects in applications, including NEMS, phonon lasing and chaos-based devices.


title: A new hybrid genetic algorithm for protein structure prediction on the 2D triangular lattice

abstract: The flawless functioning of a protein is essentially linked to its own three-dimensional structure. Therefore, the prediction of a protein structure from its amino acid sequence is a fundamental problem in many fields that draws researchers attention. This problem can be formulated as a combinatorial optimization problem based on simplified lattice models such as the hydrophobic-polar model. In this paper, we propose a new hybrid algorithm combining three different well-known heuristic algorithms: genetic algorithm, tabu search strategy and local search algorithm in order to solve the PSP problem. Regarding the assessment of suggested algorithm, an experimental study is included, where we considered the quality of the produced solution as the main quality criterion. Furthermore, we compared the suggested algorithm with state-of-the-art algorithms using a selection of well-studied benchmark instances.


title: Belief places and spaces: Mapping cognitive environments

abstract: Beliefs are not facts, but they are factive - they feel like facts. This property is what can make misinformation dangerous. Being able to deliberately navigate through a landscape of often conflicting factive statements is difficult when there is no way to show the relationships between them without incorporating the information in linear, narrative forms. In this paper, we present a mechanism to produce maps of belief places, where populations agree on salient features of fictional environments, and belief spaces, where subgroups have related but distinct perspectives. Using a model developed using agent-based simulation, we show that by observing the repeated behaviors of human participants in the same social context, it is possible to build maps that show the shared narrative environment overlaid with traces that show unique, individual or subgroup perspectives. Our contribution is a proof-of-concept system, based on the affordances of fantasy tabletop role-playing games, which support multiple groups interacting with the same dungeon in a controlled, online environment. The techniques used in this process are mathematically straightforward, and should be generalizable to auto-generating larger-scale maps of belief spaces from other corpora, such as discussions on social media.


title: Bulk viscosity of baryonic matter with trapped neutrinos

abstract: We study bulk viscosity arising from weak current Urca processes in dense baryonic matter at and beyond nuclear saturation density. We consider the temperature regime where neutrinos are trapped and therefore have nonzero chemical potential. We model the nuclear matter in a relativistic density functional approach, taking into account the trapped neutrino component. We find that the resonant maximum of the bulk viscosity would occur at or below the neutrino trapping temperature, so in the neutrino trapped regime the bulk viscosity decreases with temperature as $T^{-2}$, this decrease being interrupted by a drop to zero at a special temperature where the proton fraction becomes density-independent and the material scale-invariant. The bulk viscosity is larger for matter with lower lepton fraction, i.e., larger isospin asymmetry. We find that bulk viscosity in the neutrino-trapped regime is smaller by several orders of magnitude than in the neutrino-transparent regime, which implies that bulk viscosity in neutrino-trapped matter is probably not strong enough to affect significantly the evolution of neutron star mergers. This also implies weak damping of gravitational waves emitted by the oscillations of the postmerger remnant in the high-temperature, neutrino-trapped phase of evolution.


title: Modelling Levy space-time white noises

abstract: Based on the theory of independently scattered random measures, we introduce a natural generalisation of Gaussian space-time white noise to a Levy-type setting, which we call Levy-valued random measures. We determine the subclass of cylindrical Levy processes which correspond to Levy-valued random measures, and describe the elements of this subclass uniquely by their characteristic function. We embed the Levy-valued random measure, or the corresponding cylindrical Levy process, in the space of general and tempered distributions. For the latter case, we show that this embedding is possible if and only if a certain integrability condition is satisfied. Similar to existing definitions, we introduce Levy-valued additive sheets, and show that integrating a Levy-valued random measure in space defines a Levy-valued additive sheet. This relation is manifested by the result, that a Levy-valued random measure can be viewed as the weak derivative of a Levy-valued additive sheet in the space of distributions.


title: Adaptive Exploration for Unsupervised Person Re-Identification

abstract: Due to domain bias, directly deploying a deep person re-identification (re-ID) model trained on one dataset often achieves considerably poor accuracy on another dataset. In this paper, we propose an Adaptive Exploration (AE) method to address the domain-shift problem for re-ID in an unsupervised manner. Specifically, in the target domain, the re-ID model is inducted to 1) maximize distances between all person images and 2) minimize distances between similar person images. In the first case, by treating each person image as an individual class, a non-parametric classifier with a feature memory is exploited to encourage person images to move far away from each other. In the second case, according to a similarity threshold, our method adaptively selects neighborhoods for each person image in the feature space. By treating these similar person images as the same class, the non-parametric classifier forces them to stay closer. However, a problem of the adaptive selection is that, when an image has too many neighborhoods, it is more likely to attract other images as its neighborhoods. As a result, a minority of images may select a large number of neighborhoods while a majority of images have only a few neighborhoods. To address this issue, we additionally integrate a balance strategy into the adaptive selection. We evaluate our methods with two protocols. The first one is called "target-only re-ID", in which only the unlabeled target data is used for training. The second one is called "domain adaptive re-ID", in which both the source data and the target data are used during training. Experimental results on large-scale re-ID datasets demonstrate the effectiveness of our method. Our code has been released at https://github.com/dyh127/Adaptive-Exploration-for-Unsupervised-Person-Re-Identification.


title: Solution Landscapes in the Landau-de Gennes Theory on Rectangles

abstract: We study nematic equilibria on rectangular domains, in a reduced two-dimensional Landau-de Gennes framework. These reduced equilibria carry over to the three-dimensional framework at a special temperature. There is one essential model variable---$\epsilon$ which is a geometry-dependent and material-dependent variable. We compute the limiting profiles exactly in two distinguished limits---the $\epsilon \to 0$ limit relevant for macroscopic domains and the $\epsilon \to \infty$ limit relevant for nano-scale domains. The limiting profile has line defects near the shorter edges in the $\epsilon \to \infty$ limit whereas we observe fractional point defects in the $\epsilon \to 0$ limit. The analytical studies are complemented by bifurcation diagrams for these reduced equilibria as a function of $\epsilon$ and the rectangular aspect ratio. We also introduce the concept of `non-trivial' topologies and the relaxation of non-trivial topologies to trivial topologies mediated via point and line defects, with potential consequences for non-equilibrium phenomena and switching dynamics.


title: Finite temperature fermionic condensate in a conical space with a circular boundary and magnetic flux

abstract: We investigate the edge effects on the finite temperature fermionic condensate (FC) for a massive fermionic field in a (2+1)-dimensional conical spacetime with a magnetic flux located at the cone apex. The field obeys the bag boundary condition on a circle concentric with the apex. The analysis is presented for both the fields realizing two irreducible representations of the Clifford algebra and for general case of the chemical potential. In both the regions outside and inside the circular boundary, the FC is decomposed into the boundary-free and boundary-induced contributions. They are periodic functions of the magnetic flux with the period equal to the flux quantum and even functions under the simultaneous change of the signs for the magentic flux and the chemical potential. The dependence of the FC on the magnetic flux becomes weaker with decreasing planar angle deficit. For points near the boundary, the effects of finite temperature, of planar angle deficit and of magnetic flux are weak. For a fixed distance from the boundary and at high temperatures the FC is dominated by the Minkowskian part. The FC in parity and time-reversal symmetric (2+1)-dimensional fermionic models is discussed and applications are given to graphitic cones.


title: Attending to Emotional Narratives

abstract: Attention mechanisms in deep neural networks have achieved excellent performance on sequence-prediction tasks. Here, we show that these recently-proposed attention-based mechanisms---in particular, the Transformer with its parallelizable self-attention layers, and the Memory Fusion Network with attention across modalities and time---also generalize well to multimodal time-series emotion recognition. Using a recently-introduced dataset of emotional autobiographical narratives, we adapt and apply these two attention mechanisms to predict emotional valence over time. Our models perform extremely well, in some cases reaching a performance comparable with human raters. We end with a discussion of the implications of attention mechanisms to affective computing.


title: Sequence-to-Sequence Natural Language to Humanoid Robot Sign Language

abstract: This paper presents a study on natural language to sign language translation with human-robot interaction application purposes. By means of the presented methodology, the humanoid robot TEO is expected to represent Spanish sign language automatically by converting text into movements, thanks to the performance of neural networks. Natural language to sign language translation presents several challenges to developers, such as the discordance between the length of input and output data and the use of non-manual markers. Therefore, neural networks and, consequently, sequence-to-sequence models, are selected as a data-driven system to avoid traditional expert system approaches or temporal dependencies limitations that lead to limited or too complex translation systems. To achieve these objectives, it is necessary to find a way to perform human skeleton acquisition in order to collect the signing input data. OpenPose and skeletonRetriever are proposed for this purpose and a 3D sensor specification study is developed to select the best acquisition hardware.


title: Josephson Detection of Time Reversal Symmetry Broken Superconductivity in SnTe Nanowires

abstract: Exotic superconductors, such as high T$_C$, topological, and heavy-fermion superconductors, require phase sensitive measurements to determine the underlying pairing. Here we investigate the proximity-induced superconductivity in nanowires of SnTe, where an $s\pm is^{\prime}$ superconducting state is produced that lacks the time-reversal and valley-exchange symmetry of the parent SnTe. This effect, in conjunction with a ferroelectric distortion of the lattice at low temperatures, results in a marked alteration of the properties of Josephson junctions fabricated using SnTe nanowires. This work establishes the existence of a ferroelectric transition in SnTe nanowires and elucidates the role of ferroelectric domain walls on the flow of supercurrent through SnTe weak links. We detail two unique characteristics of these junctions: an asymmetric critical current in the DC Josephson effect and a prominent second harmonic in the AC Josephson effect. Each reveals the broken time-reversal symmetry in the junction. The novel $s\pm is^{\prime}$ superconductivity and the new Josephson effects can be used to investigate fractional vortices [1,2], topological superconductivity in multiband materials [3-5], and new types of Josephson-based devices in proximity-induced multiband and ferroelectric superconductors [6,7].


title: The admissibility theorem for the spatial X-ray transform over the two element field

abstract: We consider the Radon transform along lines in an $n$ dimensional vector space over the two element field. It is well known that this transform is injective and highly overdetermined. We classify the minimal collections of lines for which the restricted Radon transform is also injective. This is an instance of I.M.~Gelfand's {\it admissibility problem}. The solution is in stark contrast to the more uniform cases of the affine hyperplane transform and the projective line transform, which are addressed in other papers, \cite{Feld-G,Gr1}. The presentation here is intended to be widely accessible, requiring minimum background.


title: Thompson Sampling for Combinatorial Network Optimization in Unknown Environments

abstract: Influence maximization, adaptive routing, and dynamic spectrum allocation all require choosing the right action from a large set of alternatives. Thanks to the advances in combinatorial optimization, these and many similar problems can be efficiently solved given an environment with known stochasticity. In this paper, we take this one step further and focus on combinatorial optimization in unknown environments. We consider a very general learning framework called combinatorial multi-armed bandit with probabilistically triggered arms and a very powerful Bayesian algorithm called Combinatorial Thompson Sampling (CTS). Under the semi-bandit feedback model and assuming access to an oracle without knowing the expected base arm outcomes beforehand, we show that when the expected reward is Lipschitz continuous in the expected base arm outcomes CTS achieves $O(\sum_{i =1}^m\log T/(p_i\Delta_i))$ regret and $O(\max\{\mathbb{E}[m\sqrt{T\log T/p^*}],\mathbb{E}[m^2/p^*]\})$ Bayesian regret, where $m$ denotes the number of base arms, $p_i$ and $\Delta_i$ denote the minimum non-zero triggering probability and the minimum suboptimality gap of base arm $i$ respectively, $T$ denotes the time horizon, and $p^*$ denotes the overall minimum non-zero triggering probability. We also show that when the expected reward satisfies the triggering probability modulated Lipschitz continuity, CTS achieves $O(\max\{m\sqrt{T\log T},m^2\})$ Bayesian regret, and when triggering probabilities are non-zero for all base arms, CTS achieves $O(1/p^*\log(1/p^*))$ regret independent of the time horizon. Finally, we numerically compare CTS with algorithms based on upper confidence bounds in several networking problems and show that CTS outperforms these algorithms by at least an order of magnitude in majority of the cases.


title: Variational Inference MPC for Bayesian Model-based Reinforcement Learning

abstract: In recent studies on model-based reinforcement learning (MBRL), incorporating uncertainty in forward dynamics is a state-of-the-art strategy to enhance learning performance, making MBRLs competitive to cutting-edge model free methods, especially in simulated robotics tasks. Probabilistic ensembles with trajectory sampling (PETS) is a leading type of MBRL, which employs Bayesian inference to dynamics modeling and model predictive control (MPC) with stochastic optimization via the cross entropy method (CEM). In this paper, we propose a novel extension to the uncertainty-aware MBRL. Our main contributions are twofold: Firstly, we introduce a variational inference MPC, which reformulates various stochastic methods, including CEM, in a Bayesian fashion. Secondly, we propose a novel instance of the framework, called probabilistic action ensembles with trajectory sampling (PaETS). As a result, our Bayesian MBRL can involve multimodal uncertainties both in dynamics and optimal trajectories. In comparison to PETS, our method consistently improves asymptotic performance on several challenging locomotion tasks.


title: Lepton mixing patterns from $PSL_2(7)$ with a generalised CP symmetry

abstract: Lepton mixing patterns from the modular group $PSL_2(7)$ with generalised CP symmetries are studied. The residual symmetries in both charged leptons and neutrinos sector are $Z_{2}\times CP$. Seven types of mixing patterns at the $3\sigma$ level of the new global fit data are obtained. Among these patterns, three types of patterns can give the Dirac CP phase which is in the $1\sigma$ range of the global fit data. The effective mass of neutrinoless double-beta decay for these patterns are also examined.


title: Linking the pseudo-gap in the cuprates with local symmetry breaking: a commentary

abstract: Perspective on S. Mukhopadhyay, et al. Evidence for a vestigial nematic state in the cuprate pseudogap phase, PNAS 116, 13249 (2019); arXiv:1904.00915.


title: On the space of ultradistributions vanishing at infinity

abstract: We study the structural and linear topological properties of the space $\dot{\mathcal{B}}^{\prime \ast}_{\omega}$ of ultradistributions vanishing at infinity (with respect to a weight function $\omega$). Particularly, we show the first structure theorem for $\dot{\mathcal{B}}^{\prime \ast}_{\omega}$ under weaker hypotheses than were known so far. As an application, we determine the structure of the S-asymptotic behavior of ultradistributions.


abstract: The Game of Poker Chips, Dominoes and Survival fosters team building and high level cooperation in large groups, and is a tool applied in management training exercises. Each player, initially given two colored poker chips, is allowed to make exchanges with the game coordinator according to two rules, and must secure a domino before time is called in order to `survive'. Though the rules are simple, it is not evident by their form that the survival of the entire group requires that they cooperate at a high level. From the point of view of the game coordinator, the difficulty of the game for the group can be controlled not only by the time limit, but also by the initial distribution of chips, in a way we make precise by a time complexity type argument. That analysis also provides insight into good strategies for group survival, those taking the least amount of time. In addition, coordinators may also want to be aware of when the game is `solvable', that is, when their initial distribution of chips permits the survival of all group members if given sufficient time to make exchanges. It turns out that the game is solvable if and only if the initial distribution contains seven chips that have one of two particular color distributions. In addition to being a lively game to play in management training or classroom settings, the analysis of the game after play can make for an engaging exercise in any basic discrete mathematics course to give a basic introduction to elements of game theory, logical reasoning, number theory and the computation of algorithmic complexities.


title: Resonant charge-transfer in grazing collisions of H$^-$ with vicinal nanosurfaces on Cu(111), Au(100) and Pd(111) substrates: A comparative study

abstract: We compare the electron dynamics at monocrystalline Cu(111), Au(100) and Pd(111) precursor substrates with vicinal nanosteps. The unoccupied bands of a surface superlattice are populated \textit{via} the resonant charge transfer (RCT) between the surface and a H$^-$ ion that flies by at grazing angles. A quantum mechanical wave packet propagation approach is utilized to simulate the motion of the active electron where time-evolved wave packet densities are used to visualize the dynamics through the superlattice. The survived ion fraction in the reflected beam generally exhibits modulations as a function of the vicinal terrace size and shows peaks at those energies that access the image state subband dispersions. However, differences in magnitudes of the ion-survival as a function of the particular substrate selection as well as the ion-surface interaction time based on the choice of two ion-trajectories are examined. A square well model producing standing waves between the steps on the surface explains well the energies of the maxima in the ion survival probability for all the metals considered, indicating that the primary process of confinement induced subband formation is rather robust. The work may motivate measurements and applications of shallow-angle ion-scattering spectroscopy to access electronic substructures in periodically nanostructured surfaces.


title: A Direct Construction of Primitive Formally Dual Pairs Having Subsets with Unequal Sizes

abstract: The concept of formal duality was proposed by Cohn, Kumar and Sch\"urmann, which reflects a remarkable symmetry among energy-minimizing periodic configurations. This formal duality was later translated into a purely combinatorial property by Cohn, Kumar, Reiher and Sch\"urmann, where the corresponding combinatorial objects were called formally dual pairs. So far, except the results presented in Li and Pott (arXiv:1810.05433v3), we have little information about primitive formally dual pairs having subsets with unequal sizes. In this paper, we propose a direct construction of primitive formally dual pairs having subsets with unequal sizes in $\mathbb{Z}_2 \times \mathbb{Z}_4^{2m}$, where $m \ge 1$. This construction recovers an infinite family obtained in Li and Pott (arXiv:1810.05433v3), which was derived by employing a recursive approach. Although the resulting infinite family was known before, the idea of the direct construction is new and provides more insights which were not known from the recursive approach.


title: Maximum principle for stochastic optimal control problem of finite state forward-backward stochastic difference systems

abstract: In this paper, we study the maximum principle for stochastic optimal control problems of forward-backward stochastic difference systems (FBS{\Delta}Ss) where the uncertainty is modeled by a discrete time, finite state process, rather than white noises. Two types of FBS{\Delta}Ss are investigated. The first one is described by a partially coupled forward-backward stochastic difference equation (FBS{\Delta}E) and the second one is described by a fully coupled FBS{\Delta}E. By adopting an appropriate representation of the product rule and an appropriate formulation of the backward stochastic difference equation (BS{\Delta}E), we deduce the adjoint difference equation. Finally, the maximum principle for this optimal control problem with the control domain being convex is established.


title: Spherical nano-inhomogeneity with the Steigmann-Ogden interface model under general uniform far-field stress loading

abstract: An explicit solution, considering the interface bending resistance as described by the Steigmann-Ogden interface model, is derived for the problem of a spherical nano-inhomogeneity (nanoscale void/inclusion) embedded in an infinite linear-elastic matrix under a general uniform far-field-stress (including tensile and shear stresses). The Papkovich-Neuber (P-N) general solutions, which are expressed in terms of spherical harmonics, are used to derive the analytical solution. A superposition technique is used to overcome the mathematical complexity brought on by the assumed interfacial residual stress in the Steigmann-Ogden interface model. Numerical examples show that the stress field, considering the interface bending resistance as with the Steigmann-Ogden interface model, differs significantly from that considering only the interface stretching resistance as with the Gurtin-Murdoch interface model. In addition to the size-dependency, another interesting phenomenon is observed: some stress components are invariant to interface bending stiffness parameters along a certain circle in the inclusion/matrix. Moreover, a characteristic line for the interface bending stiffness parameters is presented, near which the stress concentration becomes quite severe. Finally, the derived analytical solution with the Steigmann-Ogden interface model is provided in the supplemental MATLAB code, which can be easily executed, and used as a benchmark for semi-analytical solutions and numerical solutions in future studies.


title: Universal One-Dimensional Cellular Automata Derived for Turing Machines and its Dynamical Behaviour

abstract: Universality in cellular automata theory is a central problem studied and developed from their origins by John von Neumann. In this paper, we present an algorithm where any Turing machine can be converted to one-dimensional cellular automaton with a 2-linear time and display its spatial dynamics. Three particular Turing machines are converted in three universal one-dimensional cellular automata, they are: binary sum, rule 110 and a universal reversible Turing machine.


title: On a method to construct exponential families by representation theory

abstract: Exponential family plays an important role in information geometry. In arXiv:1811.01394, we introduced a method to construct an exponential family $\mathcal{P}=\{p_\theta\}_{\theta\in\Theta}$ on a homogeneous space $G/H$ from a pair $(V,v_0)$. Here $V$ is a representation of $G$ and $v_0$ is an $H$-fixed vector in $V$. Then the following questions naturally arise: (Q1) when is the correspondence $\theta\mapsto p_\theta$ injective? (Q2) when do distinct pairs $(V,v_0)$ and $(V',v_0')$ generate the same family? In this paper, we answer these two questions (Theorems 1 and 2). Moreover, in Section 3, we consider the case $(G,H)=(\mathbb{R}_{>0}, \{1\})$ with a certain representation on $\mathbb{R}^2$. Then we see the family obtained by our method is essentially generalized inverse Gaussian distribution (GIG).


title: Dynamic Real-time Optimization of Batch Processes using Pontryagin's Minimum Principle and Set-membership Adaptation

abstract: This paper studies a dynamic real-time optimization in the context of model-based time-optimal operation of batch processes under parametric model mismatch. In order to tackle the model-mismatch issue, a receding-horizon policy is usually followed with frequent re-optimization. The main problem addressed in this study is the high computational burden that is usually required by such schemes. We propose an approach that uses parameterized conditions of optimality in the adaptive predictive-control fashion. The uncertainty in the model predictions is treated explicitly using reachable sets that are projected into the optimality conditions. Adaptation of model parameters is performed online using set-membership estimation. A class of batch membrane separation processes is in the scope of the presented applications, where the benefits of the presented approach are outlined.


title: Entropic Regularization of Markov Decision Processes

abstract: An optimal feedback controller for a given Markov decision process (MDP) can in principle be synthesized by value or policy iteration. However, if the system dynamics and the reward function are unknown, a learning agent must discover an optimal controller via direct interaction with the environment. Such interactive data gathering commonly leads to divergence towards dangerous or uninformative regions of the state space unless additional regularization measures are taken. Prior works proposed bounding the information loss measured by the Kullback-Leibler (KL) divergence at every policy improvement step to eliminate instability in the learning dynamics. In this paper, we consider a broader family of $f$-divergences, and more concretely $\alpha$-divergences, which inherit the beneficial property of providing the policy improvement step in closed form at the same time yielding a corresponding dual objective for policy evaluation. Such entropic proximal policy optimization view gives a unified perspective on compatible actor-critic architectures. In particular, common least-squares value function estimation coupled with advantage-weighted maximum likelihood policy improvement is shown to correspond to the Pearson $\chi^2$-divergence penalty. Other actor-critic pairs arise for various choices of the penalty-generating function $f$. On a concrete instantiation of our framework with the $\alpha$-divergence, we carry out asymptotic analysis of the solutions for different values of $\alpha$ and demonstrate the effects of the divergence function choice on common standard reinforcement learning problems.


title: Chimera states for a globally coupled sine circle map lattice: spatiotemporal intermittency and hyperchaos

abstract: We study the existence of chimera states, i.e. mixed states, in a globally coupled sine circle map lattice, with different strengths of inter-group and intra-group coupling. We find that at specific values of the parameters of the CML, a completely random initial condition evolves to chimera states, having a phase synchronised and a phase desynchronised group, where the space time variation of the phases of the maps in the desynchronised group shows structures similar to spatiotemporally intermittent regions. Using the complex order parameter we obtain a phase diagram that identifies the region in the parameter space which supports chimera states of this type, as well as other types of phase configurations such as globally phase synchronised states, two phase clustered states and fully phase desynchronised states. We estimate the volume of the basin of attraction of each kind of solution. The STI chimera region is studied in further detail via numerical and analytic stability analysis, and the Lyapunov spectrum is calculated. This state is identified to be hyperchaotic as the two largest Lyapunov exponents are found to be positive. The distributions of laminar and burst lengths in the incoherent region of the chimera show exponential behaviour. The average fraction of laminar/burst sites is identified to be the important quantity which governs the dynamics of the chimera. After an initial transient, these settle to steady values which can be used to reproduce the phase diagram in the chimera regime.


title: Efficient quasi-phase-matched frequency conversion in a lithium niobate racetrack microresonator

abstract: We demonstrate efficient second harmonic generation in a quasi-phase-matched, high quality factor ($Q_0 \approx 5.3\times 10^5$) racetrack microresonator. The observed normalized conversion efficiency is about $3.8\%~mW^{-1}$.


title: Streaming 1.9 Billion Hypersparse Network Updates per Second with D4M

abstract: The Dynamic Distributed Dimensional Data Model (D4M) library implements associative arrays in a variety of languages (Python, Julia, and Matlab/Octave) and provides a lightweight in-memory database implementation of hypersparse arrays that are ideal for analyzing many types of network data. D4M relies on associative arrays which combine properties of spreadsheets, databases, matrices, graphs, and networks, while providing rigorous mathematical guarantees, such as linearity. Streaming updates of D4M associative arrays put enormous pressure on the memory hierarchy. This work describes the design and performance optimization of an implementation of hierarchical associative arrays that reduces memory pressure and dramatically increases the update rate into an associative array. The parameters of hierarchical associative arrays rely on controlling the number of entries in each level in the hierarchy before an update is cascaded. The parameters are easily tunable to achieve optimal performance for a variety of applications. Hierarchical arrays achieve over 40,000 updates per second in a single instance. Scaling to 34,000 instances of hierarchical D4M associative arrays on 1,100 server nodes on the MIT SuperCloud achieved a sustained update rate of 1,900,000,000 updates per second. This capability allows the MIT SuperCloud to analyze extremely large streaming network data sets.


title: Spins of black holes in coalescing compact binaries

abstract: Modern astrophysical methods of determination of spins of rotating stellar-mass black hole in close binaries and of supermassive black holes in active galactic nuclei are briefly discussed. Effective spins of coalescing binary black holes derived from LIGO/Virgo gravitational wave observations are specially addressed. We consider three types of coalescing binaries: double black holes, black hole-neutron star binaries and primordial double black holes. The effective spins of coalescing astrophysical binary black holes and black holes with neutron stars are calculated for two plausible models of black hole formations from stellar core collapses (without or with additional fallback from the stellar envelope) taking into account the stellar metallicity and star formation rate evolution in the Universe. The calculated distributions do not contradict the reported LIGO/Virgo observations. The effective spins of primordial coalescing stellar-mass black holes can reach a few per cent due to accretion spin-up in the cold external medium.


title: Centrally generated primitive ideals of $U(\mathfrak{n})$ for exceptional types

abstract: Let $\mathfrak{g}$ be a complex semisimple Lie algebra, $\mathfrak{b}$ be a Borel subalgebra of $\mathfrak{g}$, $\mathfrak{n}$ be the nilradical of $\mathfrak{b}$, and $U(\mathfrak{n})$ be the universal enveloping algebra of $\mathfrak{n}$. We study primitive ideals of $U(\mathfrak{n})$. Almost all primitive ideals are centrally generated, i.e., are generated by their intersections with the center $Z(\mathfrak{n})$ of $U(\mathfrak{n})$. We present an explicit characterization of the centrally generated primitive ideals of $U(\mathfrak{n})$ in terms of the Dixmier map and the Kostant cascade in the case when $\mathfrak{g}$ is a simple algebra of exceptional type. (For classical simple Lie algebras, a similar characterization was obtained by Ivan Penkov and the first author.) As a corollary, we establish a classification of centrally generated primitive ideals of $U(\mathfrak{n})$ for an arbitrary semisimple algebra $\mathfrak{g}$.