title: Maxwell Parallel Imaging

abstract: Purpose: To develop a general framework for Parallel Imaging (PI) with the use of Maxwell regularization for the estimation of the sensitivity maps (SMs) and constrained optimization for the parameter-free image reconstruction. Theory and Methods: Certain characteristics of both the SMs and the images are routinely used to regularize the otherwise ill-posed optimization-based joint reconstruction from highly accelerated PI data. In this paper we rely on a fundamental property of SMs--they are solutions of Maxwell equations-- we construct the subspace of all possible SM distributions supported in a given field-of-view, and we promote solutions of SMs that belong in this subspace. In addition, we propose a constrained optimization scheme for the image reconstruction, as a second step, once an accurate estimation of the SMs is available. The resulting method, dubbed Maxwell Parallel Imaging (MPI), works seamlessly for arbitrary sequences (both 2D and 3D) with any trajectory and minimal calibration signals. Results: The effectiveness of MPI is illustrated for a wide range of datasets with various undersampling schemes, including radial, variable-density Poisson-disc, and Cartesian, and is compared against the state-of-the-art PI methods. Finally, we include some numerical experiments that demonstrate the memory footprint reduction of the constructed Maxwell basis with the help of tensor decomposition, thus allowing the use of MPI for full 3D image reconstructions. Conclusions: The MPI framework provides a physics-inspired optimization method for the accurate and efficient image reconstruction from arbitrary accelerated scans.


title: Considerations, Good Practices, Risks and Pitfalls in Developing AI Solutions Against COVID-19

abstract: The COVID-19 pandemic has been a major challenge to humanity, with 12.7 million confirmed cases as of July 13th, 2020 [1]. In previous work, we described how Artificial Intelligence can be used to tackle the pandemic with applications at the molecular, clinical, and societal scales [2]. In the present follow-up article, we review these three research directions, and assess the level of maturity and feasibility of the approaches used, as well as their potential for operationalization. We also summarize some commonly encountered risks and practical pitfalls, as well as guidelines and best practices for formulating and deploying AI applications at different scales.


title: Modelling multi-period carbon markets using singular forward backward SDEs

abstract: We introduce a model for the evolution of emissions and the price of emissions allowances in a carbon market such as the European Union Emissions Trading System (EU ETS). The model accounts for multiple trading periods, or phases, with multiple times at which compliance can occur. At the end of each trading period, the participating firms must surrender allowances for their emissions made during that period, and additional allowances can be used for compliance in the following periods. We show that the multi-period allowance pricing problem is well-posed for various mechanisms (such as banking, borrowing and withdrawal of allowances) linking the trading periods. The results are based on the analysis of a forward-backward stochastic differential equation with coupled forward and backward components, a discontinuous terminal condition and a forward component that is degenerate. We also introduce an infinite period model, for a carbon market with a sequence of compliance times and with no end date. We show that, under appropriate conditions, the value function for the multi-period pricing problem converges, as the number of periods increases, to a value function for this infinite period model, and that such functions are unique.


title: Modeling flexoelectricity in soft dielectrics at finite deformation

abstract: This paper develops the equilibrium equations describing the flexoelectric effect in soft dielectrics under large deformations. Previous works have developed related theories using a flexoelectric coupling tensor of mixed material-spatial character. Here, we formulate the model in terms of a flexoelectric tensor completely defined in the material frame, with the same symmetries of the small-strain flexocoupling tensor and leading naturally to objective flexoelectric polarization fields. The energy potential and equilibrium equations are first expressed in terms of deformation and polarization, and then rewritten in terms of deformation and electric potential, yielding an unconstrained system of fourth order partial differential equations (PDEs). We further develop a theory of geometrically nonlinear extensible flexoelectric rods under open and closed circuit conditions, with which we examine analytically cantilever bending and buckling under mechanical and electrical actuation. Besides being a simple and explicit model pertinent to slender structures, this rod theory also allows us to test our general theory and its numerical implementation using B-splines. This numerical implementation is robust as it handles the electromechanical instabilities in soft flexoelectric materials.


title: Aharonov-Bohm effect for bound states in the cosmic string spacetime in the context of rainbow gravity

abstract: We explore a scenario of the general relativity determined by the framework of the rainbow gravity with the purpose of searching for analogues of the Aharonov-Bohm effect. By focusing on the confinement of the Dirac field and the scalar field to a hard-wall confining potential in a modified background of the cosmic string spacetime, we examine the effects of the rainbow gravity and the topology of the cosmic string spacetime. Then, we compare each spectrum of energy with the cases where the rainbow gravity is absent.


title: Pose2Mesh: Graph Convolutional Network for 3D Human Pose and Mesh Recovery from a 2D Human Pose

abstract: Most of the recent deep learning-based 3D human pose and mesh estimation methods regress the pose and shape parameters of human mesh models, such as SMPL and MANO, from an input image. The first weakness of these methods is an appearance domain gap problem, due to different image appearance between train data from controlled environments, such as a laboratory, and test data from in-the-wild environments. The second weakness is that the estimation of the pose parameters is quite challenging owing to the representation issues of 3D rotations. To overcome the above weaknesses, we propose Pose2Mesh, a novel graph convolutional neural network (GraphCNN)-based system that estimates the 3D coordinates of human mesh vertices directly from the 2D human pose. The 2D human pose as input provides essential human body articulation information, while having a relatively homogeneous geometric property between the two domains. Also, the proposed system avoids the representation issues, while fully exploiting the mesh topology using a GraphCNN in a coarse-to-fine manner. We show that our Pose2Mesh outperforms the previous 3D human pose and mesh estimation methods on various benchmark datasets. For the codes, see https://github.com/hongsukchoi/Pose2Mesh_RELEASE.


title: A Survey of Coded Distributed Computing

abstract: Distributed computing has become a common approach for large-scale computation of tasks due to benefits such as high reliability, scalability, computation speed, and costeffectiveness. However, distributed computing faces critical issues related to communication load and straggler effects. In particular, computing nodes need to exchange intermediate results with each other in order to calculate the final result, and this significantly increases communication overheads. Furthermore, a distributed computing network may include straggling nodes that run intermittently slower. This results in a longer overall time needed to execute the computation tasks, thereby limiting the performance of distributed computing. To address these issues, coded distributed computing (CDC), i.e., a combination of coding theoretic techniques and distributed computing, has been recently proposed as a promising solution. Coding theoretic techniques have proved effective in WiFi and cellular systems to deal with channel noise. Therefore, CDC may significantly reduce communication load, alleviate the effects of stragglers, provide fault-tolerance, privacy and security. In this survey, we first introduce the fundamentals of CDC, followed by basic CDC schemes. Then, we review and analyze a number of CDC approaches proposed to reduce the communication costs, mitigate the straggler effects, and guarantee privacy and security. Furthermore, we present and discuss applications of CDC in modern computer networks. Finally, we highlight important challenges and promising research directions related to CDC


title: Discovering Useful Sentence Representations from Large Pretrained Language Models

abstract: Despite the extensive success of pretrained language models as encoders for building NLP systems, they haven't seen prominence as decoders for sequence generation tasks. We explore the question of whether these models can be adapted to be used as universal decoders. To be considered "universal," a decoder must have an implicit representation for any target sentence $s$, such that it can recover that sentence exactly when conditioned on its representation. For large transformer-based language models trained on vast amounts of English text, we investigate whether such representations can be easily discovered using standard optimization methods. We present and compare three representation injection techniques for transformer-based models and three accompanying methods which map sentences to and from this representation space. Experiments show that not only do representations exist for sentences from a variety of genres. More importantly, without needing complex optimization algorithms, our methods recover these sentences almost perfectly without fine-tuning the underlying language model at all.


title: Markov Chain-Based Stochastic Strategies for Robotic Surveillance

abstract: This article surveys recent advancements of strategy designs for persistent robotic surveillance tasks with the focus on stochastic approaches. The problem describes how mobile robots stochastically patrol a graph in an efficient way where the efficiency is defined with respect to relevant underlying performance metrics. We first start by reviewing the basics of Markov chains, which is the primary motion model for stochastic robotic surveillance. Then two main criteria regarding the speed and unpredictability of surveillance strategies are discussed. The central objects that appear throughout the treatment is the hitting times of Markov chains, their distributions and expectations. We formulate various optimization problems based on the concerned metrics in different scenarios and establish their respective properties.


title: Graphs whose Kronecker covers are bipartite Kneser graphs

abstract: We show that there are $k$ simple graphs whose Kronecker covers are isomorphic to the bipartite Kneser graph $H(n,k)$, and that their chromatic numbers coincide with $\chi(K(n,k)) = n - 2k + 2$. We also determine the automorphism groups of these graphs.


title: On transversality of bent hyperplane arrangements and the topological expressiveness of ReLU neural networks

abstract: Let F:R^n -> R be a feedforward ReLU neural network. It is well-known that for any choice of parameters, F is continuous and piecewise (affine) linear. We lay some foundations for a systematic investigation of how the architecture of F impacts the geometry and topology of its possible decision regions for binary classification tasks. Following the classical progression for smooth functions in differential topology, we first define the notion of a generic, transversal ReLU neural network and show that almost all ReLU networks are generic and transversal. We then define a partially-oriented linear 1-complex in the domain of F and identify properties of this complex that yield an obstruction to the existence of bounded connected components of a decision region. We use this obstruction to prove that a decision region of a generic, transversal ReLU network F: R^n -> R with a single hidden layer of dimension (n + 1) can have no more than one bounded connected component.


title: Fresh perspective on gauging the conformal group

