title: A backend-agnostic, quantum-classical framework for simulations of chemistry in C++

abstract: As quantum computing hardware systems continue to advance, the research and development of performant, scalable, and extensible software architectures, languages, models, and compilers is equally as important in order to bring this novel coprocessing capability to a diverse group of domain computational scientists. For the field of quantum chemistry, applications and frameworks exists for modeling and simulation tasks that scale on heterogeneous classical architectures, and we envision the need for similar frameworks on heterogeneous quantum-classical platforms. Here we present the XACC system-level quantum computing framework as a platform for prototyping, developing, and deploying quantum-classical software that specifically targets chemistry applications. We review the fundamental design features in XACC, with special attention to its extensibility and modularity for key quantum programming workflow interfaces, and provide an overview of the interfaces most relevant to simulations of chemistry. A series of examples demonstrating some of the state-of-the-art chemistry algorithms currently implemented in XACC are presented, while also illustrating the various APIs that would enable the community to extend, modify, and devise new algorithms and applications in the realm of chemistry.


title: Data-Efficient Reinforcement Learning for Malaria Control

abstract: Sequential decision-making under cost-sensitive tasks is prohibitively daunting, especially for the problem that has a significant impact on people's daily lives, such as malaria control, treatment recommendation. The main challenge faced by policymakers is to learn a policy from scratch by interacting with a complex environment in a few trials. This work introduces a practical, data-efficient policy learning method, named Variance-Bonus Monte Carlo Tree Search~(VB-MCTS), which can copy with very little data and facilitate learning from scratch in only a few trials. Specifically, the solution is a model-based reinforcement learning method. To avoid model bias, we apply Gaussian Process~(GP) regression to estimate the transitions explicitly. With the GP world model, we propose a variance-bonus reward to measure the uncertainty about the world. Adding the reward to the planning with MCTS can result in more efficient and effective exploration. Furthermore, the derived polynomial sample complexity indicates that VB-MCTS is sample efficient. Finally, outstanding performance on a competitive world-level RL competition and extensive experimental results verify its advantage over the state-of-the-art on the challenging malaria control task.


title: A One-Line Proof of Leversha's "Quartet of Isogonal Conjugates" Theorem

abstract: We give a short and insightful proof of Gerry Leversha's elegant theorem regarding the isogonal conjugates of each of the vertices of a non-cyclic quadrilateral with respect to the triangle formed by the other three. It uses the Maple package RENE.txt, available from . http://www.math.rutgers.edu/~zeilberg/tokhniot/RENE.txt


title: Poisoning the Unlabeled Dataset of Semi-Supervised Learning

abstract: Semi-supervised machine learning models learn from a (small) set of labeled training examples, and a (large) set of unlabeled training examples. State-of-the-art models can reach within a few percentage points of fully-supervised training, while requiring 100x less labeled data. We study a new class of vulnerabilities: poisoning attacks that modify the unlabeled dataset. In order to be useful, unlabeled datasets are given strictly less review than labeled datasets, and adversaries can therefore poison them easily. By inserting maliciously-crafted unlabeled examples totaling just 0.1% of the dataset size, we can manipulate a model trained on this poisoned dataset to misclassify arbitrary examples at test time (as any desired label). Our attacks are highly effective across datasets and semi-supervised learning methods. We find that more accurate methods (thus more likely to be used) are significantly more vulnerable to poisoning attacks, and as such better training methods are unlikely to prevent this attack. To counter this we explore the space of defenses, and propose two methods that mitigate our attack.


title: An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

abstract: We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.


title: A general model of MnSi-like spiral magnets

abstract: A general, symmetry-allowed model of nearest-neighbour interactions for MnSi-like magnets is presented. A left-handed helical magnet phase occurs within a large parameter space of the model, which is explored via numerical simulation. The relations between microscopic features of the spiral structure and various model parameters, including an external magnetic field, are determined and show good agreement with predictions from free energy considerations. A skyrmion structure is stabilized near the boundary.


title: Homoclinic RG flows, or when relevant operators become irrelevant

abstract: We study an $\mathcal{N}=1$ supersymmetric quantum field theory with $O(M)\times O(N)$ symmetry. Working in $3-\epsilon$ dimensions, we calculate the beta functions up to second loop order and analyze in detail the Renormalization Group (RG) flow and its fixed points. We allow $N$ and $M$ to assume general real values, which results in them functioning as bifurcation parameters. In studying the behaviour of the model in the space of $M$ and $N$, we demarcate the region where the RG flow is non-monotonic and determine curves along which Hopf bifurcations take place. At a number of points in the space of $M$ and $N$ we find that the model exhibits an interesting phenomenon: at these points the RG flow possesses a fixed point located at real values of the coupling constants $g_i$ but with a stability matrix $\left(\frac{\partial \beta_i}{\partial g_j}\right)$ that is not diagonalizable and has a Jordan block of size two with zero eigenvalue. Such points correspond to logarithmic CFTs and represent Bogdanov-Takens bifurcations, a type of bifurcation known to give rise to a nearby homoclinic orbit - an RG flow that originates and terminates at the same fixed point. In the present example, we are able to employ analytic and numeric evidence to display the existence of the homoclinic RG flow.


title: Microscopic formulas for thermoelectric transport coefficients in lattice systems

abstract: A macroscopic description of thermoelectric phenomena involves several tensorial transport coefficients. Textbook microscopic Kubo formulas for them are plagued with ambiguities in the definitions of the current operators and the magnetization. We derive a version of these formulas for lattice systems which is free from ambiguities but contains additional terms compared to the textbook results. For symmetric components of thermoelectric tensors, we identify a large class of lattice systems for which the additional terms vanish with a natural choice of the energy current. To eliminate ambiguities in the skew-symmetric components, one needs to interpret them as relative quantities: only their differences for pairs of materials are well-defined.


title: Black hole remnants are not too fast to be dark matter

abstract: We comment on recent claims that recoil in the final stages of Hawking evaporation gives black hole remnants large velocities, rendering them inviable as a dark matter candidate. We point out that due to cosmic expansion, such large velocities at the final stages of evaporation are not in tension with the cold dark matter paradigm so long as they are attained at sufficiently early times. In particular, the predicted recoil velocities are robustly compatible with observations if the remnants form before the epoch of big bang nucleosynthesis, a requirement which is already imposed by the physics of nucleosynthesis itself.


title: Sharp smoothing properties of averages over curves

abstract: We prove sharp smoothing properties of the averaging operator defined by convolution with a measure on a smooth nondegenerate curve $\gamma$ in $\mathbb R^d$, $d\ge 3$. Despite the simple geometric structure of such curves, the sharp smoothing estimates have remained largely unknown except for those in low dimensions. Devising a novel inductive strategy, we obtain the optimal $L^p$ Sobolev regularity estimates, which settle the conjecture raised by Beltran-Guo-Hickman-Seeger. Besides, we show the sharp local smoothing estimates for every $d$. As a result, we establish, for the first time, nontrivial $L^p$ boundedness of the maximal average over dilations of $\gamma$ for $d\ge 4$.


title: Searches for supersymmetry in CMS

abstract: The results from the CMS search for supersymmetric particles based on Run-2 data recorded at a center-of-mass energy of 13 TeV are summarized. Strong and weak production of SUSY scenarios are considered. Results presented include the searches for squarks and gluinos, direct production of charginos, neutralinos, and sleptons. These searches involve final state objects comprising jets, missing transverse momentum, electrons or muons, taus or photons, as well as long-lived particles. The data in these searches are found to be consistent with standard model predictions and no significant excess is observed. Upper limits have been set on the masses of supersymmetric particles from a variety of search channels.


title: Stochastic Process Optimization of an Integrated Biorefinery

abstract: Planning of biorefinery operations is complicated by the stochastic nature of physical and chemical characteristics of biomass feedstock, such as, moisture level and carbohydrate content. Biomass characteristics affect the performance of the equipment which feed the reactor and the efficiency of the conversion process in a biorefinery. We propose a stochastic optimization model to identify a blend of feedstocks, inventory levels, and operating conditions of equipment to ensure a continuous flowing of biomass to the reactor while meeting the requirements of the biochemical conversion process. We propose a sample average approximation (SAA) of the model, and develop an efficient algorithm to solve the SAA model. A feedstock preprocessing process consists of two-stage grinding and pelleting is used to develop a case study. Extensive numerical analysis are conducted which lead to a number of observations. Our main observation is that sequencing bales based on moisture level and carbohydrate content leads to robust solutions that improve processing time and processing rate of the reactor. We provide a number of managerial insights that facilitate the implementation of the model proposed.


title: The Hubble Constant in the Axi-Higgs Universe

abstract: The $\Lambda$CDM model provides an excellent fit to the CMB data. However, a statistically significant tension emerges when its determination of the Hubble constant $H_0$ is compared to the local distance-redshift measurements. The axi-Higgs model, which couples an ultralight axion to the Higgs field, offers a specific variation of the $\Lambda$CDM model. It relaxes the $H_0$ tension as well as explains the $^7$Li puzzle in Big-Bang nucleosynthesis, the clustering $S_8$ tension with the weak-lensing data, and the observed isotropic cosmic birefringence in CMB. In this paper, we demonstrate how the $H_0$ and $S_8$ tensions can be relaxed simultaneously, by correlating the axion impacts on the early and late universe. In a benchmark scenario ($m=2 \times 10^{-30}$ eV) selected for experimental tests soon, the analysis combining the CMB+BAO+WL+SN data yields $H_0 = 69.9 \pm 1.5$ km/s/Mpc and $S_8 = 0.8045 \pm 0.0096$. Combining this (excluding the SN (supernovae) part) with the local distance-redshift measurements yields $H_0 = 72.42 \pm 0.76$ km/s/Mpc, while $S_8$ is slightly more suppressed.


title: Solo: A Lightweight Static Analysis for Differential Privacy

