title: Kinematics of filaments in cooling flow clusters and heating by mixing

abstract: We compare a recent study of the kinematics of optical filaments in three cooling flow clusters of galaxies with previous numerical simulations of jet-inflated hot bubbles, and conclude that the velocity structure functions of the filaments better fit direct excitation by the jets than by turbulent cascade from the largest turbulent eddies. The observed velocity structure functions of the optical filaments in the three clusters are steeper than that expected from a classical cascade in turbulent dissipation. Our three-dimensional (3D) hydrodynamical simulations show that as the jets inflate bubbles in the intracluster medium (ICM), they form vortexes in a large range of scales. These vortexes might drive the ICM turbulence with eddies of over more than an order of magnitude in size. A direct excitation of turbulence by the vortexes that the jets form and the slow turbulent dissipation imply that heating the ICM by mixing with hot bubbles is more efficient than heating by turbulent dissipation.


title: Poly-freeness of Artin groups and the Farrell-Jones Conjecture

abstract: We provide two simple proofs of the fact that even Artin groups of FC-type are poly-free which was recently established by R. Blasco-Garcia, C. Mart\'inez-P\'erez and L. Paris. More generally, let $\Gamma$ be a finite simplicial graph with all edges labelled by positive even integers and $A_\Gamma$ be its associated Artin group; our new proof implies that if $A_T$ is poly-free (resp. normally poly-free) for every clique $T$ in $\Gamma$, then $A_\Gamma$ is poly-free (resp. normally poly-free). We prove similar results regarding the Farrell--Jones Conjecture for even Artin groups. In particular, we show that if $A_\Gamma$ is an even Artin group such that each clique in $\Gamma$ either has at most 3 vertices, has all of its labels at least $6$, or is the join of these two types of cliques (the edges connecting the cliques are all labelled by $2$), then $A_\Gamma$ satisfies the Farrell--Jones Conjecture. In addition, our methods enables us to obtain results for general Artin groups.


title: Efficient congruencing in ellipsephic sets: the general case

abstract: In this paper, we bound the number of solutions to a general Vinogradov system of equations $x_1^j+\dots+x_s^j=y_1^j+\dots+y_s^j$, $(1\leq j\leq k)$, as well as other related systems, in which the variables are required to satisfy digital restrictions in a given base. Specifically, our sets of permitted digits have the property that there are few representations of a natural number as sums of elements of the digit set -- the set of squares serving as a key example. We obtain better bounds using this additive structure than could be deduced purely from the size of the set of variables. In particular, when the digits are required to be squares, we obtain diagonal behaviour with $2k(k+1)$ variables.


title: Using asynchronous simulation approach for interactive simulation

abstract: This paper discusses about the advantage of using asynchronous simulation in the case of interactive simulation in which user can steer and control parameters during a simulation in progress. synchronous models allow to compute each iteration faster to address the issues of performance needed in an highly interactive context, and our hypothesis is that get partial results faster is better than getting synchronized and final results to take a decision, in a interactive simulation context.


title: A Branch-and-Price Algorithm for a Team Orienteering Problem with Fixed-Wing Drones

abstract: This paper formulates a team orienteering problem with multiple fixed-wing drones and develops a branch-and-price algorithm to solve the problem to optimality. Fixed-wing drones, unlike rotary drones, have kinematic constraints associated with them, thereby preventing them to make on-the-spot turns and restricting them to a minimum turn radius. This paper presents the implications of these constraints on the drone routing problem formulation and proposes a systematic technique to address them in the context of the team orienteering problem. Furthermore, a novel branch-and-price algorithm with branching techniques specific to the constraints imposed due to fixed-wing drones are proposed. Extensive computational experiments on benchmark instances corroborating the effectiveness of the algorithms are also presented.


title: A Systematic Search for Reddest Far-infrared and Sub-millimeter Galaxies: revealing dust-embedded starbursts at high redshifts

abstract: We present the results of our systematic search for the reddest far-infrared (FIR) and submillimeter (sub-mm) galaxies using the data from the Herschel Multi-tiered Extragalactic Survey (HerMES) and the SCUBA2 Cosmological Legacy Survey (S2CLS). The red FIR galaxies are "500~$\mu$m risers," whose spectral energy distributions increase with wavelength across the three FIR passbands of the Spectral and Photometric Imaging REceiver (SPIRE) of Herschel. Within 106.5 deg$^2$ of the HerMES fields, we have selected 629 500 $\mu$m risers. The red sub-mm galaxies are "SPIRE dropouts," which are prominent detections in the S2CLS 850 $\mu$m data but are extremely weak or invisible in the SPIRE bands. Within the 2.98 deg$^2$ common area of HerMES and S2CLS, we have selected 95 such objects. These very red sources could be dusty starbursts at high redshifts ($z\gtrsim 4$-6) because the peak of their cold-dust emission heated by star formation is shifted to the reddest FIR/sub-mm bands. The surface density of 500 $\mu$m risers is $\sim$8.2 deg$^{-2}$ at the $\geq 20$ mJy level in 500 $\mu$m, while that of SPIRE dropouts is $\sim$19.3 deg$^{-2}$ at the $\geq 5$ mJy level in 850 $\mu$m. Both type of objects could span a wide range of redshifts, however. Using deep radio data in these fields to further select the ones likely at the highest redshifts, we find that the surface density of $z>6$ candidates is 5.5 deg$^{-2}$ among 500 $\mu$m risers and is 0.8-13.6 deg$^{-2}$ among SPIRE dropouts. If this is correct, the dust-embedded star formation processes in such objects could contribute comparably as Lyman-break galaxies to the global SFR density at $z>6$.


title: Debris Disks in Multi-Planet Systems: Are Our Inferences Compromised by Unseen Planets?

abstract: Resolved debris disk features (e.g., warps, offsets, edges and gaps, azimuthal asymmetries, radially thickened rings, scale heights) contain valuable information about the underlying planetary systems, such as the posited planet's mass, semi-major axis, and other orbital parameters. Most existing models assume a single planet is sculpting the disk feature, but recent observations of mature planetary systems (e.g., by radial velocity surveys or \textit{Kepler}) have revealed that many planets reside in multi-planet systems. Here we investigate if/how planet properties inferred from single-planet models are compromised when multiple planets reside in the system. For each disk feature, we build a two-planet model that includes a planet b with fixed parameters and a planet c with a full range of possible parameters. We investigate these two-planet systems and summarize the configurations for which assuming a single planet (i.e., planet b) leads to significantly flawed inferences of that planet's properties. We find that although disk features are usually primarily dominated by a single planet, when using single-planet models we are at risk of misinterpreting planet properties by orders of magnitude in extreme cases. Specifically, we are at high risk of misinterpreting planet properties from disk warps; at moderate risk from disk edges and gaps, radially thickened rings, and scale height features; and at low risk from host star-disk center offsets and azimuthal asymmetries. We summarize situations where we can infer the need to use a multi-planet model instead of a single-planet one from disk morphology dissimilarities.


title: Interactive 3D fluid simulation: steering the simulation in progress using Lattice Boltzmann Method

abstract: This paper describes a work in progress about software and hardware architecture to steer and control an ongoing fluid simulation in a context of a serious game application. We propose to use the Lattice Boltzmann Method as the simulation approach considering that it can provide fully parallel algorithms to reach interactive time and because it is easier to change parameters while the simulation is in progress remaining physically relevant than more classical simulation approaches. We describe which parameters we can modify and how we solve technical issues of interactive steering and we finally show an application of our interactive fluid simulation approach of water dam phenomena.


title: DeepMUSIC: Multiple Signal Classification via Deep Learning

abstract: This letter introduces a deep learning (DL) framework for direction-of-arrival (DOA) estimation. Previous works in DL context mostly consider a single or two target scenario which is a strong limitation in practice. Hence, in this work, we propose a DL framework for multiple signal classification (DeepMUSIC). We design multiple deep convolutional neural networks (CNNs), each of which is dedicated to a subregion of the angular spectrum. In particular, each CNN is fed with the array covariance matrix and it learns the MUSIC spectra of the corresponding angular subregion. We have shown, through simulations, that the proposed DeepMUSIC framework has superior estimation accuracy and exhibits less computational complexity in comparison with both DL and non-DL based techniques.


title: Electron Cyclotron Resonance (ECR) Magnetometry with a Plasma Reservoir

abstract: The local magnetic field in a Penning-Malmberg trap is found by measuring the temperatures that result when electron plasmas are illuminated by microwaves pulses. Multiple heating resonances are observed as the pulse frequencies are swept. The many resonances are due to electron bounce and plasma rotation sidebands. The heating peak corresponding to the cyclotron frequency resonance is identified to determine the magnetic field. A new method for quickly preparing low density electron plasmas for destructive temperature measurements enables a rapid and automated scan of microwave frequencies. This technique can determine the magnetic field to high precision, obtaining an absolute accuracy better than $1\,\mathrm{ppm}$, and a relative precision of $26\,\mathrm{ppb}$. One important application is in situ magnetometry for antihydrogen-based tests of charge-parity-time symmetry and of the weak equivalence principle


title: Biochemically altered human erythrocytes as a carrier for targeted delivery of primaquine: an in vitro study

abstract: The aim of this study was to investigate human erythrocytes as a carrier for targeted drug delivery of primaquine (PQ). The process of PQ loading in human erythrocytes, as well as the effect of PQ loading on the oxidative status of erythrocytes, was also studied. At PQ concentrations of 2, 4, 6, and 8 mg/mL and an incubation time of 2 h, the ratios of the concentrations of PQ entrapped in erythrocytes to that in the incubation medium were 0.515, 0.688, 0.697 and 0.788, respectively. The maximal decline of erythrocyte reduced glutathione content was observed at 8 mg/mL of PQ compared with native erythrocytes p < 0.001. In contrast, malondialdehyde and protein carbonyl were significantly increased in cells loaded with PQ (p < 0.001). Furthermore, osmotic fragility of PQ carrier erythrocytes was increased in comparison with unloaded cells. Electron microscopy revealed spherocyte formation with PQ carrier erythrocytes. PQ-loaded cells showed sustained drug release over a 48 h period. Erythrocytes were loaded with PQ successfully, but there were some biochemical as well as physiological changes that resulted from the effect of PQ on the oxidative status of drug-loaded erythrocytes. These changes may result in favorable targeting of PQ-loaded cells to reticulo-endothelial organs. The relative impact of these changes remains to be explored in ongoing animal studies.