abstract: We consider the construction of gauge theories of gravity that are invariant under local conformal transformations. We first clarify the geometric nature of global conformal transformations, in both their infinitesimal and finite forms, and the consequences of global conformal invariance for field theories, before reconsidering existing approaches for gauging the conformal group, namely auxiliary conformal gauge theory and biconformal gauge theory, neither of which is generally accepted as a complete solution. We then demonstrate that, provided any matter fields belong to an irreducible representation of the Lorentz group, the recently proposed extended Weyl gauge theory (eWGT) may be considered as an alternative method for gauging the conformal group, since eWGT is invariant under the full set of local conformal transformations, including inversions, as well as possessing conservation laws that provide a natural local generalisation of those satisfied by field theories with global conformal invariance, and also having an `ungauged' limit that corresponds to global conformal transformations. By contrast, although standard Weyl gauge theory also enjoys the first of these properties, it does not share the other two, and so cannot be considered a valid gauge theory of the conformal group.


title: Ultra-thin corrugated metamaterial film as large-area transmission dynode

abstract: Large-area transmission dynodes were fabricated by depositing an ultra-thin continuous film on a silicon wafer with a 3-dimensional pattern. After removing the silicon, a corrugated membrane with enhanced mechanical properties was formed. Mechanical materials, such as this corrugated membrane, are engineered to improve its strength and robustness, which allows it to span a larger surface in comparison to flat membranes while the film thickness remains constant. The ultra-thin film consists of three layers (Al$_2$O$_3$ /TiN/Al$_2$O$_3$) and is deposited by atomic layer deposition (ALD). The encapsulated TiN layer provides in-plane conductivity, which is needed to sustain secondary electron emission. Two types of corrugated membranes were fabricated: a hexagonal honeycomb and an octagonal pattern. The latter was designed to match the square pitch of a CMOS pixel chip. The transmission secondary electron yield was determined with a collector-based method using a scanning electron microscope. The highest transmission electron yield was measured on a membrane with an octagonal pattern. A yield of 2.15 was achieved for 3.15 keV incident electrons for an Al$_2$O$_3$ /TiN/Al$_2$O$_3$ tri-layer film with layer thicknesses of 10/5/15 nm. The variation in yield across the surface of the corrugated membrane was determined by constructing a yield map. The active surface for transmission secondary electron emission is near 100%, i.e. a primary electron generates transmission secondary electrons regardless of the point of impact on the corrugated membrane.


title: An Optimal Hybrid Variance-Reduced Algorithm for Stochastic Composite Nonconvex Optimization

abstract: In this note we propose a new variant of the hybrid variance-reduced proximal gradient method in [7] to solve a common stochastic composite nonconvex optimization problem under standard assumptions. We simply replace the independent unbiased estimator in our hybrid- SARAH estimator introduced in [7] by the stochastic gradient evaluated at the same sample, leading to the identical momentum-SARAH estimator introduced in [2]. This allows us to save one stochastic gradient per iteration compared to [7], and only requires two samples per iteration. Our algorithm is very simple and achieves optimal stochastic oracle complexity bound in terms of stochastic gradient evaluations (up to a constant factor). Our analysis is essentially inspired by [7], but we do not use two different step-sizes.


title: Dislocation impediment by the grain boundaries in polycrystals

abstract: Thermodynamic dislocation theory incorporating dislocation impediment by the grain boundaries is developed to analyze the shear test of polycrystals. With a small set of physics based material parameters, we are able to simulate the stress-strain curves for the load and its reversal, which are consistent with the experimental curves of Thuillier and Manach [2009]. Representative distributions of plastic slip under load and its reversal are presented, and their evolution explains the extended length of the transition stage during load reversal.


title: The Use of Near-Axis Magnetic Fields for Stellarator Turbulence Simulations

abstract: The design of turbulence optimized stellarators has so far relied on three-dimensional equilibrium codes such as VMEC in order to find the minimum of a given objective function. In this work, we propose a complimentary approach based on the near-axis expansion to compute the geometry parameters of neoclassicaly optimized stellarators used in turbulence studies. As shown here, the near-axis expansion can be a reasonable approximation of the geometric parameters relevant for turbulence and stability simulations of the core of existing optimized stellarator designs. In particular, we examine the geometry coefficients that appear in the gyrokinetic equation, the drift-reduced fluid equations and the ideal ballooning equation. This approach may allow for the development of new stellarator optimization techniques significantly faster than conventional methods.


title: How carrier memory enters the Haus master equation of mode-locking

abstract: We present a generalization of the Haus master equation in which a dynamical boundary condition allows to describe complex pulse trains such as the Q-switched and harmonic transitions of passive mode-locking as well as the weak interactions between localized states. As an example, we investigate the influence of group velocity dispersion on the stability boundaries of the Q-switched regime. We compare our results with that of a time-delayed system.


title: Single-file diffusion in a bi-stable potential: signatures of memory in the barrier-crossing of a tagged-particle

abstract: We investigate memory effects in barrier-crossing in the overdamped setting. We focus on the scenario where the hidden degrees of freedom relax on exactly the same time scale as the observable. As a prototypical model, we analyze tagged-particle diffusion in a single file confined to a bi-stable potential. We identify the signatures of memory and explain their origin. The emerging memory is a result of the projection of collective many-body eigenmodes onto the motion of a tagged-particle. We are interested in the "confining" (all background particles in front of the tagged-particle) and "pushing" (all background particles behind the tagged-particle) scenarios for which we find non-trivial and qualitatively different relaxation behaviors. Notably and somewhat unexpectedly, at a fixed particle number, we find that the higher the barrier, the stronger the memory effects are. The fact that the external potential alters the memory is important more generally and should be taken into account in applications of generalized Langevin equations. Our results can readily be tested experimentally and may be relevant for understanding transport in biological ion-channels.


title: Bounding ramification by covers and curves

abstract: We prove that $\bar {\mathbb Q}_\ell$-local systems of bounded rank and ramification on a smooth variety $X$ defined over an algebraically closed field $k$ of characteristic $p\neq \ell$ are tamified outside of codimension $2$ by a finite separable cover of bounded degree. In rank one, there is a curve which preserves their monodromy. There is a curve defined over the algebraic closure of a purely transcendental extension of $k$ of finite degree which fulfills the Lefschetz theorem. Last version: minor typos corrected.


title: Analysis of Multivariate Scoring Functions for Automatic Unbiased Learning to Rank

abstract: Leveraging biased click data for optimizing learning to rank systems has been a popular approach in information retrieval. Because click data is often noisy and biased, a variety of methods have been proposed to construct unbiased learning to rank (ULTR) algorithms for the learning of unbiased ranking models. Among them, automatic unbiased learning to rank (AutoULTR) algorithms that jointly learn user bias models (i.e., propensity models) with unbiased rankers have received a lot of attention due to their superior performance and low deployment cost in practice. Despite their differences in theories and algorithm design, existing studies on ULTR usually use uni-variate ranking functions to score each document or result independently. On the other hand, recent advances in context-aware learning-to-rank models have shown that multivariate scoring functions, which read multiple documents together and predict their ranking scores jointly, are more powerful than uni-variate ranking functions in ranking tasks with human-annotated relevance labels. Whether such superior performance would hold in ULTR with noisy data, however, is mostly unknown. In this paper, we investigate existing multivariate scoring functions and AutoULTR algorithms in theory and prove that permutation invariance is a crucial factor that determines whether a context-aware learning-to-rank model could be applied to existing AutoULTR framework. Our experiments with synthetic clicks on two large-scale benchmark datasets show that AutoULTR models with permutation-invariant multivariate scoring functions significantly outperform those with uni-variate scoring functions and permutation-variant multivariate scoring functions.


title: Monocular Expressive Body Regression through Body-Driven Attention

abstract: To understand how people look, interact, or perform tasks, we need to quickly and accurately capture their 3D body, face, and hands together from an RGB image. Most existing methods focus only on parts of the body. A few recent approaches reconstruct full expressive 3D humans from images using 3D body models that include the face and hands. These methods are optimization-based and thus slow, prone to local optima, and require 2D keypoints as input. We address these limitations by introducing ExPose (EXpressive POse and Shape rEgression), which directly regresses the body, face, and hands, in SMPL-X format, from an RGB image. This is a hard problem due to the high dimensionality of the body and the lack of expressive training data. Additionally, hands and faces are much smaller than the body, occupying very few image pixels. This makes hand and face estimation hard when body images are downscaled for neural networks. We make three main contributions. First, we account for the lack of training data by curating a dataset of SMPL-X fits on in-the-wild images. Second, we observe that body estimation localizes the face and hands reasonably well. We introduce body-driven attention for face and hand regions in the original image to extract higher-resolution crops that are fed to dedicated refinement modules. Third, these modules exploit part-specific knowledge from existing face- and hand-only datasets. ExPose estimates expressive 3D humans more accurately than existing optimization methods at a small fraction of the computational cost. Our data, model and code are available for research at https://expose.is.tue.mpg.de .


title: A review on the mobile applications developed for COVID-19: An exploratory analysis

abstract: The objective of this research is to explore the existing mobile applications developed for the COVID-19 pandemic. To obtain this research objective, firstly the related applications were selected through the systematic search technique in the popular application stores. Secondly, data related to the app objectives, functionalities provided by the app, user ratings, and user reviews were extracted. Thirdly, the extracted data were analyzed through the affinity diagram, noticing-collecting-thinking, and descriptive analysis. As outcomes, the review provides a state-of-the-art view of mobile apps developed for COVID-19 by revealing nine functionalities or features. It revealed ten factors related to information systems design characteristics that can guide future app design. The review outcome highlights the need for new development and further refinement of the existing applications considering not only the revealed objectives and their associated functionalities, but also revealed design characteristics such as reliability, performance, usefulness, supportive, security, privacy, flexibility, responsiveness, ease of use, and cultural sensitivity.


title: Enhanced CP asymmetries in $B \to K \mu^+ \mu^-$

abstract: We show that the current values of $R_K^\mathrm{exp}$ and $R_{K^\ast}^\mathrm{exp}$ can be accommodated by allowing a nonzero New Physics coupling $\delta C_9^{\mu\mu}$ to be complex, both in the scenario in which only $\delta C_9^{\mu\mu}$ is affected, and in the scenario with complex $\delta C_{9,10}^{\mu\mu}$ satisfying $\delta C_{9}^{\mu\mu}=-\delta C_{10}^{\mu\mu}$. A presence of the weak CP-violating phase can then be tested by measuring the CP-asymmetry, $\mathcal{A}_\mathrm{CP}$. We show that this asymmetry is enhanced around the peak of each $c\bar c$-resonance, and in fact more pronounced in the close vicinity of $J/\psi$ and $\psi(2S)$. Therefore, measuring $\mathcal{A}_\mathrm{CP}$ before and after the resonances' peak could be revelatory of the CP-violation that originates from beyond the Standard Model, or to be a significant constrain when building a realistic scenario of New Physics.


title: Channels, measurements and post-selection in quantum thermodynamics

abstract: We analyse the benefit, in terms of extracting work, of having a single use of a quantum channel or measurement in quantum thermodynamics. This highlights a connection between unital and catalytic channels, and some subtleties concerning the conditional work cost of implementing a measurement given that a certain result was obtained. We also consider post-selected measurements, and show that any non-trivial post-selection leads to an unbounded work benefit.


title: Investigation of the rotation-activity relation for a sample of SuperWASP and ASAS-SN field stars

abstract: It is well established that late-type main-sequence stars display a relationship between X-ray activity and the Rossby number, $Ro$, the ratio of rotation period to the convective turnover time. This manifests itself as a saturated regime (where X-ray activity is constant) and an unsaturated regime (where X-ray activity anti-correlates with the Rossby number). However, this relationship breaks down for the fastest rotators. We cross-correlated SuperWASP visually classified photometric light curves and ASAS-SN automatically classified photometric light curves with XMM-Newton X-ray observations to identify 3,178 stars displaying a photometrically defined rotational modulation in their light curve and corresponding X-ray observations. We fitted a power-law to characterise the rotation-activity relation of 900 main-sequence stars. We identified that automatically classified rotationally modulated light curves are not as reliable as visually classified light curves for this work. We found a power-law index in the unsaturated regime of G- to M-type stars of $\beta=-1.84\pm0.18$ for the SuperWASP catalogue, in line with the canonical value of $\beta=-2$. We find evidence of supersaturation in the fastest rotating K-type stars, with a power-law index of $\beta_{s}=1.42\pm0.26$.


title: Adventures in Mathematical Reasoning

abstract: "Mathematics is not a careful march down a well-cleared highway, but a journey into a strange wilderness, where the explorers often get lost. Rigour should be a signal to the historian that the maps have been made, and the real explorers have gone elsewhere." W.S. Anglin, the Mathematical Intelligencer, 4 (4), 1982.


title: Fractional flow equations. A model for pressure deficit in an oil well

abstract: This article presents a novel system of flow equations that models the pressure deficit of a reservoir considered as a triple continuous medium formed by the rock matrix, vugular medium and fracture. In non-conventional reservoirs, the velocity of the fluid particles is altered due to physical and chemical phenomena caused by the interaction of the fluid with the medium, this behavior is defined as anomalous. A more exact model can be obtained with the inclusion of the memory formalism concept that can be expressed through the use of fractional derivatives. Using Laplace transform of the Caputo fractional derivative and Bessel functions, a semi-analytical solution is reached in the Laplace space.

categories: math.AP math.CA math.FA physics.app-ph physics.comp-ph


title: On algebras associated with invariant means on the subnormal subgroups of an amenable group

abstract: Let $G$ be an amenable group. We define and study an algebra $\mathcal{A}_{sn}(G)$, which is related to invariant means on the subnormal subgroups of $G$. For a just infinite amenable group $G$, we show that $\mathcal{A}_{sn}(G)$ is nilpotent if and only if $G$ is not a branch group, and in the case that it is nilpotent we determine the index of nilpotence. We next study $\operatorname{rad} \ell^1(G)^{**}$ for an amenable branch group $G$, and show that it always contains nilpotent left ideals of arbitrarily large index, as well as non-nilpotent elements. This provides infinitely many finitely-generated counterexamples to a question of Dales and Lau, first resolved by the author in a previous article, which asks whether we always have $(\operatorname{rad} \ell^1(G)^{**})^{\Box 2} = \{ 0 \}$. We further study this question by showing that $(\operatorname{rad} \ell^1(G)^{**})^{\Box 2} = \{ 0 \}$ imposes certain structural constraints on the group $G$.


title: Effect of spatial distribution of mesoscale heterogeneities on the shock-to-detonation transition in liquid nitromethane

abstract: The sensitizing effect of cavities in the form of microbubbles on the shock initiation of a homogeneous liquid explosive is studied computationally. While the presence of voids in an explosive has long been known to induce so-called hot spots that greatly accelerate the global reaction rate, the ability to computationally resolve the details of the interaction of the shock front with heterogeneities existing on the scale of the detonation reaction zone has only recently become feasible. In this study, the influence of the spatial distribution of air-filled cavities has been examined, enabled by the use of graphic processing unit (GPU) accelerated computations that can resolve shock initiation and detonation propagation through an explosive while fully resolving features at the mesoscale. Different spatial distributions of cavities are examined in two-dimensional simulations, including regular arrays of cavities, slightly perturbed arrays, random arrays (with varying minimum spacing being imposed on the cavities), and randomly distributed clusters of cavities. The presence of the cavities is able to reduce the time required to initiate detonation---for a given input shock strength---by greater than 50%, in agreement with previous experimental results. Randomly distributing the cavities results in a 15-20% decrease in detonation initiation time in comparison to a regular array of cavities. Clustering the cavities---as would occur in the case of agglomeration---results in an additional 10% decrease in detonation initiation time in comparison to random arrays. The effect of clustering is shown not to be a result of the clusters forming an effectively larger cavity, but rather due to interactions between clusters upon shock loading occurring on the microscale.


title: Implementation of model predictive control for tracking in embedded systems using a sparse extended ADMM algorithm

abstract: This article presents a sparse, low-memory footprint optimization algorithm for the implementation of the model predictive control (MPC) for tracking formulation in embedded systems. This MPC formulation has several advantages over standard MPC formulations, such as an increased domain of attraction and guaranteed recursive feasibility even in the event of a sudden reference change. However, this comes at the expense of the addition of a small amount of decision variables to the MPC's optimization problem that complicates the structure of its matrices. We propose a sparse optimization algorithm, based on an extension of the alternating direction method of multipliers, that exploits the structure of this particular MPC formulation. We describe the controller formulation and detail how its structure is exploited by means of the aforementioned optimization algorithm. We show closed-loop simulations comparing the proposed solver against other solvers and approaches from the literature.


title: Utilizing Explainable AI for Quantization and Pruning of Deep Neural Networks

abstract: For many applications, utilizing DNNs (Deep Neural Networks) requires their implementation on a target architecture in an optimized manner concerning energy consumption, memory requirement, throughput, etc. DNN compression is used to reduce the memory footprint and complexity of a DNN before its deployment on hardware. Recent efforts to understand and explain AI (Artificial Intelligence) methods have led to a new research area, termed as explainable AI. Explainable AI methods allow us to understand better the inner working of DNNs, such as the importance of different neurons and features. The concepts from explainable AI provide an opportunity to improve DNN compression methods such as quantization and pruning in several ways that have not been sufficiently explored so far. In this paper, we utilize explainable AI methods: mainly DeepLIFT method. We use these methods for (1) pruning of DNNs; this includes structured and unstructured pruning of \ac{CNN} filters pruning as well as pruning weights of fully connected layers, (2) non-uniform quantization of DNN weights using clustering algorithm; this is also referred to as Weight Sharing, and (3) integer-based mixed-precision quantization; this is where each layer of a DNN may use a different number of integer bits. We use typical image classification datasets with common deep learning image classification models for evaluation. In all these three cases, we demonstrate significant improvements as well as new insights and opportunities from the use of explainable AI in DNN compression.


title: On the weak scaling of the contact distance between two fluctuating interfaces with system size

abstract: A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers of mass remain separated by a nonzero distance $\ell$. Here we examine the statistics of $\ell$ at the time of first contact for surfaces that evolve in time according to the Edwards-Wilkinson equation. We present a general approach to calculate its probability distribution and determine how its most likely value $\ell^*$ depends on the surfaces' lateral size $L$. We are motivated by an interest in the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Computer simulations of this scenario verify the predicted scaling behavior in two and three dimensions. In the latter case, slow growth where $\ell^\ast$ is an algebraic function of $\log L$ implies that slab-shaped domains remain topologically intact until $\ell$ becomes very small, contradicting expectations from equilibrium thermodynamics.


title: Flight Software Development for the EIRSAT-1 Mission

abstract: The Educational Irish Research Satellite, known as EIRSAT-1, is a student-led project to design, build, test and launch Ireland's first satellite. The on-board software for this mission is being developed using Bright Ascension's GenerationOne Flight Software Development Kit. This paper provides an overview of this kit and of EIRSAT-1's on-board software design. Drawing on the team's contrasting experience with writing entirely custom firmware for the mission's science payloads, this work discusses the impact of using a kit on the software development process. The challenges associated with the educational nature of this project are the focus of this discussion. The objective of this paper is to provide useful information for other CubeSat teams assessing software development options.


title: Controlling Dialogue Generation with Semantic Exemplars

abstract: Dialogue systems pretrained with large language models generate locally coherent responses, but lack the fine-grained control over responses necessary to achieve specific goals. A promising method to control response generation is exemplar-based generation, in which models edit exemplar responses that are retrieved from training data, or hand-written to strategically address discourse-level goals, to fit new dialogue contexts. But, current exemplar-based approaches often excessively copy words from the exemplar responses, leading to incoherent replies. We present an Exemplar-based Dialogue Generation model, EDGE, that uses the semantic frames present in exemplar responses to guide generation. We show that controlling dialogue generation based on the semantic frames of exemplars, rather than words in the exemplar itself, improves the coherence of generated responses, while preserving semantic meaning and conversation goals present in exemplar responses.


title: Topological Indices Over Nonzero Component Graph of a Finite Dimensional Vector Space

abstract: The study of graphs associated with of various algebraic structures is an emerging topic in algebraic graph theory. Recently, the concept of nonzero component graph of a finite dimensional vector space $\Gamma(\mathbb{V})$ was put forward by Das \cite{5}. In this paper, we study some degree based topological indices over $\Gamma(\mathbb{V})$ the derived graphs of $\Gamma(\mathbb{V})$.


title: Sample-efficient benchmarking of multi-photon interference on a boson sampler in the sparse regime

abstract: Verification of a quantum advantage in the presence of noise is a key open problem in the study of near-term quantum devices. In this work, we show how to assess the quality of photonic interference in a linear optical quantum device (boson sampler) by using a maximum likelihood method to measure the strength at which various noise sources are present in the experiment. This allows us to use a sparse set of samples to test whether a given boson sampling experiment meets known upper bounds on the level of noise permissible to demonstrate a quantum advantage. Furthermore, this method allows us monitor the evolution of noise in real time, creating a valuable diagnostic tool. Finally, we observe that sources of noise in the experiment compound, meaning that the observed value of the mutual photon indistinguishability, which is the main imperfection in our study, is an effective value taking into account all sources of error in the experiment.


title: Testing the isomorph invariance of the bridge functions of Yukawa one-component plasmas. I. Intermediate and long range

abstract: It has been recently conjectured that bridge functions remain nearly invariant along phase diagram lines of constant excess entropy for the broad class of R-simple liquids. To test this hypothesis, the bridge functions of Yukawa systems are computed outside the correlation void with the Ornstein-Zernike inversion method and structural input from ultra-accurate molecular dynamics simulations. The effect of statistical, grid, finite-size, tail and isomorphic errors is quantified. Uncertainty propagation analysis is complemented with a detailed investigation of the sensitivity of the bridge function to periodic and aperiodic multiplicative perturbations in the radial distribution function. In the long and intermediate range, bridge functions are demonstrated to be approximately isomorph invariant.


title: Pure State Tomography with Fourier Transformation

abstract: Extracting information from quantum devices has long been a crucial problem in the field of quantum mechanics. By performing elaborate measurements, quantum state tomography, an important and fundamental tool in quantum science and technology, can be used to determine unknown quantum states completely. In this study, we explore methods to determine multi-qubit pure quantum states uniquely and directly. Two adaptive protocols are proposed, with their respective quantum circuits. Herein, two or three observables are sufficient, while the number of measurement outcomes is either the same as or fewer than those in existing methods. Additionally, experiments on the IBM 5-qubit quantum computer, as well as numerical investigations, demonstrate the feasibility of the proposed protocols.


title: First-principles studies of the electronic and magnetic structures and bonding properties of boron subnitride B$_{13}$N$_2$

abstract: Rhombohedral B$_{12}$ unit is viewed as a host matrix embedding linear tri-atomic arrangements of elements (E) resulting in a relatively large family of boron-rich compounds with B$_{12}${E-E-E} generic formulation. The present work focuses on boron subnitride, B$_{13}$N$_2$ that we express in present context as B$_{12}${N-B-N}. Within well established quantum density functional theory (DFT) a full study of its electronic properties is provided. Also linear triatomic arrangements in view of the existence in simple compounds such as sodium azide NaN$_3$, i.e., Na$^I${N-N-N} and calcium cyanamide, Ca$^{II}${N-C-N}, we devised Sc$^{III}${N-B-N} to establish comparison with B$_{12}${N-B-N}. ScBN$_2$ is calculated to be cohesive and possessing N-B-N isolated from ScIII with dB-N = 1.33 {\AA}. In B$_{12}${N-B-N} an elongated dB-N=1.43 {\AA} is identified due to the bonding of N with one of the two B12 boron substructures, B1 with the formation of "3B...N-B-N...3B"-like complex accompanied by a magnetic instability. Spin polarized (SP) calculations led to the onset of magnetization on central boron with M=1 $\mu_B$ in a stable half-ferromagnetic ground state observed from the electronic density of states (DOS). The results are backed with total energy and calculations in both non-spin-polarized (NSP) and spin-polarized stabilizing the latter configuration over a broad range of volumes from M(V) plots. Further illustrative results are given with the charge densities (total and magnetic) and electron localization function (ELF).


title: A categorical approach to dynamical quantum groups

abstract: We present a categorical point of view on dynamical quantum groups in terms of categories of Harish-Chandra bimodules. We prove Tannaka duality theorems for forgetful functors into the monoidal category of Harish-Chandra bimodules in terms of a slight modification of the notion of a bialgebroid. Moreover, we show that the standard dynamical quantum groups $F(G)$ and $F_q(G)$ are related to parabolic restriction functors for classical and quantum Harish-Chandra bimodules. Finally, we exhibit a natural Weyl symmetry of the parabolic restriction functor using Zhelobenko operators and show that it gives rise to the action of the dynamical Weyl group.


title: Exploring possible triangle singularities in the $\Xi^-_{b} \to K^- J/\psi \Lambda$ decay

abstract: We analyze possible singularities in the $J/\psi \Lambda$ invariant mass distribution of the $\Xi^-_{b}~\to~K^- J/\psi \Lambda$ process via triangle loop diagrams. Triangle singularities in the physical region are found in 18 different triangle loop diagrams. Among those with $\Xi^*$-charmonium-$\Lambda$ intermediate states, the one from the $\chi_{c1} \Xi(2120) \Lambda$ loop, which is located around 4628 MeV, is found the most likely to cause observable effects. One needs $S$- and $P$-waves in $\chi_{c1} \Lambda$ and $J/\psi \Lambda$ systems, respectively, when the quantum numbers of these systems are $1/2^+$ or $3/2^+$. When the quantum numbers of the $\Xi(2120)$ are $J^P=1/2^+$, $1/2^-$ or $3/2^+$, the peak structure should be sharper than the other $J^P$ choices. This suggests that although the whole strength is unknown, we should pay attention to the contributions from the $\Xi^*$-charmonium-$\Lambda$ triangle diagram if structures are observed in the $J/\psi \Lambda$ invariant mass spectrum experimentally. In addition, a few triangle diagrams with the $D_{s1}^*(2700)$ as one of the intermediate particles can also produce singularities in the $J/\psi\Lambda$ distribution, but at higher energies above 4.9 GeV.


title: Exact Tests for Offline Changepoint Detection in Multichannel Binary and Count Data with Application to Networks

abstract: We consider offline detection of a single changepoint in binary and count time-series. We compare exact tests based on the cumulative sum (CUSUM) and the likelihood ratio (LR) statistics, and a new proposal that combines exact two-sample conditional tests with multiplicity correction, against standard asymptotic tests based on the Brownian bridge approximation to the CUSUM statistic. We see empirically that the exact tests are much more powerful in situations where normal approximations driving asymptotic tests are not trustworthy: (i) small sample settings; (ii) sparse parametric settings; (iii) time-series with changepoint near the boundary. We also consider a multichannel version of the problem, where channels can have different changepoints. Controlling the False Discovery Rate (FDR), we simultaneously detect changes in multiple channels. This "local" approach is shown to be more advantageous than multivariate global testing approaches when the number of channels with changepoints is much smaller than the total number of channels. As a natural application, we consider network-valued time-series and use our approach with (a) edges as binary channels and (b) node-degrees or other local subgraph statistics as count channels. The local testing approach is seen to be much more informative than global network changepoint algorithms.


title: Do Syntax Trees Help Pre-trained Transformers Extract Information?

abstract: Much recent work suggests that incorporating syntax information from dependency trees can improve task-specific transformer models. However, the effect of incorporating dependency tree information into pre-trained transformer models (e.g., BERT) remains unclear, especially given recent studies highlighting how these models implicitly encode syntax. In this work, we systematically study the utility of incorporating dependency trees into pre-trained transformers on three representative information extraction tasks: semantic role labeling (SRL), named entity recognition, and relation extraction. We propose and investigate two distinct strategies for incorporating dependency structure: a late fusion approach, which applies a graph neural network on the output of a transformer, and a joint fusion approach, which infuses syntax structure into the transformer attention layers. These strategies are representative of prior work, but we introduce additional model design elements that are necessary for obtaining improved performance. Our empirical analysis demonstrates that these syntax-infused transformers obtain state-of-the-art results on SRL and relation extraction tasks. However, our analysis also reveals a critical shortcoming of these models: we find that their performance gains are highly contingent on the availability of human-annotated dependency parses, which raises important questions regarding the viability of syntax-augmented transformers in real-world applications.


title: Conway and aperiodic tilings

abstract: This is a brief introduction to the geometric aspects of aperiodic tiling and the collaboration of John Conway and the author in the decade 1990-2000.


title: Scaling and local limits of Baxter permutations and bipolar orientations through coalescent-walk processes

abstract: Baxter permutations, plane bipolar orientations, and a specific family of walks in the non-negative quadrant, called tandem walks, are well-known to be related to each other through several bijections. We introduce a further new family of discrete objects, called coalescent-walk processes and we relate it to the three families mentioned above. We prove joint Benjamini--Schramm convergence (both in the annealed and quenched sense) for uniform objects in the four families. Furthermore, we explicitly construct a new random measure on the unit square, called the Baxter permuton and we show that it is the scaling limit (in the permuton sense) of uniform Baxter permutations. In addition, we relate the limiting objects of the four families to each other, both in the local and scaling limit case. The scaling limit result is based on the convergence of the trajectories of the coalescent-walk process to the coalescing flow -- in the terminology of Le Jan and Raimond (2004) -- of a perturbed version of the Tanaka stochastic differential equation. Our scaling result entails joint convergence of the tandem walks of a plane bipolar orientation and its dual, extending the main result of Gwynne, Holden, Sun (2016), and giving an alternative answer to Conjecture 4.4 of Kenyon, Miller, Sheffield, Wilson (2019) compared to the one of Gwynne, Holden, Sun (2016).


title: Elliptic curves with non-abelian entanglements

abstract: We consider the problem of classifying quadruples $(K,E,m_1,m_2)$ where $K$ is a number field, $E$ is an elliptic curve defined over $K$ and $(m_1,m_2)$ is a pair of relatively prime positive integers for which the intersection $K(E[m_1]) \cap K(E[m_2])$ is a non-abelian extension of $K$. There is an infinite set $\mathcal{S}$ of modular curves whose $K$-rational points capture all elliptic curves over $K$ without complex multiplication that have this property. Our main theorem explicitly describes the (finite) subset of $\mathcal{S}$ consisting of those modular curves having genus zero. In the case $K = \mathbb{Q}$, this has applications to the problem of determining when the Galois representation on the torsion of $E$ is as large as possible modulo a prescribed obstruction; we illustrate this application with a specific example.


title: DeepGMR: Learning Latent Gaussian Mixture Models for Registration

abstract: Point cloud registration is a fundamental problem in 3D computer vision, graphics and robotics. For the last few decades, existing registration algorithms have struggled in situations with large transformations, noise, and time constraints. In this paper, we introduce Deep Gaussian Mixture Registration (DeepGMR), the first learning-based registration method that explicitly leverages a probabilistic registration paradigm by formulating registration as the minimization of KL-divergence between two probability distributions modeled as mixtures of Gaussians. We design a neural network that extracts pose-invariant correspondences between raw point clouds and Gaussian Mixture Model (GMM) parameters and two differentiable compute blocks that recover the optimal transformation from matched GMM parameters. This construction allows the network learn an SE(3)-invariant feature space, producing a global registration method that is real-time, generalizable, and robust to noise. Across synthetic and real-world data, our proposed method shows favorable performance when compared with state-of-the-art geometry-based and learning-based registration methods.


title: Primal-dual evolutionary dynamics for constrained population games

abstract: Population games can be regarded as a tool to study the strategic interaction of a population of players. Although several attention has been given to such field, most of the available works have focused only on the unconstrained case. That is, the allowed equilibrium of the game is not constrained. To further extend the capabilities of population games, in this paper we propose a novel class of primal-dual evolutionary dynamics that allow the consideration of constraints that must be satisfied at the equilibrium of the game. Using duality theory and Lyapunov stability theory, we provide sufficient conditions to guarantee the asymptotic stability and feasibility of the equilibria set of the game under the considered constraints. Furthermore, we illustrate the application of the developed theory to some classical population games with the addition of constraints.


title: TRU-NET: A Deep Learning Approach to High Resolution Prediction of Rainfall

abstract: Climate models (CM) are used to evaluate the impact of climate change on the risk of floods and strong precipitation events. However, these numerical simulators have difficulties representing precipitation events accurately, mainly due to limited spatial resolution when simulating multi-scale dynamics in the atmosphere. To improve the prediction of high resolution precipitation we apply a Deep Learning (DL) approach using an input of CM simulations of the model fields (weather variables) that are more predictable than local precipitation. To this end, we present TRU-NET (Temporal Recurrent U-Net), an encoder-decoder model featuring a novel 2D cross attention mechanism between contiguous convolutional-recurrent layers to effectively model multi-scale spatio-temporal weather processes. We use a conditional-continuous loss function to capture the zero-skewed %extreme event patterns of rainfall. Experiments show that our model consistently attains lower RMSE and MAE scores than a DL model prevalent in short term precipitation prediction and improves upon the rainfall predictions of a state-of-the-art dynamical weather model. Moreover, by evaluating the performance of our model under various, training and testing, data formulation strategies, we show that there is enough data for our deep learning approach to output robust, high-quality results across seasons and varying regions.


title: Maximum likelihood estimation of parameters of spherical particle size distributions from profile size measurements and its application for small samples

abstract: Microscopy research often requires recovering particle-size distributions in three dimensions from only a few (10 - 200) profile measurements in the section. This problem is especially relevant for petrographic and mineralogical studies, where parametric assumptions are reasonable and finding distribution parameters from the microscopic study of small sections is essential. This paper deals with the specific case where particles are approximately spherical (i.e. Wicksell's problem). The paper presents a novel approximation of the probability density of spherical particle profile sizes. This approximation uses the actual non-smoothness of mineral particles rather than perfect spheres. The new approximation facilitates the numerically efficient use of the maximum likelihood method, a generally powerful method that provides the distribution parameter estimates of the minimal variance in most practical cases. The variance and bias of the estimates by the maximum likelihood method were compared numerically for several typical particle-size distributions with those by alternative parametric methods (method of moments and minimum distance estimation), and the maximum likelihood estimation was found to be preferable for both small and large samples. The maximum likelihood method, along with the suggested approximation, may also be used for selecting a model, for constructing narrow confidence intervals for distribution parameters using all the profiles without random sampling and for including the measurements of the profiles intersected by section boundaries. The utility of the approach is illustrated using an example from glacier ice petrography.


title: Meta-Sim2: Unsupervised Learning of Scene Structure for Synthetic Data Generation

abstract: Procedural models are being widely used to synthesize scenes for graphics, gaming, and to create (labeled) synthetic datasets for ML. In order to produce realistic and diverse scenes, a number of parameters governing the procedural models have to be carefully tuned by experts. These parameters control both the structure of scenes being generated (e.g. how many cars in the scene), as well as parameters which place objects in valid configurations. Meta-Sim aimed at automatically tuning parameters given a target collection of real images in an unsupervised way. In Meta-Sim2, we aim to learn the scene structure in addition to parameters, which is a challenging problem due to its discrete nature. Meta-Sim2 proceeds by learning to sequentially sample rule expansions from a given probabilistic scene grammar. Due to the discrete nature of the problem, we use Reinforcement Learning to train our model, and design a feature space divergence between our synthesized and target images that is key to successful training. Experiments on a real driving dataset show that, without any supervision, we can successfully learn to generate data that captures discrete structural statistics of objects, such as their frequency, in real images. We also show that this leads to downstream improvement in the performance of an object detector trained on our generated dataset as opposed to other baseline simulation methods. Project page: https://nv-tlabs.github.io/meta-sim-structure/.


title: PARADE: Passage Representation Aggregation for Document Reranking

abstract: Pretrained transformer models, such as BERT and T5, have shown to be highly effective at ad-hoc passage and document ranking. Due to inherent sequence length limits of these models, they need to be run over a document's passages, rather than processing the entire document sequence at once. Although several approaches for aggregating passage-level signals have been proposed, there has yet to be an extensive comparison of these techniques. In this work, we explore strategies for aggregating relevance signals from a document's passages into a final ranking score. We find that passage representation aggregation techniques can significantly improve over techniques proposed in prior work, such as taking the maximum passage score. We call this new approach PARADE. In particular, PARADE can significantly improve results on collections with broad information needs where relevance signals can be spread throughout the document (such as TREC Robust04 and GOV2). Meanwhile, less complex aggregation techniques may work better on collections with an information need that can often be pinpointed to a single passage (such as TREC DL and TREC Genomics). We also conduct efficiency analyses, and highlight several strategies for improving transformer-based aggregation.


title: Scruples: A Corpus of Community Ethical Judgments on 32,000 Real-Life Anecdotes

abstract: As AI systems become an increasing part of people's everyday lives, it becomes ever more important that they understand people's ethical norms. Motivated by descriptive ethics, a field of study that focuses on people's descriptive judgments rather than theoretical prescriptions on morality, we investigate a novel, data-driven approach to machine ethics. We introduce Scruples, the first large-scale dataset with 625,000 ethical judgments over 32,000 real-life anecdotes. Each anecdote recounts a complex ethical situation, often posing moral dilemmas, paired with a distribution of judgments contributed by the community members. Our dataset presents a major challenge to state-of-the-art neural language models, leaving significant room for improvement. However, when presented with simplified moral situations, the results are considerably more promising, suggesting that neural models can effectively learn simpler ethical building blocks. A key take-away of our empirical analysis is that norms are not always clean-cut; many situations are naturally divisive. We present a new method to estimate the best possible performance on such tasks with inherently diverse label distributions, and explore likelihood functions that separate intrinsic from model uncertainty.


title: A multiplicative Tate spectral sequence for compact Lie group actions

abstract: Given a compact Lie group $G$ and a commutative orthogonal ring spectrum $R$ such that $R[G]_* = \pi_*(R \wedge G_+)$ is finitely generated and projective over $\pi_*(R)$, we construct a multiplicative $G$-Tate spectral sequence for each $R$-module $X$ in orthogonal $G$-spectra, with $E^2$-page given by the Hopf algebra Tate cohomology of $R[G]_*$ with coefficients in $\pi_*(X)$. Under mild hypotheses, such as $X$ being bounded below and the derived page $RE^\infty$ vanishing, this spectral sequence converges strongly to the homotopy $\pi_*(X^{tG})$ of the $G$-Tate construction $X^{tG} = [\widetilde{EG} \wedge F(EG_+, X)]^G$.


title: Characterising the Gaia Radial Velocity sample selection function in its native photometry

abstract: The Gaia DR2 radial velocity sample (GDR2RVS), which provides six-dimensional phase-space information on 7.2 million stars, is of great value for inferring properties of the Milky Way. Yet a quantitative and accurate modelling of this sample is hindered without knowledge and inclusion of a well-characterized selection function. Here we derive the selection function through estimates of the internal completeness, i.e. the ratio of GDR2RVS sources compared to all Gaia DR2 sources (GDR2all). We show that this selection function or "completeness" depends on basic observables, in particular the apparent magnitude GRVS and colour G-GRP, but also on the surrounding source density and on sky position, where the completeness exhibits distinct small-scale structure. We identify a region of magnitude and colour that has high completeness, providing an approximate but simple way of implementing the selection function. For a more rigorous and detailed description we provide python code to query our selection function, as well as tools and ADQL queries that produce custom selection functions with additional quality cuts.


title: Tuning the magnetism and band topology through antisite defects in Sb doped MnBi4Te7

abstract: The fine control of magnetism and electronic structure is crucial since the interplay between magnetism and band topology can lead to various novel magnetic topological states including axion insulators, magnetic Weyl semimetals and Chern insulators etc. Through crystal growth, transport, thermodynamic, neutron diffraction measurements, we show that with Sb-doping, the newly-discovered intrinsic antiferromagnetic topological insulator MnBi4Te7 evolves from antiferro-magnetic to ferromagnetic and then ferrimagnetic. We attribute this to the formation of Mn(Bi,Sb) antisites upon doping, which result in additional Mn sublattices that modify the delicate interlayer magnetic interactions and cause the dominant Mn sublattice to go from antiferromagnetic to ferro-magnetic. We further investigate the effect of antisites on the band topology using the first-principles calculations. Without considering antisites, the series evolves from antiferromagnetic topological insulator (x = 0) to ferromagnetic axion insulators. In the exaggerated case of 16.7% of periodic antisites, the band topology is modified and type-I magnetic Weyl semimetal phase can be realized at intermediate dopings. Therefore, this doping series provides a fruitful platform with continuously tunable magnetism and topology for investigating emergent phenomena, including quantum anomalous Hall effect, Fermi arc states, etc.


title: Modified Emergent Dark Energy and its Astronomical Constraints

abstract: We introduce a modified form of the Phenomenologically Emergent Dark Energy (PEDE) model by showing a very elegant approach. The model is named as Modified Emergent Dark Energy (MEDE) to distinguish from PEDE model and it includes $\Lambda$CDM, PEDE model, Chevallier-Polarski-Linder model and other cosmological models of interest. We show that the present article offers a very fantastic route to construct other PEDE models in a simple but very elegant way. The model has seven free parameters where the six parameters are the same as in $\Lambda$CDM or PEDE model and the remaining one parameter `$\alpha$' quantifies the generalization. The present model predicts that dark energy equation of state at present assumes, $w_{\rm DE}\; (z=0) = -1 - \frac{\alpha}{3 \ln (10)}$ and in the far future (i.e., for $z \longrightarrow -1$), it will evolve asymptotically to $w_{\rm DE} \longrightarrow -1$. We perform a very robust observational analysis of the model using various observational datasets including cosmic microwave background radiation, baryon acoustic oscillations distance measurements, a local value of the Hubble constant and the pantheon sample of Supernovae Type Ia. We find that the $H_0$ tension here is alleviated but not solved, while the only free parameter, $\alpha$ could recover the $\Lambda$CDM or PEDE for different datasets. In summary, the present article describes a novel route to construct other modified versions of the PEDE model in a very simplified manner.


