title: Asymptotically matched quasi-circular inspiral and transition-to-plunge in the small mass ratio expansion

abstract: In the small mass ratio expansion and on the equatorial plane, the two-body problem for point particles in general relativity admits a quasi-circular inspiral motion followed by a transition-to-plunge motion. We first derive the equations governing the quasi-circular inspiral in the Kerr background at adiabatic, post-adiabatic and post-post-adiabatic orders in the slow-timescale expansion in terms of the self-force and we highlight the structure of the equations of motion at higher subleading orders. We derive in parallel the equations governing the transition-to-plunge motion to any subleading order, and demonstrate that they are governed by sourced linearized Painlev\'e transcendental equations of the first kind. The first ten perturbative orders do not require any further developments in self-force theory, as they are determined by the second-order self-force. We propose a scheme that matches the slow-timescale expansion of the inspiral with the transition-to-plunge motion to all perturbative orders in the overlapping region exterior to the last stable orbit where both expansions are valid. We explicitly verify the validity of the matching conditions for a large set of coefficients involved, on the one hand, in the adiabatic or post-adiabatic inspiral and, on the other hand, in the leading, subleading or higher subleading transition-to-plunge motion. This result is instrumental for deriving gravitational waveforms within the self-force formalism beyond the innermost stable circular orbit.


title: Astraeus V: The emergence and evolution of metallicity scaling relations during the Epoch of Reionization

abstract: In this work, we have implemented a detailed physical model of galaxy chemical enrichment into the ${\it Astraeus}$ (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dark matter simUlationS) framework which couples galaxy formation and reionization in the first billion years. Simulating galaxies spanning over 2.5 orders of magnitude in halo mass with $M_h \sim 10^{8.9-11.5} M_\odot$ ($M_h \sim 10^{8.9-12.8} M_\odot$) at $z \sim 10 ~ (5)$, we find: (i) smooth-accretion of metal-poor gas from the intergalactic medium (IGM) plays a key role in diluting the interstellar medium (ISM) metallicity which is effectively restored due to self-enrichment from star formation; (ii) a redshift averaged gas-mass loading factor that depends on the stellar mass as $\eta_g \approx 1.38 ({M_*}/{10^{10} M_\odot})^{-0.43}$; (iii) the mass-metallicity relation is already in place at $z \sim 10$ and shows effectively no redshift evolution down to $z \sim 5$; (iv) for a given stellar mass, the metallicity decreases with an increase in the star formation rate (SFR); (v) the key properties of the gas-phase metallicity (in units of 12+log(O/H), stellar mass, SFR and redshift are linked through a high-redshift fundamental plane of metallicity (HFPZ) for which we provide a functional form; (vi) the mass-metallicity-SFR relations are effectively independent of the reionization radiative feedback model for $M_* \geq 10^{6.5} M_\odot$ galaxies; (vii) while low-mass galaxies ($M_h \leq 10^9 M_\odot$) are the key contributors to the metal budget of the IGM at early times, higher mass halos provide about 50% of the metal budget at lower-redshifts.


title: Reconstruction of bottomonium spectral functions in thermal QCD using Kernel Ridge Regression

abstract: We discuss results for bottomonium at nonzero temperature obtained using NRQCD on Fastsum Generation 2L ensembles, as part of the Fastsum collaboration's programme to determine the spectrum of the bottomonium system as a function of temperature using a variety of approaches. Here we give an update on results for spectral functions obtained using Kernel Ridge Regression. We pay in particular attention to the generation of training data and introduce the notion of using lattice QCD ensembles to learn how to improve the generation of training data. A practical implementation is given.


title: APOGEE Detection of N-rich stars in the tidal tails of Palomar 5

abstract: Recent results from chemical tagging studies using APOGEE data suggest a strong link between the chemical abundance patterns of stars found within globular clusters, and chemically peculiar populations in the Galactic halo field. In this paper we analyse the chemical compositions of stars within the cluster body and tidal streams of Palomar 5, a globular cluster that is being tidally disrupted by interaction with the Galactic gravitational potential. We report the identification of nitrogen-rich (N-rich) stars both within and beyond the tidal radius of Palomar 5, with the latter being clearly aligned with the cluster tidal streams; this acts as confirmation that N-rich stars are lost to the Galactic halo from globular clusters, and provides support to the hypothesis that field N-rich stars identified by various groups have a globular cluster origin.


title: Relational observables in Asymptotically safe gravity

abstract: We introduce an approach to compute the renormalisation group flow of relational observables in quantum gravity which evolve from their microscopic expressions towards the full quantum expectation value. This is achieved by using the composite operator formalism of the functional renormalisation group. These methods can be applied to a large class of relational observables within a derivative expansion for different physical coordinate systems. As a first application we consider four scalar fields coupled to gravity to represent the physical coordinate frame from which relational observables can be constructed. At leading order of the derivative expansion the observables are the inverse relational metric and the relational scalar curvature. We evaluate their scaling dimensions at the fixed point, both in the standard renormalisation group scheme and in the essential scheme. This represents the first steps to describe running observables within asymptotic safety; this treatment can be generalised to other observables constructed from different tensors and in different physical coordinate systems.


title: Global Natural Orbital Functional: Towards the Complete Description of the Electron Correlation

abstract: The current work presents a natural orbital functional (NOF) for electronic systems with any spin value independent of the external potential being considered, that is, a global NOF (GNOF). It is based on a new two-index reconstruction of the two-particle reduced density matrix for spin multiplets. The emergent functional describes the complete intrapair electron correlation, and the correlation between orbitals that make up both the pairs and the individual electrons. The interorbital correlation is composed of static and dynamic terms. The concept of dynamic part of the occupation numbers is introduced. To evaluate the accuracy achieved with GNOF, calculation of a variety of properties is presented. They include the total energies and energy differences between the ground state and the lowest-lying excited state with different spin of atoms from H to Ne, ionization potentials of the first-row transition-metal atoms (Sc-Zn), and the total energies of a selected set of 55 molecular systems in different spin states. GNOF is also applied to the homolytic dissociation of selected diatomic molecules in different spin states and to the rotation barrier of ethylene, both paradigmatic cases of systems with significant multi-configurational character. The values obtained agree with those reported at high level of theory and experimental data.


title: Prediction and compression of lattice QCD data using machine learning algorithms on quantum annealer

abstract: We present regression and compression algorithms for lattice QCD data utilizing the efficient binary optimization ability of quantum annealers. In the regression algorithm, we encode the correlation between the input and output variables into a sparse coding machine learning algorithm. The trained correlation pattern is used to predict lattice QCD observables of unseen lattice configurations from other observables measured on the lattice. In the compression algorithm, we define a mapping from lattice QCD data of floating-point numbers to the binary coefficients that closely reconstruct the input data from a set of basis vectors. Since the reconstruction is not exact, the mapping defines a lossy compression, but, a reasonably small number of binary coefficients are able to reconstruct the input vector of lattice QCD data with the reconstruction error much smaller than the statistical fluctuation. In both applications, we use D-Wave quantum annealers to solve the NP-hard binary optimization problems of the machine learning algorithms.


title: Generalized geodesic deviation in de Sitter spacetime

abstract: The geodesic deviation equation (GDE) describes the tendency of objects to accelerate towards or away from each other due to spacetime curvature. The GDE assumes that nearby geodesics have a small rate of separation, which is formally treated as the same order in smallness as the separation itself. This assumption is discussed in various papers, but is typically not recognized in textbooks. Relaxing this assumption leads to the generalized geodesic deviation equation (GGDE). We demonstrate the distinction between the GDE and the GGDE by computing the relative acceleration between timelike geodesics in two-dimensional de Sitter spacetime. We do this by considering a fiducial geodesic and a secondary geodesic (both timelike) that cross with nonzero speed. These geodesics are spanned by a spacelike geodesic, whose tangent evaluated at the fiducial geodesic defines the separation. The second derivative of the separation describes the relative acceleration between the fiducial and secondary geodesics. Near the crossing point, where the separation between the timelike geodesics is small but their rates of separation can be large, we show that the GGDE holds but the GDE fails to apply.


title: Moduli spaces of BPS lumps with holomorphic impurities

abstract: A self-dual generalization of the lump-impurity system is introduced. This model possesses lump-antilump-like pairs as static solutions of the pertinent Bogomolny equations. This allows for a moduli space approximation analysis of the BPS solutions which are identified as lump-antilump configurations. Some geometrical properties of the resulting moduli are analyzed. In addition, it is argued that, this type of impurity models can be interpreted as a limit of certain non-impurity theories.


title: Phonon anharmonicity in Cu-based layered thiophosphates

abstract: In this work, we performed temperature-dependent Raman spectroscopy studies on three layered metal thiophosphates: CuInP2S6, CuCrP2S6, and CuInP2S6-In4/3P2S6 heterostructures. These materials are emerging multiferroics, with CuInP2S6 and CuInP2S6-In4/3P2S6 exhibiting room temperature ferrielectric polarization while CuCrP2S6 is both anti-ferroelectric and antiferromagnetic at low temperatures (<145 and <30 K, respectively). We studied two prominent peaks in the Raman spectra of these compounds, namely their P-P and P-S stretching modes. First order linear fits to the mode softening yielded coefficients similar to other layered materials. Furthermore, analysis of the phonon decay revealed the dominant role of anharmonic effects rather than thermal expansion. In particular, the decay of the optical phonons was found to be governed by three-phonon processes in CuInP2S6 and CuCrP2S6 and by the four-phonon process in the CuInP2S6-In4/3P2S6 heterostructures. These phonon lifetimes were found to be less than 1 ps, indicating high scattering rates and showing the potential for these materials to be used in thermoelectric applications.


title: Dissipation-Induced Order: The $S=1/2$ Quantum Spin Chain Coupled to an Ohmic Bath

abstract: We consider an $S=1/2$ antiferromagnetic quantum Heisenberg chain where each site is coupled to an independent bosonic bath with ohmic dissipation. The coupling to the bath preserves the global SO(3) spin symmetry. Using large-scale, approximation-free quantum Monte Carlo simulations, we show that any finite coupling to the bath suffices to stabilize long-range antiferromagnetic order. This is in stark contrast to the isolated Heisenberg chain where spontaneous breaking of the SO(3) symmetry is forbidden by the Mermin-Wagner theorem. A linear spin-wave theory analysis confirms that the memory of the bath and the concomitant retarded interaction stabilize the order. For the Heisenberg chain, the ohmic bath is a marginal perturbation so that exponentially large system sizes are required to observe long-range order at small couplings. Below this length scale, our numerics is dominated by a crossover regime where spin correlations show different power-law behaviors in space and time. We discuss the experimental relevance of this crossover phenomena.


title: Examining Zero-Shot Vulnerability Repair with Large Language Models

abstract: Human developers can produce code with cybersecurity bugs. Can emerging 'smart' code completion tools help repair those bugs? In this work, we examine the use of large language models (LLMs) for code (such as OpenAI's Codex and AI21's Jurassic J-1) for zero-shot vulnerability repair. We investigate challenges in the design of prompts that coax LLMs into generating repaired versions of insecure code. This is difficult due to the numerous ways to phrase key information - both semantically and syntactically - with natural languages. We perform a large scale study of five commercially available, black-box, "off-the-shelf" LLMs, as well as an open-source model and our own locally-trained model, on a mix of synthetic, hand-crafted, and real-world security bug scenarios. Our experiments demonstrate that while the approach has promise (the LLMs could collectively repair 100% of our synthetically generated and hand-crafted scenarios), a qualitative evaluation of the model's performance over a corpus of historical real-world examples highlights challenges in generating functionally correct code.


title: Quo Vadis, Mathematical General Relativity?

abstract: A review of selected topics in mathematical general relativity


title: Abelian monopoles of the Dirac type and color confinement in QCD

abstract: We present results of $SU(3)$ Monte-Carlo studies of a new color confinement scheme due to Abelian-like monopoles of the Dirac type without any gauge-fixing. We get (1) perfect Abelian dominance with respect to the static potentials on $12^4\sim 16^4$ lattice at $\beta=5.6-5.8$ using the multilevel method, (2) monopole as well as Abelian dominances with respect to the static potentials by evaluating Polyakov loop correlators on $24^3\times4$ lattice at $\beta=5.6$. (3) Abelian dual Meissner effects are studied directly by measuring Abelian color fields and monopole currents around the static source. The vacuum in pure $SU(3)$ seems to be of the type 1 near the border between both types, although scaling is not studied yet.


title: MIMO Networks with One-Bit ADCs: Receiver Design and Communication Strategies

abstract: High resolution analog to digital converters (ADCs) are conventionally used at the receiver terminals to store an accurate digital representation of the received signal, thereby allowing for reliable decoding of transmitted messages. However, in a wide range of applications, such as communication over millimeter wave and massive multiple-input multiple-output (MIMO) systems, the use of high resolution ADCs is not feasible due to power budget limitations. In the conventional fully digital receiver design, where each receiver antenna is connected to a distinct ADC, reducing the ADC resolution leads to performance loss in terms of achievable rates. One proposed method to mitigate the rate-loss is to use analog linear combiners leading to design of hybrid receivers. Here, the hybrid framework is augmented by the addition of delay elements to allow for temporal analog processing. Two new classes of receivers consisting of delay elements, analog linear combiners, and one-bit ADCs are proposed. The fundamental limits of communication in single and multi-user (uplink and downlink) MIMO systems employing the proposed receivers are investigated. In the high signal to noise ratio regime, it is shown that the proposed receivers achieve the maximum achievable rates among all receivers with the same number of one-bit ADCs.


title: Distributed Adaptive Learning Under Communication Constraints

abstract: This work examines adaptive distributed learning strategies designed to operate under communication constraints. We consider a network of agents that must solve an online optimization problem from continual observation of streaming data. The agents implement a distributed cooperative strategy where each agent is allowed to perform local exchange of information with its neighbors. In order to cope with communication constraints, the exchanged information must be unavoidably compressed. We propose a diffusion strategy nicknamed as ACTC (Adapt-Compress-Then-Combine), which relies on the following steps: i) an adaptation step where each agent performs an individual stochastic-gradient update with constant step-size; ii) a compression step that leverages a recently introduced class of stochastic compression operators; and iii) a combination step where each agent combines the compressed updates received from its neighbors. The distinguishing elements of this work are as follows. First, we focus on adaptive strategies, where constant (as opposed to diminishing) step-sizes are critical to respond in real time to nonstationary variations. Second, we consider the general class of directed graphs and left-stochastic combination policies, which allow us to enhance the interplay between topology and learning. Third, in contrast with related works that assume strong convexity for all individual agents' cost functions, we require strong convexity only at a network level, a condition satisfied even if a single agent has a strongly-convex cost and the remaining agents have non-convex costs. Fourth, we focus on a diffusion (as opposed to consensus) strategy. Under the demanding setting of compressed information, we establish that the ACTC iterates fluctuate around the desired optimizer, achieving remarkable savings in terms of bits exchanged between neighboring agents.