title: On Character Varieties with Non-Connected Structure Groups

abstract: We define character varieties with non-connected structure groups of finitely presented discrete groups and study some basic aspects, such as generic conjugacy classes and relation with fixed points in character varieties with connected structure groups. In the process, we give the classification of the isomorphism classes of the semi-direct products $G\rtimes\Gamma$, with fixed $G$ and $\Gamma$.


title: Origin of Earth's water: sources and constraints

abstract: We review the state of knowledge on the origin of Earth's water. Empirical constraints come from chemical and isotopic measurements of solar system bodies and of Earth itself. Dynamical models have revealed pathways for water delivery to Earth during its formation; most are anchored to specific models for terrestrial planet formation. Meanwhile, disk chemical models focus on determining how the isotopic ratios of the building blocks of planets varied as a function of radial distance and time, defining markers of material transported along those pathways. Carbonaceous chondrite meteorites -- representative of the outer asteroid belt -- match Earth's water isotopes (although mantle plumes have been measured at lower D/H). But how was this connection established -- did Earth's water originate among the asteroids (as in the classical model of terrestrial planet formation)? Or, more likely, was Earth's water delivered from the same parent population as the hydrated asteroids (e.g., external pollution, as in the Grand Tack model)? We argue that the outer asteroid belt -- the boundary between the inner and outer solar system -- is the next frontier for new discoveries. The outer asteroid belt contains a population of icy bodies with volatile-driven activity seen on twelve main belt comets (MBCs); seven of which exhibit sublimation-driven activity on repeated perihelion passages. Measurements of the isotopic characteristics of MBCs would provide essential missing links in the chain between disk models and dynamical models. Finally, we extrapolate to rocky exoplanets. Migration is the only mechanism likely to produce very water-rich planets with more than a few percent water by mass (and even with migration, some planets are pure rock). While water loss mechanisms remain to be studied in more detail, we expect that water should be delivered to the vast majority of rocky exoplanets.


title: Sampling in Thermoacoustic Tomography

abstract: We explore the effect of sampling rates when measuring data given by $Mf$ for special operators $M$ arising in Thermoacoustic Tomography. We start with sampling requirements on $Mf$ given $f$ satisfying certain conditions. After this we discuss the resolution limit on $f$ posed by the sampling rate of $Mf$ without assuming any conditions on these sampling rates. Next we discuss aliasing artifacts when $Mf$ is known to be under sampled in one or more of its variables. Finally, we discuss averaging of measurement data and resulting aliasing and artifacts, along with a scheme for anti-aliasing.


title: Car Pose in Context: Accurate Pose Estimation with Ground Plane Constraints

abstract: Scene context is a powerful constraint on the geometry of objects within the scene in cases, such as surveillance, where the camera geometry is unknown and image quality may be poor. In this paper, we describe a method for estimating the pose of cars in a scene jointly with the ground plane that supports them. We formulate this as a joint optimization that accounts for varying car shape using a statistical atlas, and which simultaneously computes geometry and internal camera parameters. We demonstrate that this method produces significant improvements for car pose estimation, and we show that the resulting 3D geometry, when computed over a video sequence, makes it possible to improve on state of the art classification of car behavior. We also show that introducing the planar constraint allows us to estimate camera focal length in a reliable manner.


title: On the well-posedness of Galbrun's equation

abstract: Galbrun's equation, which is a second order partial differential equation describing the evolution of a so-called Lagrangian displacement vector field, can be used to study acoustics in background flows as well as perturbations of astrophysical flows. Our starting point for deriving Galbrun's equation is linearized Euler's equations, which is a first order system of partial differential equations that describe the evolution of the so-called Eulerian flow perturbations. Given a solution to linearized Euler's equations, we introduce the Lagrangian displacement as the solution to a linear first order partial differential equation, driven by the Eulerian perturbation of the fluid velocity. Our Lagrangian displacement solves Galbrun's equation, provided it is regular enough and that the so-called "no resonance" assumption holds. In the case that the background flow is steady and tangential to the domain boundary, we prove existence, uniqueness, and continuous dependence on data of solutions to an initial--boundary value problem for linearized Euler's equations. For such background flows, we demonstrate that the Lagrangian displacement is well-defined, that the initial datum of the Lagrangian displacement can be chosen in order to fulfill the "no resonance" assumption, and derive a classical energy estimate for (sufficiently regular solutions to) Galbrun's equation. Due to the presence of zeroth order terms of indefinite signs in the equations, the energy estimate allows solutions that grow exponentially with time.


title: Trust-Region Newton-CG with Strong Second-Order Complexity Guarantees for Nonconvex Optimization

abstract: Worst-case complexity guarantees for nonconvex optimization algorithms have been a topic of growing interest. Multiple frameworks that achieve the best known complexity bounds among a broad class of first- and second-order strategies have been proposed. These methods have often been designed primarily with complexity guarantees in mind and, as a result, represent a departure from the algorithms that have proved to be the most effective in practice. In this paper, we consider trust-region Newton methods, one of the most popular classes of algorithms for solving nonconvex optimization problems. By introducing slight modifications to the original scheme, we obtain two methods -- one based on exact subproblem solves and one exploiting inexact subproblem solves as in the popular "trust-region Newton-Conjugate-Gradient" (trust-region Newton-CG) method -- with iteration and operation complexity bounds that match the best known bounds for the aforementioned class of first- and second-order methods. The resulting trust-region Newton-CG method also retains the attractive practical behavior of classical trust-region Newton-CG, which we demonstrate with numerical comparisons on a standard benchmark test set.


title: Interleaving by Parts: Join Decompositions of Interleavings and Join-Assemblage of Geodesics

abstract: Metrics of interest in topological data analysis (TDA) are often explicitly or implicitly in the form of an interleaving distance $d_{\mathrm{I}}$ between poset maps (i.e. order-preserving maps), e.g. the Gromov-Hausdorff distance between metric spaces can be reformulated in this way. We propose a representation of a poset map $\mathbf{F}:\mathcal{P}\to\mathcal{Q}$ as a join (i.e. supremum) $\bigvee_{b\in B} \mathbf{F}_b$ of simpler poset maps $\mathbf{F}_b$ (for a join dense subset $B\subset \mathcal{Q}$) which in turn yields a decomposition of $d_{\mathrm{I}}$ into a product metric. The decomposition of $d_{\mathrm{I}}$ is simple, but its ramifications are manifold: (1) We can construct a geodesic path between any poset maps $\mathbf{F}$ and $\mathbf{G}$ with $d_{\mathrm{I}}(\mathbf{F},\mathbf{G})<\infty$ by assembling geodesics between all $\mathbf{F}_b$s and $\mathbf{G}_b$s via the join operation. This construction generalizes at least three constructions of geodesic paths that have appeared in the literature. (2) We can extend the Gromov-Hausdorff distance to a distance between simplicial filtrations over an arbitrary poset with a flow, preserving its universality and geodesicity. (3) We can clarify equivalence between several known metrics on multiparameter hierarchical clusterings. (4) We can illuminate the relationship between the erosion distance by Patel and the graded rank function by Betthauser, Bubenik, and Edwards, which in turn takes us to an interpretation on the representation $\bigvee_b \mathbf{F}_b$ as a generalization of persistence landscapes and graded rank functions.


title: Derived equivalences between skew-gentle algebras using orbifolds

abstract: Skew-gentle algebras are skew-group algebras of gentle algebras equipped with a certain $\Z_2$-action. Building on the bijective correspondence between gentle algebras and dissected surfaces, we obtain in this paper a bijection between skew-gentle algebras and certain dissected orbifolds that admit a double cover. We prove the compatibility of the $\Z_2$-action on the double cover with the skew-group algebra construction. This allows us to investigate the derived equivalence relation between skew-gentle algebras in geometric terms: We associate to each skew-gentle algebra a line field on the orbifold, and on its double cover, and interpret different kinds of derived equivalences of skew-gentle algebras in terms of diffeomorphisms respecting the homotopy class of the line fields associated to the algebras.


title: Learning Non-Markovian Quantum Noise from Moir\'{e}-Enhanced Swap Spectroscopy with Deep Evolutionary Algorithm

abstract: Two-level-system (TLS) defects in amorphous dielectrics are a major source of noise and decoherence in solid-state qubits. Gate-dependent non-Markovian errors caused by TLS-qubit coupling are detrimental to fault-tolerant quantum computation and have not been rigorously treated in the existing literature. In this work, we derive the non-Markovian dynamics between TLS and qubits during a SWAP-like two-qubit gate and the associated average gate fidelity for frequency-tunable Transmon qubits. This gate dependent error model facilitates using qubits as sensors to simultaneously learn practical imperfections in both the qubit's environment and control waveforms. We combine the-state-of-art machine learning algorithm with Moir\'{e}-enhanced swap spectroscopy to achieve robust learning using noisy experimental data. Deep neural networks are used to represent the functional map from experimental data to TLS parameters and are trained through an evolutionary algorithm. Our method achieves the highest learning efficiency and robustness against experimental imperfections to-date, representing an important step towards in-situ quantum control optimization over environmental and control defects.


title: Classifying sections of del Pezzo fibrations, I

abstract: We develop a strategy to classify the components of the space of sections of a del Pezzo fibration over $\mathbb{P}^{1}$. In particular, we prove the Movable Bend and Break lemma for del Pezzo fibrations. Our approach is motivated by Geometric Manin's Conjecture and proves upper bounds on the associated counting function. We also give applications to enumerativity of Gromov-Witten invariants and to the study of the Abel-Jacobi map.


title: Cross-Language Aphasia Detection using Optimal Transport Domain Adaptation