title: Non-Oscillatory No-Scale Inflation

abstract: We propose a non-oscillatory no-scale supergravity model of inflation (NO-NO inflation) in which the inflaton does not oscillate at the end of the inflationary era. Instead, the Universe is then dominated by the inflaton kinetic energy density (kination). During the transition from inflation to kination, the Universe preheats instantly through a coupling to Higgs-like fields. These rapidly annihilate and scatter into ultra-relativistic matter particles, which subsequently dominate the energy density, and reheating occurs at a temperature far above that of Big Bang Nucleosynthesis. After the electroweak transition, the inflaton enters a tracking phase as in some models of quintessential inflation. The model predictions for cosmic microwave background observables are consistent with Planck 2018 data, and the density of gravitational waves is below the upper bound from Big Bang Nucleosynthesis. We also find that the density of supersymmetric cold dark matter produced by gravitino decay is consistent with Planck 2018 data over the expected range of supersymmetric particle masses.


title: The Role of Domain Expertise in User Trust and the Impact of First Impressions with Intelligent Systems

abstract: Domain-specific intelligent systems are meant to help system users in their decision-making process. Many systems aim to simultaneously support different users with varying levels of domain expertise, but prior domain knowledge can affect user trust and confidence in detecting system errors. While it is also known that user trust can be influenced by first impressions with intelligent systems, our research explores the relationship between ordering bias and domain expertise when encountering errors in intelligent systems. In this paper, we present a controlled user study to explore the role of domain knowledge in establishing trust and susceptibility to the influence of first impressions on user trust. Participants reviewed an explainable image classifier with a constant accuracy and two different orders of observing system errors (observing errors in the beginning of usage vs. in the end). Our findings indicate that encountering errors early-on can cause negative first impressions for domain experts, negatively impacting their trust over the course of interactions. However, encountering correct outputs early helps more knowledgable users to dynamically adjust their trust based on their observations of system performance. In contrast, novice users suffer from over-reliance due to their lack of proper knowledge to detect errors.