categories: cs.LG cs.IT cs.MA math.IT math.OC stat.ML


title: NN Based Active Disturbance Rejection Controller for a Multi-Axis Gimbal System

abstract: The increasing demand for target tracking, environmental surveys, surveillance and mapping requires multi-axis gimbal systems with high tracking and stabilization performance. In this paper, first, computed torque model is generated to estimate the complex disturbances acting on the system. Then, two different control strategies based on active disturbance rejection control (ADRC) and computed torque model are implemented on a two-axis gimbal system. The purpose is to improve the robustness, environmental adaptability and tracking accuracy of the system and reduce the tuning effort of ADRC by integrating a neural network (NN) based disturbance compensator (NN assisted ADRC). In the second control strategy, NN is replaced with a computed torque model (CTM assisted ADRC), whose inputs come from plant outputs. The simulation results show that, NN and CTM assisted ADRC structures can decrease mean tracking errors up to 85.4% and 40.8%, respectively.


title: In-silico synthesis of lowest-pressure high-$T_\text{c}$ ternary superhydrides

abstract: We report the theoretical prediction of two high-performing hydride superconductors BaSiH$_8$ and SrSiH$_8$. They are thermodynamically stable above pressures of $130$ and $174\,\text{GPa}$, respectively, and metastable below that. Employing anharmonic phonon calculations, we determine the minimum pressures of dynamical stability to be around $3\,\text{GPa}$ for BaSiH$_8$ and $27\,\text{GPa}$ for SrSiH$_8$, and using the fully anisotropic Migdal-Eliashberg theory, we predict $T_\text{c}$'s around $71$ and $126\,\text{K}$, respectively. We also introduce a method to estimate the lowest pressure of synthesis, based on the calculation of the enthalpy barriers protecting the BaSiH$_8$ $Fm\bar{3}m$ structure from decomposition at various pressures. This $kinetic$ pressure threshold is sensibly higher than the one based on $dynamic$ stability, but gives a much more rigorous limit for synthesizability.


title: Optimal Vehicle Path Planning Using Quadratic Optimization for Baidu Apollo Open Platform

abstract: Path planning is a key component in motion planning for autonomous vehicles. A path specifies the geometrical shape that the vehicle will travel, thus, it is critical to safe and comfortable vehicle motions. For urban driving scenarios, autonomous vehicles need the ability to navigate in cluttered environment, e.g., roads partially blocked by a number of vehicles/obstacles on the sides. How to generate a kinematically feasible and smooth path, that can avoid collision in complex environment, makes path planning a challenging problem. In this paper, we present a novel quadratic programming approach that generates optimal paths with resolution-complete collision avoidance capability.


title: Optimal Trajectory Generation for Autonomous Vehicles Under Centripetal Acceleration Constraints for In-lane Driving Scenarios

abstract: This paper presents a noval method that generates optimal trajectories for autonomous vehicles for in-lane driving scenarios. The method computes a trajectory using a two-phase optimization procedure. In the first phase, the optimization procedure generates a close-form driving guide line with differetiable curvatures. In the second phase, the procedure takes the driving guide line as input, and outputs dynamically feasible, jerk and time optimal trajectories for vehicles driving along the guide line. This method is especially useful for generating trajectories at curvy road where the vehicles need to apply frequent accelerations and decelerations to accommodate centripetal acceleration limits.


title: New approach to lattice QCD at finite density: reweighting without an overlap problem

abstract: Approaches to finite baryon density lattice QCD usually suffer from uncontrolled systematic uncertainties in addition to the well-known sign problem. We test a method - sign reweighting - that works directly at finite chemical potential and is yet free from any such uncontrolled systematics: with this approach the only problem is the sign problem itself. In practice the approach involves the generation of configurations with the positive fermionic weights given by the absolute value of the real part of the quark determinant, and a reweighting by a sign. There are only two sectors, +1 and -1 and as long as the average $\left\langle \pm \right\rangle \neq 0$ (with respect to the positive weight) this discrete reweighting has no overlap problem - unlike reweighting from $\mu=0$ - and the results are reliable. We also present results based on this algorithm on the phase diagram of lattice QCD with two different actions: as a first test, we apply the method to calculate the position of the critical endpoint with unimproved staggered fermions at $N_\tau=4$; as a second application, we study the phase diagram with 2stout improved staggered fermions at $N_\tau=6$. This second one is already a reasonably fine lattice - relevant for phenomenology. We demonstrate that the method penetrates the region of the phase diagram where the Taylor and imaginary chemical potential methods lose predictive power.


title: H{\alpha} Dots: Direct-Method Metal Abundances of Low-Luminosity Star-Forming Systems