abstract: All current approaches for statically enforcing differential privacy in higher order languages make use of either linear or relational refinement types. A barrier to adoption for these approaches is the lack of support for expressing these "fancy types" in mainstream programming languages. For example, no mainstream language supports relational refinement types, and although Rust and modern versions of Haskell both employ some linear typing techniques, they are inadequate for embedding enforcement of differential privacy, which requires "full" linear types a la Girard. We propose a new type system that enforces differential privacy, avoids the use of linear and relational refinement types, and can be easily embedded in mainstream richly typed programming languages such as Scala, OCaml and Haskell. We demonstrate such an embedding in Haskell, demonstrate its expressiveness on case studies, and prove that our type-based enforcement of differential privacy is sound.


title: An Estimation of Online Video User Engagement from Features of Continuous Emotions

abstract: Portraying emotion and trustworthiness is known to increase the appeal of video content. However, the causal relationship between these signals and online user engagement is not well understood. This limited understanding is partly due to a scarcity in emotionally annotated data and the varied modalities which express user engagement online. In this contribution, we utilise a large dataset of YouTube review videos which includes ca. 600 hours of dimensional arousal, valence and trustworthiness annotations. We investigate features extracted from these signals against various user engagement indicators including views, like/dislike ratio, as well as the sentiment of comments. In doing so, we identify the positive and negative influences which single features have, as well as interpretable patterns in each dimension which relate to user engagement. Our results demonstrate that smaller boundary ranges and fluctuations for arousal lead to an increase in user engagement. Furthermore, the extracted time-series features reveal significant (p<0.05) correlations for each dimension, such as, count below signal mean (arousal), number of peaks (valence), and absolute energy (trustworthiness). From this, an effective combination of features is outlined for approaches aiming to automatically predict several user engagement indicators. In a user engagement prediction paradigm we compare all features against semi-automatic (cross-task), and automatic (task-specific) feature selection methods. These selected feature sets appear to outperform the usage of all features, e.g., using all features achieves 1.55 likes per day (Lp/d) mean absolute error from valence; this improves through semi-automatic and automatic selection to 1.33 and 1.23 Lp/d, respectively (data mean 9.72 Lp/d with a std. 28.75 Lp/d).


title: Remote Pathological Gait Classification System

abstract: Several pathologies can alter the way people walk, i.e. their gait. Gait analysis can therefore be used to detect impairments and help diagnose illnesses and assess patient recovery. Using vision-based systems, diagnoses could be done at home or in a clinic, with the needed computation being done remotely. State-of-the-art vision-based gait analysis systems use deep learning, requiring large datasets for training. However, to our best knowledge, the biggest publicly available pathological gait dataset contains only 10 subjects, simulating 4 gait pathologies. This paper presents a new dataset called GAIT-IT, captured from 21 subjects simulating 4 gait pathologies, with 2 severity levels, besides normal gait, being considerably larger than publicly available gait pathology datasets, allowing to train a deep learning model for gait pathology classification. Moreover, it was recorded in a professional studio, making it possible to obtain nearly perfect silhouettes, free of segmentation errors. Recognizing the importance of remote healthcare, this paper proposes a prototype of a web application allowing to upload a walking person's video, possibly acquired using a smartphone camera, and execute a web service that classifies the person's gait as normal or across different pathologies. The web application has a user friendly interface and could be used by healthcare professionals or other end users. An automatic gait analysis system is also developed and integrated with the web application for pathology classification. Compared to state-of-the-art solutions, it achieves a drastic reduction in the number of model parameters, which means significantly lower memory requirements, as well as lower training and execution times. Classification accuracy is on par with the state-of-the-art.


title: Time evolution of vortex rings with large radius and very concentrated vorticity

abstract: We study the time evolution of an incompressible fluid with axial symmetry without swirl when the vorticity is sharply concentrated on $N$ annuli of radii $\approx$ $r_0$ and thickness $\epsilon$. We prove that when $r_0= |\log \epsilon|^\alpha, \,\, \alpha>2$, the vorticity field of the fluid converges as $\epsilon \to 0$ to the point vortex model, at least for a small but positive time. This result generalizes a previous paper that assumed a power law for the relation between $r_0$ and $\epsilon$.


title: Learning 3D Granular Flow Simulations

abstract: Recently, the application of machine learning models has gained momentum in natural sciences and engineering, which is a natural fit due to the abundance of data in these fields. However, the modeling of physical processes from simulation data without first principle solutions remains difficult. Here, we present a Graph Neural Networks approach towards accurate modeling of complex 3D granular flow simulation processes created by the discrete element method LIGGGHTS and concentrate on simulations of physical systems found in real world applications like rotating drums and hoppers. We discuss how to implement Graph Neural Networks that deal with 3D objects, boundary conditions, particle - particle, and particle - boundary interactions such that an accurate modeling of relevant physical quantities is made possible. Finally, we compare the machine learning based trajectories to LIGGGHTS trajectories in terms of particle flows and mixing entropies.


title: Implicit differentiation for fast hyperparameter selection in non-smooth convex learning

abstract: Finding the optimal hyperparameters of a model can be cast as a bilevel optimization problem, typically solved using zero-order techniques. In this work we study first-order methods when the inner optimization problem is convex but non-smooth. We show that the forward-mode differentiation of proximal gradient descent and proximal coordinate descent yield sequences of Jacobians converging toward the exact Jacobian. Using implicit differentiation, we show it is possible to leverage the non-smoothness of the inner problem to speed up the computation. Finally, we provide a bound on the error made on the hypergradient when the inner optimization problem is solved approximately. Results on regression and classification problems reveal computational benefits for hyperparameter optimization, especially when multiple hyperparameters are required.


title: The impact of nuclear deformation on relativistic heavy-ion collisions: assessing consistency in nuclear physics across energy scales

abstract: In the hydrodynamic framework of heavy-ion collisions, elliptic flow, $v_2$, is sensitive to the quadrupole deformation, $\beta$, of the colliding ions. This enables one to test whether the established knowledge on the low-energy structure of nuclei is consistent with collider data from high-energy experiments. We derive a formula based on generic scaling laws of hydrodynamics to relate the difference in $v_2$ measured between collision systems that are close in size to the value of $\beta$ of the respective species. We validate our formula in simulations of 238U+238U and 197Au+197Au collisions at top Relativistic Heavy Ion Collider (RHIC) energy, and subsequently apply it to experimental data. Using the deformation of 238U from low-energy experiments, we find that RHIC $v_2$ data implies $0.16 \lesssim |\beta| \lesssim 0.20$ for 197Au nuclei, i.e., significantly more deformed than reported in the literature, posing an interesting puzzle in nuclear phenomenology.


title: Practical Verification of Quantum Properties in Quantum Approximate Optimization Runs

abstract: In order to assess whether quantum resources can provide an advantage over classical computation, it is necessary to characterize and benchmark the non-classical properties of quantum algorithms in a practical manner. In this paper, we show that using measurements in no more than 3 out of the possible $3^N$ bases, one can not only reconstruct the single-qubit reduced density matrices and measure the ability to create coherent superpositions, but also possibly verify entanglement across all $N$ qubits participating in the algorithm. We introduce a family of generalized Bell-type observables for which we establish an upper bound to the expectation values in fully separable states by proving a generalization of the Cauchy-Schwarz inequality, which may serve of independent interest. We demonstrate that a subset of such observables can serve as entanglement witnesses for QAOA-MaxCut states, and further argue that they are especially well tailored for this purpose by defining and computing an entanglement potency metric on witnesses. A subset of these observables also certify, in a weaker sense, the entanglement in GHZ states, which share the $\mathbb{Z}_2$ symmetry of QAOA-MaxCut. The construction of such witnesses follows directly from the cost Hamiltonian to be optimized, and not through the standard technique of using the projector of the state being certified. It may thus provide insights to construct similar witnesses for other variational algorithms prevalent in the NISQ era. We demonstrate our ideas with proof-of-concept experiments on the Rigetti Aspen-9 chip for ansatze containing up to 24 qubits.


title: Reduction of order, resummation and radiation reaction