title: The intrinsic topology of linear maps

abstract: Based on many experts' former work in the Jacobian conjecture and an essential analysis of intrinsic topology of linear maps, I completely prove the Jacobian conjecture by demonstrating the injectivity of real Keller map of any $n$-dimensions.


title: Finite presentation, the local lifting property, and local approximation properties of operator modules

abstract: We introduce notions of finite presentation and co-exactness which serve as qualitative and quantitative analogues of finite-dimensionality for operator modules over completely contractive Banach algebras. With these notions we begin the development of a local theory of operator modules by introducing analogues of the local lifting property, nuclearity, and semi-discreteness. For a large class of operator modules we prove that the local lifting property is equivalent to flatness, generalizing the operator space result of Kye and Ruan. We pursue applications to abstract harmonic analysis, where, for a locally compact quantum group $\mathbb{G}$, we show that $L^1(\mathbb{G})$-nuclearity of $\mathrm{LUC}(\mathbb{G})$ and $L^1(\mathbb{G})$-semi-discreteness of $L^\infty(\mathbb{G})$ are both equivalent to co-amenability of $\mathbb{G}$. We establish the equivalence between $A(G)$-injectivity of $G\bar{\ltimes}M$, $A(G)$-semi-discreteness of $G\bar{\ltimes} M$, and amenability of $W^*$-dynamical systems $(M,G,\alpha)$ with $M$ injective. We end with remarks on future directions.