abstract: Multi-language speech datasets are scarce and often have small sample sizes in the medical domain. Robust transfer of linguistic features across languages could improve rates of early diagnosis and therapy for speakers of low-resource languages when detecting health conditions from speech. We utilize out-of-domain, unpaired, single-speaker, healthy speech data for training multiple Optimal Transport (OT) domain adaptation systems. We learn mappings from other languages to English and detect aphasia from linguistic characteristics of speech, and show that OT domain adaptation improves aphasia detection over unilingual baselines for French (6% increased F1) and Mandarin (5% increased F1). Further, we show that adding aphasic data to the domain adaptation system significantly increases performance for both French and Mandarin, increasing the F1 scores further (10% and 8% increase in F1 scores for French and Mandarin, respectively, over unilingual baselines).


title: Galaxy interactions in IllustrisTNG-100, I: The power and limitations of visual identification

abstract: We present a sample of 446 galaxy pairs constructed using the cosmological simulation IllustrisTNG-100 at z = 0, with M$_{FoF, dm}$ = 10$^{11}$-10$^{13.5}$ M$_{\odot}$. We produce ideal mock SDSS g-band images of all pairs to test the reliability of visual classification schema employed to produce samples of interacting galaxies. We visually classify each image as interacting or not based on the presence of a close neighbour, the presence of stellar debris fields, disturbed discs, and/or tidal features. By inspecting the trajectories of the pairs, we determine that these indicators correctly identify interacting galaxies $\sim45\%$ of the time. We subsequently split the sample into the visually identified interacting pairs (VIP; 38 pairs) and those which are interacting but are not visually identified (nonVIP; 47 pairs). We find that VIP have undergone a close passage nearly twice as recently as the nonVIP, and typically have higher stellar masses. Further, the VIP sit in dark matter haloes that are approximately 2.5 times as massive, in environments nearly 2 times as dense, and are almost a factor of 10 more affected by the tidal forces of their surroundings than the nonVIP. These factors conspire to increase the observability of tidal features and disturbed morphologies, making the VIP more likely to be identified. Thus, merger rate calculations which rely on stellar morphologies are likely to be significantly biased toward massive galaxy pairs which have recently undergone a close passage.


title: On Post-Starburst Galaxies Dominating Tidal Disruption Events

abstract: A starburst induced by a galaxy merger may create a relatively thin central stellar disk at radius $\le 100$pc. We calculate the rate of tidal disruption events (TDEs) by the inspiraling secondary supermassive black (SMBH) through the disk. With a small enough stellar velocity dispersion ($\sigma/v_c \le 0.1$) in the disk, it is shown that $10^5-10^6$ TDEs of solar-type main sequence stars per post-starburst galaxy (PSB) can be produced to explain their dominance in producing observed TDEs. Although the time it takes to bring the secondary SMBH to the disk apparently varies in the range of $\sim 0.1-1$Gyr since the starburst, depending on its landing location and subsequently due to dynamical friction with stars exterior to the central stellar disk in question, the vast majority of TDEs by the secondary SMBH in any individual PSB occurs within a space of time shorter than $\sim 30$Myr. Five unique testable predictions of this model are suggested.


title: Module constructions for certain subgroups of the largest Mathieu group

abstract: For certain subgroups of $M_{24}$, we give vertex operator algebraic module constructions whose associated trace functions are meromorphic Jacobi forms. These meromorphic Jacobi forms are canonically associated to the mock modular forms of Mathieu moonshine. The construction is related to the Conway moonshine module and employs a technique introduced by Anagiannis--Cheng--Harrison. With this construction we are able to give concrete vertex algebraic realizations of certain cuspidal Hecke eigenforms of weight two. In particular, we give explicit realizations of trace functions whose integralities are equivalent to divisibility conditions on the number of $\mathbb{F}_p$ points on the Jacobians of modular curves.


title: Geometry of multigraded rings and embeddings of toric varieties

abstract: We use homogeneous spectra of multigraded rings to construct toric embeddings of a large family of projective varieties which preserve some of the birational geometry of the underlying variety, generalizing the well-known construction associated to Mori Dream Spaces.


title: Sequential generation of linear cluster states from a single photon emitter

abstract: Light states composed of multiple entangled photons - such as cluster states - are essential for developing and scaling-up quantum computing networks. Photonic cluster states with discrete variables can be obtained from single-photon sources and entangling gates, but so far this has only been done with probabilistic sources constrained to intrinsically-low efficiencies, and an increasing hardware overhead. Here, we report the resource-efficient generation of polarization-encoded, individually-addressable, photons in linear cluster states occupying a single spatial mode. We employ a single entangling-gate in a fiber loop configuration to sequentially entangle an ever-growing stream of photons originating from the currently most efficient single-photon source technology - a semiconductor quantum dot. With this apparatus, we demonstrate the generation of linear cluster states up to four photons in a single-mode fiber. The reported architecture can be programmed to generate linear-cluster states of any number of photons with record scaling ratios, potentially enabling practical implementation of photonic quantum computing schemes.


title: Modular Multimodal Architecture for Document Classification

abstract: Page classification is a crucial component to any document analysis system, allowing for complex branching control flows for different components of a given document. Utilizing both the visual and textual content of a page, the proposed method exceeds the current state-of-the-art performance on the RVL-CDIP benchmark at 93.03% test accuracy.


title: LSTM Neural Networks: Input to State Stability and Probabilistic Safety Verification

abstract: The goal of this paper is to analyze Long Short Term Memory (LSTM) neural networks from a dynamical system perspective. The classical recursive equations describing the evolution of LSTM can be recast in state space form, resulting in a time-invariant nonlinear dynamical system. A sufficient condition guaranteeing the Input-to-State (ISS) stability property of this class of systems is provided. The ISS property entails the boundedness of the output reachable set of the LSTM. In light of this result, a novel approach for the safety verification of the network, based on the Scenario Approach, is devised. The proposed method is eventually tested on a pH neutralization process.


title: Depth-Width Trade-offs for ReLU Networks via Sharkovsky's Theorem

abstract: Understanding the representational power of Deep Neural Networks (DNNs) and how their structural properties (e.g., depth, width, type of activation unit) affect the functions they can compute, has been an important yet challenging question in deep learning and approximation theory. In a seminal paper, Telgarsky highlighted the benefits of depth by presenting a family of functions (based on simple triangular waves) for which DNNs achieve zero classification error, whereas shallow networks with fewer than exponentially many nodes incur constant error. Even though Telgarsky's work reveals the limitations of shallow neural networks, it does not inform us on why these functions are difficult to represent and in fact he states it as a tantalizing open question to characterize those functions that cannot be well-approximated by smaller depths. In this work, we point to a new connection between DNNs expressivity and Sharkovsky's Theorem from dynamical systems, that enables us to characterize the depth-width trade-offs of ReLU networks for representing functions based on the presence of generalized notion of fixed points, called periodic points (a fixed point is a point of period 1). Motivated by our observation that the triangle waves used in Telgarsky's work contain points of period 3 - a period that is special in that it implies chaotic behavior based on the celebrated result by Li-Yorke - we proceed to give general lower bounds for the width needed to represent periodic functions as a function of the depth. Technically, the crux of our approach is based on an eigenvalue analysis of the dynamical system associated with such functions.


title: General Matrix-Matrix Multiplication Using SIMD features of the PIII

abstract: Generalised matrix-matrix multiplication forms the kernel of many mathematical algorithms. A faster matrix-matrix multiply immediately benefits these algorithms. In this paper we implement efficient matrix multiplication for large matrices using the floating point Intel Pentium SIMD (Single Instruction Multiple Data) architecture. A description of the issues and our solution is presented, paying attention to all levels of the memory hierarchy. Our results demonstrate an average performance of 2.09 times faster than the leading public domain matrix-matrix multiply routines.


title: Thermodynamic theory of crystal plasticity: formulation and application to polycrystal fcc copper

abstract: We present a thermodynamic description of crystal plasticity. Our formulation is based on the Langer-Bouchbinder-Lookman thermodynamic dislocation theory (TDT), which asserts the fundamental importance of an effective temperature that describes the state of configurational disorder and therefore the dislocation density of the crystalline material. We extend the TDT description from isotropic plasticity to crystal plasticity with many slip systems. Finite-element simulations show favourable comparison with experiments on polycrystal fcc copper under uniaxial compression, tension, and simple shear. The thermodynamic theory of crystal plasticity thus provides a thermodynamically consistent and physically rigorous description of dislocation motion in crystals. We also discuss new insights about the interaction of dislocations belonging to different slip systems.

categories: cond-mat.mtrl-sci cond-mat.mes-hall cond-mat.stat-mech physics.comp-ph


title: A Dataset for measuring reading levels in India at scale

abstract: One out of four children in India are leaving grade eight without basic reading skills. Measuring the reading levels in a vast country like India poses significant hurdles. Recent advances in machine learning opens up the possibility of automating this task. However, the datasets of children's speech are not only rare but are primarily in English. To solve this assessment problem and advance deep learning research in regional Indian languages, we present the ASER dataset of children in the age group of 6-14. The dataset consists of 5,301 subjects generating 81,330 labeled audio clips in Hindi, Marathi and English. These labels represent expert opinions on the child's ability to read at a specified level. Using this dataset, we built a simple ASR-based classifier. Early results indicate that we can achieve a prediction accuracy of 86% for the English language. Considering the ASER survey spans half a million subjects, this dataset can grow to those scales.


title: The hadronic vacuum polarization of the muon from four-flavor lattice QCD

abstract: We present an update on the ongoing calculations by the Fermilab Lattice, HPQCD, and MILC Collaboration of the leading-order (in electromagnetism) hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon. Our project employs ensembles with four flavors of highly improved staggered fermions, physical light-quark masses, and four lattice spacings ranging from $a \approx 0.06$ to 0.15 fm for most of the results thus far.


title: On the Topology of Fano Smoothings

abstract: Suppose that X is a Fano manifold that corresponds under Mirror Symmetry to a Laurent polynomial f, and that P is the Newton polytope of f. In this setting it is expected that there is a family of algebraic varieties over the unit disc with general fiber X and special fiber the toric variety defined by the spanning fan of P. Building on recent work and conjectures by Corti--Hacking--Petracci, who construct such families of varieties, we determine the topology of the general fiber from combinatorial data on P. This provides evidence for the Corti--Hacking--Petracci conjectures, and verifies that their construction is compatible with expectations from Mirror Symmetry.