abstract: The Landau-Lifshitz equation is the first in an infinite series of approximations to the Lorentz-Abraham-Dirac equation obtained from `reduction of order'. We show that this series is divergent, predicting wildly different dynamics at successive perturbative orders. Iterating reduction of order ad infinitum in a constant crossed field, we obtain an equation of motion which is free of the erratic behaviour of perturbation theory. We show that Borel-Pad\'e resummation of the divergent series accurately reproduces the dynamics of this equation, using as little as two perturbative coefficients. Comparing with the Lorentz-Abraham-Dirac equation, our results show that for large times the optimal order of truncation typically amounts to using the Landau-Lifshitz equation, but that this fails to capture the resummed dynamics over short times.


title: Real-Time Guarantees for Critical Traffic in IEEE 802.1Qbv TSN Networks with Unscheduled and Unsynchronized End-Systems

abstract: Time-Sensitive Networking (TSN) aims to extend the IEEE 802.1Q Ethernet standard with real-time and time-aware capabilities. Each device's transmission of time-critical frames is done according to a so-called Gate Control List (GCL) schedule via the timed-gate mechanism described in IEEE 802.1Qbv. Most schedule generation mechanisms for TSN have a constraining assumption that both switches and end-systems in the network must have at least the TSN capabilities related to scheduled gates and time synchronization. However, many TSN networks use off-the-shelf end-systems, e.g., for providing sensor data, which are not scheduled and/or synchronized. In this paper, we propose a more flexible scheduling strategy that considers a worst-case delay analysis within the scheduling synthesis step, leveraging the solution's optimality to support TSN networks with unscheduled and unsynchronized end-systems while still being able to guarantee bounded latency for critical messages. Our method enables real-world systems that feature off-the-shelf microcontrollers and sensor nodes without TSN capabilities connected to state-of-the-art TSN networks to communicate critical messages in a real-time fashion. We evaluate our approach using both synthetic and real-world test cases, comparing it with existing scheduling mechanisms. Furthermore, we use OMNET++ to validate the generated GCL schedules.


title: Overlooked centrifugal and Coriolis forces in the Atmosphere due to the Earth's motion in the solar system

abstract: The slow revolution of the Earth and Moon around their barycentrum does not induce Coriolis accelerations. On the other hand, the motion of Sun and Earth is a rotation with Coriolis forces which appear not to have been calculated yet, nor have the inertial accelerations within the system of motion of all three celestial bodies. It is the purpose of this contribution to evaluate the related Coriolis and centrifugal terms and to compare them to the available atmospheric standard terms. It is a main result that the revolution is of central importance in the combined dynamics of Earth, Moon and Sun. Covariant flow equations are well known tools for dealing with such complicated flow settings. They are used here to quantify the effects of the Earth's revolution around the Earth-Moon barycenter and its rotation around the Sun on the atmospheric circulation. It is found that the motion around the Sun adds time dependent terms to the standard Coriolis forces. The related centrifugal accelerations are presented. A major part of these accelerations is balanced by the gravitational attraction by Moon and Sun, but important unbalanced contributions remain. New light on the consequences of the Earth's revolution is shed by repeating the calculations for a rotating Earth-Moon pair. It is found that the revolution complicates the atmospheric dynamics.


title: Electrochemical removal of amphoteric ions

abstract: Several harmful or valuable ionic species present in sea, brackish and wastewaters are amphoteric, and thus their properties depend on the local water pH. Effective removal of these species can be challenging by conventional membrane technologies, necessitating chemical dosing of the feedwater to adjust its pH. Capacitive deionization (CDI) is an emerging membraneless technique for water treatment and desalination, based on electrosorption of salt ions into charging microporous electrodes. CDI cells show strong internally-generated pH variations during operation, and thus CDI can potentially remove amphoteric species without chemical dosing. However, development of this technique is inhibited by the complexities inherent to coupling of pH dynamics and amphoteric ion properties in a charging CDI cell. Here, we present a novel theoretical framework predicting the electrosorption of amphoteric species in flow-through electrode CDI cells. We demonstrate that such a model enables deep insight into factors affecting amphoteric species electrosorption, and conclude that important design rules for such systems are highly counter-intuitive. For example, we show both theoretically and experimentally that for boron removal the anode should be placed upstream, which runs counter to accepted wisdom in the CDI field. Overall, we show that to achieve target separations relying on coupled, complex phenomena, such as in the removal of amphoteric species, a theoretical CDI model is essential.


title: Market Potential for CO$_2$ Removal and Sequestration from Renewable Natural Gas Production in California

abstract: Bioenergy with Carbon Capture and Sequestration (BECCS) is critical for stringent climate change mitigation, but is commercially and technologically immature and resource-intensive. In California, state and federal fuel and climate policies can drive first-markets for BECCS. We develop a spatially explicit optimization model to assess niche markets for renewable natural gas (RNG) production with carbon capture and sequestration (CCS) from waste biomass in California. Existing biomass residues produce biogas and RNG and enable low-cost CCS through the upgrading process and CO$_2$ truck transport. Under current state and federal policy incentives, we could capture and sequester 2.9 million MT CO$_2$/year (0.7% of California's 2018 CO$_2$ emissions) and produce 93 PJ RNG/year (4% of California's 2018 natural gas demand) with a profit maximizing objective. Existing federal and state policies produce profits of \$11/GJ. Distributed RNG production with CCS potentially catalyzes markets and technologies for CO$_2$ capture, transport, and storage in California.


title: Refection Matrices for the (n+1)(2n+1)- Vertex Models: Diagonal Solutions

abstract: We have find the diagonal K matrix solutions of the reflection equations for a class of vertex models. These models have (n+1)(2n+1) vertices and are defined as two set of (n + 1) R matrices, solutions of the equations of Yang-Baxter equations. For a given value of \text{n} we find n!-\frac{1}{2}(n-2)(n-1) K diagonal matrices.


title: Motion-Augmented Self-Training for Video Recognition at Smaller Scale

abstract: The goal of this paper is to self-train a 3D convolutional neural network on an unlabeled video collection for deployment on small-scale video collections. As smaller video datasets benefit more from motion than appearance, we strive to train our network using optical flow, but avoid its computation during inference. We propose the first motion-augmented self-training regime, we call MotionFit. We start with supervised training of a motion model on a small, and labeled, video collection. With the motion model we generate pseudo-labels for a large unlabeled video collection, which enables us to transfer knowledge by learning to predict these pseudo-labels with an appearance model. Moreover, we introduce a multi-clip loss as a simple yet efficient way to improve the quality of the pseudo-labeling, even without additional auxiliary tasks. We also take into consideration the temporal granularity of videos during self-training of the appearance model, which was missed in previous works. As a result we obtain a strong motion-augmented representation model suited for video downstream tasks like action recognition and clip retrieval. On small-scale video datasets, MotionFit outperforms alternatives for knowledge transfer by 5%-8%, video-only self-supervision by 1%-7% and semi-supervised learning by 9%-18% using the same amount of class labels.


title: Particle creation in nonstationary large N quantum mechanics

abstract: We consider an analog of particle production in a quartic $O(N)$ quantum oscillator with time-dependent frequency, which is a toy model of particle production in the dynamical Casimir effect and de Sitter space. We calculate exact quantum averages, Keldysh propagator, and particle number using two different methods. First, we employ a kind of rotating wave approximation to estimate these quantities for small deviations from stationarity. Second, we extend these results to arbitrarily large deviations using the Schwinger-Keldysh diagrammatic technique. We show that in strongly nonstationary situations, including resonant oscillations, loop corrections to the tree-level expressions effectively result in an additional degree of freedom, $N \to N + \frac{3}{2}$, which modifies the average number and energy of created particles.


title: On Lottery Tickets and Minimal Task Representations in Deep Reinforcement Learning

abstract: The lottery ticket hypothesis questions the role of overparameterization in supervised deep learning. But how is the performance of winning lottery tickets affected by the distributional shift inherent to reinforcement learning problems? In this work, we address this question by comparing sparse agents who have to address the non-stationarity of the exploration-exploitation problem with supervised agents trained to imitate an expert. We show that feed-forward networks trained with behavioural cloning compared to reinforcement learning can be pruned to higher levels of sparsity without performance degradation. This suggests that in order to solve the RL-specific distributional shift agents require more degrees of freedom. Using a set of carefully designed baseline conditions, we find that the majority of the lottery ticket effect in both learning paradigms can be attributed to the identified mask rather than the weight initialization. The input layer mask selectively prunes entire input dimensions that turn out to be irrelevant for the task at hand. At a moderate level of sparsity the mask identified by iterative magnitude pruning yields minimal task-relevant representations, i.e., an interpretable inductive bias. Finally, we propose a simple initialization rescaling which promotes the robust identification of sparse task representations in low-dimensional control tasks.


title: Chemical post-processing of magneto-hydrodynamical simulations of star-forming regions: robustness and pitfalls

abstract: A common approach to model complex chemistry in numerical simulations is via post-processing of existing magneto-hydrodynamic simulations, relying on computing the evolution of chemistry over the dynamic history of a subset of particles from within the raw simulation. Here, we validate such a technique, assessing its ability to recover the abundances of chemical species, using the chemistry package KROME. We also assess, for the first time, the importance of the main free input parameters, by means of a direct comparison with a self-consistent state-of-the-art simulation in which chemistry was directly coupled to hydrodynamics. We have found that the post-processing is highly reliable, with an accuracy at the percent level, even when the most relaxed input parameters are employed. In particular, our results show that the number of particles used does not affect significantly the average properties, although it suppresses the appearance of possibly important spatial features. On the other hand, the choice of the integration time-step plays a crucial role. Longer integration time-steps can produce large errors, as the post-processing solution will be forced towards chemical equilibrium, a condition that does not always necessarily apply. When the interpolation-based reconstruction of chemical properties is performed, the errors further increase up to a factor of $\sim2$. Concluding, our results suggest that this technique is extremely useful when exploring the relative quantitative effect of different chemical parameters and/or networks, without the need of re-running simulations multiple times, but some care should be taken in the choice of particles sub-sample and integration time-step.


title: Stochastic gradient descent with noise of machine learning type. Part I: Discrete time analysis

abstract: Stochastic gradient descent (SGD) is one of the most popular algorithms in modern machine learning. The noise encountered in these applications is different from that in many theoretical analyses of stochastic gradient algorithms. In this article, we discuss some of the common properties of energy landscapes and stochastic noise encountered in machine learning problems, and how they affect SGD-based optimization. In particular, we show that the learning rate in SGD with machine learning noise can be chosen to be small, but uniformly positive for all times if the energy landscape resembles that of overparametrized deep learning problems. If the objective function satisfies a Lojasiewicz inequality, SGD converges to the global minimum exponentially fast, and even for functions which may have local minima, we establish almost sure convergence to the global minimum at an exponential rate from any finite energy initialization. The assumptions that we make in this result concern the behavior where the objective function is either small or large and the nature of the gradient noise, but the energy landscape is fairly unconstrained on the domain where the objective function takes values in an intermediate regime.


title: Pricing Private Data with Personalized Differential Privacy and Partial Arbitrage Freeness

abstract: There is a growing trend regarding perceiving personal data as a commodity. Existing studies have built frameworks and theories about how to determine an arbitrage-free price of a given query according to the privacy loss quantified by differential privacy. However, those studies have assumed that data buyers can purchase query answers with the arbitrary privacy loss of data owners, which may not be valid under strict privacy regulations and data owners' increasing privacy concerns. In this paper, we study how to empower data owners to control privacy loss in data trading. First, we propose a framework for trading personal data that enables data owners to bound their personalized privacy losses. Second, since bounded privacy losses indicate bounded utilities of query answers, we propose a reasonable relaxation of arbitrage freeness named partial arbitrage freeness, i.e., the guarantee of arbitrage-free pricing only for a limited range of utilities, which provides more possibilities for our market design. Third, to avoid arbitrage, we propose a general method for ensuring arbitrage freeness under personalized differential privacy. Fourth, to fully utilize data owners' personalized privacy loss bounds, we propose privacy budget allocation techniques to allocate privacy losses for queries under arbitrage freeness. Finally, we conduct experiments to verify the effectiveness of our proposed trading protocols.


title: The Pursuit of Knowledge: Discovering and Localizing Novel Categories using Dual Memory

abstract: We tackle object category discovery, which is the problem of discovering and localizing novel objects in a large unlabeled dataset. While existing methods show results on datasets with less cluttered scenes and fewer object instances per image, we present our results on the challenging COCO dataset. Moreover, we argue that, rather than discovering new categories from scratch, discovery algorithms can benefit from identifying what is already known and focusing their attention on the unknown. We propose a method that exploits prior knowledge about certain object types to discover new categories by leveraging two memory modules, namely Working and Semantic memory. We show the performance of our detector on the COCO minival dataset to demonstrate its in-the-wild capabilities.


title: Topologically-driven three-spin chiral exchange interactions treated from first principles

abstract: The mechanism behind the three-spin chiral interaction (TCI) included in the extended Heisenberg Hamiltonian and represented by an expression worked out recently (Phys.\ Rev.\ B, {\bf 101}, 174401 (2020)) is discussed. It is stressed that this approach provides a unique set of the multispin exchange parameters which are independent of each other either due to their different order of perturbation or due to different symmetry. This ensures in particular the specific properties of the TCI that were demonstrated previously via fully relativistic first principles calculations, and that result from the common influence of several issues not explicitly seen from the expression for the TCI parameters. Therefore, an interpretation of the TCI is suggested, showing explicitly its dependence on the relativistic spin-orbit coupling and on the topological orbital susceptibility (TOS). This is based on an expression for the TOS that is worked out on the same footing as the expression for the TCI. Using first-principles calculations we demonstrate in addition numerically the common topological properties of the TCI and TOS. To demonstrate the role of the relativistic spin-orbit coupling (SOC) for the TCI, a so-called topological' spin susceptibility (TSS) is introduced. This quantity characterizes the SOC induced spin magnetic moment on the atom in the presence of non-collinear magnetic structure, giving a connection between the TOS and TCI. Numerical results again support our conclusions.


title: Resampling Methods for Detecting Anisotropic Correlation Structure

abstract: This paper proposes parametric and non-parametric hypothesis testing algorithms for detecting anisotropy -- rotational variance of the covariance function in random fields. Both algorithms are based on resampling mechanisms, which enable avoiding relying on asymptotic assumptions, as is common in previous algorithms. The algorithms' performance is illustrated numerically in simulation experiments and on real datasets representing a variety of potential challenges.


title: Hamiltonian structure of 2D fluid dynamics with broken parity

abstract: Isotropic fluids in two spatial dimensions can break parity symmetry and sustain transverse stresses which do not lead to dissipation. Corresponding transport coefficients include odd viscosity, odd torque, and odd pressure. We consider an isotropic Galilean invariant fluid dynamics in the adiabatic regime with momentum and particle density conservation. We find conditions on transport coefficients that correspond to dissipationless and separately to Hamiltonian fluid dynamics. The restriction on the transport coefficients will help identify what kind of hydrodynamics can be obtained by coarse-graining a microscopic Hamiltonian system. Interestingly, not all parity-breaking transport coefficients lead to energy conservation and, generally, the fluid dynamics is energy conserving but not Hamiltonian. We show how this dynamics can be realized by imposing a nonholonomic constraint on the Hamiltonian system.


title: Vibrations of Dimers of Mechanically Coupled Nanostructures: Analytical and Numerical Modeling

abstract: The coupling effects affecting the vibrations of two close nanostructures (e.g., two metal nanoplates or nanospheres separated by a thin dielectric layer) may considerably alter their vibrational eigenfrequencies, as demonstrated by several recent experimental studies. In this work, we present theoretical investigations of these coupling processes based on a continuum mechanics approach, considering various systems composed by two identical nanostructures mechanically coupled by a spacer made of a different material and computing their eigenfrequencies as a function of the spacer thickness. We first discuss the vibrations of stacked slabs, a one-dimensional problem which can be treated analytically. The more complex configurations of dimers of rods or spheres coupled by a finite cylindrical spacer are then treated numerically. In all cases, the frequency shifts occurring for thin spacers can be simply interpreted as a modification of the boundary conditions of the problem as compared to the single nanostructure case, while those predicted near specific spacer thicknesses are ascribed to an avoided crossing effect, happening when the individual building blocks of the dimers (nanostructures and spacer) present common eigenfrequencies.


title: QuantumCumulants.jl: A Julia framework for generalized mean-field equations in open quantum systems

abstract: A full quantum mechanical treatment of open quantum systems via a Master equation is often limited by the size of the underlying Hilbert space. As an alternative, the dynamics can also be formulated in terms of systems of coupled differential equations for operators in the Heisenberg picture. This typically leads to an infinite hierarchy of equations for products of operators. A well-established approach to truncate this infinite set at the level of expectation values is to neglect quantum correlations of high order. This is systematically realized with a so-called cumulant expansion, which decomposes expectation values of operator products into products of a given lower order, leading to a closed set of equations. Here we present an open-source framework that fully automizes this approach: first, the equations of motion of operators up to a desired order are derived symbolically using predefined canonical commutation relations. Next, the resulting equations for the expectation values are expanded employing the cumulant expansion approach, where moments up to a chosen order specified by the user are included. Finally, a numerical solution can be directly obtained from the symbolic equations. After reviewing the theory we present the framework and showcase its usefulness in a few example problems.


title: Resolved Dwarf Galaxy Searches within ~5 Mpc with the Vera Rubin Observatory and Subaru Hyper Suprime-Cam

abstract: We present a preview of the faint dwarf galaxy discoveries that will be possible with the Vera C. Rubin Observatory and Subaru Hyper Suprime-Cam in the next decade. In this work, we combine deep ground-based images from the Panoramic Imaging Survey of Centaurus and Sculptor (PISCeS) and extensive image simulations to investigate the recovery of faint, resolved dwarf galaxies in the Local Volume with a matched-filter technique. We adopt three fiducial distances - 1.5, 3.5, 5 Mpc, and quantitatively evaluate the effects on dwarf detection of varied stellar backgrounds, ellipticity, and Milky Way foreground contamination and extinction. We show that our matched-filter method is powerful for identifying both compact and extended systems, and near-future surveys will be able to probe at least ~4.5 mag below the tip of the red giant branch (TRGB) for a distance of up to 1.5 Mpc, and ~2 mag below the TRGB at 5 Mpc. This will push the discovery frontier for resolved dwarf galaxies to fainter magnitudes, lower surface brightnesses, and larger distances. Our simulations show the secure census of dwarf galaxies down to $M_{V}$$\approx$-5, -7, -8, will be soon within reach, out to 1.5 Mpc, 3.5 Mpc, and 5 Mpc, respectively, allowing us to quantify the statistical fluctuations in satellite abundances around hosts, and parse environmental effects as a function of host properties.


title: Non-Gaussianity and the induced gravitational wave background

abstract: Scalar metric fluctuations generically source a spectrum of gravitational waves at second order in perturbation theory, poising gravitational wave experiments as potentially powerful probes of the small-scale curvature power spectrum. We perform a detailed study of the imprint of primordial non-Gaussianity on these induced gravitational waves, emphasizing the role of both the disconnected and connected components of the primoridal trispectrum. Specializing to local-type non-Gaussianity, we numerically compute all contributions and present results for a variety of enhanced primordial curvature power spectra.


title: MaNGA galaxy properties -- I. An extensive optical, mid-infrared photometric, and environmental catalogue

abstract: We present an extensive catalog of non-parametric structural properties derived from optical and mid-infrared imaging for 4585 galaxies from the MaNGA survey. DESI and WISE imaging are used to extract surface brightness profiles in the g, r, z, W1, W2 photometric bands. Our optical photometry takes advantage of the automated algorithm AutoProf and probes surface brightnesses that typically reach below 29 mag/arcsec^2 in the r band, while our WISE photometry achieves 28 mag/arcsec^2 in the W1 band. Neighbour density measures and central/satellite classifications are also provided for a large sub-sample of the MaNGA galaxies. Highlights of our analysis of galaxy light profiles include: (i) an extensive comparison of galaxian structural properties that illustrates the robustness of non-parametric extraction of light profiles over parametric methods; (ii) the ubiquity of bimodal structural properties suggesting the existence of galaxy families in multiple dimensions; and, (iii) an appreciation that structural properties measured relative to total light, regardless of the fractional level, are uncertain. We study galaxy scaling relations based on photometric parameters, and present detailed comparisons with literature and theory. Salient features of this analysis include the near-constancy of the slope and scatter of the size-luminosity and size-stellar mass relations for late-type galaxies with wavelength, and the saturation of the central surface density, measured within 1 kpc, for elliptical galaxies with M* > 10.7 Msol (corresponding to Sigma_1 ~ 10^{10} Msol/kpc^2). The multi-band photometry, environmental parameters, and structural scaling relations presented are useful constraints for stellar population and galaxy formation models.


title: The chemical make-up of the Sun: A 2020 vision

abstract: The chemical composition of the Sun is a fundamental yardstick in astronomy, relative to which essentially all cosmic objects are referenced. We reassess the solar abundances of all 83 long-lived elements, using highly realistic solar modelling and state-of-the-art spectroscopic analysis techniques coupled with the best available atomic data and observations. Our new improved analysis confirms the relatively low solar abundances of C, N, and O obtained in our previous 3D-based studies: $\log\epsilon_{\text{C}}=8.46\pm0.04$, $\log\epsilon_{\text{N}}=7.83\pm0.07$, and $\log\epsilon_{\text{O}}=8.69\pm0.04$. The revised solar abundances for the other elements also typically agree well with our previously recommended values with just Li, F, Ne, Mg, Cl, Kr, Rb, Rh, Ba, W, Ir, and Pb differing by more than $0.05$ dex. The here advocated present-day photospheric metal mass fraction is only slightly higher than our previous value, mainly due to the revised Ne abundance from Genesis solar wind measurements: $X_{\rm surface}=0.7438\pm0.0054$, $Y_{\rm surface}=0.2423\pm 0.0054$, $Z_{\rm surface}=0.0139\pm 0.0006$, and $Z_{\rm surface}/X_{\rm surface}=0.0187\pm 0.0009$. Overall the solar abundances agree well with those of CI chondritic meteorites but we identify a correlation with condensation temperature such that moderately volatile elements are enhanced by $\approx 0.04$ dex in the CI chondrites and refractory elements possibly depleted by $\approx 0.02$ dex, conflicting with conventional wisdom of the past half-century. Instead the solar chemical composition resembles more closely that of the fine-grained matrix of CM chondrites. The so-called solar modelling problem remains intact with our revised solar abundances, suggesting shortcomings with the computed opacities and/or treatment of mixing below the convection zone in existing standard solar models.


title: A diffuse tidal dwarf galaxy destined to fade out as a "dark galaxy"

abstract: We have explored the properties of a peculiar object detected in deep optical imaging and located at the tip of an HI tail emerging from Hickson Compact Group 16. Using multiband photometry from infrared to ultraviolet, we were able to constrain its stellar age to 58$^{+22}_{-9}$ Myr with a rather high metallicity of [Fe/H] = $-$0.16$^{+0.43}_{-0.41}$ for its stellar mass of M$_\star$ = 4.2$\times$10$^6$ M$_\odot$, a typical signature of tidal dwarf galaxies. The structural properties of this object are similar to those of diffuse galaxies, with a round and featureless morphology, a large effective radius (r$_{eff}$ = 1.5 kpc), and a low surface brightness (<$\mu_{g}$>$_{eff}$ = 25.6 mag arcsec$^{-2}$). Assuming that the object is dynamically stable and able to survive in the future, its fading in time via the aging of its stellar component will make it undetectable in optical observations in just $\sim$2 Gyr of evolution, even in the deepest current or future optical surveys. Its high HI mass, M(HI) = 3.9$\times$10$^8$ M$_\odot$, and future undetectable stellar component will make the object match the observational properties of dark galaxies, that is, dark matter halos that failed to turn gas into stars. Our work presents further observational evidence of the feasibility of HI tidal features becoming fake dark galaxies; it also shows the impact of stellar fading, particularly in high metallicity systems such as tidal dwarfs, in hiding aged stellar components beyond detection limits in optical observations.


title: Lyman-alpha opacities at z=4-6 require low mass, radiatively-suppressed galaxies to drive cosmic reionization

abstract: The high redshift Lyman-alpha forest, in particular the Gunn-Peterson trough, is the most unambiguous signature of the neutral to ionized transition of the intergalactic medium (IGM) taking place during the Epoch of Reionization (EoR). Recent studies, e.g. Kulkarni et al. (2019a) and Keating et al. (2019), showed that reproducing the observed Lyman-alpha opacities after overlap required a non-monotonous evolution of cosmic emissivity: rising, peaking at z=6, and then decreasing onwards to z=4. Such an evolution is puzzling considering galaxy build-up and the cosmic star formation rate are still continously on the rise at these epochs. Here, we use new RAMSES-CUDATON simulations to show that such a peaked evolution may occur naturally in a fully coupled radiation-hydrodynamical framework. In our fiducial run, cosmic emissivity at z>6 is dominated by a low mass (M$_{\rm DM}<2.10^9$ M$_{\odot}$), high escape fraction halo population, driving reionization, up to overlap. Approaching z=6, this population is radiatively suppressed due to the rising ionizing UV background, and its emissivity drops. In the meantime, the high mass halo population builds up and its emissivity rises, but not fast enough to compensate the dimming of the low mass haloes, because of low escape fractions. The combined ionizing emissivity of these two populations therefore naturally results in a rise and fall of the cosmic emissivity, from z=12 to z=4, with a peak at z=6. An alternative run, which features higher escape fractions for the high mass haloes and later suppression at low mass, leads to overshooting the ionizing rate, over-ionizing the IGM and therefore too low Lyman-alpha opacities.


title: Performance of Reservoir Discretizations in Quantum Transport Simulations

abstract: Quantum transport simulations often use explicit, yet finite, electronic reservoirs. These should converge to the correct continuum limit, albeit with a trade-off between discretization and computational cost. Here, we study this interplay for extended reservoir simulations, where relaxation maintains a bias or temperature drop across the system. Our analysis begins in the non-interacting limit, where we parameterize different discretizations to compare them on an even footing. For many-body systems, we develop a method to estimate the relaxation that best approximates the continuum by controlling virtual transitions in Kramers' turnover for the current. While some discretizations are more efficient for calculating currents, there is little benefit with regard to the overall state of the system. Any gains become marginal for many-body, tensor network simulations, where the relative performance of discretizations varies when sweeping other numerical controls. These results indicate that a given reservoir discretization may have little impact on numerical efficiency for certain computational tools. The choice of a relaxation parameter, however, is crucial, and the method we develop provides a reliable estimate of the optimal relaxation for finite reservoirs.


title: Fracton-elasticity duality in twisted moir\'e superlattices

abstract: We formulate a fracton-elasticity duality for twisted moir\'e superlattices, taking into account that they are incommensurate crystals with dissipative phason dynamics. From a dual tensor-gauge formulation, as compared to standard crystals, we identify twice the number of conserved charges that describe topological lattice defects, namely, disclinations and a new type of defect that we dub discompressions. The key implication of these conservation laws is that both glide and climb motions of lattice dislocations are suppressed, indicating that dislocation networks may become exceptionally stable. Our results also apply to other planar incommensurate crystals and quasicrystals.


title: A dispersive estimate of scalar contributions to hadronic light-by-light scattering

abstract: We consider the contribution of scalar resonances to hadronic light-by-light scattering in the anomalous magnetic moment of the muon. While the $f_0(500)$ has already been addressed in previous work using dispersion relations, heavier scalar resonances have only been estimated in hadronic models so far. Here, we compare an implementation of the $f_0(980)$ resonance in terms of the coupled-channel $S$-waves for $\gamma^*\gamma^*\to \pi\pi/\bar K K$ to a narrow-width approximation, which indicates $a_\mu^{\text{HLbL}}[f_0(980)]=-0.2(2)\times 10^{-11}$. With a similar estimate for the $a_0(980)$, the combined effect is thus well below $1\times 10^{-11}$ in absolute value. We also estimate the contribution of heavier scalar resonances. In view of the very uncertain situation concerning their two-photon couplings we suggest to treat them together with other resonances of similar mass when imposing the matching to short-distance constraints. Our final result is a refined estimate of the $S$-wave rescattering effects in the $\pi \pi$ and $\bar K K$ channel up to about $1.3$ GeV and including a narrow-width evaluation of the $a_0(980)$: $a_\mu^\text{HLbL}[\text{scalars}]=-9(1)\times 10^{-11}$.


title: A metallicity study of F, G, K and M dwarfs in the Coma Berenices open cluster from the APOGEE survey

abstract: We present a study of metallicities in a sample of main sequence stars with spectral types M, K, G and F ($T_{\rm eff}$ $\sim$ 3200 -- 6500K and log $g$ $\sim$ 4.3 -- 5.0 dex) belonging to the solar neighborhood young open cluster Coma Berenices. Metallicities were determined using the high-resolution (R=$\lambda$/$\Delta$ $\lambda$ $\sim$ 22,500) NIR spectra ($\lambda$1.51 -- $\lambda$1.69 $\mu$m) of the SDSS-IV APOGEE survey. Membership to the cluster was confirmed using previous studies in the literature along with APOGEE radial velocities and Gaia DR2. An LTE analysis using plane-parallel MARCS model atmospheres and the APOGEE DR16 line list was adopted to compute synthetic spectra and derive atmospheric parameters ($T_{\rm eff}$ and log $g$) for the M dwarfs and metallicities for the sample. The derived metallicities are near solar and are homogeneous at the level of the expected uncertainties, in particular when considering stars from a given stellar class. The mean metallicity computed for the sample of G, K, and M dwarfs is $\langle$[Fe/H]$\rangle$ = +0.04 $\pm$ 0.02 dex; however, the metallicities of the F-type stars are slightly lower, by about 0.04 dex, when compared to cooler and less massive members. Models of atomic diffusion can explain this modest abundance dip for the F dwarfs, indicating that atomic diffusion operates in Coma Berenices stars. The [Fe/H] dip occurs in nearly the same effective temperature range as that found in previous analyses of the lithium and beryllium abundances in Coma Berenices.


title: The ANITA Anomalous Events and Axion Quark Nuggets

abstract: The Antarctic Impulse Transient Antenna (\textsc{ANITA}) collaboration [1-3] has reported two anomalous events with noninverted polarity. These events are hard to explain in terms of conventional cosmic rays (CRs). We explore a new possible explanation for these anomalous events by suggesting that these events can be related to the dark matter (DM) annihilations within the so-called axion quark nugget (AQN) DM model. This model was initially invented for a completely different purpose to explain the observed similarity between the dark and the visible components in the Universe, i.e. $\Omega_{\rm DM}\sim \Omega_{\rm visible}$ without any fitting parameters. We investigate the signal properties of the upward-going AQN events, including the event rate, the pulse duration, and the electric field strength, and find them consistent with the observations. We list several features of the upward-going AQN events distinct from conventional CR air showers. The observations (or nonobservation) of these features may substantiate (or refute) our proposal.


title: Dynamics of binary black holes in low-mass young star clusters

abstract: Young star clusters are dynamically active stellar systems and are a common birthplace for massive stars. Low-mass star clusters ($\sim{}300-10^3$ M$_\odot$) are more numerous than massive systems and are characterized by a two-body relaxation time scale of a few Myr: the most massive stars sink to the cluster core and dynamically interact with each other even before they give birth to compact objects. Here, we explore the properties of black holes (BHs) and binary black holes (BBHs) formed in low-mass young star clusters, by means of a suite of $10^5$ direct $N$-body simulations with a high original binary fraction (100 % for stars with mass $>5$ M$_\odot$). Most BHs are ejected in the first $\sim{}20$ Myr by dynamical interactions. Dynamical exchanges are the main formation channel of BBHs, accounting for $\sim{}40-80$ % of all the systems. Most BBH mergers in low-mass young star clusters involve primary BHs with mass $<40$ M$_\odot$ and low mass ratios are extremely more common than in the field. Comparing our data with those of more massive star clusters ($10^3-3\times{}10^4$ M$_\odot$), we find a strong dependence of the percentage of exchanged BBHs on the mass of the host star cluster. In contrast, our results show just a mild correlation between the mass of the host star cluster and the efficiency of BBH mergers.


title: Interrelation of the environment of Ly$\alpha$ emitters and massive galaxies at $2<z<4.5$

abstract: We present a comparison of the spatial distributions of Ly$\alpha$ emitters (LAEs) and massive star-forming and quiescent galaxies (SFGs and QGs) at $2<z<4.5$. We use the photometric redshift catalog to select SFGs and QGs and a LAE catalog from intermediate/narrow bands obtained from the Subaru Telescope and Isaac-Newton Telescope in Cosmic Evolution Survey (COSMOS). We derive the auto-/cross- correlation signals of SFGs, QGs, and LAEs, and the galaxy overdensity distributions at the position of them. Whereas the cross-correlation signals of SFGs and QGs are explained solely by their halo mass differences, those of SFGs and LAEs are significantly lower than those expected from their auto-correlation signals, suggesting that some additional physical processes are segregating these two populations. Such segregation of SFGs and LAEs becomes stronger for rest-frame ultraviolet faint LAEs ($M_{\rm UV}>-20$). From the overdensity distributions, LAEs are located in less dense regions than SFGs and QGs, whereas SFGs and QGs tend to be in the same overdensity distributions. The different spatial distributions of LAEs compared to those of massive galaxies may be attributed to assembly bias or large amounts of neutral hydrogen gas associated with massive halos. These results reinforce the importance of exploring multiple galaxy populations in quantifying the intrinsic galaxy environment of the high-$z$ universe.


title: On the Mass Ratio Distribution of Black Hole Mergers in Triple Systems

abstract: Observations have shown that the majority of massive stars, progenitors of black holes (BHs), have on average more than one stellar companion. In triple systems, wide inner binaries can be driven to a merger by the third body due to long-term secular interactions, most notably by the eccentric Lidov-Kozai effect. In this study, we explore the properties of BH mergers in triple systems and compare their population properties to those of binaries produced in isolation and assembled in dense star clusters. Using the same stellar physics and identical assumptions for the initial populations of binaries and triples, we show that stellar triples yield a significantly flatter mass ratio distribution from $q=1$ down to $q\sim0.3$ than either binary stars or dense stellar clusters, similar to the population properties inferred from the most recent catalog of gravitational-wave events, though we do not claim that all the observed events can be accounted for with triples. While hierarchical mergers in clusters can also produce asymmetric mass ratios, the unique spins of such mergers can be used to distinguished them from those produced from stellar triples. All three channels occupy distinct regions in total mass-mass ratio space, which may allow them to be disentangled as more BH mergers are detected by LIGO, Virgo, and KAGRA.


title: Global 21-cm signal extraction from foreground and instrumental effects IV: Accounting for realistic instrument uncertainties and their overlap with foreground and signal models

abstract: All 21-cm signal experiments rely on electronic receivers that affect the data via both multiplicative and additive biases through the receiver's gain and noise temperature. While experiments attempt to remove these biases, the residuals of their imperfect calibration techniques can still confuse signal extraction algorithms. In this paper, the fourth and final installment of our pipeline series, we present a technique for fitting out receiver effects as efficiently as possible. The fact that the gain and global signal, which are multiplied in the observation equation, must both be modeled implies that the model of the data is nonlinear in its parameters, making numerical sampling the only way to explore the parameter distribution rigorously. However, multi-spectra fits, which are necessary to extract the signal confidently as demonstrated in the third paper of the series, often require large numbers of foreground parameters, increasing the dimension of the posterior distribution that must be explored and therefore causing numerical sampling inefficiencies. Building upon techniques in the second paper of the series, we outline a method to explore the full parameter distribution by numerically sampling a small subset of the parameters and analytically marginalizing over the others. We test this method in simulation using a type-I Chebyshev band-pass filter gain model and a fast signal model based on a spline between local extrema. The method works efficiently, converging quickly to the posterior signal parameter distribution. The final signal uncertainties are of the same order as the noise in the data.


title: Telltale signs of metal recycling in the circumgalactic medium of a $z \sim 0.77$ galaxy

abstract: We present gravitational-arc tomography of the cool-warm enriched circumgalactic medium (CGM) of an isolated galaxy (``G1'') at $z \approx 0.77$. Combining VLT/MUSE adaptive-optics and Magellan/MagE echelle spectroscopy we obtain partially-resolved kinematics of MgII in absorption and [OII] in emission. The unique arc configuration allows us to probe 42 spatially independent arc positions transverse to G1, plus 4 positions in front of it. The transverse positions cover G1's minor and major axes at impact parameters of $\approx 10-30$ kpc and $\approx 60$ kpc, respectively. We observe a direct kinematic connection between the cool-warm enriched CGM (traced by MgII) and the interstellar medium (traced by [OII]). This provides strong evidence for the existence of an extended disc that co-rotates with the galaxy out to tens of kiloparsecs. The MgII velocity dispersion ($\sigma \approx 30-100$ km s$^{-1}$, depending on position) is of the same order as the modeled galaxy rotational velocity ($v_{\rm rot} \approx 80$ km s$^{-1}$), providing evidence for the presence of a turbulent and pressure-supported CGM component. We regard the absorption to be modulated by a galactic-scale outflow, as it offers a natural scenario for the observed line-of-sight dispersion and asymmetric profiles observed against both the arcs and the galaxy. An extended enriched co-rotating disc together with the signatures of a galactic outflow, are telltale signs of metal recycling in the $z\sim 1$ CGM.


title: Hubble Space Telescope imaging of the extremely metal-poor globular cluster EXT8 in Messier 31

abstract: We recently found the globular cluster (GC) EXT8 in M31 to have an extremely low metallicity of [Fe/H]=-2.91+/-0.04 using high-resolution spectroscopy. Here we present a colour-magnitude diagram (CMD) for EXT8, obtained with the Wide Field Camera 3 on board the Hubble Space Telescope. Compared with the CMDs of metal-poor Galactic GCs, we find that the upper red giant branch (RGB) of EXT8 is about 0.03 mag bluer in F606W-F814W and slightly steeper, as expected from the low spectroscopic metallicity. The observed colour spread on the upper RGB is consistent with being caused entirely by the measurement uncertainties, and we place an upper limit of sigma(F606W-F814W)=0.015 mag on any intrinsic colour spread. The corresponding metallicity spread can be up to sigma([Fe/H])=0.2 dex or >0.7 dex, depending on the isochrone library adopted. The horizontal branch (HB) is located mostly on the blue side of the instability strip and has a tail extending to at least M(F606W)=+3, as in the Galactic GC M15. We identify two candidate RR Lyrae variables and several UV-luminous post-HB/post AGB star candidates, including one very bright (M(F300X)=-3.2) source near the centre of EXT8. The surface brightness of EXT8 out to a radius of 25 arcsec is well fitted by a Wilson-type profile with an ellipticity of epsilon=0.20, a semi-major axis core radius of 0.25", and a central surface brightness of 15.2 mag per square arcsec in the F606W band, with no evidence of extra-tidal structure. Overall, EXT8 has properties consistent with it being a "normal", but very metal-poor GC, and its combination of relatively high mass and very low metallicity thus remains challenging to explain in the context of GC formation theories operating within the hierarchical galaxy assembly paradigm.


title: Galaxy Evolution in all Five CANDELS Fields and IllustrisTNG: Morphological, Structural, and the Major Merger Evolution to $z \sim 3$

abstract: A fundamental feature of galaxies is their structure, yet we are just now understanding the evolution of structural properties in quantitative ways. As such, we explore the quantitative non-parametric structural evolution of 16,778 galaxies up to $z\sim3$ in all five CANDELS fields, the largest collection of high resolution images of distant galaxies to date. Our goal is to investigate how the structure, as opposed to size, surface brightness, or mass, changes with time. In particular, we investigate how the concentration and asymmetry of light evolve in the rest-frame optical. To interpret our galaxy structure measurements, we also run and analyse 300 simulation realisations from IllustrisTNG to determine the timescale of mergers for the CAS system. We measure that from $z=0-3$, the median asymmetry merger timescale is $0.56^{+0.23}_{-0.18}$Gyr, and find it does not vary with redshift. Using this data, we find that galaxies become progressively asymmetric at a given mass at higher redshifts and we derive merger rates which scale as $\sim(1+z)^{1.87\pm0.04}$Gyr$^{-1}$, which agrees well with recent machine learning and galaxy pair approaches, removing previous inconsistencies. We also show that far-infrared selected galaxies that are invisible to \textit{HST} have a negligible effect on our measurements. We also find that galaxies are more concentrated at higher redshifts. We interpret this as a sign of how their formation occurs from a smaller initial galaxy that later grows into a larger one through mergers, consistent with the size growth of galaxies from `inside-out', suggesting that the centres are the oldest parts of most galaxies.