title: A Plug-and-play Scheme to Adapt Image Saliency Deep Model for Video Data

abstract: With the rapid development of deep learning techniques, image saliency deep models trained solely by spatial information have occasionally achieved detection performance for video data comparable to that of the models trained by both spatial and temporal information. However, due to the lesser consideration of temporal information, the image saliency deep models may become fragile in the video sequences dominated by temporal information. Thus, the most recent video saliency detection approaches have adopted the network architecture starting with a spatial deep model that is followed by an elaborately designed temporal deep model. However, such methods easily encounter the performance bottleneck arising from the single stream learning methodology, so the overall detection performance is largely determined by the spatial deep model. In sharp contrast to the current mainstream methods, this paper proposes a novel plug-and-play scheme to weakly retrain a pretrained image saliency deep model for video data by using the newly sensed and coded temporal information. Thus, the retrained image saliency deep model will be able to maintain temporal saliency awareness, achieving much improved detection performance. Moreover, our method is simple yet effective for adapting any off-the-shelf pre-trained image saliency deep model to obtain high-quality video saliency detection. Additionally, both the data and source code of our method are publicly available.


title: A review of deep learning in medical imaging: Imaging traits, technology trends, case studies with progress highlights, and future promises

abstract: Since its renaissance, deep learning has been widely used in various medical imaging tasks and has achieved remarkable success in many medical imaging applications, thereby propelling us into the so-called artificial intelligence (AI) era. It is known that the success of AI is mostly attributed to the availability of big data with annotations for a single task and the advances in high performance computing. However, medical imaging presents unique challenges that confront deep learning approaches. In this survey paper, we first present traits of medical imaging, highlight both clinical needs and technical challenges in medical imaging, and describe how emerging trends in deep learning are addressing these issues. We cover the topics of network architecture, sparse and noisy labels, federating learning, interpretability, uncertainty quantification, etc. Then, we present several case studies that are commonly found in clinical practice, including digital pathology and chest, brain, cardiovascular, and abdominal imaging. Rather than presenting an exhaustive literature survey, we instead describe some prominent research highlights related to these case study applications. We conclude with a discussion and presentation of promising future directions.