abstract: Utilizing low-luminosity star-forming systems discovered in the H$\alpha$ Dots survey, we present spectroscopic observations undertaken using the KPNO 4m telescope for twenty-six sources. With determinations of robust, "direct"-method metal abundances, we examine the properties of these dwarf systems, exploring their utility in characterizing starburst galaxies at low luminosities and stellar masses. We find that the H$\alpha$ Dots survey provides an effective new avenue for identifying star-forming galaxies in these regimes. In addition, we examine abundance characteristics and metallicity scaling relations with these sources, highlighting a flattening of both the luminosity-metallicity ($L$-$Z$) and stellar mass-metallicity ($M_{*}$-$Z$) relation slopes in these regimes as compared with those utilizing samples covering wider respective dynamic ranges. These local, accessible analogues to the kinds of star-forming dwarfs common at high redshift will help shed light on the building blocks which assembled into the massive galaxies common today.


title: Physics-based Model for Nonuniform Thermionic Electron Emission from Polycrystalline Cathodes

abstract: Experimental observations of thermionic electron emission demonstrate a smooth transition between TL and FSCL regions of the emitted-current-density-versus-temperature (J-T) (Miram) curve and the emitted-current-density-versus-voltage (J-V) curve. Knowledge of the temperature and shape of the TL-FSCL transition is important in evaluating the thermionic electron emission performance of cathodes, including predicting the lifetime. However, there have been no first-principles physics-based models that can predict the smooth TL-FSCL transition region for real thermionic cathodes without applying physically difficult to justify a priori assumptions or empirical phenomenological equations. Previous work detailing the nonuniform thermionic emission found that the effects of 3-D space charge, patch fields, and Schottky barrier lowering can lead to a smooth TL-FSCL transition region from a model thermionic cathode surface with a checkerboard spatial distribution of work function values. In this work, we construct a physics-based nonuniform emission model for commercial dispenser cathodes for the first time. This emission model is obtained by incorporating the cathode surface grain orientation via electron backscatter diffraction (EBSD) and the facet-orientation-specific work function values from density functional theory (DFT) calculations. The model enables construction of two-dimensional emitted current density maps of the cathode surface and corresponding J-T and J-V curves. The predicted emission curves show excellent agreement with experiment, not only in TL and FSCL regions but, crucially, also in the TL-FSCL transition region. This model improves the understanding on the relationship between thermionic emission and cathode microstructure, which is beneficial to the design of vacuum electronic devices.


title: Loops in 4+1d Topological Phases

abstract: 2+1d topological phases are well characterized by the fusion rules and braiding/exchange statistics of fractional point excitations. In 4+1d, some topological phases contain only fractional loop excitations. What kind of loop statistics exist? We study the 4+1d gauge theory with 2-form $\mathbb{Z}_2$ gauge field (the loop only toric code) and find that while braiding statistics between two different types of loops can be nontrivial, the self `exchange' statistics are all trivial. In particular, we show that the electric, magnetic, and dyonic loop excitations in the 4+1d toric code are not distinguished by their self-statistics. They tunnel into each other across 3+1d invertible domain walls which in turn give explicit unitary circuits that map the loop excitations into each other. The $SL(2,\mathbb{Z}_2)$ symmetry that permutes the loops, however, cannot be consistently gauged and we discuss the associated obstruction in the process. Moreover, we discuss a gapless boundary condition dubbed the `fractional Maxwell theory' and show how it can be Higgsed into gapped boundary conditions. We also discuss the generalization of these results from the $\mathbb{Z}_2$ gauge group to $\mathbb{Z}_N$.


title: Cavity-Modulated Proton Transfer Reactions

abstract: Proton transfer is ubiquitous in many fundamental chemical and biological processes, and the ability to modulate and control the proton transfer rate would have a major impact on numerous quantum technological advances. One possibility to modulate the reaction rate of proton transfer processes is given by exploiting the strong light-matter coupling of chemical systems inside optical or nanoplasmonic cavities. In this work, we investigate the proton transfer reactions in the prototype malonaldehyde and Z-3-amino-propenal (aminopropenal) molecules using different quantum electrodynamics methods, in particular quantum electrodynamics coupled cluster theory (QED-CC) and quantum electrodynamical density functional theory (QEDFT). Depending on the cavity mode polarization direction, we show that the optical cavity can increase the reaction energy barrier by 10--20$\%$ or decrease the reaction barrier by $\sim$5$\%$. By using first principles methods, this work establishes strong light-matter coupling as a viable and practical route to alter and catalyze proton transfer reactions.


title: Face Reconstruction with Variational Autoencoder and Face Masks

abstract: Variational AutoEncoders (VAE) employ deep learning models to learn a continuous latent z-space that is subjacent to a high-dimensional observed dataset. With that, many tasks are made possible, including face reconstruction and face synthesis. In this work, we investigated how face masks can help the training of VAEs for face reconstruction, by restricting the learning to the pixels selected by the face mask. An evaluation of the proposal using the celebA dataset shows that the reconstructed images are enhanced with the face masks, especially when SSIM loss is used either with l1 or l2 loss functions. We noticed that the inclusion of a decoder for face mask prediction in the architecture affected the performance for l1 or l2 loss functions, while this was not the case for the SSIM loss. Besides, SSIM perceptual loss yielded the crispest samples between all hypotheses tested, although it shifts the original color of the image, making the usage of the l1 or l2 losses together with SSIM helpful to solve this issue.


title: Combining Embeddings and Fuzzy Time Series for High-Dimensional Time Series Forecasting in Internet of Energy Applications

abstract: The prediction of residential power usage is essential in assisting a smart grid to manage and preserve energy to ensure efficient use. An accurate energy forecasting at the customer level will reflect directly into efficiency improvements across the power grid system, however forecasting building energy use is a complex task due to many influencing factors, such as meteorological and occupancy patterns. In addiction, high-dimensional time series increasingly arise in the Internet of Energy (IoE), given the emergence of multi-sensor environments and the two way communication between energy consumers and the smart grid. Therefore, methods that are capable of computing high-dimensional time series are of great value in smart building and IoE applications. Fuzzy Time Series (FTS) models stand out as data-driven non-parametric models of easy implementation and high accuracy. Unfortunately, the existing FTS models can be unfeasible if all features were used to train the model. We present a new methodology for handling high-dimensional time series, by projecting the original high-dimensional data into a low dimensional embedding space and using multivariate FTS approach in this low dimensional representation. Combining these techniques enables a better representation of the complex content of multivariate time series and more accurate forecasts.


title: Tuning the coherent propagation of organic exciton-polaritons through dark state delocalization

abstract: While there have been numerous reports of long-range polariton transport at room-temperature in organic cavities, the spatio-temporal evolution of the propagation is scarcely reported, particularly in the initial coherent sub-ps regime, where photon and exciton wavefunctions are inextricably mixed. Hence the detailed process of coherent organic exciton-polariton transport and in particular the role of dark states has remained poorly understood. Here, we use femtosecond transient absorption microscopy to directly image coherent polariton motion in microcavities of varying quality factor. We find the transport to be well-described by a model of band-like propagation of an initially Gaussian distribution of exciton-polaritons in real space. The velocity of the polaritons reaches values of ~0.65x10^6 m s-1, substantially lower than expected from the polariton dispersion. Further, we find that the velocity is proportional to the quality factor of the microcavity. We suggest this unexpected link between the quality-factor and polariton velocity and slow coherent transport to be a result of varying admixing between delocalised dark and polariton states.


title: An Impossible Asylum

abstract: In 1982, Raymond Smullyan published an article, "The Asylum of Doctor Tarr and Professor Fether," that consists of a series of puzzles. These were later reprinted in the anthology, "The Lady or The Tiger? and Other Logic Puzzles." The last puzzle, which describes the asylum alluded to in the title, was designed to be especially difficult. With the help of automated reasoning, we show that the puzzle's hypotheses are, in fact, inconsistent, which is to say, no such asylum can possibly exist.


title: CTIN: Robust Contextual Transformer Network for Inertial Navigation

abstract: Recently, data-driven inertial navigation approaches have demonstrated their capability of using well-trained neural networks to obtain accurate position estimates from inertial measurement units (IMU) measurements. In this paper, we propose a novel robust Contextual Transformer-based network for Inertial Navigation~(CTIN) to accurately predict velocity and trajectory. To this end, we first design a ResNet-based encoder enhanced by local and global multi-head self-attention to capture spatial contextual information from IMU measurements. Then we fuse these spatial representations with temporal knowledge by leveraging multi-head attention in the Transformer decoder. Finally, multi-task learning with uncertainty reduction is leveraged to improve learning efficiency and prediction accuracy of velocity and trajectory. Through extensive experiments over a wide range of inertial datasets~(e.g. RIDI, OxIOD, RoNIN, IDOL, and our own), CTIN is very robust and outperforms state-of-the-art models.


title: PyBryt: auto-assessment and auto-grading for computational thinking

abstract: We continuously interact with computerized systems to achieve goals and perform tasks in our personal and professional lives. Therefore, the ability to program such systems is a skill needed by everyone. Consequently, computational thinking skills are essential for everyone, which creates a challenge for the educational system to teach these skills at scale and allow students to practice these skills. To address this challenge, we present a novel approach to providing formative feedback to students on programming assignments. Our approach uses dynamic evaluation to trace intermediate results generated by student's code and compares them to the reference implementation provided by their teachers. We have implemented this method as a Python library and demonstrate its use to give students relevant feedback on their work while allowing teachers to challenge their students' computational thinking skills.


title: Exploratory Data Analysis of Urdu Poetry