title: Training Deep Neural Networks to Detect Repeatable 2D Features Using Large Amounts of 3D World Capture Data

abstract: Image space feature detection is the act of selecting points or parts of an image that are easy to distinguish from the surrounding image region. By combining a repeatable point detection with a descriptor, parts of an image can be matched with one another, which is useful in applications like estimating pose from camera input or rectifying images. Recently, precise indoor tracking has started to become important for Augmented and Virtual reality as it is necessary to allow positioning of a headset in 3D space without the need for external tracking devices. Several modern feature detectors use homographies to simulate different viewpoints, not only to train feature detection and description, but test them as well. The problem is that, often, views of indoor spaces contain high depth disparity. This makes the approximation that a homography applied to an image represents a viewpoint change inaccurate. We claim that in order to train detectors to work well in indoor environments, they must be robust to this type of geometry, and repeatable under true viewpoint change instead of homographies. Here we focus on the problem of detecting repeatable feature locations under true viewpoint change. To this end, we generate labeled 2D images from a photo-realistic 3D dataset. These images are used for training a neural network based feature detector. We further present an algorithm for automatically generating labels of repeatable 2D features, and present a fast, easy to use test algorithm for evaluating a detector in an 3D environment.


title: Multi-Stage Preconditioners for Thermal-Compositional-Reactive Flow in Porous Media

abstract: We present a family of multi-stage preconditioners for coupled thermal-compositional-reactive reservoir simulation problems. The most common preconditioner used in industrial practice, the Constrained Pressure Residual (CPR) method, was designed for isothermal models and does not offer a specific strategy for the energy equation. For thermal simulations, inadequate treatment of the temperature unknown can cause severe convergence degradation. When strong thermal diffusion is present, the energy equation exhibits significant elliptic behavior that cannot be accurately corrected by CPR's second stage. In this work, we use Schur-complement decompositions to extract a temperature subsystem and apply an Algebraic MultiGrid (AMG) approximation as an additional preconditioning stage to improve the treatment of the energy equation. We present results for several two-dimensional hot air injection problems using an extra heavy oil, including challenging reactive In-Situ Combustion (ISC) cases. We show improved performance and robustness across different thermal regimes, from advection dominated (high Peclet number) to diffusion dominated (low Peclet number). The number of linear iterations is reduced by 40-85% compared to standard CPR for both homogeneous and heterogeneous media, and the new methods exhibit almost no sensitivity to the thermal regime.


title: Carbon and Oxygen Isotopic Ratios. II. Semiregular Variable M Giants

abstract: Carbon and oxygen isotopic ratios are reported for a sample of 51 SRb- and Lb-type variable asymptotic giant branch stars. Vibration-rotation first- and second-overtone CO lines in 1.5-2.5 mum spectra were measured to derive isotopic ratios for 12C/13C, 16O/17O, and 16O/18O. Comparisons with previous measurements for individual stars and with various samples of evolved stars, as available in the extant literature, are discussed. Using the oxygen isotopic ratios, the masses of the SRb stars can be derived. Combining the masses with Gaia luminosities, the SRb stars are shown to be antecedents of the Mira variables. The limiting parameters where plane-parallel, hydrostatic equilibrium model atmospheres can be used for abundance analysis of M giants are explored.


title: Determination of the strong coupling constant $\alpha_{S}(m_\mathrm{Z})$ from measurements of inclusive W$^\pm$ and Z boson production cross sections in proton-proton collisions at $\sqrt{s} =$ 7 and 8 TeV

abstract: Twelve measurements of inclusive cross sections of W$^\pm$ and Z boson production, performed in proton-proton collisions at centre-of-mass energies of 7 and 8 TeV, are compared with perturbative quantum chromodynamics calculations at next-to-next-to-leading order (NNLO) accuracy obtained with the CT14, HERAPDF2.0, MMHT14, and NNPDF3.0 parton distribution functions (PDFs). Data and theory agree well for all PDF sets, taking into account the experimental and theoretical uncertainties. A novel procedure is employed to extract the strong coupling constant at the Z pole mass from a detailed comparison of all the experimental fiducial cross sections to the corresponding NNLO theoretical predictions, yielding $\alpha_{S}(m_\mathrm{Z}) =$ 0.1163 $^{+0.0024}_{-0.0031}$ (CT14), 0.1072 $^{+0.0043}_{-0.0040}$ (HERAPDF2.0), 0.1186 $\pm$ 0.0025 (MMHT14), and 0.1147 $\pm$ 0.0023 (NNPDF3.0). Using the results obtained with the CT14 and MMHT14 PDFs, which yield the most robust and stable $\alpha_{S}(m_\mathrm{Z})$ extractions, a value $\alpha_{S}(m_\mathrm{Z}) =$ 0.1175 $^{+0.0025}_{-0.0028}$ is determined.


title: Convergence of the method of reflections for particle suspensions in Stokes flows

abstract: We study the convergence of the method of reflections for the Stokes equations in domains perforated by countably many spherical particles with boundary conditions typical for the suspension of rigid particles. We prove that a relaxed version of the method is always convergent in $\dot H^1$ under a mild separation condition on the particles. Moreover, we prove optimal convergence rates of the method in $\dot W^{1,q}$, $1 < q < \infty$ and in $L^\infty$ in terms of the particle volume fraction under a stronger separation condition of the particles.


title: The Influence of Stellar Contamination on the Interpretation of Near Infrared Transmission Spectra of Sub-Neptune Worlds around M-Dwarfs

abstract: The impact of unocculted stellar surface heterogeneities in the form of cool spots and hot faculae on the spectrum of a transiting planet has been a daunting problem for the characterization of exoplanet atmospheres. The wavelength-dependent nature of stellar surface heterogeneities imprinting their signatures on planetary transmission spectra are of concern particularly for systems of sub-Neptunes orbiting M-dwarfs. Here we present a systematic exploration of the impact of this spot-contamination on simulated near infrared transmission spectra of sub-Neptune planets. From our analysis, we find that improper correction of stellar surface heterogeneities on transmission spectra can lead to significant bias when inferring planetary atmospheric properties. However, this bias is negligible for lower fractions of heterogeneities (<1%). Additionally, we find that acquiring a priori knowledge of stellar heterogeneities does not improve precision in constraining planetary parameters if the heterogeneities are appropriately marginalized within a retrieval---however these are conditional on our confidence of stellar atmospheric models being accurate representations of the true photosphere. In sum, to acquire unbiased constraints when characterizing planetary atmospheres with the James Webb Space Telescope, we recommend performing joint retrievals of both the disk-integrated spectrum of the star and the stellar contamination corrected transmission spectrum.


title: A refined machinery to calculate large moments from coupled systems of linear differential equations

abstract: The large moment method can be used to compute a large number of moments of physical quantities that are described by coupled systems of linear differential equations. Besides these systems the algorithm requires a certain number of initial values as input, that are often hard to derive in a preprocessing step.Thus a major challenge is to keep the number of initial values as small as possible. We present the basic ideas of the underlying large moment method and present refined versions that reduce significantly the number of required initial values.


title: Semi-supervised Learning Approach to Generate Neuroimaging Modalities with Adversarial Training

abstract: Magnetic Resonance Imaging (MRI) of the brain can come in the form of different modalities such as T1-weighted and Fluid Attenuated Inversion Recovery (FLAIR) which has been used to investigate a wide range of neurological disorders. Current state-of-the-art models for brain tissue segmentation and disease classification require multiple modalities for training and inference. However, the acquisition of all of these modalities are expensive, time-consuming, inconvenient and the required modalities are often not available. As a result, these datasets contain large amounts of \emph{unpaired} data, where examples in the dataset do not contain all modalities. On the other hand, there is smaller fraction of examples that contain all modalities (\emph{paired} data) and furthermore each modality is high dimensional when compared to number of datapoints. In this work, we develop a method to address these issues with semi-supervised learning in translating between two neuroimaging modalities. Our proposed model, Semi-Supervised Adversarial CycleGAN (SSA-CGAN), uses an adversarial loss to learn from \emph{unpaired} data points, cycle loss to enforce consistent reconstructions of the mappings and another adversarial loss to take advantage of \emph{paired} data points. Our experiments demonstrate that our proposed framework produces an improvement in reconstruction error and reduced variance for the pairwise translation of multiple modalities and is more robust to thermal noise when compared to existing methods.


title: Multistage Graph Problems on a Global Budget

abstract: Time-evolving or temporal graphs gain more and more popularity when studying the behavior of complex networks. In this context, the multistage view on computational problems is among the most natural frameworks. Roughly speaking, herein one studies the different (time) layers of a temporal graph (effectively meaning that the edge set may change over time, but the vertex set remains unchanged), and one searches for a solution of a given graph problem for each layer. The twist in the multistage setting is that the solutions found must not differ too much between subsequent layers. We relax on this already established notion by introducing a global instead of the local budget view studied so far. More specifically, we allow for few disruptive changes between subsequent layers but request that overall, that is, summing over all layers, the degree of change is moderate. Studying several classical graph problems (both NP-hard and polynomial-time solvable ones) from a parameterized complexity angle, we encounter both fixed-parameter tractability and parameterized hardness results. Somewhat surprisingly, we find that sometimes the global multistage versions of NP-hard problems such as Vertex Cover turn out to be computationally more tractable than the ones of polynomial-time solvable problems such as Matching.


title: Direct Numerical Simulation of Rotating Ellipsoidal Particles using Moving Nonconforming Schwarz-Spectral Element Method

abstract: We present application of a highly-scalable overlapping grid-based nonconforming Schwarz-spectral element method (Schwarz-SEM) to study the dynamics of rotating ellipsoidal particles. The current study is one of the first to explore the effect of rotation on ellipsoidal particles using fully resolved simulations (direct numerical simulation). The rotating ellipsoidal particles show substantial difference in the dynamics of the flow, when compared against non-rotating particles. The difference is primarily due to periodic attachment and separation of the flow to the surface of the particle for the rotating cases, which results in a higher drag on the particles when compared to the corresponding non-rotating cases. The dynamics is also different from a rotating spherical particle, where a steady shear layer develops near the surface of the sphere. For the rotating ellipsoidal particles, this mechanism results in a phase-difference between the position of observed maximum and minimum drag, and the position of expected maximum and minimum drag (i.e., maximum and minimum projected area). A similar phase-difference is also observed for the lift acting on the rotating ellipsoidal particles. The results presented here demonstrate the importance of explicitly modeling the shape and rotation of particles when we study the dynamics of non-spherical particles. Finally, the study also validates the use of non-conforming Schwarz-SEM for tackling problems in fully resolved particulate flow dynamics.