title: Location-aware Graph Convolutional Networks for Video Question Answering

abstract: We addressed the challenging task of video question answering, which requires machines to answer questions about videos in a natural language form. Previous state-of-the-art methods attempt to apply spatio-temporal attention mechanism on video frame features without explicitly modeling the location and relations among object interaction occurred in videos. However, the relations between object interaction and their location information are very critical for both action recognition and question reasoning. In this work, we propose to represent the contents in the video as a location-aware graph by incorporating the location information of an object into the graph construction. Here, each node is associated with an object represented by its appearance and location features. Based on the constructed graph, we propose to use graph convolution to infer both the category and temporal locations of an action. As the graph is built on objects, our method is able to focus on the foreground action contents for better video question answering. Lastly, we leverage an attention mechanism to combine the output of graph convolution and encoded question features for final answer reasoning. Extensive experiments demonstrate the effectiveness of the proposed methods. Specifically, our method significantly outperforms state-of-the-art methods on TGIF-QA, Youtube2Text-QA, and MSVD-QA datasets. Code and pre-trained models are publicly available at: https://github.com/SunDoge/L-GCN


title: Generative View Synthesis: From Single-view Semantics to Novel-view Images

abstract: Content creation, central to applications such as virtual reality, can be a tedious and time-consuming. Recent image synthesis methods simplify this task by offering tools to generate new views from as little as a single input image, or by converting a semantic map into a photorealistic image. We propose to push the envelope further, and introduce Generative View Synthesis (GVS), which can synthesize multiple photorealistic views of a scene given a single semantic map. We show that the sequential application of existing techniques, e.g., semantics-to-image translation followed by monocular view synthesis, fail at capturing the scene's structure. In contrast, we solve the semantics-to-image translation in concert with the estimation of the 3D layout of the scene, thus producing geometrically consistent novel views that preserve semantic structures. We first lift the input 2D semantic map onto a 3D layered representation of the scene in feature space, thereby preserving the semantic labels of 3D geometric structures. We then project the layered features onto the target views to generate the final novel-view images. We verify the strengths of our method and compare it with several advanced baselines on three different datasets. Our approach also allows for style manipulation and image editing operations, such as the addition or removal of objects, with simple manipulations of the input style images and semantic maps respectively. Visit the project page at https://gvsnet.github.io.