abstract: The study presented here provides numerical insight into ghazal -- the most appreciated genre in Urdu poetry. Using 48,761 poetic works from 4,754 poets produced over a period of 800 years, this study explores the main features of Urdu ghazal that make it popular and admired more than other forms. A detailed explanation is provided as to the types of words used for expressing love, nature, birds, and flowers etc. Also considered is the way in which the poets addressed their loved ones in their poetry. The style of poetry is numerically analyzed using Multi Dimensional Scaling to reveal the lexical diversity and similarities/differences between the different poetic works that have drawn the attention of critics, such as Iqbal and Ghalib, Mir Taqi Mir and Mir Dard. The analysis produced here is particularly helpful for research in computational stylistics, neurocognitive poetics, and sentiment analysis.


title: Skyrmion-skyrmion interaction induced by itinerant electrons in a ferromagnetic strip

abstract: Magnetic skyrmions are promising spin textures for building next-generation magnetic memories and spintronic devices. Nevertheless, one of the major challenges in realizing skyrmion-based devices is the stabilization of ordered arrays of these spin textures in different geometries. Here we numerically study the skyrmion-skyrmion interaction potential that arises due to the dynamics of itinerant electrons coupled to the magnetic texture in a ferromagnetic background with racetrack geometry. We consider different topological textures (ferromagnetic (FM) and antiferromagnetic (AFM)), namely: skyrmions, antiskyrmions and biskyrmions. We show that at low electron filling, for sufficiently short separation, the skyrmions strongly couple each other yielding a bound-state bound by electronic dynamics. However, when the filling is increased, the interaction potential energy presents local minima at specific values of the skyrmion-skyrmion distance. Each of these local minima correspond to energetically stable positions of skyrmions which are "protected" by well defined energy barriers. By inspecting the local charge density, we find that in the case of AFM skyrmions, the local antiferromagnetic nature prevents electronic penetration into the core, allowing the AFM skyrmions to be seen as infinite potential barriers for electrons.


title: Hall-Littlewood polynomials, boundaries, and $p$-adic random matrices

abstract: We prove that the boundary of the Hall-Littlewood $t$-deformation of the Gelfand-Tsetlin graph is parametrized by infinite integer signatures, extending results of Gorin and Cuenca on boundaries of related deformed Gelfand-Tsetlin graphs. In the special case when $1/t$ is a prime $p$ we use this to recover results of Bufetov-Qiu and Assiotis on infinite $p$-adic random matrices, placing them in the general context of branching graphs derived from symmetric functions. Our methods rely on explicit formulas for certain skew Hall-Littlewood polynomials. As a separate corollary to these, we obtain a simple expression for the joint distribution of the cokernels of products $A_1, A_2A_1, A_3A_2A_1,\ldots$ of independent Haar-distributed matrices $A_i$ over the $p$-adic integers $\mathbb{Z}_p$. This expression generalizes the explicit formula for the classical Cohen-Lenstra measure on abelian $p$-groups.


title: On Rees algebras of ideals and modules over hypersurface rings

abstract: The acquisition of the defining equations of Rees algebras is a natural way to study these algebras and allows certain invariants and properties to be deduced. In this paper, we consider Rees algebras of codimension 2 perfect ideals of hypersurface rings and produce a minimal generating set for their defining ideals. Then, using generic Bourbaki ideals, we study Rees algebras of modules with projective dimension one over hypersurface rings. We describe the defining ideal of such algebras and determine Cohen-Macaulayness and other invariants.


title: Strategies and Resources for teaching your physics course online on short notice

abstract: This article summarizes strategies and resources for teaching physics courses online on short notice.


title: Connecting Chromospheric Condensation Signatures to Reconnection Driven Heating Rates in an Observed Flare

abstract: Observations of solar flare reconnection at very high spatial and temporal resolution can be made indirectly at the footpoints of reconnected loops into which flare energy is deposited. The response of the lower atmosphere to this energy input includes a downward-propagating shock called chromospheric condensation, which can be observed in the UV and visible. In order to characterize reconnection using high-resolution observations of this response, one must develop a quantitative relationship between the two. Such a relation was recently developed and here we test it on observations of chromospheric condensation in a single footpoint from a flare ribbon of the X1.0 flare on 25 Oct. 2014 (SOL2014-10-25T16:56:36). Measurements taken of Si iv 1402.77 {\AA} emission spectra using the Interface Region Imaging Spectrograph (IRIS) in a single pixel show red-shifted component undergoing characteristic condensation evolution. We apply the technique called the Ultraviolet Footpoint Calorimeter (UFC) to infer energy deposition into the one footpoint. This energy profile, persisting much longer than the observed condensation, is input into a one-dimensional, hydrodynamic simulation to compute the chromospheric response, which contains a very brief condensation episode. From this simulation we synthesize Si iv spectra and compute the time-evolving Doppler velocity. The synthetic velocity evolution is found to compare reasonably well with the IRIS observation, thus corroborating our reconnection-condensation relationship. The exercise reveals that the chromospheric condensation characterizes a particular portion of the reconnection energy release rather than its entirety, and that the time scale of condensation does not necessarily reflect the time scale of energy input.


title: Symbolic dynamics of planar piecewise smooth vector fields

abstract: Recently, the theory concerning piecewise smooth vector fields (PSVFs for short) have been undergoing important improvements. In fact, many results obtained do not have an analogous for smooth vector fields. For example, the chaoticity of planar PSVFs, which is impossible for the smooth ones. These differences are generated by the non-uniqueness of trajectory passing through a point. Inspired by the classical fact that one-dimensional discrete dynamic systems can produce chaotic behavior, we construct a conjugation between the shift map and PSVFs. By means of the results obtained and the techniques employed, a new perspective on the study of PSVFs is brought to light and, through already established results for discrete dynamic systems, we will be able to obtain results regarding PSVFs.


title: A reliable Turing machine

abstract: We consider computations of a Turing machine subjected to noise. In every step, the action (the new state and the new content of the observed cell, the direction of the head movement) can differ from that prescribed by the transition function with a small probability (independently of previous such events). We construct a universal 1-tape Turing machine that for a low enough (constant) noise probability performs arbitrarily large computations. For this unavoidably, the input needs to be encoded -- by a simple code depending on its size. The work uses a technique familiar from reliable cellular automata, complemented by some new ones.


title: Standard Self-Confinement and Extrinsic Turbulence Models for Cosmic Ray Transport are Fundamentally Incompatible with Observations

abstract: Models for cosmic ray (CR) dynamics fundamentally depend on the rate of CR scattering from magnetic fluctuations. In the ISM, for CRs with energies ~MeV-TeV, these fluctuations are usually attributed either to 'extrinsic turbulence' (ET) - a cascade from larger scales - or 'self-confinement' (SC) - self-generated fluctuations from CR streaming. Using simple analytic arguments and detailed live numerical CR transport calculations in galaxy simulations, we show that both of these, in standard form, cannot explain even basic qualitative features of observed CR spectra. For ET, any spectrum that obeys critical balance or features realistic anisotropy, or any spectrum that accounts for finite damping below the dissipation scale, predicts qualitatively incorrect spectral shapes and scalings of B/C and other species. Even if somehow one ignored both anisotropy and damping, observationally-required scattering rates disagree with ET predictions by orders-of-magnitude. For SC, the dependence of driving on CR energy density means that it is nearly impossible to recover observed CR spectral shapes and scalings, and again there is an orders-of-magnitude normalization problem. But more severely, SC solutions with super-Alfvenic streaming are unstable. In live simulations, they revert to either arbitrarily-rapid CR escape with zero secondary production, or to bottleneck solutions with far-too-strong CR confinement and secondary production. Resolving these fundamental issues without discarding basic plasma processes requires invoking different drivers for scattering fluctuations. These must act on a broad range of scales with a power spectrum obeying several specific (but plausible) constraints.

categories: astro-ph.HE astro-ph.CO astro-ph.GA physics.plasm-ph physics.space-ph


title: Wireless Power Transmission on Martian Surface for Zero-Energy Devices

abstract: Exploration of the Red Planet is essential on the way through both human colonization and establishing a habitat on the planet. Due to the high costs of space missions, the use of distributed sensor networks has been investigated to make in situ explorations affordable. Along with this, the devices with ultra-low-power receivers, which are called zero-energy (ZE) devices, can pave the way to further discoveries about the environment of Mars. This study focuses on wireless power transmission to provide the power required by ZE devices on the Martian surface. The main motivation of this study is to investigate whether conventional harvesters and communication units can supply the required power for a long distance. The numerical results show that it is possible to deliver power to ZE devices without utilizing any sophisticated hardware. In addition, the effects of pointing error and dust storms on harvesting performance are investigated. Comprehensive simulation results reveal that harvester selection and design should be done by considering propagation channel and transmitter characteristics.


title: What do abelian categories form?

abstract: This is mostly an overview. Given finitely presentable abelian categories $A$ and $B$, we sketch the construction of an abelian category of continuous functors from $A$ to $B$ that has nice $2$-categorical behaviour and gives an explicit model for the stable category of stable functors between the derived categories of $A$ and $B$. The construction is absolute, thus allows one to recover not only Hochschild but also Mac Lane Cohomology.


title: New topological observables in a model of Causal Dynamical Triangulations on a torus

abstract: The structure of simplicial manifolds in a model of Causal Dynamical Triangulations in 3+1 dimensions with the spatial topology of a 3-torus is analyzed with the help of topological observables, such as loops with nonzero winding numbers and coordinates based on scalar fields with jumps at the boundaries of the elementary cell of the torus. The results are given an interpretation and used in the measurements of a more local observable that is the quantum Ricci curvature. We moreover analyze the influence of scalar matter fields on the geometry of CDT spacetimes.


title: A Web of Confocal Parabolas in a Grid of Hexagons

abstract: If one erects regular hexagons upon the sides of a triangle $T$, several surprising properties emerge, including: (i) the triangles which flank said hexagons have an isodynamic point common with $T$, (ii) the construction can be extended iteratively, forming an infinite grid of regular hexagons and flank triangles, (iii) a web of confocal parabolas with only three distinct foci interweaves the vertices of hexagons in the grid. Finally, (iv) said foci are the vertices of an equilateral triangle.