title: Multiregeneration for polynomial system solving

abstract: We demonstrate our implementation of a continuation method as described in \cite{HR2015} for solving polynomials systems. Given a sequence of (multi)homogeneous polynomials, the software "multiregeneration" outputs the respective (multi)degree in a wide range of cases and partial multidegree in all others. We use Python for the file processing, while Bertini is needed for the continuation. Moreover, parallelization options and several strategies for solving structured polynomial systems are available.


title: Formation of ultracold weakly bound dimers of bosonic $^{23}\textrm{Na}^{39}\textrm{K}$

abstract: We create weakly bound bosonic $^{23}\textrm{Na}^{39}\textrm{K}$ molecules in a mixture of ultracold $^{23}\textrm{Na}$ and $^{39}\textrm{K}$. The creation is done in the vicinity of a so far undetected Feshbach resonance at about $196\,\text{G}$ which we identify in this work by atom-loss spectroscopy. We investigate the involved molecular state by performing destructive radio frequency binding energy measurements. For the constructive molecule creation we use radio frequency pulses with which we assemble up to $6000$ molecules. We analyze the molecule creation efficiency as a function of the radio frequency pulse duration and the atom number ratio between $^{23}\textrm{Na}$ and $^{39}\textrm{K}$. We find an overall optimal efficiency of $6\,\%$ referring to the $^{39}\textrm{K}$ atom number. The measured lifetime of the molecules in the bath of trapped atoms is about $0.3\,\textrm{ms}$.


title: A New Task for the Philosophy of Science

abstract: This paper argues that philosophers of science have before them an important new task that they urgently need to take up. It is to convince the scientific community to adopt and implement a new philosophy of science that does better justice to the deeply problematic basic intellectual aims of science than that which we have at present. Problematic aims evolve with evolving knowledge, that part of philosophy of science concerned with aims and methods thus becoming an integral part of science itself. The outcome of putting this new philosophy into scientific practice would be a new kind of science, both more intellectually rigorous and one that does better justice to the best interests of humanity.


title: Excess noise in Al${}_\text{x}$Ga${}_\text{1-x}$As/GaAs based quantum rings

abstract: Cross-correlated noise measurements are performed in etched Al${}_\text{x}$Ga${}_\text{1-x}$As/GaAs based quantum rings in equilibrium at bath temperature of $T_\text{bath}=4.2\text{ K}$. The measured white noise exceeds the thermal (Johnson-Nyquist) noise expected from the measured electron temperature $T_\text{e}$ and the electrical resistance $R$. This excess part of the white noise decreases as $T_\text{bath}$ increases and vanishes for $T_\text{bath}\geq 12\text{ K}$. Excess noise is neither observed if one arm of a quantum ring is depleted of electrons nor in 1D-constrictions that have a length and width comparable to the quantum rings. A model is presented that suggests that the excess noise originates from the correlation of noise sources, mediated by phase-coherent propagation of electrons.


title: System-bath entanglement theorem with Gaussian environments

abstract: In this work, we establish a so-called "system-bath entanglement theorem", for arbitrary systems coupled with Gaussian environments. This theorem connects the entangled system-bath response functions in the total composite space to those of local systems, as long as the interacting bath spectral densities are given. We validate the theorem with the direct evaluation via the exact dissipaton-equation-of-motion approach. Therefore, this work enables various quantum dissipation theories, which originally describe only the reduced system dynamics, for their evaluations on the system-bath entanglement properties. Numerical demonstrations are carried out on the Fano interference spectroscopies of spin-boson systems.


title: Entanglement, violation of Kramers-Kronig relation and curvature in spacetime

abstract: Recent studies show that a maximally entangled Schwinger pair creates a nontraversable Einstein-Rosen~(ER) bridge (wormhole) in the gravitational theory, which bridges causally not connected parts of the AdS spacetime~[PRL 111, 211603 and 211602]. Same authors raise the possibility of a general correspondence between a weaker entanglement and curvature. Here, we provide some clues for such a generalized relation. First, we show that (i) entanglement and (ii) violation of Kramers-Kronig~(KK) relations appear at the same critical parameters for a standard two-mode (squeezing) entanglement interaction. We also include the dampings. Second, we bring a study into attention: Presence of a spacetime-curvature in vacuum polarization, electron-positron Schwinger pairs, makes QED also violate the KK relations without violating the causality. Then, we discuss that these findings are possible to provide new clues for a generalized relation between entanglement and spacetime-curvature. Such an interpretation may also save violation of KK relations from implying the violation of the causality.


title: Evolution of 3D Printing Methods and Materials for Electrochemical Energy Storage

abstract: Additive manufacturing has revolutionized the building of materials direct from design, allowing high resolution rapid prototyping in complex 3D designs with many materials. 3D printing hasenabled high strength damage-tolerant structures, bioprinted artificial organs and tissues, ultralight metals, medicine, education, prosthetics, architecture, consumer electronics,and as a prototyping tool for engineers and hobbyists alike. 3D printing has emerged as a useful tool for complex electrode and material assembly method for batteries and supercapacitors in recent years. The field initially grew from extrusion-based methods such as fused deposition modelling, and evolved to photopolymerization printing of intricate composites, while supercapacitor technologies less sensitive to solvents more often involved material jetting processes. Underpinning every part of a 3D printable battery and many other devices is the printing method and the nature of the feed material. Material purity, printing fidelity, accuracy, complexity, and the ability to form conductive, ceramic, glassy, or solvent-stable plastics relies on the nature of the feed material or composite to such an extent, that the future of 3D printable batteries and electrochemical energy storage devices will depend on materials and printing methods that are co-operatively informed by the requirements of the device and how it is fabricated. In this Perspective, we address the materials and methods requirements in 3D printable batteries and supercapacitors and outline requirements for the future of the field by linking existing performance limitations to the requirements of printable energy storage materials, casing materials and the direct printing of electrodes and electrolytes. We also look to the future by taking inspiration from additive manufacturing, to posit links between materials and printing methods to allow new form factor cells.


title: An Introduction to Mordell Weil Theorem

abstract: This article is an introduction to Mordell-Weil theorem. In this article, I introduced some basic properties about ellptic curves and proved the theorem in two different ways.


title: Optimal management of impaired self-avoiding random walks for minimizing spatial coverage

abstract: Self-avoidance is a common mechanism to improve the efficiency of a random walker for covering a spatial domain. However, how this efficiency decreases when self-avoidance is impaired or limited by other processes has remained largely unexplored. Here we use simulations to study the case when the self-avoiding signal left by a walker both (i) saturates after successive revisits to a site, and (ii) evaporates, or dissappears, after some characteristic time. We surprisingly reveal that the mean cover time becomes minimum for intermediate values of the evaporation time, leading to the existence of a nontrivial optimum management of the self-avoiding signal. We argue that this is a consequence of complex blocking effects caused by the interplay with the signal saturation and, remarkably, we show that the optimum becomes more and more significant as the domain size increases.


title: Analysis of dysautonomia in patients with Chagas Cardiomyopathy

abstract: Chagas disease American trypanosomiasis is caused by a flagellated parasite: trypanosoma cruzi, transmitted by an insect of the genus Triatoma and also by blood transfusions. In Latin America the number of infected people is approximately 6 million, with a population exposed to the risk of infection of 550000. It is our interest to develop a non-invasive, low-cost methodology, capable of detecting any alteration early on cardiaca produced by T. cruzi. We analyzed the 24 hour RR records in patients with ECG abnormalities (CH2), patients without ECG alterations (CH1) who had positive serological findings for Chagas disease and healthy (Control) matched by sex and age. We found significant differences between the Control, CH1 and CH2 groups that show dysautonomy and enervation of the autonomic nervous system.


title: Characterization of the new free-air primary standard for medium-energy X-rays at CMI

abstract: In 2011 a decision was made by Czech Metrology Institute to build a free-air ionization chamber (FAC) intended to be used as a primary standard of air kerma rate for medium-energy X-rays (photon energy from 40 to 300 keV, including mammography X-ray qualities) in order to replace currently used secondary ionization chamber and to decrease the uncertainty of air kerma reference value. In period of 2013-2017, the FAC has been designed, manufactured and put into operation. Correction factors were measured or calculated by a Monte Carlo method. FAC performance was preliminary tested using a calibrated secondary chamber. Physical characteristics of the FAC are described and a summary of the correction factors with the uncertainty budget is presented.


title: Explosive higher-order Kuramoto dynamics on simplicial complexes

abstract: The higher-order interactions of complex systems, such as the brain are captured by their simplicial complex structure and have a significant effect on dynamics. However, the existing dynamical models defined on simplicial complexes make the strong assumption that the dynamics resides exclusively on the nodes. Here we formulate the higher-order Kuramoto model which describes the interactions between oscillators placed not only on nodes but also on links, triangles, and so on. We show that higher-order Kuramoto dynamics can lead to an explosive synchronization transition by using an adaptive coupling dependent on the solenoidal and the irrotational component of the dynamics.

categories: nlin.AO cond-mat.dis-nn cond-mat.stat-mech cs.SI physics.bio-ph physics.soc-ph


title: Semiparametric Regression for Dual Population Mortality

abstract: Parameter shrinkage applied optimally can always reduce error and projection variances from those of maximum likelihood estimation. Many variables that actuaries use are on numerical scales, like age or year, which require parameters at each point. Rather than shrinking these towards zero, nearby parameters are better shrunk towards each other. Semiparametric regression is a statistical discipline for building curves across parameter classes using shrinkage methodology. It is similar to but more parsimonious than cubic splines. We introduce it in the context of Bayesian shrinkage and apply it to joint mortality modeling for related populations. Bayesian shrinkage of slope changes of linear splines is an approach to semiparametric modeling that evolved in the actuarial literature. It has some theoretical and practical advantages, like closed-form curves, direct and transparent determination of degree of shrinkage and of placing knots for the splines, and quantifying goodness of fit. It is also relatively easy to apply to the many nonlinear models that arise in actuarial work. We find that it compares well to a more complex state-of-the-art statistical spline shrinkage approach on a popular example from that literature.