title: Greedoids from flames

abstract: A digraph $ D $ with $ r\in V(D) $ is an $ r $-flame if for every $ {v\in V(D)-r} $, the in-degree of $ v $ is equal to the local edge-connectivity $ \lambda_D(r,v) $. We show that for every digraph $ D $ and $ r\in V(D) $, the edge sets of the $ r $-flame subgraphs of $ D $ form a greedoid. Our method yields a new proof of Lov\'asz' theorem stating: for every digraph $ D $ and $ r\in V(D) $, there is an $ r $-flame subdigraph $ F $ of $ D $ such that $ \lambda_F(r,v) =\lambda_D(r,v) $ for $ v\in V(D)-r $. We also give a strongly polynomial algorithm to find such an $ F $ working with a fractional generalization of Lov\'asz' theorem.


title: Analytic Calibration in Andreasen-Huge SABR Model

abstract: We derive analytic formulae which link $\alpha$, $\nu$ and $\rho$ parameters in Andreasen-Huge style SABR model to the ATM price and option prices at four strikes close to ATM. Based on these formulae we give a characterisation for the SABR parameters in terms of derivatives of the swap rate forward probability density function. We test the analytic result in the application to the interest rate futures option market.

categories: q-fin.CP math.PR q-fin.MF q-fin.PR q-fin.RM


title: Modeling the effects of prosocial awareness on COVID-19 dynamics: A case study on Colombia

abstract: The ongoing COVID-19 pandemic has affected most of the countries on Earth. It has become a pandemic outbreak with more than 24 million confirmed infections and above 840 thousand deaths worldwide. In this study, we consider a mathematical model on COVID-19 transmission with the prosocial awareness effect. The proposed model can have four equilibrium states based on different parametric conditions. The local and global stability conditions for awareness free, disease-free equilibrium is studied. Using Lyapunov function theory and LaSalle Invariance Principle, the disease-free equilibrium is shown globally asymptotically stable under some parametric constraints. The existence of unique awareness free, endemic equilibrium and unique endemic equilibrium is presented. We calibrate our proposed model parameters to fit daily cases and deaths from Colombia. Sensitivity analysis indicates that the transmission rate and learning factor related to awareness of susceptibles are very crucial for reduction in disease related deaths. Finally, we assess the impact of prosocial awareness during the outbreak and compare this strategy with popular control measures. Results indicate that prosocial awareness has competitive potential to flatten the curve.


title: Weakly-supervised 3D Shape Completion in the Wild

abstract: 3D shape completion for real data is important but challenging, since partial point clouds acquired by real-world sensors are usually sparse, noisy and unaligned. Different from previous methods, we address the problem of learning 3D complete shape from unaligned and real-world partial point clouds. To this end, we propose a weakly-supervised method to estimate both 3D canonical shape and 6-DoF pose for alignment, given multiple partial observations associated with the same instance. The network jointly optimizes canonical shapes and poses with multi-view geometry constraints during training, and can infer the complete shape given a single partial point cloud. Moreover, learned pose estimation can facilitate partial point cloud registration. Experiments on both synthetic and real data show that it is feasible and promising to learn 3D shape completion through large-scale data without shape and pose supervision.


title: Varying-coefficient stochastic differential equations with applications in ecology

abstract: Stochastic differential equations (SDEs) are popular tools to analyse time series data in many areas, such as mathematical finance, physics, and biology. They provide a mechanistic description of the phenomeon of interest, and their parameters often have a clear interpretation. These advantages come at the cost of requiring a relatively simple model specification. We propose a flexible model for SDEs with time-varying dynamics where the parameters of the process are non-parametric functions of covariates, similar to generalized additive models. Combining the SDEs and non-parametric approaches allows for the SDE to capture more detailed, non-stationary, features of the data-generating process. We present a computationally efficient method of approximate inference, where the SDE parameters can vary according to fixed covariate effects, random effects, or basis-penalty smoothing splines. We demonstrate the versatility and utility of this approach with three applications in ecology, where there is often a modelling trade-off between interpretability and flexibility.


title: Inducing Language-Agnostic Multilingual Representations

abstract: Cross-lingual representations have the potential to make NLP techniques available to the vast majority of languages in the world. However, they currently require large pretraining corpora or access to typologically similar languages. In this work, we address these obstacles by removing language identity signals from multilingual embeddings. We examine three approaches for this: (i) re-aligning the vector spaces of target languages (all together) to a pivot source language; (ii) removing language-specific means and variances, which yields better discriminativeness of embeddings as a by-product; and (iii) increasing input similarity across languages by removing morphological contractions and sentence reordering. We evaluate on XNLI and reference-free MT across 19 typologically diverse languages. Our findings expose the limitations of these approaches -- unlike vector normalization, vector space re-alignment and text normalization do not achieve consistent gains across encoders and languages. Due to the approaches' additive effects, their combination decreases the cross-lingual transfer gap by 8.9 points (m-BERT) and 18.2 points (XLM-R) on average across all tasks and languages, however. Our code and models are publicly available.


title: Perturbative instability towards delocalization at phase transitions between MBL phases

abstract: We examine the stability of marginally Anderson localized phase transitions between localized phases to the addition of many-body interactions, focusing in particular on the spin-glass to paramagnet transition in a disordered transverse field Ising model in one dimension. We find evidence for a perturbative instability of localization at finite energy densities once interactions are added, i.e. evidence for the relevance of interactions - in a renormalization group sense - to the non-interacting critical point governed by infinite randomness scaling. We introduce a novel diagnostic, the "susceptibility of entanglement", which allows us to perturbatively probe the effect of adding interactions on the entanglement properties of eigenstates, and helps us elucidate the resonant processes that can cause thermalization. The susceptibility serves as a much more sensitive probe, and its divergence can detect the perturbative beginnings of an incipient instability even in regimes and system sizes for which conventional diagnostics point towards localization. We expect this new measure to be of independent interest for analyzing the stability of localization in a variety of different settings.


title: Thermodynamics and and phase transition of topological dilatonic Lifshitz-like black holes

abstract: It is known that scalar-tensor gravity models can be studied in Einstein and Jordan frames. In this paper, we consider a model of scalar-tensor gravity in Einstein's frame to calculate the Lifshitz-like black hole solutions with different horizon topologies. We study thermodynamic properties and first order van der Waals like phase transition, and find that the Lifshitz parameter affects the phase structure. In addition, we investigate thermal stability by using the behavior of heat capacity and various methods of geometrical thermodynamics.


title: Constraining the Generalized Uncertainty Principle Through Black Hole Shadow and Quasiperiodic Oscillations

abstract: In this paper we study the effect of the Generalized Uncertainty Principle (GUP) on the shadow of GUP-modified Kerr black hole and the correspondence between the shadow radius and the real part of the quasinormal modes (QNMs). We find that the shadow curvature radius of the GUP-modfied Kerr black hole is bigger compared to the Kerr vacuum solution and increases linearly monotonically with the increase of the GUP parameter. We then investigate the characteristic points of intrinsic curvature of the shadow from a topological point of view to calculate the the angular size for these curvature radii of the shadow. To this end, we have used the EHT data for the M87* black hole to constrain the upper limits of the GUP parameter red and our result is $\beta<10^{95}$. Finally, we have explored the connection between the shadow radius and the scalar/electromagnetic/gravitational QNMs. The GUP-modified Kerr black hole is also used to provide perfect curve fitting of the particle oscillation upper and lower frequencies to the observed frequencies for three microquasars and to restrict the values of the correction parameter in the metric of the modified black hole to very reasonable bound $\beta<10^{77}$.


title: Electric Dipole Moments, New Forces and Dark Matter

abstract: New sources of CP violation beyond the Standard Model are crucial to explain the baryon asymmetry in the Universe. We discuss the impact of new CP violating interactions in theories where a dark matter candidate is predicted by the cancellation of gauge anomalies. In these theories, the constraint on the dark matter relic density implies an upper bound on the new symmetry breaking scale from which all new states acquire their masses. We investigate in detail the predictions for electric dipole moments and show that if the relevant CP-violating phase is large, experiments such as the ACME collaboration will be able to fully probe the theory.


title: The Fate of Discrete 1-Form Symmetries in 6d

abstract: Recently introduced generalized global symmetries have been useful in order to understand non-perturbative aspects of quantum field theories in four and lower dimensions. In this paper we focus on 1-form symmetries of weakly coupled 6d supersymmetric gauge theories coupled to dynamical tensor multiplets. We study the consistency of global 1-form symmetries corresponding to the center of the gauge groups, or subgroups thereof, by activating their background fields, which makes the instanton density fractional. In 6d, an instanton background for a given gauge theory sources BPS strings via tadpole cancellation. The non-trivial 1-form symmetry background configurations contribute to the charge of the BPS strings. However, Dirac quantization imposes restrictions on the consistent 1-form backgrounds, since they can in general lead to and induce fractional charges, thus making (part of) the putative higher-form symmetry inconsistent. This gives explicit criteria to determine whether the discrete 1-form symmetries are realized. We implement these criteria in concrete examples originating from string compactifications. We also corroborate this by finding that a non-trivial fractional contribution is related to states which explicitly break the global 1-form symmetry appearing as massive excitations of the 6d BPS strings. For 6d theories consistently coupled to gravity, this hints at a symmetry breaking tower of states. When the fractional contributions are absent, the F-theory realization of the theories points to the gauging of the 1-form symmetry via the presence of non-trivial Mordell--Weil torsion.