title: On Topological Entropy of Piecewise Smooth Vector Fields 2

abstract: Non-smooth vector fields does not have necessarily the property of uniqueness of solution passing through a point and this is responsible to enrich the behavior of the system. Even on the plane non-smooth vector fields can be chaotic, a feature impossible for the smooth or continuous case. We propose a new approach towards a better understanding of chaos for non-smooth vector fields and this is done by studying the entropy of the system. In this work we set the ground for one to begin the study of entropy for non-smooth vector fields. We construct a metric space of all possible trajectories of a non-smooth vector field, where we define a flow inherited by the vector field and then define the topological entropy in this scenario. As a consequence, we are able to obtain some general results of this theory and give some examples of planar non-smooth vector fields with positive (finite and infinite) entropy.


title: Optimizing Performance and Satisfaction in Matching and Movement Tasks in Virtual Reality with Interventions Using the Data Visualization Literacy Framework

abstract: Virtual reality (VR) has seen increased use for training and instruction. Designers can enable VR users to gain insights into their own performance by visualizing telemetry data from their actions in VR. Our ability to detect patterns and trends visually suggests the use of data visualization as a tool for users to identify strategies for improved performance. Typical tasks in VR training scenarios are manipulation of 3D objects (e.g., for learning how to maintain a jet engine) and navigation (e.g., to learn the geography of a building or landscape before traveling on-site). In this paper, we present the results of the RUI VR (84 subjects) and Luddy VR studies (68 subjects), where participants were divided into experiment and control cohorts. All subjects performed a series of tasks: 44 cube-matching tasks in RUI VR and 48 navigation tasks through a virtual building in Luddy VR (all divided into two sets). All Luddy VR subjects used VR gear; RUI VR subjects were divided across three setups: 2D Desktop (with laptop and mouse), VR Tabletop (in VR, sitting at a table), and VR Standup (in VR, standing). In an intervention called "Reflective phase," the experiment cohorts were presented with data visualizations, designed with the Data Visualization Literacy Framework (DVL-FW), of the data they generated during the first set of tasks before continuing to the second part of the study. For Luddy VR, we found that experiment users had significantly faster completion times in their second trial (p = 0.014) while scoring higher in a mid-questionnaire about the virtual building (p = 0.009). For RUI VR, we found no significant differences for completion time and accuracy between the two cohorts in the VR setups; however, 2D Desktop subjects in the experiment cohort had significantly higher rotation accuracy as well as satisfaction (p(rotation) = 0.031, p(satisfaction) = 0.040).


title: A Detection of Red Noise in PSR J1824$-$2452A and Projections for PSR B1937+21 using NICER X-ray Timing Data

abstract: We have used X-ray data from the Neutron Star Interior Composition Explorer (NICER) to search for long time-scale, correlated variations ("red noise") in the pulse times of arrival from the millisecond pulsars PSR J1824$-$2452A and PSR B1937+21. These data more closely track intrinsic noise because X-rays are unaffected by the radio-frequency dependent propagation effects of the interstellar medium. Our Bayesian search methodology yields strong evidence (natural log Bayes factor of $9.634 \pm 0.016$) for red noise in PSR J1824$-$2452A, but is inconclusive for PSR B1937+21. In the interest of future X-ray missions, we devise and implement a method to simulate longer and higher precision X-ray datasets to determine the timing baseline necessary to detect red noise. We find that the red noise in PSR B1937+21 can be reliably detected in a 5-year mission with a time-of-arrival (TOA) error of 2 microseconds and an observing cadence of 20 observations per month compared to the 5 microsecond TOA error and 11 observations per month that NICER currently achieves in PSR B1937+21. We investigate detecting red noise in PSR B1937+21 with other combinations of observing cadences and TOA errors. We also find that an injected stochastic gravitational wave background (GWB) with an amplitude of $A_{\rm GWB}=2\times10^{-15}$ and spectral index of $\gamma_{\rm GWB}=13/3$ can be detected in a pulsar with similar TOA precision to PSR B1937+21, but with no additional red noise, in a 10-year mission that observes the pulsar 15 times per month and has an average TOA error of 1 microsecond.


title: MAC address randomization tolerant crowd monitoring system using Wi-Fi packets

abstract: Media access control (MAC) addresses inside Wi-Fi packets can be used for beneficial activities such as crowdedness estimation, marketing, and hazard maps. However, the MAC address randomization systems introduced around 2014 make all conventional MAC-address-based crowd monitoring systems count the same device more than once. Therefore, there is a need to create a new crowd monitoring system tolerant to MAC address randomization to estimate the number of devices accurately. In this paper, Vision and TrueSight, two new crowd monitoring algorithms that estimate the number of devices, are proposed to prove that MAC-address-based crowd monitoring is still possible. In addition to probe requests, Vision uses data packets and beacon packets to mitigate the influence of randomization. Moreover, TrueSight uses sequence numbers and hierarchical clustering to estimate the number of devices. The experimental results of this study show that even without installing any special software, Vision can gather 440 randomly generated MAC addresses into one group and count only once, and TrueSight can estimate the number of devices with an accuracy of more than 75% with an acceptable error range of 1.


title: A Low-cost Robot with Autonomous Recharge and Navigation for Weed Control in Fields with Narrow Row Spacing

abstract: Modern herbicide application in agricultural settings typically relies on either large scale sprayers that dispense herbicide over crops and weeds alike or portable sprayers that require labor intensive manual operation. The former method results in overuse of herbicide and reduction in crop yield while the latter is often untenable in large scale operations. This paper presents the first fully autonomous robot for weed management for row crops capable of computer vision based navigation, weed detection, complete field coverage, and automatic recharge for under \$400. The target application is autonomous inter-row weed control in crop fields, e.g. flax and canola, where the spacing between croplines is as small as one foot. The proposed robot is small enough to pass between croplines at all stages of plant growth while detecting weeds and spraying herbicide. A recharging system incorporates newly designed robotic hardware, a ramp, a robotic charging arm, and a mobile charging station. An integrated vision algorithm is employed to assist with charger alignment effectively. Combined, they enable the robot to work continuously in the field without access to electricity. In addition, a color-based contour algorithm combined with preprocessing techniques is applied for robust navigation relying on the input from the onboard monocular camera. Incorporating such compact robots into farms could help automate weed control, even during late stages of growth, and reduce herbicide use by targeting weeds with precision. The robotic platform is field-tested in the flaxseed fields of North Dakota.


title: WebRTC-based measurement tool for peer-to-peer applications and preliminary findings with real users

abstract: Direct peer-to-peer (P2P) communication is often used to minimize the end-to-end latency for real-time applications that require accurate synchronization, such as remote musical ensembles. However, there are few studies on the performance of P2P communication between home network environments, thus hindering the deployment of services that require synchronization. In this study, we developed a P2P performance measurement tool using the Web Real-Time Communication (WebRTC) statistics application programming interface. Using this tool, we can easily measure P2P performance between home network environments on a web browser without downloading client applications. We also verified the reliability of round-trip time (RTT) measurements using WebRTC and confirmed that our system could provide the necessary measurement accuracy for RTT and jitter measurements for real-time applications. In addition, we measured the performance of a full mesh topology connection with 10 users in an actual environment in Japan. Consequently, we found that only 66% of the peer connections had a latency of 30 ms or less, which is the minimum requirement for high synchronization applications, such as musical ensembles.


title: Bridging the gap between prostate radiology and pathology through machine learning

abstract: Prostate cancer is the second deadliest cancer for American men. While Magnetic Resonance Imaging (MRI) is increasingly used to guide targeted biopsies for prostate cancer diagnosis, its utility remains limited due to high rates of false positives and false negatives as well as low inter-reader agreements. Machine learning methods to detect and localize cancer on prostate MRI can help standardize radiologist interpretations. However, existing machine learning methods vary not only in model architecture, but also in the ground truth labeling strategies used for model training. In this study, we compare different labeling strategies, namely, pathology-confirmed radiologist labels, pathologist labels on whole-mount histopathology images, and lesion-level and pixel-level digital pathologist labels (previously validated deep learning algorithm on histopathology images to predict pixel-level Gleason patterns) on whole-mount histopathology images. We analyse the effects these labels have on the performance of the trained machine learning models. Our experiments show that (1) radiologist labels and models trained with them can miss cancers, or underestimate cancer extent, (2) digital pathologist labels and models trained with them have high concordance with pathologist labels, and (3) models trained with digital pathologist labels achieve the best performance in prostate cancer detection in two different cohorts with different disease distributions, irrespective of the model architecture used. Digital pathologist labels can reduce challenges associated with human annotations, including labor, time, inter- and intra-reader variability, and can help bridge the gap between prostate radiology and pathology by enabling the training of reliable machine learning models to detect and localize prostate cancer on MRI.


title: On Submodular Contextual Bandits

abstract: We consider the problem of contextual bandits where actions are subsets of a ground set and mean rewards are modeled by an unknown monotone submodular function that belongs to a class $\mathcal{F}$. We allow time-varying matroid constraints to be placed on the feasible sets. Assuming access to an online regression oracle with regret $\mathsf{Reg}(\mathcal{F})$, our algorithm efficiently randomizes around local optima of estimated functions according to the Inverse Gap Weighting strategy. We show that cumulative regret of this procedure with time horizon $n$ scales as $O(\sqrt{n \mathsf{Reg}(\mathcal{F})})$ against a benchmark with a multiplicative factor $1/2$. On the other hand, using the techniques of (Filmus and Ward 2014), we show that an $\epsilon$-Greedy procedure with local randomization attains regret of $O(n^{2/3} \mathsf{Reg}(\mathcal{F})^{1/3})$ against a stronger $(1-e^{-1})$ benchmark.


title: Mertens' Third Theorem for Number Fields: A New Proof, Cram\'er's Inequality, Oscillations, and Bias