title: $\widehat{\mathfrak{sl}}(n)_N$ WZW conformal blocks from $SU(N)$ instanton partition functions on ${\mathbb {C}}^2/{\mathbb {Z}}_n$

abstract: Generalizations of the AGT correspondence between 4D $\mathcal{N}=2$ $SU(2)$ supersymmetric gauge theory on ${\mathbb {C}}^2$ with $\Omega$-deformation and 2D Liouville conformal field theory include a correspondence between 4D $\mathcal{N}=2$ $SU(N)$ supersymmetric gauge theories, $N = 2, 3, \ldots$, on ${\mathbb {C}}^2/{\mathbb {Z}}_n$, $n = 2, 3, \ldots$, with $\Omega$-deformation and 2D conformal field theories with $\mathcal{W}^{\, para}_{N, n}$ ($n$-th parafermion $\mathcal{W}_N$) symmetry and $\widehat{\mathfrak{sl}}(n)_N$ symmetry. In this work, we trivialize the factor with $\mathcal{W}^{\, para}_{N, n}$ symmetry in the 4D $SU(N)$ instanton partition functions on ${\mathbb {C}}^2/{\mathbb {Z}}_n$ (by using specific choices of parameters and imposing specific conditions on the $N$-tuples of Young diagrams that label the states), and extract the 2D $\widehat{\mathfrak{sl}}(n)_N$ WZW conformal blocks, $n = 2, 3, \ldots$, $N = 1, 2, \ldots\, .$


title: Exploiting Model Sparsity in Adaptive MPC: A Compressed Sensing Viewpoint

abstract: This paper proposes an Adaptive Stochastic Model Predictive Control (MPC) strategy for stable linear time-invariant systems in the presence of bounded disturbances. We consider multi-input, multi-output systems that can be expressed by a Finite Impulse Response (FIR) model. The parameters of the FIR model corresponding to each output are unknown but assumed sparse. We estimate these parameters using the Recursive Least Squares algorithm. The estimates are then improved using set-based bounds obtained by solving the Basis Pursuit Denoising [1] problem. Our approach is able to handle hard input constraints and probabilistic output constraints. Using tools from distributionally robust optimization, we reformulate the probabilistic output constraints as tractable convex second-order cone constraints, which enables us to pose our MPC design task as a convex optimization problem. The efficacy of the developed algorithm is highlighted with a thorough numerical example, where we demonstrate performance gain over the counterpart algorithm of [2], which does not utilize the sparsity information of the system impulse response parameters during control design.


title: Electromagnetic calculations for multiscale and multiphysics simulations: a new perspective

abstract: Present day electromagnetic field calculations have limitations that are due to techniques employing edge-based discretization methods. While these vector finite element methods solve the issues of tangential continuity of fields and the removal of spurious solutions, resulting fields do not have a unique directionality at nodes in the discretization mesh. This review presents electromagnetic field calculations in waveguides, cavity fields, and photonic crystals. We develop Hermite interpolation polynomials and node-based finite element methods based on variational principles. We show that the Hermite-finite element method (HFEM) delivers high accuracy and suitable for multiscale calculations with mixed physics. We use group representation theory to derive the HFEM polynomial basis set in two-dimensions. The energy level degeneracy in a cubic cavity can be denumerably large even though the symmetry of the cube. We show that the additional operators available for the problem lead to accidental degeneracy. We discuss this remarkable degeneracy and its reduction in detail. We consider photonic crystals corresponding to a 2D checkerboard superlattice structure, and the Escher drawing of the Horsemen which satisfies the nonsymmorphic group pg. We show that HFEM is able to deliver high accuracy in such spatially complex examples with far less computational effort than Fourier expansion methods. The algorithms developed here hold the promise of successful modeling of multi-physics systems. This general method is applicable to a broad class of physical systems, to semiconducting lasers which require simultaneous modeling of transitions in quantum wells or dots together with EM cavity calculations, to modeling plasmonic structures in the presence of EM field emissions, and to on-chip propagation within monolithic integrated circuits.

categories: physics.comp-ph cond-mat.mes-hall physics.app-ph physics.optics


title: Ionization correction factors for sodium, potassium, and calcium in planetary nebulae

abstract: We use a large grid of photoionization models that are representative of observed planetary nebulae (PNe) to derive ionization correction factors (ICFs) for sodium, potassium, and calcium. In addition to the analytical expressions of the ICFs, we provide the range of validity where the ICFs can be safely used and an estimate of the typical uncertainties associated with the ICFs. We improved the previous ICFs for calcium and potassium in the literature and suggest for the first time an ICF for sodium. We tested our ICFs with a sample of 39 PNe with emission lines of some ion of these elements. No obvious trend is found between the derived abundances and the degree of ionization, suggesting that our ICFs do not seem to be introducing an artificial bias in the results. The abundances found in the studied PNe range from -2.88$^{+0.21}_{-0.22}$ to -2.09$\pm$0.21 in log(Na/O), from -4.20$^{+0.31}_{-0.45}$ to -3.05$^{+0.26}_{-0.47}$ in log(K/O), and from -3.71$^{+0.41}_{-0.34}$ to -1.57$^{+0.33}_{-0.47}$ in log(Ca/O). These numbers imply that some of the studied PNe have up to 65, 75, or 95 per cent of their Na, K, and/or Ca atoms condensed into dust grains, respectively. As expected, the highest depletions are found for calcium which is the element with the highest condensation temperature.


title: A Revisitation of Low-Rate Bounds on the Reliability Function of Discrete Memoryless Channels for List Decoding

abstract: We revise the proof of low-rate upper bounds on the reliability function of discrete memoryless channels for ordinary and list-decoding schemes, in particular Berlekamp and Blinovsky's zero-rate bound, as well as Blahut's bound for low rates. The available proofs of the zero-rate bound devised by Berlekamp and Blinovsky are somehow complicated in that they contain in one form or another some cumbersome "non-standard" procedures or computations. Here we follow Blinovsky's idea of using a Ramsey-theoretic result by Komlos, and we complement it with some missing steps to present a proof which is rigorous and easier to inspect. Furthermore, we show how these techniques can be used to fix an error that invalidated the proof of Blahut's low-rate bound, which is here presented in an extended form for list decoding and for general channels.


title: Probabilistic cosmic web classification using fast-generated training data

abstract: We present a novel method of robust probabilistic cosmic web particle classification in three dimensions using a supervised machine learning algorithm. Training data was generated using a simplified $\Lambda$CDM toy model with pre-determined algorithms for generating halos, filaments, and voids. While this framework is not constrained by physical modeling, it can be generated substantially more quickly than an N-body simulation without loss in classification accuracy. For each particle in this dataset, measurements were taken of the local density field magnitude and directionality. These measurements were used to train a random forest algorithm, which was used to assign class probabilities to each particle in a $\Lambda$CDM, dark matter-only N-body simulation with $256^3$ particles, as well as on another toy model data set. By comparing the trends in the ROC curves and other statistical metrics of the classes assigned to particles in each dataset using different feature sets, we demonstrate that the combination of measurements of the local density field magnitude and directionality enables accurate and consistent classification of halo, filament, and void particles in varied environments We also show that this combination of training features ensures that the construction of our toy model does not affect classification. The use of a fully supervised algorithm allows greater control over the information deemed important for classification, preventing issues arising from hyperparameters and mode collapse in deep learning models. Due to the speed of training data generation, our method is highly scalable, making it particularly suited for classifying large datasets, including observed data.


title: The Spectral Hecke Algebra

abstract: We introduce a derived enhancement of local Galois deformation rings that we call the "spectral Hecke algebra", in analogy to a construction in the Geometric Langlands program. This is a Hecke algebra that acts on the spectral side of the Langlands correspondence, i.e. on moduli spaces of Galois representations. We verify the simplest form of derived local-global compatibility between the action of the spectral Hecke algebra on the derived Galois deformation ring of Galatius-Venkatesh, and the action of Venkatesh's derived Hecke algebra on the cohomology of arithmetic groups.


title: CCCP and MENeaCS: (updated) weak-lensing masses for 100 galaxy clusters

abstract: Large area surveys have detected significant samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at $0.05<z<0.15$, observed as part of the Multi Epoch Nearby Cluster Survey (MENeaCS), and reevaluate the mass estimates for 52 clusters from the Canadian Cluster Comparison Project (CCCP). Updated high-fidelity photometric redshift catalogues of reference deep fields are used in combination with advances in shape measurements and state-of-the-art cluster simulations, yielding an average systematic uncertainty in the lensing signal below 5%, similar to the statistical uncertainty for our cluster sample. We derive a scaling relation with Planck measurements for the full sample and find a bias in the Planck masses of $1-b=0.84 \pm 0.04$. We find no statistically significant trend of the mass bias with redshift or cluster mass, but find that different selections could change the bias by up to 1.5$\sigma$. We find a gas fraction of $0.139 \pm 0.014$ for 8 relaxed clusters in our sample, which can also be used to infer cosmological parameters.


title: Tuning spatial entanglement in interacting few-electron quantum dots

abstract: Confined geometries such as semiconductor quantum dots are promising candidates for fabricating quantum computing devices. When several quantum dots are in proximity, spatial correlation between electrons in the system becomes significant. In this article, we develop a fully variational action integral formulation for calculating accurate few-electron wavefunctions in configuration space, irrespective of potential geometry. To evaluate the Coulomb integrals with high accuracy, a novel numerical integration method using multiple Gauss quadratures is proposed. Using this approach, we investigate the confinement of two electrons in double quantum dots, and evaluate the spatial entanglement. We investigate the dependence of spatial entanglement on various geometrical parameters. We derive the two-particle wavefunctions in the asymptotic limit of the separation distance between quantum dots, and obtain universal saturation values for the spatial entanglement. Resonances in the entanglement values due to avoided level-crossings of states are observed. We also demonstrate the formation of electron clusters, and show that the entanglement value is a good indicator for the formation of such clusters. Further, we show that a precise tuning of the entanglement values is feasible with applied external electric fields.