title: Interaction of a cold cloud with a hot wind: the regimes of cloud growth and destruction and the impact of magnetic fields

abstract: Multiphase galaxy winds, the accretion of cold gas through galaxy haloes, and gas stripping from jellyfish galaxies are examples of interactions between cold and hot gaseous phases. There are two important regimes in such systems. A sufficiently small cold cloud is destroyed by the hot wind as a result of Kelvin-Helmholtz instabilities, which shatter the cloud into small pieces that eventually mix and dissolve in the hot wind. On the contrary, stripped cold gas from a large cloud mixes with the hot wind to intermediate temperatures, and then becomes thermally unstable and cools, causing a net accretion of hot gas to the cold tail. Using the magneto-hydrodynamical code AREPO, we perform cloud crushing simulations and test analytical criteria for the transition between the growth and destruction regimes to clarify a current debate in the literature. We find that the hot-wind cooling time sets the transition radius and not the cooling time of the mixed phase. Magnetic fields modify the wind-cloud interaction. Draping of wind magnetic field enhances the field upstream of the cloud and fluid instabilities are suppressed by a turbulently magnetised wind beyond what is seen for a wind with a uniform magnetic field. We furthermore predict jellyfish galaxies to have ordered magnetic fields aligned with their tails. We finally discuss how the results of idealised simulations can be used to provide input to subgrid models in cosmological (magneto-)hydrodynamical simulations, which cannot resolve the detailed small-scale structure of cold gas clouds in the circum-galactic medium.


title: Accidental SO(10) axion from gauged flavour

abstract: An accidental U(1) Peccei-Quinn (PQ) symmetry automatically arises in a class of SO(10) unified theories upon gauging the SU(3)$_{f}$ flavour group. The PQ symmetry is protected by the $\mathbb{Z}_4 \times \mathbb{Z}_3$ center of $\text{SO(10)} \times \text{SU(3)}_f$ up to effective operators of canonical dimension six. However, high-scale contributions to the axion potential posing a PQ quality problem arise only at $d=9$. In the pre-inflationary PQ breaking scenario the axion mass window is predicted to be $m_a \in [7 \times 10^{-8}, 10^{-3}]$ eV, where the lower end is bounded by the seesaw scale and the upper end by iso-curvature fluctuations. A high-quality axion, that is immune to the PQ quality problem, is obtained for $m_a \gtrsim 0.02$ eV. We finally offer a general perspective on the PQ quality problem in grand unified theories.


title: Cosmological DCBH formation sites hostile for their growth

abstract: The direct collapse (DC) is a promising mechanism that provides massive seed black holes (BHs) with $\sim 10^{5}~M_{\odot}$ in the early universe. To study a long-term accretion growth of a DCBH thus formed, we perform cosmological radiation-hydrodynamics simulations, extending our previous work where we investigated its formation stage. With a high spatial resolution down below the Bondi radius, we show that the accretion rate onto the BH is far below the Eddington value. Such slow mass growth is partly because of the strong radiative feedback from the accreting BH. Moreover, we find that the BH has a large velocity of $\sim 100~{\rm km~s^{-1}}$ relative to the gas after it falls into the first galaxy, which substantially reduces the accretion rate. The latter effect stems from the fact that the DCBHs form in metal-free environments typically at $\sim 1~$kpc from the galaxy. The BH accelerates as it approaches the galactic center due to the gravity. The relative velocity never damps after that, and the BH does not settle down to the galactic center but continues to wander around it. An analytic estimate predicts that the DCBH formation within $\sim 100$~pc around the galactic center is necessary to decelerate the BH with dynamical friction before $z=7$. Since metal enrichment with $Z \sim 10^{-5} - 10^{-3}~Z_\odot$ is expected in such a case, the formation of DCBHs in the metal-poor environments is preferable for the subsequent rapid growth.


title: Cold gas in the Milky Way's nuclear wind

abstract: The centre of the Milky Way is the site of several high-energy processes that have strongly impacted the inner regions of our Galaxy. Activity from the super-massive black hole, Sgr A*, and/or stellar feedback from the inner molecular ring expel matter and energy from the disc in the form of a galactic wind. Multiphase gas has been observed within this outflow, from hot highly-ionized, to warm ionized and cool atomic gas. To date, however, there has been no evidence of the cold and dense molecular phase. Here we report the first detection of molecular gas outflowing from the centre of our Galaxy. This cold material is associated with atomic hydrogen clouds travelling in the nuclear wind. The morphology and the kinematics of the molecular gas, resolved on ~1 pc scale, indicate that these clouds are mixing with the warmer medium and are possibly being disrupted. The data also suggest that the mass of molecular gas driven out is not negligible and could impact the rate of star formation in the central regions. The presence of this cold, dense, high-velocity gas is puzzling, as neither Sgr A* at its current level of activity, nor star formation in the inner Galaxy seem viable sources for this material.


title: Designing three-dimensional flat bands in nodal-line semimetals

abstract: Electrons with large kinetic energy have a superconducting instability for infinitesimal attractive interactions. Quenching the kinetic energy and creating a flat band renders an infinitesimal repulsive interaction the relevant perturbation. Thus, flat band systems are an ideal platform to study the competition of superconductivity and magnetism and their possible coexistence. Recent advances in the field of twisted bilayer graphene highlight this in the context of two-dimensional materials. Two dimensions, however, put severe restrictions on the stability of the low-temperature phases due to enhanced fluctuations. Only three-dimensional flat bands can solve the conundrum of combining the exotic flat-band phases with stable order existing at high temperatures. Here, we present a way to generate such flat bands through strain engineering in topological nodal-line semimetals. We present analytical and numerical evidence for this scenario and study the competition of the arising superconducting and magnetic orders as a function of externally controlled parameters. We show that the order parameter is rigid because the quantum geometry of the Bloch wave functions leads to a large superfluid stiffness. Using density-functional theory and numerical tight-binding calculations we further apply our theory to strained rhombohedral graphite and CaAgP materials.


title: Cosmological perturbations for two cold fluids in $\Lambda$CDM

abstract: The cosmic large-scale structure of our Universe is comprised of baryons and cold dark matter (CDM). Yet it is customary to treat these two components as a combined single-matter fluid with vanishing pressure, which is justified only for sufficiently large scales and late times. Here we go beyond the single-fluid approximation and develop the perturbation theory for two gravitationally coupled fluids while still assuming vanishing pressure. We mostly focus on perturbative expansions in powers of $D$ (or $D_+$), the linear structure growth of matter in a $\Lambda$CDM Universe with cosmological constant $\Lambda$. We derive in particular (1) explicit recursion relations for the two fluid densities, (2) complementary all-order results in the Lagrangian-coordinates approach, as well as (3) the associated component wavefunctions in a semi-classical approach to cosmic large-scale structure. In our companion paper (Hahn et al. 2020) we apply these new theoretical results to generate novel higher-order initial conditions for cosmological hydrodynamical simulations.


title: Higher-order initial conditions for mixed baryon-CDM simulations

abstract: We present a novel approach to generate higher-order initial conditions (ICs) for cosmological simulations that take into account the distinct evolution of baryons and dark matter. We focus on the numerical implementation and the validation of its performance, based on both collisionless N-body simulations and full hydrodynamic Eulerian and Lagrangian simulations. We improve in various ways over previous approaches that were limited to first-order Lagrangian perturbation theory (LPT). Specifically, we (1) generalize nth-order LPT to multi-fluid systems, allowing 2LPT or 3LPT ICs for two-fluid simulations, (2) employ a novel propagator perturbation theory to set up ICs for Eulerian codes that are fully consistent with 1LPT or 2LPT, (3) demonstrate that our ICs resolve previous problems of two-fluid simulations by using variations in particle masses that eliminate spurious deviations from expected perturbative results, (4) show that the improvements achieved by going to higher-order PT are comparable to those seen for single-fluid ICs, and (5) demonstrate the excellent (i.e., few per cent level) agreement between Eulerian and Lagrangian simulations, once high-quality initial conditions are used. The rigorous development of the underlying perturbation theory is presented in a companion paper. All presented algorithms are implemented in the Monofonic Music-2 package that we make publicly available.


title: Convective Differential Rotation in Stars and Planets I: Theory

abstract: We derive the scaling of differential rotation in both slowly- and rapidly-rotating convection zones using order of magnitude methods. Our calculations apply across stars and fluid planets and all rotation rates, as well as to both magnetized and purely hydrodynamic systems. We find shear $|R\nabla\Omega|$ of order the angular frequency $\Omega$ for slowly-rotating systems with $\Omega \ll |N|$, where $N$ is the \brvs\ frequency, and find that it declines as a power-law in $\Omega$ for rapidly-rotating systems with $\Omega \gg |N|$. We further calculate the meridional circulation rate and baroclinicity and examine the magnetic field strength in the rapidly rotating limit. Our results are in general agreement with simulations and observations and we perform a detailed comparison with those in a companion paper.


title: Convective Differential Rotation in Stars and Planets II: Observational and Numerical Tests

abstract: Differential rotation is central to a great many mysteries in stars and planets. In Part I we predicted the order of magnitude and scaling of the differential rotation in both hydrodynamic and magnetohydrodynamic convection zones. Our results apply to both slowly- and rapidly-rotating systems, and provide a general picture of differential rotation in stars and fluid planets. We further calculated the scalings of the meridional circulation, entropy gradient and baroclinicity. In this companion paper we compare these predictions with a variety of observations and numerical simulations. With a few exceptions we find that these are consistent in both the slowly-rotating and rapidly-rotating limits. Our results help to localize core-envelope shear in red~giant stars, suggest a rotation-dependent frequency shift in the internal gravity waves of massive stars and potentially explain observed deviations from von Zeipel's gravity darkening in late-type stars.