abstract: The first result of our article is another proof of Mertens' third theorem in the number field setting, which generalises a method of Hardy. The second result concerns the sign of the error term in Mertens' third theorem. Diamond and Pintz showed that the error term in the classical case changes sign infinitely often and in our article, we establish this result for number fields assuming a reasonable technical condition. In order to do so, we needed to prove Cram\'er's inequality for number fields, which is interesting in its own right. Lamzouri built upon Diamond and Pintz's work to prove the existence of the logarithmic density of the set of real numbers $x \ge 2$ such that the error term in Mertens' third theorem is positive, so the third result of our article generalises Lamzouri's results for number fields. We also include numerical investigations for the number fields $\mathbb{Q}(\sqrt{5})$ and $\mathbb{Q}(\sqrt{13})$, building upon similar work done by Rubinstein and Sarnak in the classical case.


title: The five largest satellites of Uranus: astrometric observations spread over 29 years at the Pico dos Dias Observatory

abstract: We present the astrometry of the five largest satellites of Uranus from observations spread over almost three decades with photographic plates and CCDs (mainly), taken at the Pico dos Dias Observatory - Brazil. All positions presented here are obtained from the reanalysis of measurements and images used in previous publications. Reference stars are those from the Gaia Early Data Release 3 (Gaia EDR3) allowing, in addition to a higher accuracy, a larger number of positions of the largest satellites as compared to our previous works. From 1982 to 1987, positions were obtained from photographic plates. From 1989 to 2011, CCDs were used. On average, we obtained $\Delta\alpha{\rm cos}\delta=-11~(\pm52)$ milli-arcseconds and $\Delta\delta=-14~(\pm43)$ milli-arcseconds for the differences in the sense observation minus ephemerides (DE435$+$ura111). Comparisons with different ephemerides (DE440, INPOP21a, INPOP19a and NOE-7-2013-MAIN) and results from stellar occultations indicate a possible offset in the (Solar System) barycentric position of the Uranian system barycenter. Overall, our results are useful to improve dynamical models of the Uranian largest satellites as well as the orbit of Uranus.


title: The Scorpion Planet Survey: Wide-Orbit Giant Planets Around Young A-type Stars

abstract: The first directly imaged exoplanets indicated that wide-orbit giant planets could be more common around A-type stars. However, the relatively small number of nearby A-stars has limited the precision of exoplanet demographics studies to $\gtrsim$10%. We aim to constrain the frequency of wide-orbit giant planets around A-stars using the VLT/SPHERE extreme adaptive optics system, which enables targeting $\gtrsim$100 A-stars between 100$-$200 pc. We present the results of a survey of 84 A-stars within the nearby $\sim$5$-$17 Myr-old Sco OB2 association. The survey detected three companions$-$one of which is a new discovery (HIP75056Ab), whereas the other two (HD 95086b and HIP65426b) are now-known planets that were included without a priori knowledge of their existence. We assessed the image sensitivity and observational biases with injection and recovery tests combined with Monte Carlo simulations to place constraints on the underlying demographics. We measure a decreasing frequency of giant planets with increasing separation, with measured values falling between 10$-$2% for separations of 30$-$100 au, and 95% confidence-level (CL) upper limits of $\lesssim$45$-$8% for planets on 30$-$100 au orbits, and $\lesssim$5% between 200$-$500 au. These values are in excellent agreement with recent surveys of A-stars in the solar neighborhood$-$supporting findings that giant planets at $\lesssim$100 au are more frequent around A-stars than around solar-type hosts. Finally, the relatively low occurrence rate of super-Jupiters on wide orbits, the positive correlation with stellar mass, and the inverse correlation with orbital separation are consistent with core accretion being their dominant formation mechanism.


title: UV-Visible Absorption Spectra of Solvated Molecules by Quantum Chemical Machine Learning

abstract: Predicting UV-visible absorption spectra is essential to understanding photochemical processes and designing energy materials. Quantum chemical methods can deliver accurate calculations of UV-visible absorption spectra, but they are computationally expensive, especially for large systems or when one computes line shapes from thermal averages. Here, we present an approach to predicting UV-visible absorption spectra of solvated aromatic molecules by quantum chemistry (QC) and machine learning (ML). We show that a ML model, trained on the high-level QC calculation of the excitation energy of a set of aromatic molecules, can accurately predict the line shape of the lowest-energy UV-visible absorption band of several related molecules with less than 0.1 eV deviation with respect to reference experimental spectra. Applying linear decomposition analysis on the excitation energies, we unveil that our ML models probe vertical excitations of these aromatic molecules primarily by learning the atomic environment of their phenyl rings, which align with the physical origin of the $\pi\rightarrow\pi^\star$ electronic transition. Our study provides an effective workflow that combines ML with quantum chemical methods to accelerate the calculations of UV-visible absorption spectra for various molecular systems.


title: Counterfactual Fairness in Mortgage Lending via Matching and Randomization

abstract: Unfairness in mortgage lending has created generational inequality among racial and ethnic groups in the US. Many studies address this problem, but most existing work focuses on correlation-based techniques. In our work, we use the framework of counterfactual fairness to train fair machine learning models. We propose a new causal graph for the variables available in the Home Mortgage Disclosure Act (HMDA) data. We use a matching-based approach instead of the latent variable modeling approach, because the former approach does not rely on any modeling assumptions. Furthermore, matching provides us with counterfactual pairs in which the race variable is isolated. We first demonstrate the unfairness in mortgage approval and interest rates between African-American and non-Hispanic White sub-populations. Then, we show that having balanced data using matching does not guarantee perfect counterfactual fairness of the machine learning models.


title: The Importance of Electron Landau Damping for the Dissipation of Turbulent Energy in Terrestrial Magnetosheath Plasma

abstract: Heliospheric plasma turbulence plays a key role in transferring the energy of large-scale magnetic field and plasma flow fluctuations to smaller scales where the energy can be dissipated, ultimately leading to plasma heating. High-quality measurements of electromagnetic fields and electron velocity distributions by the Magnetospheric Multiscale (MMS) mission in Earth's magnetosheath present a unique opportunity to characterize plasma turbulence and to determine the mechanisms responsible for its dissipation. We apply the field-particle correlation technique to a set of twenty MMS magnetosheath intervals to identify the dissipation mechanism and quantify the dissipation rate. It is found that 95% of the intervals have velocity-space signatures of electron Landau damping that are quantitatively consistent with linear kinetic theory for the collisionless damping of kinetic Alfv\'en waves. About 75% of the intervals contain asymmetric signatures, indicating a local imbalance of kinetic Alfv\'en wave energy flux in one direction along the magnetic field than the other. About one third of the intervals have an electron energization rate with the same order-of-magnitude as the estimated turbulent cascade rate, suggesting that electron Landau damping plays a significant, and sometimes dominant, role in the dissipation of the turbulent energy in these magnetosheath intervals.


title: Quantum Control of Two Critically Dressed Spin 1/2 Species in Magnetic Fluctuations

abstract: The neutron electric dipole moment experiment at the Spallation Neutron Source (nEDM@SNS experiment) proposes to measure the nEDM using the spin-dependent capture cross section of neutrons on $^3$He. The critical dressing mode of this experiment uses an oscillating magnetic field to dress the gyromagnetic ratios of neutrons and $^3$He to the same value. While this technique grants increased sensitivity to the nEDM by improving the signal-to-noise ratio, this mode of measurement also introduces additional noise from the power supply used to drive the dressing field. This can lead to randomly fluctuating magnetic fields which cause the spins of neutrons and $^3$He to drift apart over time. Here we use second-order time-dependent perturbation theory to compute relaxation and frequency shifts due to fluctuations in the dressing field in terms of the magnetic field noise power spectrum and compare these calculations to numerical solutions obtained by integrating the Bloch equations. We then use these results to develop mitigation strategies for this type of noise. Furthermore, we report on spin dressing modulation techniques that significantly amplify coherence times for the critically dressed system, and attempt to quantify the coherence time achievable.


title: The stability of 3D skyrmions under mechanical stress studied via Monte Carlo calculations

abstract: Using Monte Carlo (MC) simulations, we study the skyrmion stability/instability as a response to uniaxial mechanical stresses. Skyrmions emerge in chiral magnetic materials as a stable spin configuration under external magnetic field $\vec{B}$ with the competition of ferromagnetic interaction and Dzyaloshinskii-Moriya interaction (DMI) at low temperature $T$. Skyrmion configurations are also known to be stable (unstable) under a compressive stress applied parallel (perpendicular) to $\vec{B}$. To understand the origin of such experimentally confirmed stability/instability, we use the Finsler geometry modeling technique with a new degree of freedom for strains, which plays an essential role in DMI being anisotropic. We find from MC data that the area of the skyrmion state on the $B$-$T$ phase diagram increases (decreases) depending on the direction of applied stresses, in agreement with reported experimental results. This change in the area of the skyrmion state indicates that skyrmions become more (less) stable if the tensile strain direction is parallel (perpendicular) to $\vec{B}$. From the numerical data in this paper, we find that the so-called magneto-elastic effect is suitably implemented in the effective DMI theory with the strain degree of freedom without complex magneto-elastic coupling terms for chiral magnetic materials. This result confirms that experimentally-observed skyrmion stability and instability are caused by DMI anisotropy.


title: Unifying the geometric decompositions of full and trimmed polynomial spaces in finite element exterior calculus

abstract: Arnold, Falk, & Winther, in "Finite element exterior calculus, homological techniques, and applications" (2006), show how to geometrically decompose the full and trimmed polynomial spaces on simplicial elements into direct sums of trace-free subspaces and in "Geometric decompositions and local bases for finite element differential forms" (2009) the same authors give direct constructions of extension operators for the same spaces. The two families -- full and trimmed -- are treated separately, using differently defined isomorphisms between each and the other's trace-free subspaces and mutually incompatible extension operators. This work describes a single operator $\mathring{\star}_T$ that unifies the two isomorphisms and also defines a weighted-$L^2$ norm appropriate for defining well-conditioned basis functions and dual-basis functionals for geometric decomposition. This work also describes a single extension operator $\dot{E}_{\sigma,T}$ that implements geometric decompositions of all differential forms as well as for the full and trimmed polynomial spaces separately.