title: Structural Properties of the Stability of Jamitons

abstract: It is known that inhomogeneous second-order macroscopic traffic models can reproduce the phantom traffic jam phenomenon: whenever the sub-characteristic condition is violated, uniform traffic flow is unstable, and small perturbations grow into nonlinear traveling waves, called jamitons. In contrast, what is essentially unstudied is the question: which jamiton solutions are dynamically stable? To understand which stop-and-go traffic waves can arise through the dynamics of the model, this question is critical. This paper first presents a computational study demonstrating which types of jamitons do arise dynamically, and which do not. Then, a procedure is presented that characterizes the stability of jamitons. The study reveals that a critical component of this analysis is the proper treatment of the perturbations to the shocks, and of the neighborhood of the sonic points.


title: Bargmann transfoms associated with reproducing kernel Hilbert space and application to Dirichlet spaces

abstract: The aim of the present paper is three folds. For a reproducing kernel Hilbert space $\mathcal{A}$ (R.K.H.S) and a $\sigma-$finite measure space $(M_{1},d\mu_{1})$ for which the corresponding $L^{2}-$space is a separable Hilbert space, we first build an isometry of Bargmann type as an integral transform from $L^{2}(M_{1},d\mu_{1})$ into $\mathcal{A}$. Secondly, in the case where there exists a $\sigma-$finite measure space $(M_{2},d\mu_{2})$ such that the Hilbert space $L^{2}(M_{2},d\mu_{2})$ is separable and $\mathcal{A}\subset L^{2}(M_{2},d\mu_{2})$ the inverse isometry is also given in an explicit form as an integral transform. As consequence, we recover some classical isometries of Bargmann type. Thirdly, for the classical Dirichlet space as R.K.H.S, we elaborate a new isometry of Bargmann type. Furthermore, for this Dirichlet space, we give a new characterization, as harmonic space of a single second order elliptic partial differential operator for which, we present some spectral properties. Finally, we extend the same results to a class of generalized Bergman-Dirichlet space.


title: Deep Autoencoders with Value-at-Risk Thresholding for Unsupervised Anomaly Detection

abstract: Many real-world monitoring and surveillance applications require non-trivial anomaly detection to be run in the streaming model. We consider an incremental-learning approach, wherein a deep-autoencoding (DAE) model of what is normal is trained and used to detect anomalies at the same time. In the detection of anomalies, we utilise a novel thresholding mechanism, based on value at risk (VaR). We compare the resulting convolutional neural network (CNN) against a number of subspace methods, and present results on changedetection net.


title: Analysis of the Ethiopic Twitter Dataset for Abusive Speech in Amharic

abstract: In this paper, we present an analysis of the first Ethiopic Twitter Dataset for the Amharic language targeted for recognizing abusive speech. The dataset has been collected since 2014 that is written in Fidel script. Since several languages can be written using the Fidel script, we have used the existing Amharic, Tigrinya and Ge'ez corpora to retain only the Amharic tweets. We have analyzed the tweets for abusive speech content with the following targets: Analyze the distribution and tendency of abusive speech content over time and compare the abusive speech content between a Twitter and general reference Amharic corpus.


title: Multiplication-invariant operators and the classification of LCA group frames

abstract: In this paper we study the properties of multiplication invariant (MI) operators acting on subspaces of the vector-valued space $L^2(X;\mathcal H)$. We characterize such operators in terms of range functions by showing that there is an isomorphism between the category of MI spaces (with MI operators as morphisms) and the category of measurable range functions whose morphisms are measurable range operators. We investigate how global properties of an MI operator are reflected by local pointwise properties of its corresponding range operator. We also establish several results about frames generated by multiplications in $L^2(X;\mathcal H)$. This includes the classification of frames of multiplications with respect to unitary equivalence by measurable fields of Gramians. Finally, we show applications of our results in the study of abelian group frames and translation-invariant (TI) operators acting on subspaces of $L^2(\mathcal G)$, where $\mathcal G$ is a locally compact group.


title: Basis Prediction Networks for Effective Burst Denoising with Large Kernels

abstract: Bursts of images exhibit significant self-similarity across both time and space. This motivates a representation of the kernels as linear combinations of a small set of basis elements. To this end, we introduce a novel basis prediction network that, given an input burst, predicts a set of global basis kernels -- shared within the image -- and the corresponding mixing coefficients -- which are specific to individual pixels. Compared to state-of-the-art techniques that output a large tensor of per-pixel spatiotemporal kernels, our formulation substantially reduces the dimensionality of the network output. This allows us to effectively exploit comparatively larger denoising kernels, achieving both significant quality improvements (over 1dB PSNR) and faster run-times over state-of-the-art methods.


title: On the mistake in defining fractional derivative using a non-singular kernel

abstract: Definitions of fractional derivative of order $\alpha$ ($0 < \alpha \leq 1$) using non-singular kernels have been recently proposed. In this note we show that these definitions cannot be useful in modelling problems with a initial value condition (like, for example, the fractional diffusion equation) because the solutions obtained for these equations do not satisfy the initial condition (except for the integer case $\alpha = 1$). In order to satisfy an arbitrary initial condition the definitions of fractional derivative must necessarily involve a singular kernel.


title: Robust, Extensible, and Fast: Teamed Classifiers for Vehicle Tracking and Vehicle Re-ID in Multi-Camera Networks

abstract: As camera networks have become more ubiquitous over the past decade, the research interest in video management has shifted to analytics on multi-camera networks. This includes performing tasks such as object detection, attribute identification, and vehicle/person tracking across different cameras without overlap. Current frameworks for management are designed for multi-camera networks in a closed dataset environment where there is limited variability in cameras and characteristics of the surveillance environment are well known. Furthermore, current frameworks are designed for offline analytics with guidance from human operators for forensic applications. This paper presents a teamed classifier framework for video analytics in heterogeneous many-camera networks with adversarial conditions such as multi-scale, multi-resolution cameras capturing the environment with varying occlusion, blur, and orientations. We describe an implementation for vehicle tracking and vehicle re-identification (re-id), where we implement a zero-shot learning (ZSL) system that performs automated tracking of all vehicles all the time. Our evaluations on VeRi-776 and Cars196 show the teamed classifier framework is robust to adversarial conditions, extensible to changing video characteristics such as new vehicle types/brands and new cameras, and offers real-time performance compared to current offline video analytics approaches.


title: Implementation of the XY interaction family with calibration of a single pulse

abstract: Near-term applications of quantum information processors will rely on optimized circuit implementations to minimize gate depth and therefore mitigate the impact of gate errors in noisy intermediate-scale quantum (NISQ) computers. More expressive gate sets can significantly reduce the gate depth of generic circuits. Similarly, structured algorithms can benefit from a gate set that more directly matches the symmetries of the problem. The XY interaction generates a family of gates that provides expressiveness well tailored to quantum chemistry as well as to combinatorial optimization problems, while also offering reductions in circuit depth for more generic circuits. Here we implement the full family of XY entangling gates in a transmon-based superconducting qubit architecture. We use a composite pulse scheme that requires calibration of only a single gate pulse and maintains constant gate time for all members of the family. This allows us to maintain a high fidelity implementation of the gate across all entangling angles. The average fidelity of gates sampled from this family ranges from $95.67 \pm 0.60\%$ to $99.01 \pm 0.15\%$, with a median fidelity of $97.35 \pm 0.17\%$, which approaches the coherence-limited gate fidelity of the qubit pair. We furthermore demonstrate the utility of XY in a quantum approximation optimization algorithm in enabling circuit depth reductions as compared to the CZ only case.


title: Renormalization of a Contact Interaction on a Lattice

abstract: Contact interactions can be used to describe a system of particles at unitarity, contribute to the leading part of nuclear interactions and are numerically non-trivial because they require a proper regularization and renormalization scheme. We explain how to tune the coefficient of a contact interaction between non-relativistic particles on a discretized space in 1, 2, and 3 spatial dimensions such that we can remove all discretization artifacts. By taking advantage of a latticized L\"uscher zeta function, we can achieve a momentum-independent scattering amplitude at any finite lattice spacing.


title: Turning chiral skyrmion inside out

abstract: The stability of two-dimensional chiral skyrmions in a tilted magnetic field is studied. It is shown that by changing the direction of the field and its magnitude, one can continuously transform chiral skyrmion into a skyrmion with opposite polarity and vorticity. This turned inside out skyrmion can be considered as an antiparticle for ordinary axisymmetric skyrmion. For any tilt angle of the magnetic field, there is a range of its absolute values where two types of skyrmions may coexist. In a tilted field, the potentials for inter-skyrmion interactions are characterized by the presence of local minima suggesting attractive interaction between the particles. The potentials of inter-particle interactions also have so-called fusion channels allowing either annihilation of two particles or the emergence of a new particle. The presented results are general for a wide class of magnetic crystals with both easy-plane and easy-axis anisotropy.


title: Winning the Lottery with Continuous Sparsification

abstract: The search for efficient, sparse deep neural network models is most prominently performed by pruning: training a dense, overparameterized network and removing parameters, usually via following a manually-crafted heuristic. Additionally, the recent Lottery Ticket Hypothesis conjectures that, for a typically-sized neural network, it is possible to find small sub-networks which, when trained from scratch on a comparable budget, match the performance of the original dense counterpart. We revisit fundamental aspects of pruning algorithms, pointing out missing ingredients in previous approaches, and develop a method, Continuous Sparsification, which searches for sparse networks based on a novel approximation of an intractable $\ell_0$ regularization. We compare against dominant heuristic-based methods on pruning as well as ticket search -- finding sparse subnetworks that can be successfully re-trained from an early iterate. Empirical results show that we surpass the state-of-the-art for both objectives, across models and datasets, including VGG trained on CIFAR-10 and ResNet-50 trained on ImageNet. In addition to setting a new standard for pruning, Continuous Sparsification also offers fast parallel ticket search, opening doors to new applications of the Lottery Ticket Hypothesis.