title: The Three Hundred Project: The stellar angular momentum evolution of cluster galaxies

abstract: Using 324 numerically modelled galaxy clusters as provided by THE THREE HUNDRED project, we study the evolution of the kinematic properties of the stellar component of haloes on first infall. We select objects with M$_{\textrm{star}}>5\times10^{10} h^{-1}M_{\odot}$ within $3R_{200}$ of the main cluster halo at $z=0$ and follow their progenitors. We find that although haloes are stripped of their dark matter and gas after entering the main cluster halo, there is practically no change in their stellar kinematics. For the vast majority of our `galaxies' -- defined as the central stellar component found within the haloes that form our sample -- their kinematic properties, as described by the fraction of ordered rotation, and their position in the specific stellar angular momentum$-$stellar mass plane $j_{\rm star}$ -- M$_{\rm star}$, are mostly unchanged by the influence of the central host cluster. However, for a small number of infalling galaxies, stellar mergers and encounters with remnant stellar cores close to the centre of the main cluster, particularly during pericentre passage, are able to spin-up their stellar component by $z=0$.


title: A fluid simulation system based on the MPS method

abstract: Fluid flow simulation is a highly active area with applications in a wide range of engineering problems and interactive systems. Meshless methods like the Moving Particle Semi-implicit (MPS) are a great alternative to deal efficiently with large deformations and free-surface flow. However, mesh-based approaches can achieve higher numerical precision than particle-based techniques with a performance cost. This paper presents a numerically stable and parallelized system that benefits from advances in the literature and parallel computing to obtain an adaptable MPS method. The proposed technique can simulate liquids using different approaches, such as two ways to calculate the particles' pressure, turbulent flow, and multiphase interaction. The method is evaluated under traditional test cases presenting comparable results to recent techniques. This work integrates the previously mentioned advances into a single solution, which can switch on improvements, such as better momentum conservation and less spurious pressure oscillations, through a graphical interface. The code is entirely open-source under the GPLv3 free software license. The GPU-accelerated code reached speedups ranging from 3 to 43 times, depending on the total number of particles. The simulation runs at one fps for a case with approximately 200,000 particles. Code: https://github.com/andreluizbvs/VoxarMPS


title: A Parameter Space Exploration of the Minimal SU(5) Unification

abstract: We present phenomenological study of the most minimal realistic $SU(5)$ model that owns its predictivity solely to the gauge symmetry and the representational content. The model is built entirely out of the fields residing in the first five lowest dimensional representations that transform non-trivially under the $SU(5)$ gauge group. It has eighteen real parameters and fourteen phases, all in all, to address experimental observables of the Standard Model fermions and accomplishes that via simultaneous use of three different mass generation mechanisms. Furthermore, it inextricably links the origin of the neutrino mass to the experimentally observed difference between the down-type quark and charged lepton masses. The main predictions of the model are that $(i)$ the neutrinos are Majorana particles, $(ii)$ one neutrino is massless, $(iii)$ the neutrinos have normal mass ordering, and $(iv)$ there are four new scalar multiplets at or below a $120$ TeV mass scale. A one-loop analysis demonstrates that an improvement of the current $p \rightarrow \pi^0 e^+$ partial lifetime limit by a factor of $2$, $15$, and $96$ would require these four scalar multiplets to reside at or below the $100$ TeV, $10$ TeV, and $1$ TeV mass scales, respectively.


title: Non-perturbative determination of collisional broadening and medium induced radiation in QCD plasmas

abstract: We supply recently obtained results from lattice EQCD with the correct UV limit to construct the collisional broadening kernel $C(b_\perp)$ in a QCD plasma. We discuss the limiting behavior of $C(b_\perp)$ at small and large impact parameters $b_\perp$, and illustrate how the results can be used to compute medium-induced radiation rates.


title: On Moment Matching for Stochastic Systems

abstract: In this paper we study the problem of model reduction by moment matching for stochastic systems. We characterize the mathematical object which generalizes the notion of moment to stochastic differential equations and we find a class of models which achieve moment matching. However, differently from the deterministic case, these reduced-order models cannot be considered "simpler" because of the high computational cost paid to determine the moment. To overcome this difficulty, we relax the moment matching problem in two different ways and we present two classes of reduced-order models which, approximately matching the stochastic moment, are computationally tractable.


title: HI-to-H$_2$ Transitions in Dust-Free Interstellar Gas

abstract: We present numerical computations and analysis of atomic to molecular (HI-to-H$_2$) transitions in cool ($\sim$100 K) low-metallicity dust-free (primordial) gas, in which molecule formation occurs via cosmic-ray driven negative ion chemistry, and removal is by a combination of far-UV photodissociation and cosmic-ray ionization and dissociation. For any gas temperature, the behavior depends on the ratio of the Lyman-Werner (LW) band FUV intensity to gas density, $I_{\rm LW}/n$, and the ratio of the cosmic-ray ionization rate to the gas density, $\zeta/n$. We present sets of HI-to-H$_2$ abundance profiles for a wide range of $\zeta/n$ and $I_{\rm LW}/n$, for dust-free gas. We determine the conditions for which H$_2$ absorption line self-shielding in optically thick clouds enables a transition from atomic to molecular form for ionization-driven chemistry. We also examine the effects of cosmic-ray energy losses on the atomic and molecular density profiles and transition points. For a unit Galactic interstellar FUV field intensity ($I_{\rm LW}=1$) with LW flux $2.07\times 10^7$ photons cm$^{-2}$ s$^{-1}$, and a uniform cosmic-ray ionization rate $\zeta=10^{-16}$ s$^{-1}$, an HI-to-H$_2$ transition occurs at a total hydrogen gas column density of $4\times 10^{21}$ cm$^{-2}$, within $3\times 10^7$ yr, for a gas volume density of $n=10^6$ cm$^{-3}$ at 100 K. For these parameters, the dust-free limit obtains for a dust-to-gas ratio Z$^\prime_d \lesssim 10^{-5}$, which may be reached for overall metallicities $Z^\prime\lesssim 0.01$ relative to Galactic solar values.


title: Properties of the Obscuring Torus in NGC 1052 from Multi-epoch Broadband X-ray Spectroscopy

abstract: Obscuration of the innermost parts of active galactic nuclei (AGN) is observed in the majority of the population both in the nearby universe and at high redshift. However, the nature of the structures causing obscuration, especially in low-luminosity AGN, is poorly understood at present. We present a novel approach to multi-epoch broadband X-ray spectroscopy, anchored in the long-term average spectrum in the hard X-ray band, applied to the nearby, X-ray bright AGN in the galaxy NGC 1052. From spectral features due to X-ray reprocessing in the circumnuclear material, based on a simple, uniform-density torus X-ray reprocessing model, we find a covering factor of 80-100% and a globally averaged column density in the range (1-2) x 10^23 cm^-2. This closely matches the independently measured variable line-of-sight column density range, leading to a straightforward and self-consistent picture of the obscuring torus in NGC 1052, similar to several other AGN in recent literature. Comparing this X-ray-constrained torus model with measurements of spatially resolved sub-parsec absorption from radio observations, we find that it may be possible to account for both X-ray and radio data with a torus model featuring a steep density gradient along the axis of the relativistic jets. This provides a valuable direction for the development of improved physical models for the circumnuclear environment in NGC 1052 and potentially in a wider class of AGN.


title: A reconfigurable neural network ASIC for detector front-end data compression at the HL-LHC

abstract: Despite advances in the programmable logic capabilities of modern trigger systems, a significant bottleneck remains in the amount of data to be transported from the detector to off-detector logic where trigger decisions are made. We demonstrate that a neural network autoencoder model can be implemented in a radiation tolerant ASIC to perform lossy data compression alleviating the data transmission problem while preserving critical information of the detector energy profile. For our application, we consider the high-granularity calorimeter from the CMS experiment at the CERN Large Hadron Collider. The advantage of the machine learning approach is in the flexibility and configurability of the algorithm. By changing the neural network weights, a unique data compression algorithm can be deployed for each sensor in different detector regions, and changing detector or collider conditions. To meet area, performance, and power constraints, we perform a quantization-aware training to create an optimized neural network hardware implementation. The design is achieved through the use of high-level synthesis tools and the hls4ml framework, and was processed through synthesis and physical layout flows based on a LP CMOS 65 nm technology node. The flow anticipates 200 Mrad of ionizing radiation to select gates, and reports a total area of 3.6 mm^2 and consumes 95 mW of power. The simulated energy consumption per inference is 2.4 nJ. This is the first radiation tolerant on-detector ASIC implementation of a neural network that has been designed for particle physics applications.


title: 2-distance list $(\Delta + 3)$-coloring of sparse graphs

abstract: A 2-distance list k-coloring of a graph is a proper coloring of the vertices where each vertex has a list of at least k available colors and vertices at distance at most 2 cannot share the same color. We prove the existence of a 2-distance list $(\Delta + 3)$-coloring for graphs with maximum average degree less than $\frac83$ and maximum degree $\Delta\geq 4$ as well as graphs with maximum average degree less than $\frac{14}5$ and maximum degree $\Delta\geq 6$.


title: Witnessing Bell violations through probabilistic negativity

abstract: Bell's theorem shows that no hidden-variable model can explain the measurement statistics of a quantum system shared between two parties, thus ruling out a classical (local) understanding of nature. In this work we demonstrate that by relaxing the positivity restriction in the hidden-variable probability distribution it is possible to derive quasiprobabilistic Bell inequalities whose sharp upper bound is written in terms of a negativity witness of said distribution. This provides an analytic solution for the amount of negativity necessary to violate the CHSH inequality by an arbitrary amount, therefore revealing the amount of negativity required to emulate the quantum statistics in a Bell test.


title: Colourful Poincar\'e symmetry, gravity and particle actions

abstract: We construct a generalisation of the three-dimensional Poincar\'e algebra that also includes a colour symmetry factor. This algebra can be used to define coloured Poincar\'e gravity in three space-time dimensions as well as to study generalisations of massive and massless free particle models. We present various such generalised particle models that differ in which orbits of the coloured Poincar\'e symmetry are described. Our approach can be seen as a stepping stone towards the description of particles interacting with a non-abelian background field or as a starting point for a worldline formulation of an associated quantum field theory.


title: Solving the pooling problem at scale with extensible solver GALINI

abstract: This paper presents a Python library to model pooling problems, a class of network flow problems with many engineering applications. The library automatically generates a mixed-integer quadratically-constrained quadratic optimization problem from a given network structure. The library additionally uses the network structure to build 1) a convex linear relaxation of the non-convex quadratic program and 2) a mixed-integer linear restriction of the problem. We integrate the pooling network library with galini, an open-source extensible global solver for quadratic optimization. We demonstrate galini's extensible characteristics by using the pooling library to develop two galini plug-ins: 1) a cut generator plug-in that adds valid inequalities in the galini cut loop and 2) a primal heuristic plug-in that uses the mixed-integer linear restriction. We test galini on large scale pooling problems and show that, thanks to the good upper bound provided by the mixed-integer linear restriction and the good lower bounds provided by the convex relaxation, we obtain optimality gaps that are competitive with Gurobi 9.1 on the largest problem instances.


title: Height Estimation of Children under Five Years using Depth Images

abstract: Malnutrition is a global health crisis and is the leading cause of death among children under five. Detecting malnutrition requires anthropometric measurements of weight, height, and middle-upper arm circumference. However, measuring them accurately is a challenge, especially in the global south, due to limited resources. In this work, we propose a CNN-based approach to estimate the height of standing children under five years from depth images collected using a smart-phone. According to the SMART Methodology Manual [5], the acceptable accuracy for height is less than 1.4 cm. On training our deep learning model on 87131 depth images, our model achieved an average mean absolute error of 1.64% on 57064 test images. For 70.3% test images, we estimated height accurately within the acceptable 1.4 cm range. Thus, our proposed solution can accurately detect stunting (low height-for-age) in standing children below five years of age.


title: Twofold van Hove singularity and origin of charge order in topological kagome superconductor CsV3Sb5

abstract: The layered vanadium antimonides AV3Sb5 (A = K, Rb, Cs) are a recently discovered family of topological kagome metals with a rich phenomenology of strongly correlated electronic phases including charge order and superconductivity. Understanding how the singularities inherent to the kagome electronic structure are linked to the observed many-body phases is a topic of great interest and relevance. Here, we combine angle-resolved photoemission spectroscopy and density functional theory to reveal multiple kagome-derived van Hove singularities (vHs) coexisting near the Fermi level of CsV3Sb5 and analyze their contribution to electronic symmetry breaking. Intriguingly, the vHs in CsV3Sb5 have two distinct flavors - p-type and m-type - which originate from their pure and mixed sublattice characters, respectively. This twofold vHs is unique property of the kagome lattice, and its flavor critically determines the pairing symmetry and ground states emerging in AV3Sb5 series. We establish that, among the multiple vHs in CsV3Sb5, the m-type vHs of the dxz/dyz kagome band and the p-type vHs of the dxy/dx2-y2 kagome band cross the Fermi level to set the stage for electronic symmetry breaking. The former band exhibits pronounced Fermi surface nesting, while the latter contributes via higher-order vHs. Our work reveals the essential role of kagome-derived vHs for the collective phenomena realized in the AV3Sb5 family, paving the way to a deeper understanding of strongly correlated topological kagome systems.


title: Metric Comparisons of Relations

abstract: This paper defines a new pseudometric for binary relations between finite sets that measures consensus among subsets. The main results are (1) a concise restatement of this pseudometric with an intuitively appealing interpretation via a full and faithful functor from the category of relations to a category of relation multisets and (2) that the pseudometric can be bounded without an expensive search of possible mappings, based solely on the dimensions of the relations themselves. Additionally, (3) an algorithm is described to calculate this bound with time and memory complexity at worst proportional to the product of those dimensions: $\mathcal{O}(m \times n).$. The tools developed in this paper should find broad application in exploratory data analysis. We provide one such application by briefly exploring ad hoc consensus specifications for the well-known PDF file format.


title: Data Augmentation by Concatenation for Low-Resource Translation: A Mystery and a Solution

abstract: In this paper, we investigate the driving factors behind concatenation, a simple but effective data augmentation method for low-resource neural machine translation. Our experiments suggest that discourse context is unlikely the cause for the improvement of about +1 BLEU across four language pairs. Instead, we demonstrate that the improvement comes from three other factors unrelated to discourse: context diversity, length diversity, and (to a lesser extent) position shifting.


title: Unconditional energy dissipation and error estimates of the SAV Fourier spectral method for nonlinear fractional generalized wave equation

abstract: In this paper, we consider a second-order scalar auxiliary variable (SAV) Fourier spectral method to solve the nonlinear fractional generalized wave equation. Unconditional energy conservation or dissipation properties of the fully discrete scheme are first established. Next, we utilize the temporal-spatial error splitting argument to obtain unconditional optimal error estimate of the fully discrete scheme, which overcomes time-step restrictions caused by strongly nonlinear system, or the restrictions that the nonlinear term needs to satisfy the assumption of global Lipschitz condition in all previous works for fractional undamped or damped wave equations. Finally, some numerical experiments are presented to confirm our theoretical analysis.


title: On the Stability of Citation Networks

abstract: Citation networks can reveal many important information regarding the development of science and the relationship between different areas of knowledge. Thus, many studies have analyzed the topological properties of such networks. Frequently, citation networks are created using articles acquired from a set of relevant keywords or queries. Here, we study the robustness of citation networks with regards to the keywords that were used for collecting the respective articles. A perturbation approach is proposed, in which the influence of missing keywords on the topology and community structure of citation networks is quantified. In addition, the relationship between keywords and the community structure of citation networks is studied using networks generated from a simple model. We find that, owing to its highly modular structure, the community structure of citation networks tends to be preserved even when many relevant keywords are left out. Furthermore, the proposed model can reflect the impact of missing keywords on different situations.


title: Non-local magnon transport in a magnetic domain wall wave guide

abstract: Magnetic domain walls function as a wave guide for low energy magnons. In this paper we develop the theory for the non-local transport of these bound magnons through a ferromagnetic insulator that are injected and detected electrically in adjacent normal metal leads by spin-flip scatting processes and the (inverse) spin-Hall effect. Our set-up requires a twofold degeneracy of the magnetic ground state, which we realize by an easy axis and hard axis anisotropy, in the ferromagnetic insulator. This is readily provided by a broad range of materials. The domain wall is a a topologically protected feature of the system and we obtain the non-local spin transport through it. Thereby we provide a framework for reconfigureable magnonic devices.


title: Effectively Leveraging Attributes for Visual Similarity

abstract: Measuring similarity between two images often requires performing complex reasoning along different axes (e.g., color, texture, or shape). Insights into what might be important for measuring similarity can can be provided by annotated attributes, but prior work tends to view these annotations as complete, resulting in them using a simplistic approach of predicting attributes on single images, which are, in turn, used to measure similarity. However, it is impractical for a dataset to fully annotate every attribute that may be important. Thus, only representing images based on these incomplete annotations may miss out on key information. To address this issue, we propose the Pairwise Attribute-informed similarity Network (PAN), which breaks similarity learning into capturing similarity conditions and relevance scores from a joint representation of two images. This enables our model to identify that two images contain the same attribute, but can have it deemed irrelevant (e.g., due to fine-grained differences between them) and ignored for measuring similarity between the two images. Notably, while prior methods of using attribute annotations are often unable to outperform prior art, PAN obtains a 4-9% improvement on compatibility prediction between clothing items on Polyvore Outfits, a 5% gain on few shot classification of images using Caltech-UCSD Birds (CUB), and over 1% boost to Recall@1 on In-Shop Clothes Retrieval. Implementation available at https://github.com/samarth4149/PAN


title: Learning to Continuously Optimize Wireless Resource in a Dynamic Environment: A Bilevel Optimization Perspective

abstract: There has been a growing interest in developing data-driven, and in particular deep neural network (DNN) based methods for modern communication tasks. For a few popular tasks such as power control, beamforming, and MIMO detection, these methods achieve state-of-the-art performance while requiring less computational efforts, less resources for acquiring channel state information (CSI), etc. However, it is often challenging for these approaches to learn in a dynamic environment. This work develops a new approach that enables data-driven methods to continuously learn and optimize resource allocation strategies in a dynamic environment. Specifically, we consider an ``episodically dynamic" setting where the environment statistics change in ``episodes", and in each episode the environment is stationary. We propose to build the notion of continual learning (CL) into wireless system design, so that the learning model can incrementally adapt to the new episodes, {\it without forgetting} knowledge learned from the previous episodes. Our design is based on a novel bilevel optimization formulation which ensures certain ``fairness" across different data samples. We demonstrate the effectiveness of the CL approach by integrating it with two popular DNN based models for power control and beamforming, respectively, and testing using both synthetic and ray-tracing based data sets. These numerical results show that the proposed CL approach is not only able to adapt to the new scenarios quickly and seamlessly, but importantly, it also maintains high performance over the previously encountered scenarios as well.


title: Episodic Learning for Safe Bipedal Locomotion with Control Barrier Functions and Projection-to-State Safety

abstract: This paper combines episodic learning and control barrier functions in the setting of bipedal locomotion. The safety guarantees that control barrier functions provide are only valid with perfect model knowledge; however, this assumption cannot be met on hardware platforms. To address this, we utilize the notion of projection-to-state safety paired with a machine learning framework in an attempt to learn the model uncertainty as it affects the barrier functions. The proposed approach is demonstrated both in simulation and on hardware for the AMBER-3M bipedal robot in the context of the stepping-stone problem, which requires precise foot placement while walking dynamically.


title: Model-Free Incremental Adaptive Dynamic Programming Based Approximate Robust Optimal Regulation

abstract: This paper presents a new formulation for model-free robust optimal regulation of continuous-time nonlinear systems. The proposed reinforcement learning based approach, referred to as incremental adaptive dynamic programming (IADP), exploits measured data to allow the design of the approximate optimal incremental control strategy, which stabilizes the controlled system incrementally under model uncertainties, environmental disturbances, and input saturation. By leveraging the time delay estimation (TDE) technique, we first exploit sensory data to reduce the requirement of a complete dynamics, where measured data are adopted to construct an incremental dynamics that reflects the system evolution in an incremental form. Then, the resulting incremental dynamics serves to design the approximate optimal incremental control strategy based on adaptive dynamic programming, which is implemented as a simplified single critic structure to get the approximate solution to the value function of the Hamilton-Jacobi-Bellman equation. Furthermore, for the critic artificial neural network, experience data are used to design an off-policy weight update law with guaranteed weight convergence. Rather importantly, to address the unintentionally introduced TDE error, we incorporate a TDE error bound related term into the cost function, whereby the TDE error is attenuated during the optimization process. The system stability proof and the weight convergence proof are provided. Numerical simulations are conducted to validate the effectiveness and superiority of our proposed IADP, especially regarding the reduced control energy expenditure and the enhanced robustness.


title: Determining 4-edge-connected components in linear time

abstract: In this work, we present the first linear time deterministic algorithm computing the 4-edge-connected components of an undirected graph. First, we show an algorithm listing all 3-edge-cuts in a given 3-edge-connected graph, and then we use the output of this algorithm in order to determine the 4-edge-connected components of the graph.