title: The h-index

abstract: The h-index is a mainstream bibliometric indicator, since it is widely used in academia, research management and research policy. While its advantages have been highlighted, such as its simple calculation, it has also received widespread criticism. The criticism is mainly based on the negative effects it may have on scholars, when the index is used to describe the quality of a scholar. The "h" means "highly-cited" and "high achievement", and should not be confused with the last name of its inventor, Hirsch. Put simply, the h-index combines a measure of quantity and impact in a single indicator. Several initiatives try to provide alternatives to the h-index to counter some of its shortcomings.


title: The TESS Faint Star Search: 1,617 TOIs from the TESS Primary Mission

abstract: We present the detection of 1,617 new transiting planet candidates, identified in the Transiting Exoplanet Survey Satellite (TESS) full-frame images (FFIs) observed during the Primary Mission (Sectors 1 - 26). These candidates were initially detected by the Quick-Look Pipeline (QLP), which extracts FFI lightcurves for and searches all stars brighter than TESS magnitude T = 13.5 mag in each sector. However, QLP heavily relies on manual inspection for the identification of planet candidates, limiting vetting efforts to planet-hosting stars brighter than T = 10.5 mag and leaving millions of potential transit signals un-vetted. We describe an independent vetting pipeline applied to QLP transit search results, incorporating both automated vetting tests and manual inspection to identify promising planet candidates around these fainter stars. The new candidates discovered by this ongoing project will allow TESS to significantly improve the statistical power of demographics studies of giant, close-in exoplanets.


title: On-line Policy Iteration with Policy Switching for Markov Decision Processes

abstract: Motivated from Bertsekas' recent study on policy iteration (PI) for solving the problems of infinite-horizon discounted Markov decision processes (MDPs) in an on-line setting, we develop an off-line PI integrated with a multi-policy improvement method of policy switching and then adapt its asynchronous variant into on-line PI algorithm that generates a sequence of policies over time. The current policy is updated into the next policy by switching the action only at the current state while ensuring the monotonicity of the value functions of the policies in the sequence. Depending on MDP's state-transition structure, the sequence converges in a finite time to an optimal policy for an associated local MDP. When MDP is communicating, the sequence converges to an optimal policy for the original MDP.


title: On some hyperelliptic Hurwitz-Hodge integrals

abstract: This short note addresses Hodge integrals over the hyperelliptic locus. Recently Afandi computed, via localisation techniques, such one-descendant integrals and showed that they are Stirling numbers. We give another proof of the same statement by a very short argument, exploiting Chern classes of spin structures and relations arising from Topological Recursion in the sense of Eynard and Orantin. These techniques seem also suitable to deal with three orthogonal generalisations: 1. the extension to the r-hyperelliptic locus, 2. the extension to an arbitrary number of non-Weierstrass pairs of points, 3. the extension to multiple descendants.


title: Projective Clustering Product Quantization

abstract: This paper suggests the use of projective clustering based product quantization for improving nearest neighbor and max-inner-product vector search (MIPS) algorithms. We provide anisotropic and quantized variants of projective clustering which outperform previous clustering methods used for this problem such as ScaNN. We show that even with comparable running time complexity, in terms of lookup-multiply-adds, projective clustering produces more quantization centers resulting in more accurate dot-product estimates. We provide thorough experimentation to support our claims.


title: Generalized Transitional Markov Chain Monte Carlo Sampling Technique for Bayesian Inversion

abstract: In the context of Bayesian inversion for scientific and engineering modeling, Markov chain Monte Carlo sampling strategies are the benchmark due to their flexibility and robustness in dealing with arbitrary posterior probability density functions (PDFs). However, these algorithms been shown to be inefficient when sampling from posterior distributions that are high-dimensional or exhibit multi-modality and/or strong parameter correlations. In such contexts, the sequential Monte Carlo technique of transitional Markov chain Monte Carlo (TMCMC) provides a more efficient alternative. Despite the recent applicability for Bayesian updating and model selection across a variety of disciplines, TMCMC may require a prohibitive number of tempering stages when the prior PDF is significantly different from the target posterior. Furthermore, the need to start with an initial set of samples from the prior distribution may present a challenge when dealing with implicit priors, e.g. based on feasible regions. Finally, TMCMC can not be used for inverse problems with improper prior PDFs that represent lack of prior knowledge on all or a subset of parameters. In this investigation, a generalization of TMCMC that alleviates such challenges and limitations is proposed, resulting in a tempering sampling strategy of enhanced robustness and computational efficiency. Convergence analysis of the proposed sequential Monte Carlo algorithm is presented, proving that the distance between the intermediate distributions and the target posterior distribution monotonically decreases as the algorithm proceeds. The enhanced efficiency associated with the proposed generalization is highlighted through a series of test inverse problems and an engineering application in the oil and gas industry.


title: Projections onto hyperbolas or bilinear constraint sets in Hilbert spaces

abstract: Sets of bilinear constraints are important in various machine learning models. Mathematically, they are hyperbolas in a product space. In this paper, we give a complete formula for projections onto sets of bilinear constraints or hyperbolas in a general Hilbert space.


title: High return level estimates of daily ERA-5 precipitation in Europe estimated using regionalised extreme value distributions

abstract: Accurate estimation of daily rainfall return levels associated with large return periods is needed for a number of hydrological planning purposes, including protective infrastructure, dams, and retention basins. This is especially relevant at small spatial scales. The ERA-5 reanalysis product provides seasonal daily precipitation over Europe on a 0.25 x 0.25 grid (about 27 x 27 km). This translates more than 20,000 land grid points and leads to models with a large number of parameters when estimating return levels. To bypass this abundance of parameters, we build on the regional frequency analysis (RFA), a well-known strategy in statistical hydrology. This approach consists in identifying homogeneous regions, by gathering locations with similar distributions of extremes up to a normalizing factor and developing sparse regional models. In particular, we propose a step-by-step blueprint that leverages a recently developed and fast clustering algorithm to infer return level estimates over large spatial domains. This enables us to produce maps of return level estimates of ERA-5 reanalysis daily precipitation over continental Europe for various return periods and seasons. We discuss limitations and practical challenges and also provide a git hub repository. We show that a relatively parsimonious model with only a spatially varying scale parameter can compete well against statistical models of higher complexity.


title: International Conferences of Bibliometrics

abstract: Conferences are deeply connected to research fields, in this case bibliometrics. As such, they are a venue to present and discuss current and innovative research, and play an important role for the scholarly community. In this article, we provide an overview on the history of conferences in bibliometrics. We conduct an analysis to list the most prominent conferences that were announced in the newsletter by ISSI, the International Society for Scientometrics and Informetrics. Furthermore, we describe how conferences are connected to learned societies and journals. Finally, we provide an outlook on how conferences might change in future.


title: V2X Misbehavior and Collective Perception Service: Considerations for Standardization

abstract: Connected and Automated Vehicles use sensors and wireless communication to improve road safety and efficiency. However, attackers may target Vehicle-to-Everything communication. Indeed, an attacker may send authenticated but wrong data to send false location information, alert incorrect events or report a bogus object endangering the safety of other CAVs. Currently, Standardization Development Organizations are working on developing security standards against such attacks. Unfortunately, current standardization efforts do not include misbehavior specifications for advanced V2X services such as Collective Perception yet. This work assesses the security of Collective Perception Messages and proposes inputs for consideration in existing standards.


title: Neural Pseudo-Label Optimism for the Bank Loan Problem

abstract: We study a class of classification problems best exemplified by the \emph{bank loan} problem, where a lender decides whether or not to issue a loan. The lender only observes whether a customer will repay a loan if the loan is issued to begin with, and thus modeled decisions affect what data is available to the lender for future decisions. As a result, it is possible for the lender's algorithm to ``get stuck'' with a self-fulfilling model. This model never corrects its false negatives, since it never sees the true label for rejected data, thus accumulating infinite regret. In the case of linear models, this issue can be addressed by adding optimism directly into the model predictions. However, there are few methods that extend to the function approximation case using Deep Neural Networks. We present Pseudo-Label Optimism (PLOT), a conceptually and computationally simple method for this setting applicable to DNNs. \PLOT{} adds an optimistic label to the subset of decision points the current model is deciding on, trains the model on all data so far (including these points along with their optimistic labels), and finally uses the resulting \emph{optimistic} model for decision making. \PLOT{} achieves competitive performance on a set of three challenging benchmark problems, requiring minimal hyperparameter tuning. We also show that \PLOT{} satisfies a logarithmic regret guarantee, under a Lipschitz and logistic mean label model, and under a separability condition on the data.


title: Solar System Dynamics and Multiyear Droughts of the Western USA

abstract: The recent addition of orbit-spin coupling torques to atmospheric global circulation models has enabled successful years-in-advance forecasts of global and regional-scale dust storms on Mars. Here we explore the applicability of the orbit-spin coupling mechanism for understanding and forecasting anomalous weather and climate events on Earth. We calculate the time history of orbit-spin coupling torques on the Earth system for the interval from 1860-2040. The torque exhibits substantial variability on decadal to bidecadal timescales. Deep minima recur at intervals from 15-26 years; eight such episodes are documented within the study period prior to 2020. Each of the identified torque minima corresponds in time to an episode of widespread drought in the Western USA extending over several years. The multiyear droughts of the 1930s, the 1950s, the mid-1970s, the early 1990s, and of 2011-2015 were each coincident in time with orbit-spin coupling torque minima. The upcoming torque minimum of 2030 is the deepest such minimum of the 180-yr study interval. A multiyear episode of widespread drought in the Western USA is likely to be underway by 2028 plus or minus 4 years (2 standard deviations). The potential benefits to societies of improved drought predictions justify an immediate high-priority effort to include forcing by orbit-spin coupling within state-of-the-art Earth system GCMs. Future targeted numerical modeling investigations are likely to yield forecasts with considerably lower uncertainties and with much improved temporal resolution in comparison to that obtained here.