title: Generic Rotation Sets

abstract: Let $(X,T)$ be a topological dynamical system. Given a continuous vector-valued function $F \in C(X, \mathbb{R}^{d})$ called a potential we define its rotation set $R(F)$ as the set of integrals of $F$ with respect to all $T$-invariant probability measures, which is a convex body of $\mathbb{R}^{d}$. In this paper, we study the geometry of rotation sets. We prove that if $T$ is a non-uniquely ergodic topological dynamical system with a dense set of periodic measures, then the map $R(\cdot)$ is open with respect to the uniform topologies. As a consequence, we obtain that the rotation set of a generic potential is strictly convex and has $C^{1}$ boundary. Furthermore, we prove that the map $R(\cdot)$ is surjective, extending a result of Kucherenko and Wolf.


title: Formation of point shocks for 3D compressible Euler

abstract: We consider the 3D isentropic compressible Euler equations with the ideal gas law. We provide a constructive proof of shock formation from smooth initial datum of finite energy, with no vacuum regions, with nontrivial vorticity present at the shock, and under no symmetry assumptions. We prove that for an open set of Sobolev-class initial data which are a small $L^ \infty $ perturbation of a constant state, there exist smooth solutions to the Euler equations which form a generic stable shock in finite time. The blow up time and location can be explicitly computed, and solutions at the blow up time are smooth except for a single point, where they are of cusp-type with H\"{o}lder $C^ {\frac{1}{3}}$ regularity. Our proof is based on the use of modulated self-similar variables that are used to enforce a number of constraints on the blow up profile, necessary to establish the stability in self-similar variables of the generic shock profile.


title: HalluciNet-ing Spatiotemporal Representations Using a 2D-CNN

abstract: Spatiotemporal representations learned using 3D convolutional neural networks (CNN) are currently used in state-of-the-art approaches for action related tasks. However, 3D-CNN are notorious for being memory and compute resource intensive as compared with more simple 2D-CNN architectures. We propose to hallucinate spatiotemporal representations from a 3D-CNN teacher with a 2D-CNN student. By requiring the 2D-CNN to predict the future and intuit upcoming activity, it is encouraged to gain a deeper understanding of actions and how they evolve. The hallucination task is treated as an auxiliary task, which can be used with any other action related task in a multitask learning setting. Thorough experimental evaluation shows that the hallucination task indeed helps improve performance on action recognition, action quality assessment, and dynamic scene recognition tasks. From a practical standpoint, being able to hallucinate spatiotemporal representations without an actual 3D-CNN can enable deployment in resource-constrained scenarios, such as with limited computing power and/or lower bandwidth. Codebase is available here: https://github.com/ParitoshParmar/HalluciNet.


title: Swift monitoring of M51: A 38-day super-orbital period for the pulsar ULX7 and a new transient ULX

abstract: We present the results from a monitoring campaign made with the Neil Gehrels Swift Observatory of the M51 galaxies, which contain several variable ultraluminous X-ray sources (ULXs). The ongoing campaign started in May 2018, and we report here on $\sim1.5$ years of observations. The campaign, which consists of 105 observations, has a typical cadence of 3--6 days, and has the goal of determining the long-term X-ray variability of the ULXs. Two of the most variable sources were ULX7 and ULX8, both of which are known to be powered by neutron stars that are exceeding their isotropic Eddington luminosities by factors of up to 100. This is further evidence that neutron star powered ULXs are the most variable. Our two main results are, first, that ULX7 exhibits a periodic flux modulation with a period of 38 days varying over a magnitude and a half in flux from peak to trough. Since the orbital period of the system is known to be 2 days, the modulation is super-orbital, which is a near-ubiquitous property of ULX pulsars. Secondly we identify a new transient ULX, M51 XT-1, the onset of which occurred during our campaign, reaching a peak luminosity of $\sim10^{40}$ erg s$^{-1}$, before gradually fading over the next $\sim200$ days until it slipped below the detection limit of our observations. Combined with the high-quality Swift/XRT lightcurve of the transient, serendipitous observations made with Chandra and XMM-Newton provide insights into the onset and evolution of a likely super-Eddington event.


title: Variable selection for transportability

abstract: Transportability provides a principled framework to address the problem of applying study results to new populations. Here, we consider the problem of selecting variables to include in transport estimators. We provide a brief overview of the transportability framework and illustrate that while selection diagrams are a vital first step in variable selection, these graphs alone identify a sufficient but not strictly necessary set of variables for generating an unbiased transport estimate. Next, we conduct a simulation experiment assessing the impact of including unnecessary variables on the performance of the parametric g-computation transport estimator. Our results highlight that the types of variables included can affect the bias, variance, and mean squared error of the estimates. We find that addition of variables that are not causes of the outcome but whose distributions differ between the source and target populations can increase the variance and mean squared error of the transported estimates. On the other hand, inclusion of variables that are causes of the outcome (regardless of whether they modify the causal contrast of interest or differ in distribution between the populations) reduces the variance of the estimates without increasing the bias. Finally, exclusion of variables that cause the outcome but do not modify the causal contrast of interest does not increase bias. These findings suggest that variable selection approaches for transport should prioritize identifying and including all causes of the outcome in the study population rather than focusing on variables whose distribution may differ between the study sample and target population.


title: Landau Levels of Topologically-Protected Surface States Probed by Dual-Gated Quantum Capacitance

abstract: Spectroscopy of discrete Landau levels (LLs) in bulk-insulating three-dimensional topological insulators (3D TIs) in perpendicular magnetic field characterizes the Dirac nature of their surface states. Despite a number of studies demonstrating the quantum Hall effect (QHE) of topological surface states, quantitative evaluation of the LL energies, which serve as fundamental electronic quantities for study of the quantum states, is still limited. In this work, we explore the density of states of LLs by measuring quantum capacitance (CQ) in a truly bulk insulating 3D TI via a van der Waals heterostructure configuration. By applying dual gate voltages, we access the individual surface states' LLs and extract their chemical potentials to quantify the LL spacings of each surface. We evaluate the LLs' energies at two distinguished QH states, namely dissipationless ({\nu}= +/-1) and dissipative ({\nu}= 0) states in the 3D TI.


title: Generating nonequilibrium stationary state from ground state condensate through an almost-adiabatic cycle

abstract: It is shown that the ground state of weakly interacting Bose particles in a quasi one-dimensional box trap can be converted into an excited stationary state by an adiabatic cyclic operation that involves a quench of interaction strength: A sharp impurity potential is applied and its strength is varied during the cycle, which induces a nonequilibrium stationary state exhibiting the inversion of population. This process is robust in the sense that the resultant stationary state is almost independent of the details of the cycle, such as the position of the impurity as long as the cycle is far enough from critical regions. The case of the failure of the population inversion due to the strong interparticle interactions is also examined.


title: Stylised Choropleth Maps for New Zealand Regions and District Health Boards

abstract: New Zealand has two top-level sets of administrative divisions: the District Health Boards and the Regions. In this note I describe a hexagonal layout for creating stylised maps of these divisions, and using colour, size, and triangular subdivisions to compare data between divisions and across multiple variables. I present an implementation in the DHBins package for R using both base graphics and ggplot2; the concepts and specific hexagonal layout could be used in any software.


title: A new bound on the acyclic edge chromatic index

abstract: In this note we obtain a new bound for the acyclic edge chromatic number $a'(G)$ of a graph $G$ with maximum degree $D$ proving that $a'(G)\leq 3.569(D-1)$. To get this result we revisit and slightly modify the method described in [Giotis, Kirousis, Psaromiligkos and Thilikos, Theoretical Computer Science, 66: 40-50, 2017].


title: Design Trade-offs for Decentralized Baseband Processing in Massive MU-MIMO Systems

abstract: Massive multi-user (MU) multiple-input multiple-output (MIMO) provides high spectral efficiency by means of spatial multiplexing and fine-grained beamforming. However, conventional base-station (BS) architectures for systems with hundreds of antennas that rely on centralized baseband processing inevitably suffer from (i) excessive interconnect data rates between radio-frequency circuitry and processing fabrics, and (ii) prohibitive complexity at the centralized baseband processor. Recently, decentralized baseband processing (DBP) architectures and algorithms have been proposed, which mitigate the interconnect bandwidth and complexity bottlenecks. This paper systematically explores the design trade-offs between error-rate performance, computational complexity, and data transfer latency of DBP architectures under different system configurations and channel conditions. Considering architecture, algorithm, and numerical precision aspects, we provide practical guidelines to select the DBP architecture and algorithm that are able to realize the full benefits of massive MU-MIMO in the uplink and downlink.


title: Accurate determination of the translational correlation function of two-dimensional solids

abstract: The identification of the different phases of a two-dimensional (2d) system, which might be in solid, hexatic, or liquid, requires the accurate determination of the correlation function of the translational and of the bond-orientational order parameters. According to the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) theory, in the solid phase the translational correlation function decays algebraically, as a consequence of the Mermin-Wagner long-wavelength fluctuations. Recent results have however reported an exponential-like decay. By revisiting different definitions of the translational correlation function commonly used in the literature, here we clarify that the observed exponential-like decay in the solid phase results from an inaccurate determination of the symmetry axis of the solid; the expected power-law behaviour is recovered when the symmetry axis is properly identified. We show that, contrary to the common assumption, the symmetry axis of a 2d solid is not fixed by the direction of its global bond-orientational parameter, and introduce an approach allowing to determine the symmetry axis from a real space analysis of the sample.


title: Privacy-preserving data sharing via probabilistic modelling

abstract: Differential privacy allows quantifying privacy loss resulting from accessing sensitive personal data. Repeated accesses to underlying data incur increasing loss. Releasing data as privacy-preserving synthetic data would avoid this limitation, but would leave open the problem of designing what kind of synthetic data. We propose formulating the problem of private data release through probabilistic modelling. This approach transforms the problem of designing the synthetic data into choosing a model for the data, allowing also including prior knowledge, which improves the quality of the synthetic data. We demonstrate empirically, in an epidemiological study, that statistical discoveries can be reliably reproduced from the synthetic data. We expect the method to have broad use in creating high-quality anonymized data twins of key data sets for research.