title: Stripes and spin-density waves in the doped two-dimensional Hubbard model: ground state phase diagram

abstract: We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an anti-ferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods, established the stripe order as the ground state near $1/8$ doping at strong interactions. In this work, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about $1/5$ doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes in the pure Hubbard model.


title: Notes on peculiarities of quantum fields in space-times with horizons

abstract: We consider massive scalar field theory on static four-dimensional space-times with horizons. We study the near horizon behaviour of the quantum expectation values of the stress-energy tensor operator for thermal state with generic temperatures. It turns out that the dependence of the expectation values on the temperature and tensor structure of the stress-energy tensor differs from the usual one in the Minkowski space-time. Moreover, for non-canonical temperatures these expectation values are divergent on the horizons. We also show that the Wightman functions have additional infrared peculiarities near the horizons.


title: Homogeneous codimension-one foliations on reducible symmetric spaces of noncompact type

abstract: We extend the classification of homogeneous codimension-one foliations on irreducible Riemannian symmetric spaces of noncompact type obtained by Berndt and Tamaru to the reducible case, thus completing it for all noncompact symmetric spaces.


title: Markov Chain Monte-Carlo Enhanced Variational Quantum Algorithms

abstract: Variational quantum algorithms are poised to have significant impact on high-dimensional optimization, with applications in classical combinatorics, quantum chemistry, and condensed matter. Nevertheless, the optimization landscape of these algorithms is generally nonconvex, causing suboptimal solutions due to convergence to local, rather than global, minima. In this work, we introduce a variational quantum algorithm that uses classical Markov chain Monte Carlo techniques to provably converge to global minima. These performance gaurantees are derived from the ergodicity of our algorithm's state space and enable us to place analytic bounds on its time-complexity. We demonstrate both the effectiveness of our technique and the validity of our analysis through quantum circuit simulations for MaxCut instances, solving these problems deterministically and with perfect accuracy. Our technique stands to broadly enrich the field of variational quantum algorithms, improving and gauranteeing the performance of these promising, yet often heuristic, methods.


title: NN-LUT: Neural Approximation of Non-Linear Operations for Efficient Transformer Inference

abstract: Non-linear operations such as GELU, Layer normalization, and Softmax are essential yet costly building blocks of Transformer models. Several prior works simplified these operations with look-up tables or integer computations, but such approximations suffer inferior accuracy or considerable hardware cost with long latency. This paper proposes an accurate and hardware-friendly approximation framework for efficient Transformer inference. Our framework employs a simple neural network as a universal approximator with its structure equivalently transformed into a LUT. The proposed framework called NN-LUT can accurately replace all the non-linear operations in popular BERT models with significant reductions in area, power consumption, and latency.


title: Self-energy corrected band-gap tuning induced by strain in the hexagonal boron phosphide monolayer

abstract: We performed a first-principles study of the electronic behavior of a 2D hexagonal boron phosphide monolayer (2D-h-BP). The system was deformed isotropically by applying a simultaneous tensile strain along the a and b crystal axes. We analyzed the band-gap evolution as function of the deformation percentage, ranging from 1% to 8%. Results show that the system behaves as a direct band-gap semiconductor, with the valence band maximum and conduction band minimum located at the K point (1/3, 1/3, 0) of the Brillouin zone. This behavior is unchanged despite the strain application. The band gap underestimation, as computed within the standard DFT, was corrected by applying the G0W0 approach. Trends in the band-gap behavior are the same within both approaches: for low deformation percentages, the band-gap grows linearly with a small slope, and at higher values, it grows very slow with a tendency to achieve saturation. Therefore, the band-gap is less sensitive to tensile strain for deformations near 8%. The origin of this band gap behavior is explained in terms of the projected density of states and charge densities, and it can be attributed to Coulomb interactions, and charge redistributions due to the applied tensile strain. To study the carrier mobility, we computed the electron and hole effective masses, finding high mobility for both carriers. Finally, the stability analysis of each strained system includes the calculation of phonon spectra, to assure the dynamical stability, the computation of elastic constants to evaluate the mechanical stability, and cohesive energies for exploring the thermodynamical stability. Results indicate that the boron phosphide monolayer is stable under the calculated tensile strains.


title: Minimal informational requirements for fitness

abstract: The existing concept of the "fitness value of information" provides a theoretical upper bound on the fitness advantage of using information concerning a fluctuating environment. Using concepts from rate-distortion theory, we develop a theoretical framework to answer a different pair of questions: What is the minimal amount of information needed for a population to achieve a certain growth rate? What is the minimal amount of information gain needed for one sub-population to achieve a certain average selection coefficient over another? We introduce a correspondence between fitness and distortion and solve for the rate-distortion functions of several systems using analytical and numerical methods. Because accurate information processing is energetically costly, our approach provides a theoretical basis for understanding evolutionary "design principles" underlying information-cost trade-offs.


title: ALX: Large Scale Matrix Factorization on TPUs

abstract: We present ALX, an open-source library for distributed matrix factorization using Alternating Least Squares, written in JAX. Our design allows for efficient use of the TPU architecture and scales well to matrix factorization problems of O(B) rows/columns by scaling the number of available TPU cores. In order to spur future research on large scale matrix factorization methods and to illustrate the scalability properties of our own implementation, we also built a real world web link prediction dataset called WebGraph. This dataset can be easily modeled as a matrix factorization problem. We created several variants of this dataset based on locality and sparsity properties of sub-graphs. The largest variant of WebGraph has around 365M nodes and training a single epoch finishes in about 20 minutes with 256 TPU cores. We include speed and performance numbers of ALX on all variants of WebGraph. Both the framework code and the dataset is open-sourced.


title: Revisiting local branching with a machine learning lens

abstract: Finding high-quality solutions to mixed-integer linear programming problems (MILPs) is of great importance for many practical applications. In this respect, the refinement heuristic local branching (LB) has been proposed to produce improving solutions and has been highly influential for the development of local search methods in MILP. The algorithm iteratively explores a sequence of solution neighborhoods defined by the so-called local branching constraint, namely, a linear inequality limiting the distance from a reference solution. For a LB algorithm, the choice of the neighborhood size is critical to performance. In this work, we study the relation between the size of the search neighborhood and the behavior of the underlying LB algorithm, and we devise a leaning based framework for predicting the best size for the specific instance to be solved. Furthermore, we have also investigated the relation between the time limit for exploring the LB neighborhood and the actual performance of LB scheme, and devised a strategy for adapting the time limit. We computationally show that the neighborhood size and time limit can indeed be learned, leading to improved performances and that the overall algorithm generalizes well both with respect to the instance size and, remarkably, across instances.


title: Exploring the s-process history in the Galactic disk: Cerium abundances and gradients in Open Clusters from the OCCAM/APOGEE sample

abstract: The APOGEE Open Cluster Chemical Abundances and Mapping (OCCAM) survey is used to probe the chemical evolution of the s-process element cerium in the Galactic disk. Cerium abundances were derived from measurements of Ce II lines in the APOGEE spectra using the Brussels Automatic Code for Characterizing High Accuracy Spectra (BACCHUS) in 218 stars belonging to 42 open clusters. Our results indicate that, in general, for Ages $<$ 4 Gyr, younger open clusters have higher [Ce/Fe] and [Ce/$\alpha$-element] ratios than older clusters. In addition, metallicity segregates open clusters in the [Ce/X]-Age plane (where X can be H, Fe, and the $\alpha$-elements O, Mg, Si, or Ca). These metallicity-dependant relations result in [Ce/Fe] and [Ce/$\alpha$] ratios with age that are not universal clocks. Radial gradients of [Ce/H] and [Ce/Fe] ratios in open clusters, binned by age, were derived for the first time, with d[Ce/H]dR$_{GC}$ being negative, while d[Ce/Fe]/dR$_{GC}$ is positive. [Ce/H] and [Ce/Fe] gradients are approximately constant over time, with the [Ce/Fe] gradient becoming slightly steeper, changing by $\sim$+0.009 dex-kpc$^{-1}$-Gyr$^{-1}$. Both the [Ce/H] and [Ce/Fe] gradients are shifted to lower values of [Ce/H] and [Ce/Fe] for older open clusters. The chemical pattern of Ce in open clusters across the Galactic disk is discussed within the context of s-process yields from AGB stars, $\sim$Gyr time delays in Ce enrichment of the interstellar medium, and the strong dependence of Ce nucleosynthesis on the metallicity of its AGB stellar sources.


title: A Divide-and-Conquer Algorithm for Distributed Optimization on Networks

abstract: In this paper, we consider networks with topologies described by some connected undirected graph ${\mathcal{G}}=(V, E)$ and with some agents (fusion centers) equipped with processing power and local peer-to-peer communication, and optimization problem $\min_{{\boldsymbol x}}\big\{F({\boldsymbol x})=\sum_{i\in V}f_i({\boldsymbol x})\big\}$ with local objective functions $f_i$ depending only on neighboring variables of the vertex $i\in V$. We introduce a divide-and-conquer algorithm to solve the above optimization problem in a distributed and decentralized manner. The proposed divide-and-conquer algorithm has exponential convergence, its computational cost is almost linear with respect to the size of the network, and it can be fully implemented at fusion centers of the network. Our numerical demonstrations also indicate that the proposed divide-and-conquer algorithm has superior performance than popular decentralized optimization methods do for the least squares problem with/without $\ell^1$ penalty.