title: RIAV-MVS: Recurrent-Indexing an Asymmetric Volume for Multi-View Stereo

abstract: This paper presents a learning-based method for multi-view depth estimation from posed images. Our core idea is a "learning-to-optimize" paradigm that iteratively indexes a plane-sweeping cost volume and regresses the depth map via a convolutional Gated Recurrent Unit (GRU). Since the cost volume plays a paramount role in encoding the multi-view geometry, we aim to improve its construction both at pixel- and frame- levels. At the pixel level, we propose to break the symmetry of the Siamese network (which is typically used in MVS to extract image features) by introducing a transformer block to the reference image (but not to the source images). Such an asymmetric volume allows the network to extract global features from the reference image to predict its depth map. Given potential inaccuracies in the poses between reference and source images, we propose to incorporate a residual pose network to correct the relative poses. This essentially rectifies the cost volume at the frame level. We conduct extensive experiments on real-world MVS datasets and show that our method achieves state-of-the-art performance in terms of both within-dataset evaluation and cross-dataset generalization.


title: Automatic Expert Selection for Multi-Scenario and Multi-Task Search

abstract: Multi-scenario learning (MSL) enables a service provider to cater for users' fine-grained demands by separating services for different user sectors, e.g., by user's geographical region. Under each scenario there is a need to optimize multiple task-specific targets e.g., click through rate and conversion rate, known as multi-task learning (MTL). Recent solutions for MSL and MTL are mostly based on the multi-gate mixture-of-experts (MMoE) architecture. MMoE structure is typically static and its design requires domain-specific knowledge, making it less effective in handling both MSL and MTL. In this paper, we propose a novel Automatic Expert Selection framework for Multi-scenario and Multi-task search, named AESM^{2}. AESM^{2} integrates both MSL and MTL into a unified framework with an automatic structure learning. Specifically, AESM^{2} stacks multi-task layers over multi-scenario layers. This hierarchical design enables us to flexibly establish intrinsic connections between different scenarios, and at the same time also supports high-level feature extraction for different tasks. At each multi-scenario/multi-task layer, a novel expert selection algorithm is proposed to automatically identify scenario-/task-specific and shared experts for each input. Experiments over two real-world large-scale datasets demonstrate the effectiveness of AESM^{2} over a battery of strong baselines. Online A/B test also shows substantial performance gain on multiple metrics. Currently, AESM^{2} has been deployed online for serving major traffic.


title: Ab initio study of the crystal and electronic structure of mono- and bi-alkali antimonides: Stability, Goldschmidt-like tolerance factors, and optical properties

abstract: Mono- and bi-alkali antimonides, X$_2$YSb (X and Y from Group I), are promising for next-generation electron emitters due to their capability of producing high-quality electron beams. However, these materials are not yet well understood, in part due to the technical challenges in growing pure, ordered alkali antimonides. For example, in the current literature there is a lack of complete understanding of the mechanically stable crystal structures of these materials. As a first step towards understanding this issue, this paper presents an ab initio study of stability of single-crystal mono- and bi-alkali antimonides in the $D0_3$ structure, the structure generally assumed in the literature for these materials. Finding that many of these materials actually are unstable in the $D0_3$ structure, we formulate a new set of Goldschmidt-like tolerance factors that accurately predict $D0_3$ stability using a procedure analogous to machine-learning perceptron-based analysis. Next, we consider possible stable structures for materials that we predict to be unstable in the $D0_3$ structure. Taking as examples the mono- and bi-alkali antimonides Cs$_3$Sb and Cs$_2$KSb, which also are technologically interesting for photoemission and photoabsorption applications, respectively, we note that the most unstable phonon displacements are consistent with the cubic structure, and we therefore perform extensive ab initio searches to identify potential ground-state structures in a cubic lattice. Our X-ray diffraction experiments confirm that indeed these two materials are not stable in the $D0_3$ structure and show scattering that is consistent with our new, proposed stable structures. Finally, we explore ab initio the implications of the breaking of the $D0_3$ symmetry on the electronic structure, showing significant impact on the location of the optical absorption edge.


title: Multi-agent Databases via Independent Learning

abstract: Machine learning is rapidly being used in database research to improve the effectiveness of numerous tasks included but not limited to query optimization, workload scheduling, physical design, etc. Currently, the research focus has been on replacing a single database component responsible for one task by its learning-based counterpart. However, query performance is not simply determined by the performance of a single component, but by the cooperation of multiple ones. As such, learning based database components need to collaborate during both training and execution in order to develop policies that meet end performance goals. Thus, the paper attempts to address the question "Is it possible to design a database consisting of various learned components that cooperatively work to improve end-to-end query latency?". To answer this question, we introduce MADB (Multi-Agent DB), a proof-of-concept system that incorporates a learned query scheduler and a learned query optimizer. MADB leverages a cooperative multi-agent reinforcement learning approach that allows the two components to exchange the context of their decisions with each other and collaboratively work towards reducing the query latency. Preliminary results demonstrate that MADB can outperform the non-cooperative integration of learned components.


title: Differentially Private Covariance Revisited

abstract: In this paper, we present two new algorithms for covariance estimation under concentrated differential privacy (zCDP). The first algorithm achieves a Frobenius error of $\tilde{O}(d^{1/4}\sqrt{\mathrm{tr}}/\sqrt{n} + \sqrt{d}/n)$, where $\mathrm{tr}$ is the trace of the covariance matrix. By taking $\mathrm{tr}=1$, this also implies a worst-case error bound of $\tilde{O}(d^{1/4}/\sqrt{n})$, which improves the standard Gaussian mechanism's $\tilde{O}(d/n)$ for the regime $d>\widetilde{\Omega}(n^{2/3})$. Our second algorithm offers a tail-sensitive bound that could be much better on skewed data. The corresponding algorithms are also simple and efficient. Experimental results show that they offer significant improvements over prior work.


title: Feature subset selection for kernel SVM classification via mixed-integer optimization

abstract: We study the mixed-integer optimization (MIO) approach to feature subset selection in nonlinear kernel support vector machines (SVMs) for binary classification. First proposed for linear regression in the 1970s, this approach has recently moved into the spotlight with advances in optimization algorithms and computer hardware. The goal of this paper is to establish an MIO approach for selecting the best subset of features for kernel SVM classification. To measure the performance of subset selection, we use the kernel-target alignment, which is the distance between the centroids of two response classes in a high-dimensional feature space. We propose a mixed-integer linear optimization (MILO) formulation based on the kernel-target alignment for feature subset selection, and this MILO problem can be solved to optimality using optimization software. We also derive a reduced version of the MILO problem to accelerate our MILO computations. Experimental results show good computational efficiency for our MILO formulation with the reduced problem. Moreover, our method can often outperform the linear-SVM-based MILO formulation and recursive feature elimination in prediction performance, especially when there are relatively few data instances.


title: Adaptive Activation Network For Low Resource Multilingual Speech Recognition

abstract: Low resource automatic speech recognition (ASR) is a useful but thorny task, since deep learning ASR models usually need huge amounts of training data. The existing models mostly established a bottleneck (BN) layer by pre-training on a large source language, and transferring to the low resource target language. In this work, we introduced an adaptive activation network to the upper layers of ASR model, and applied different activation functions to different languages. We also proposed two approaches to train the model: (1) cross-lingual learning, replacing the activation function from source language to target language, (2) multilingual learning, jointly training the Connectionist Temporal Classification (CTC) loss of each language and the relevance of different languages. Our experiments on IARPA Babel datasets demonstrated that our approaches outperform the from-scratch training and traditional bottleneck feature based methods. In addition, combining the cross-lingual learning and multilingual learning together could further improve the performance of multilingual speech recognition.


title: Efficient Policy Iteration for Robust Markov Decision Processes via Regularization

abstract: Robust Markov decision processes (MDPs) provide a general framework to model decision problems where the system dynamics are changing or only partially known. Efficient methods for some \texttt{sa}-rectangular robust MDPs exist, using its equivalence with reward regularized MDPs, generalizable to online settings. In comparison to \texttt{sa}-rectangular robust MDPs, \texttt{s}-rectangular robust MDPs are less restrictive but much more difficult to deal with. Interestingly, recent works have established the equivalence between \texttt{s}-rectangular robust MDPs and policy regularized MDPs. But we don't have a clear understanding to exploit this equivalence, to do policy improvement steps to get the optimal value function or policy. We don't have a clear understanding of greedy/optimal policy except it can be stochastic. There exist no methods that can naturally be generalized to model-free settings. We show a clear and explicit equivalence between \texttt{s}-rectangular $L_p$ robust MDPs and policy regularized MDPs that resemble very much policy entropy regularized MDPs widely used in practice. Further, we dig into the policy improvement step and concretely derive optimal robust Bellman operators for \texttt{s}-rectangular $L_p$ robust MDPs. We find that the greedy/optimal policies in \texttt{s}-rectangular $L_p$ robust MDPs are threshold policies that play top $k$ actions whose $Q$ value is greater than some threshold (value), proportional to the $(p-1)$th power of its advantage. In addition, we show time complexity of (\texttt{sa} and \texttt{s}-rectangular) $L_p$ robust MDPs is the same as non-robust MDPs up to some log factors. Our work greatly extends the existing understanding of \texttt{s}-rectangular robust MDPs and naturally generalizable to online settings.


title: Point RCNN: An Angle-Free Framework for Rotated Object Detection

abstract: Rotated object detection in aerial images is still challenging due to arbitrary orientations, large scale and aspect ratio variations, and extreme density of objects. Existing state-of-the-art rotated object detection methods mainly rely on angle-based detectors. However, angle regression can easily suffer from the long-standing boundary problem. To tackle this problem, we propose a purely angle-free framework for rotated object detection, called Point RCNN, which mainly consists of PointRPN and PointReg. In particular, PointRPN generates accurate rotated RoIs (RRoIs) by converting the learned representative points with a coarse-to-fine manner, which is motivated by RepPoints. Based on the learned RRoIs, PointReg performs corner points refinement for more accurate detection. In addition, aerial images are often severely unbalanced in categories, and existing methods almost ignore this issue. In this paper, we also experimentally verify that re-sampling the images of the rare categories will stabilize training and further improve the detection performance. Experiments demonstrate that our Point RCNN achieves the new state-of-the-art detection performance on commonly used aerial datasets, including DOTA-v1.0, DOTA-v1.5, and HRSC2016.


title: Speech Augmentation Based Unsupervised Learning for Keyword Spotting

abstract: In this paper, we investigated a speech augmentation based unsupervised learning approach for keyword spotting (KWS) task. KWS is a useful speech application, yet also heavily depends on the labeled data. We designed a CNN-Attention architecture to conduct the KWS task. CNN layers focus on the local acoustic features, and attention layers model the long-time dependency. To improve the robustness of KWS model, we also proposed an unsupervised learning method. The unsupervised loss is based on the similarity between the original and augmented speech features, as well as the audio reconstructing information. Two speech augmentation methods are explored in the unsupervised learning: speed and intensity. The experiments on Google Speech Commands V2 Dataset demonstrated that our CNN-Attention model has competitive results. Moreover, the augmentation based unsupervised learning could further improve the classification accuracy of KWS task. In our experiments, with augmentation based unsupervised learning, our KWS model achieves better performance than other unsupervised methods, such as CPC, APC, and MPC.


title: Differentiable Point-Based Radiance Fields for Efficient View Synthesis

abstract: We propose a differentiable rendering algorithm for efficient novel view synthesis. By departing from volume-based representations in favor of a learned point representation, we improve on existing methods more than an order of magnitude in memory and runtime, both in training and inference. The method begins with a uniformly-sampled random point cloud and learns per-point position and view-dependent appearance, using a differentiable splat-based renderer to evolve the model to match a set of input images. Our method is up to 300x faster than NeRF in both training and inference, with only a marginal sacrifice in quality, while using less than 10~MB of memory for a static scene. For dynamic scenes, our method trains two orders of magnitude faster than STNeRF and renders at near interactive rate, while maintaining high image quality and temporal coherence even without imposing any temporal-coherency regularizers.


title: Survival Analysis on Structured Data using Deep Reinforcement Learning

abstract: Survival analysis is playing a major role in manufacturing sector by analyzing occurrence of any unwanted event based on the input data. Predictive maintenance, which is a part of survival analysis, helps to find any device failure based on the current incoming data from different sensor or any equipment. Deep learning techniques were used to automate the predictive maintenance problem to some extent, but they are not very helpful in predicting the device failure for the input data which the algorithm had not learned. Since neural network predicts the output based on previous learned input features, it cannot perform well when there is more variation in input features. Performance of the model is degraded with the occurrence of changes in input data and finally the algorithm fails in predicting the device failure. This problem can be solved by our proposed method where the algorithm can predict the device failure more precisely than the existing deep learning algorithms. The proposed solution involves implementation of Deep Reinforcement Learning algorithm called Double Deep Q Network (DDQN) for classifying the device failure based on the input features. The algorithm is capable of learning different variation of the input feature and is robust in predicting whether the device will fail or not based on the input data. The proposed DDQN model is trained with limited or lesser amount of input data. The trained model predicted larger amount of test data efficiently and performed well compared to other deep learning and machine learning models.


title: V4D: Voxel for 4D Novel View Synthesis

abstract: Neural radiance fields have made a remarkable breakthrough in the novel view synthesis task at the 3D static scene. However, for the 4D circumstance (e.g., dynamic scene), the performance of the existing method is still limited by the capacity of the neural network, typically in a multilayer perceptron network (MLP). In this paper, we utilize 3D Voxel to model the 4D neural radiance field, short as V4D, where the 3D voxel has two formats. The first one is to regularly model the 3D space and then use the sampled local 3D feature with the time index to model the density field and the texture field by a tiny MLP. The second one is in look-up tables (LUTs) format that is for the pixel-level refinement, where the pseudo-surface produced by the volume rendering is utilized as the guidance information to learn a 2D pixel-level refinement mapping. The proposed LUTs-based refinement module achieves the performance gain with little computational cost and could serve as the plug-and-play module in the novel view synthesis task. Moreover, we propose a more effective conditional positional encoding toward the 4D data that achieves performance gain with negligible computational burdens. Extensive experiments demonstrate that the proposed method achieves state-of-the-art performance at a low computational cost.


title: One Reference Is Not Enough: Diverse Distillation with Reference Selection for Non-Autoregressive Translation

abstract: Non-autoregressive neural machine translation (NAT) suffers from the multi-modality problem: the source sentence may have multiple correct translations, but the loss function is calculated only according to the reference sentence. Sequence-level knowledge distillation makes the target more deterministic by replacing the target with the output from an autoregressive model. However, the multi-modality problem in the distilled dataset is still nonnegligible. Furthermore, learning from a specific teacher limits the upper bound of the model capability, restricting the potential of NAT models. In this paper, we argue that one reference is not enough and propose diverse distillation with reference selection (DDRS) for NAT. Specifically, we first propose a method called SeedDiv for diverse machine translation, which enables us to generate a dataset containing multiple high-quality reference translations for each source sentence. During the training, we compare the NAT output with all references and select the one that best fits the NAT output to train the model. Experiments on widely-used machine translation benchmarks demonstrate the effectiveness of DDRS, which achieves 29.82 BLEU with only one decoding pass on WMT14 En-De, improving the state-of-the-art performance for NAT by over 1 BLEU. Source code: https://github.com/ictnlp/DDRS-NAT


title: Teaching Models to Express Their Uncertainty in Words

abstract: We show that a GPT-3 model can learn to express uncertainty about its own answers in natural language -- without use of model logits. When given a question, the model generates both an answer and a level of confidence (e.g. "90% confidence" or "high confidence"). These levels map to probabilities that are well calibrated. The model also remains moderately calibrated under distribution shift, and is sensitive to uncertainty in its own answers, rather than imitating human examples. To our knowledge, this is the first time a model has been shown to express calibrated uncertainty about its own answers in natural language. For testing calibration, we introduce the CalibratedMath suite of tasks. We compare the calibration of uncertainty expressed in words ("verbalized probability") to uncertainty extracted from model logits. Both kinds of uncertainty are capable of generalizing calibration under distribution shift. We also provide evidence that GPT-3's ability to generalize calibration depends on pre-trained latent representations that correlate with epistemic uncertainty over its answers.


title: Convergence and Sample Complexity of Policy Gradient Methods for Stabilizing Linear Systems

abstract: System stabilization via policy gradient (PG) methods has drawn increasing attention in both control and machine learning communities. In this paper, we study their convergence and sample complexity for stabilizing linear time-invariant systems in terms of the number of system rollouts. Our analysis is built upon a discounted linear quadratic regulator (LQR) method which alternatively updates the policy and the discount factor of the LQR problem. Firstly, we propose an explicit rule to adaptively adjust the discount factor by exploring the stability margin of a linear control policy. Then, we establish the sample complexity of PG methods for stabilization, which only adds a coefficient logarithmic in the spectral radius of the state matrix to that for solving the LQR problem with a prior stabilizing policy. Finally, we perform simulations to validate our theoretical findings and demonstrate the effectiveness of our method on a class of nonlinear systems.


title: Gating Dropout: Communication-efficient Regularization for Sparsely Activated Transformers

abstract: Sparsely activated transformers, such as Mixture of Experts (MoE), have received great interest due to their outrageous scaling capability which enables dramatical increases in model size without significant increases in computational cost. To achieve this, MoE models replace the feedforward sub-layer with Mixture-of-Experts sub-layer in transformers and use a gating network to route each token to its assigned experts. Since the common practice for efficient training of such models requires distributing experts and tokens across different machines, this routing strategy often incurs huge cross-machine communication cost because tokens and their assigned experts likely reside in different machines. In this paper, we propose \emph{Gating Dropout}, which allows tokens to ignore the gating network and stay at their local machines, thus reducing the cross-machine communication. Similar to traditional dropout, we also show that Gating Dropout has a regularization effect during training, resulting in improved generalization performance. We validate the effectiveness of Gating Dropout on multilingual machine translation tasks. Our results demonstrate that Gating Dropout improves a state-of-the-art MoE model with faster wall-clock time convergence rates and better BLEU scores for a variety of model sizes and datasets.


title: List-Decodable Sparse Mean Estimation

abstract: Robust mean estimation is one of the most important problems in statistics: given a set of samples in $\mathbb{R}^d$ where an $\alpha$ fraction are drawn from some distribution $D$ and the rest are adversarially corrupted, we aim to estimate the mean of $D$. A surge of recent research interest has been focusing on the list-decodable setting where $\alpha \in (0, \frac12]$, and the goal is to output a finite number of estimates among which at least one approximates the target mean. In this paper, we consider that the underlying distribution $D$ is Gaussian with $k$-sparse mean. Our main contribution is the first polynomial-time algorithm that enjoys sample complexity $O\big(\mathrm{poly}(k, \log d)\big)$, i.e. poly-logarithmic in the dimension. One of our core algorithmic ingredients is using low-degree sparse polynomials to filter outliers, which may find more applications.


title: Object-wise Masked Autoencoders for Fast Pre-training

abstract: Self-supervised pre-training for images without labels has recently achieved promising performance in image classification. The success of transformer-based methods, ViT and MAE, draws the community's attention to the design of backbone architecture and self-supervised task. In this work, we show that current masked image encoding models learn the underlying relationship between all objects in the whole scene, instead of a single object representation. Therefore, those methods bring a lot of compute time for self-supervised pre-training. To solve this issue, we introduce a novel object selection and division strategy to drop non-object patches for learning object-wise representations by selective reconstruction with interested region masks. We refer to this method ObjMAE. Extensive experiments on four commonly-used datasets demonstrate the effectiveness of our model in reducing the compute cost by 72% while achieving competitive performance. Furthermore, we investigate the inter-object and intra-object relationship and find that the latter is crucial for self-supervised pre-training.


title: Spectral Evidence for Unidirectional Charge Density Wave in Detwinned BaNi$_2$As$_2$

abstract: The emergence of unconventional superconductivity in proximity to intertwined electronic orders is especially relevant in the case of iron-based superconductors. Such order consists of an electronic nematic order and a spin density wave in these systems. BaNi$_2$As$_2$, like its well-known iron-based analog BaFe$_2$As$_2$, also hosts a symmetry-breaking structural transition that is coupled to a unidirectional charge density wave (CDW), providing a novel platform to study intertwined orders. Here, through a systematic angle-resolved photoemission spectroscopy study combined with a detwinning $B_1g$ uniaxial strain, we identify distinct spectral evidence of band evolution due to the structural transition as well as CDW-induced band folding. In contrast to the nematicity and spin density wave in BaFe$_2$As$_2$, the structural and CDW order parameters in BaNi$_2$As$_2$ are observed to be strongly coupled and do not separate in the presence of uniaxial strain. Our measurements point to a likely lattice origin of the CDW in BaNi$_2$As$_2$.


title: Insights from an Industrial Collaborative Assembly Project: Lessons in Research and Collaboration

abstract: Significant progress in robotics reveals new opportunities to advance manufacturing. Next-generation industrial automation will require both integration of distinct robotic technologies and their application to challenging industrial environments. This paper presents lessons from a collaborative assembly project between three academic research groups and an industry partner. The goal of the project is to develop a flexible, safe, and productive manufacturing cell for sub-centimeter precision assembly. Solving this problem in a high-mix, low-volume production line motivates multiple research thrusts in robotics. This work identifies new directions in collaborative robotics for industrial applications and offers insight toward strengthening collaborations between institutions in academia and industry on the development of new technologies.


title: Full waveform inversion by model extension: theory, design and optimization

abstract: We describe a new method, full waveform inversion by model extension (FWIME) that recovers accurate acoustic subsurface velocity models from seismic data, when conventional methods fail. We leverage the advantageous convergence properties of wave-equation migration velocity analysis (WEMVA) with the accuracy and high-resolution nature of acoustic full waveform inversion (FWI) by combining them into a robust mathematically-consistent workflow with minimal need for user inputs. The novelty of FWIME resides in the design of a new cost function using the variable projection method, and a novel optimization strategy to combine the two techniques, making our approach more efficient and powerful than applying them sequentially. We observe that FWIME mitigates the need for accurate initial models and low-frequency long-offset data, which can be challenging to acquire. We generate three cycle-skipped 2D synthetic datasets, each containing only one type of wave (transmitted, reflected, refracted), and we analyze how FWIME successfully recovers accurate solutions with the same procedure for all three cases. In a second paper, we apply FWIME to challenging realistic examples where we simultaneously invert all wave modes.


title: Helicity-independent all-optical switching of magnetization in ferrimagnetic alloys

abstract: We review and discuss the process of single-shot helicity-independent all-optical switching of magnetization by which a single suitably-ultrafast excitation, under the right conditions, toggles magnetization from one stable state to another. For almost a decade, this phenomenon was only consistently observed in specific rare-earth-transition-metal ferrimagnetic alloys of GdFeCo, but breakthrough experiments in recent years have revealed that the same behavior can be achieved in a wide range of multi-sublattice magnets including TbCo alloys doped with minute amounts of Gd, Gd/Co and Tb/Co synthetic ferrimagnets, and the rare-earth-free Heusler alloy Mn$_2$Ru$_x$Ga. Aiming to resolve the conditions that allow switching, a series of experiments have shown that the process in the ferrimagnetic alloys GdFeCo and Mn$_2$Ru$_x$Ga is highly sensitive to the pulse duration, starting temperature and the alloy composition. We argue here that the switching displayed by these two very different ferrimagnetic alloys can be generally understood within a single phenomenological framework describing the flow of angular momentum between the constituent sublattices and from the sublattices to the environment. The conditions that facilitate switching stem from the properties of these channels of angular momentum flow in combination with the size of the angular momentum reservoirs. We conclude with providing an outlook in this vibrant research field, with emphasis on the outstanding open questions pertaining to the underlying physics along with noting the advances in exploiting this switching process in technological applications.


title: Varieties of Lazy Magmas Characterized by Forbidden Substructure Theorems

abstract: A magma (or groupoid) is a set with a binary operation $(A,f)$. Roughly speaking, a magma is said to be lazy if compositions such as $f(x,f(f(y,z),u))$ depend on at most two variables. Recently, Kaprinai, Machida and Waldhauser described the lattice of all the varieties of lazy groupoids. A forbidden structure theorem is one that charcaterizes a smaller class $A$ inside a larger class $B$ as all the elements in $B$ that avoid some substructures. For example, a lattice is distributive (smaller class $A$) if and only if it is a lattice (larger class $B$) and avoids the pentagon and the diamond. In this paper we provide a characterization of all pairs of lazy groupoid varieties $A\le B$ by forbidden substructure theorems. Some of the results are straightforward, but some other are very involved. All of these results and proofs were found using a computational tool that proves theorems of this type (for many different classes of relational algebras) and that we make available to every mathematician.


title: Data-Driven Evolutionary Multi-Objective Optimization Based on Multiple-Gradient Descent for Disconnected Pareto Fronts

abstract: Data-driven evolutionary multi-objective optimization (EMO) has been recognized as an effective approach for multi-objective optimization problems with expensive objective functions. The current research is mainly developed for problems with a 'regular' triangle-like Pareto-optimal front (PF), whereas the performance can significantly deteriorate when the PF consists of disconnected segments. Furthermore, the offspring reproduction in the current data-driven EMO does not fully leverage the latent information of the surrogate model. Bearing these considerations in mind, this paper proposes a data-driven EMO algorithm based on multiple-gradient descent. By leveraging the regularity information provided by the up-to-date surrogate model, it is able to progressively probe a set of well distributed candidate solutions with a convergence guarantee. In addition, its infill criterion recommends a batch of promising candidate solutions to conduct expensive objective function evaluations. Experiments on $33$ benchmark test problem instances with disconnected PFs fully demonstrate the effectiveness of our proposed method against four selected peer algorithms.


title: Reinforcement Learning for Branch-and-Bound Optimisation using Retrospective Trajectories

abstract: Combinatorial optimisation problems framed as mixed integer linear programmes (MILPs) are ubiquitous across a range of real-world applications. The canonical branch-and-bound algorithm seeks to exactly solve MILPs by constructing a search tree of increasingly constrained sub-problems. In practice, its solving time performance is dependent on heuristics, such as the choice of the next variable to constrain ('branching'). Recently, machine learning (ML) has emerged as a promising paradigm for branching. However, prior works have struggled to apply reinforcement learning (RL), citing sparse rewards, difficult exploration, and partial observability as significant challenges. Instead, leading ML methodologies resort to approximating high quality handcrafted heuristics with imitation learning (IL), which precludes the discovery of novel policies and requires expensive data labelling. In this work, we propose retro branching; a simple yet effective approach to RL for branching. By retrospectively deconstructing the search tree into multiple paths each contained within a sub-tree, we enable the agent to learn from shorter trajectories with more predictable next states. In experiments on four combinatorial tasks, our approach enables learning-to-branch without any expert guidance or pre-training. We outperform the current state-of-the-art RL branching algorithm by 3-5x and come within 20% of the best IL method's performance on MILPs with 500 constraints and 1000 variables, with ablations verifying that our retrospectively constructed trajectories are essential to achieving these results.


title: Magnetic collapse in Fe$_3$Se$_4$ under high pressure

abstract: Electronic structure and magnetic properties of Fe$_3$Se$_4$ are calculated using the density functional approach. Due to the metallic properties, magnetic moments of the iron atoms in two nonequivalent positions in the unit cell are different from ionic values for Fe$^{3+}$ and Fe$^{2+}$ and are equal to $M_1=2.071 \mu_B$ and $M_2=-2.042 \mu_B$, making the system ferrimagnetic. The total magnetic moment for the unit cell is $2.135 \mu_B$. Under isotropic compression, the total magnetic moment decreases non-monotonically and correlates with the non-monotonic dependence of the density of states at the Fermi level $N(E_F)$. For 7% compression, the magnetic order changes from the ferrimagnetic to the ferromagnetic. At 14% compression, the magnetic order disappears and the total magnetic moment becomes zero, leaving the system in a paramagnetic state. This compression corresponds to the pressure of 114 GPa. The magnetic ordering changes faster upon application of an isotropic external pressure due to the sizeable anisotropy of the chemical bondings in Fe$_3$Se$_4$. The ferrimagnetic and paramagnetic states occur under pressures of 5.0 and 8.0 GPa, respectively. The system remains in the metallic state for all values of compression.


title: Estimation of 3D Body Shape and Clothing Measurements from Frontal- and Side-view Images

abstract: The estimation of 3D human body shape and clothing measurements is crucial for virtual try-on and size recommendation problems in the fashion industry but has always been a challenging problem due to several conditions, such as lack of publicly available realistic datasets, ambiguity in multiple camera resolutions, and the undefinable human shape space. Existing works proposed various solutions to these problems but could not succeed in the industry adaptation because of complexity and restrictions. To solve the complexity and challenges, in this paper, we propose a simple yet effective architecture to estimate both shape and measures from frontal- and side-view images. We utilize silhouette segmentation from the two multi-view images and implement an auto-encoder network to learn low-dimensional features from segmented silhouettes. Then, we adopt a kernel-based regularized regression module to estimate the body shape and measurements. The experimental results show that the proposed method provides competitive results on the synthetic dataset, NOMO-3d-400-scans Dataset, and RGB Images of humans captured in different cameras.


title: Exponential mixing and limit theorems of quasi-periodically forced 2D stochastic Navier-Stokes Equations in the hypoelliptic setting

abstract: We consider the incompressible 2D Navier-Stokes equations on the torus driven by a deterministic time quasi-periodic force and a noise that is white in time and degenerate in Fourier space. We show that the asymptotic statistical behavior is characterized by a quasi-periodic invariant measure that exponentially attracts the law of all solutions. The result is true for any value of the viscosity $\nu>0$ and does not depend on the strength of the external forces. By utilizing this quasi-periodic invariant measure, we establish a quantitative version of the strong law of large numbers and central limit theorem for the continuous time inhomogeneous solution processes with explicit convergence rates. It turns out that the convergence rate in the central limit theorem depends on the time inhomogeneity through the Diophantine approximation property on the quasi-periodic frequency of the quasi-periodic force.


title: Do We Really Need to Use Constraint Violation in Constrained Evolutionary Multi-Objective Optimization?

abstract: Constraint violation has been a building block to design evolutionary multi-objective optimization algorithms for solving constrained multi-objective optimization problems. However, it is not uncommon that the constraint violation is hardly approachable in real-world black-box optimization scenarios. It is unclear that whether the existing constrained evolutionary multi-objective optimization algorithms, whose environmental selection mechanism are built upon the constraint violation, can still work or not when the formulations of the constraint functions are unknown. Bearing this consideration in mind, this paper picks up four widely used constrained evolutionary multi-objective optimization algorithms as the baseline and develop the corresponding variants that replace the constraint violation by a crisp value. From our experiments on both synthetic and real-world benchmark test problems, we find that the performance of the selected algorithms have not been significantly influenced when the constraint violation is not used to guide the environmental selection.


title: Norm inflation for a non-linear heat equation with Gaussian initial conditions

abstract: We consider a non-linear heat equation $\partial_t u = \Delta u + B(u,Du)+P(u)$ posed on the $d$-dimensional torus, where $P$ is a polynomial of degree at most $3$ and $B$ is a bilinear map that is not a total derivative. We show that, if the initial condition $u_0$ is taken from a sequence of smooth Gaussian fields with a specified covariance, then $u$ exhibits norm inflation with high probability. A consequence of this result is that there exists no Banach space of distributions which carries the Gaussian free field on the 3D torus and to which the DeTurck-Yang-Mills heat flow extends continuously, which complements recent well-posedness results in arXiv:2111.10652 and arXiv:2201.03487. Another consequence is that the (deterministic) non-linear heat equation exhibits norm inflation, and is thus locally ill-posed, at every point in the Besov space $B^{-1/2}_{\infty,\infty}$; the space $B^{-1/2}_{\infty,\infty}$ is an endpoint since the equation is locally well-posed for $B^{\eta}_{\infty,\infty}$ for every $\eta>-\frac12$.


title: Modified L\'evy Laplacian on manifold and Yang-Mills instantons

abstract: An infinite dimensional Laplacian defined as the Ces\'aro mean of the second order directional derivatives on manifold is considered. This Laplacian is parameterized by the choice of a curve in the group of orthogonal rotations. It is shown that, under certain conditions on the curve, this operator is related to instantons on a 4-dimensional manifold.


title: Travelling Salesman Problem: Parallel Implementations & Analysis

abstract: The Traveling Salesman Problem (often called TSP) is a classic algorithmic problem in the field of computer science and operations research. It is an NP-Hard problem focused on optimization. TSP has several applications even in its purest formulation, such as planning, logistics, and the manufacture of microchips; and can be slightly modified to appear as a sub-problem in many areas, such as DNA sequencing. In this paper, a study on parallelization of the Brute Force approach (under several paradigms) of the Travelling Salesman Problem is presented. Detailed timing studies for the serial and various parallel implementations of the Travelling Salesman Problem have also been illustrated.


title: A macroscopic delayed-choice quantum eraser using a commercial laser

abstract: The heart of quantum mechanics is quantum superposition between orthogonal bases of a single particle. In the particle nature of quantum mechanics, quantum superposition is represented by probability amplitudes between mutually exclusive natures such as orthogonal polarization bases. The delayed-choice quantum eraser is for the post-determination of the photon nature, raising the cause-effect relation issue. Over the last several decades, quantum erasers have been intensively studied using nearly all kinds of photons. Here, the macroscopic delayed-choice quantum eraser is experimentally demonstrated using a continuous wave laser and discussed for quantum superposition in a macroscopic regime. For this, a noninterfering Mach-Zehnder interferometer composed of two polarizing beam splitters is chosen to manipulate polarization bases of lights and to measure them in a delayed-choice manner via polarization-basis projection.


title: Multi-Task Learning with Multi-Query Transformer for Dense Prediction

abstract: Previous multi-task dense prediction studies developed complex pipelines such as multi-modal distillations in multiple stages or searching for task relational contexts for each task. The core insight beyond these methods is to maximize the mutual effects of each task. Inspired by the recent query-based Transformers, we propose a simple pipeline named Multi-Query Transformer (MQTransformer) that is equipped with multiple queries from different tasks to facilitate the reasoning among multiple tasks and simplify the cross-task interaction pipeline. Instead of modeling the dense per-pixel context among different tasks, we seek a task-specific proxy to perform cross-task reasoning via multiple queries where each query encodes the task-related context. The MQTransformer is composed of three key components: shared encoder, cross-task query attention module and shared decoder. We first model each task with a task-relevant query. Then both the task-specific feature output by the feature extractor and the task-relevant query are fed into the shared encoder, thus encoding the task-relevant query from the task-specific feature. Secondly, we design a cross-task query attention module to reason the dependencies among multiple task-relevant queries; this enables the module to only focus on the query-level interaction. Finally, we use a shared decoder to gradually refine the image features with the reasoned query features from different tasks. Extensive experiment results on two dense prediction datasets (NYUD-v2 and PASCAL-Context) show that the proposed method is an effective approach and achieves state-of-the-art results. Code and models are available at https://github.com/yangyangxu0/MQTransformer.


title: Comparison of limit shapes for Bernoulli first-passage percolation

abstract: We consider Bernoulli first-passage percolation on the $d$-dimensional hypercubic lattice with $d \geq 2$. The passage time of edge $e$ is $0$ with probability $p$ and $1$ with probability $1-p$, independently of each other. Let $p_c$ be the critical probability for percolation of edges with passage time $0$. When $0\leq p<p_c$, there exists a nonrandom, nonempty compact convex set $\mathcal{B}_p$ such that the set of vertices to which the first-passage time from the origin is within $t$ is well-approximated by $t\mathcal{B}_p$ for all large $t$, with probability one. The aim of this paper is to prove that for $0\leq p<q<p_c$, the Hausdorff distance between $\mathcal{B}_p$ and $\mathcal{B}_q$ grows linearly in $q-p$. Moreover, we mention that the approach taken in the paper provides a lower bound for the expected size of the intersection of geodesics, that gives a nontrivial consequence for the \textit{critical} case.


title: Differences between Lyapunov exponents for the simple random walk in Bernoulli potentials

abstract: We consider the simple random walk on the $d$-dimensional lattice $\mathbb{Z}^d$ ($d \geq 1$), traveling in potentials which are Bernoulli distributed. The so-called Lyapunov exponent describes the cost of traveling for the simple random walk in the potential, and it is known that the Lyapunov exponent is strictly monotone in the parameter of the Bernoulli distribution. Hence, the aim of this paper is to investigate the effect of the potential on the Lyapunov exponent more precisely, and we derive some Lipschitz-type estimates for the difference between the Lyapunov exponents.


title: Theoretical modeling of catalytic self-propulsion

abstract: Self-propelling particles or microswimmers have opened a new field of investigation with both fundamental and practical perspectives. They represent very convenient model objects for experimental studies of active matter, and have implications in nano-robotics, drug delivery, and more. Here, we summarize recent advances in theoretical description of the self-propulsion of catalytic microswimmers that non-uniformly release ions, including its physical origins and the current switch in focus to non-linear effects, geometric tuning and more. In particular, we show that the ionic self-propulsion always includes both diffusiophoretic and electrophoretic contributions, and that non-linear effects are physical causes of a number of intriguing phenomena, such as the reverse in the direction of the particle motion in response to variations of the salt concentration or self-propulsion of electro-neutral particles. We finally suggest several remaining theoretical challenges in the field.


title: Fair Labeled Clustering

abstract: Numerous algorithms have been produced for the fundamental problem of clustering under many different notions of fairness. Perhaps the most common family of notions currently studied is group fairness, in which proportional group representation is ensured in every cluster. We extend this direction by considering the downstream application of clustering and how group fairness should be ensured for such a setting. Specifically, we consider a common setting in which a decision-maker runs a clustering algorithm, inspects the center of each cluster, and decides an appropriate outcome (label) for its corresponding cluster. In hiring for example, there could be two outcomes, positive (hire) or negative (reject), and each cluster would be assigned one of these two outcomes. To ensure group fairness in such a setting, we would desire proportional group representation in every label but not necessarily in every cluster as is done in group fair clustering. We provide algorithms for such problems and show that in contrast to their NP-hard counterparts in group fair clustering, they permit efficient solutions. We also consider a well-motivated alternative setting where the decision-maker is free to assign labels to the clusters regardless of the centers' positions in the metric space. We show that this setting exhibits interesting transitions from computationally hard to easy according to additional constraints on the problem. Moreover, when the constraint parameters take on natural values we show a randomized algorithm for this setting that always achieves an optimal clustering and satisfies the fairness constraints in expectation. Finally, we run experiments on real world datasets that validate the effectiveness of our algorithms.


title: Large order breathers of the nonlinear Schr\"odinger equation

abstract: Multi-soliton and high-order soliton solutions are two type of famous ones in the integrable focusing nonlinear Schr\"odinger equation. The dynamics of multi-soliton was well known to us since 70s of the last century by the determinant analysis. However, there is few progress on the high-order solitons. In this work, we would like to analyze the large order asymptotics for the high-order breathers, which are special cases of double high-order solitons with the same velocity to the nonlinear Schr\"odinger equation. To analyze the large order dynamics, we first convert the representation of Darboux transformation into a framework of Riemann-Hilbert problem. Then we show that there exist five distinct asymptotic regions by the Deift-Zhou nonlinear steepest descent method. It is very interesting that a novel genus-three asymptotic region is first found in the large order asymptotics to large high-order breathers, which enriches the dynamic behaviors in the field of large order solitons. All results to the asymptotic analysis are verified by the numerical method.


title: A discrete analogue of Terrell's characterization of rectangular distributions

abstract: George R. Terrell (1983, {Ann. Probab., vol. 11(3), pp. 823--826) showed that the Pearson coefficient of correlation of an ordered pair from a random sample of size two is at most one-half, and the equality is attained only for rectangular (uniform over some interval) distributions. In the present note it is proved that the same is true for the discrete case, in the sense that the correlation coefficient attains its maximal value only for discrete rectangular (uniform over some finite lattice) distributions. MSC: Primary 60E15; 62E10; Secondary 62G30. Key words and phrases: discrete rectangular distribution; order statistics; Hahn polynomials; Pearson coefficient of correlation.


title: Boosting Facial Expression Recognition by A Semi-Supervised Progressive Teacher

abstract: In this paper, we aim to improve the performance of in-the-wild Facial Expression Recognition (FER) by exploiting semi-supervised learning. Large-scale labeled data and deep learning methods have greatly improved the performance of image recognition. However, the performance of FER is still not ideal due to the lack of training data and incorrect annotations (e.g., label noises). Among existing in-the-wild FER datasets, reliable ones contain insufficient data to train robust deep models while large-scale ones are annotated in lower quality. To address this problem, we propose a semi-supervised learning algorithm named Progressive Teacher (PT) to utilize reliable FER datasets as well as large-scale unlabeled expression images for effective training. On the one hand, PT introduces semi-supervised learning method to relieve the shortage of data in FER. On the other hand, it selects useful labeled training samples automatically and progressively to alleviate label noise. PT uses selected clean labeled data for computing the supervised classification loss and unlabeled data for unsupervised consistency loss. Experiments on widely-used databases RAF-DB and FERPlus validate the effectiveness of our method, which achieves state-of-the-art performance with accuracy of 89.57% on RAF-DB. Additionally, when the synthetic noise rate reaches even 30%, the performance of our PT algorithm only degrades by 4.37%.


title: Variations of Milnor's triple linking number

abstract: Topological polymers have various topological types, and they are expressed by graphs. However, the Jones polynomial, we have a difficulty to compute it; computational time is growing exponentially with respect to the crossing number. The simplest Vassiliev invariant is the linking number and thus we will seek a next simple one is as the Milnor's triple linking number. In this paper, we introduce simple Gauss diagram formulas of Vassiliev invariants of Milnor type. These are non-torsion valued, whereas the base-point-free Milnor's triple linking number is usually torsion-valued.


title: Flickering returns as RS Oph reestablishes quiescent conditions following its 2021 nova outburst

abstract: RS Oph has persistently displayed flickering at optical wavelengths when observed away from its repeating nova outbursts. During the 2006 eruption the flickering disappeared, and this repeated during the recent 2021 event. We have been monitoring RS Oph looking for the reappearance of flickering at B-band following the 2021 outburst. The flickering was still absent (sigma(B)<0.002 mag) on day +210 (counted from nova optical maximum), appeared at sigma(B)=0.008 mag on day +224, and raised to sigma(B)=0.029 mag on day +250. On following dates the amplitude remained large, although fluctuating. The recovery of B-band quiescence brightness by RS Oph begun around day +225 and was completed by day +260. The parallel patterns followed by the rise in system brightness and the reappearance of flickering confirm the central role played in RS Oph by the return to pre-outburst conditions of the accretion disk and the refilling by the RG wind of the immediate circumstellar space.


title: Stationary broken parity states in active matter models

abstract: We demonstrate that several nonvariational continuum models commonly used to describe active matter as well as other active systems exhibit nongeneric behavior: each model supports asymmetric but stationary localized states even in the absence of pinning at heterogeneities. Moreover, such states only begin to drift following a drift-transcritical bifurcation as the activity increases. Asymmetric stationary states should only exist in variational systems, i.e., in models with gradient structure. In other words, such states are expected in passive systems, but not in active systems where the gradient structure of the model is broken by activity. We identify a "spurious" gradient dynamics structure of these models that is responsible for this nongeneric behavior, and determine the types of additional terms that render the models generic, i.e., with asymmetric states that appear via drift-pitchfork bifurcations and are generically moving. We provide detailed illustrations of our results using numerical continuation of resting and steadily drifting states in both generic and nongeneric cases.


title: Granular Generalized Variable Precision Rough Sets and Rational Approximations

abstract: Rational approximations are introduced and studied in granular graded rough sets and generalizations thereof by the first author in recent research papers. The concept of rationality is determined by related ontologies and coherence between granularity, mereology and approximations in the context. In addition, a framework for rational approximations is introduced by her in the mentioned paper(s). Granular approximations constructed as per the procedures of variable precision rough sets (VPRS) are likely to be more rational than those constructed from a classical perspective under certain conditions. This may continue to hold for some generalizations of the former. However, a formal characterization of such conditions is not available in the previously published literature. In this research, theoretical aspects of the problem are critically examined, uniform generalizations of granular VPRS are introduced, new connections with granular graded rough sets are proved, appropriate concepts of substantial parthood are introduced, their extent of compatibility with the framework is accessed, and the framework is extended. Basic assumptions are explained in detail, and additional examples are constructed for readability. Furthermore, meta applications to cluster validation, image segmentation and dynamic sorting are invented. Extensions to direct generalizations of VPRS such as probabilistic rough sets are a natural consequence of the work.


title: Topology Optimization of Graded Truss Lattices Based on On-the-Fly Homogenization

abstract: We introduce a computational framework for the topology optimization of cellular structures with spatially varying architecture, which is applied to functionally graded truss lattices under quasistatic loading. We make use of a first-order homogenization approach, which replaces the discrete truss by an effective continuum description to be treated by finite elements in a macroscale boundary value problem. By defining the local truss architecture through a set of Bravais vectors, we formulate the optimization problem with regards to the spatially varying basis vectors and demonstrate its feasibility and performance through a series of benchmark problems in 2D (though the method is sufficiently general to also apply in 3D, as discussed). Both the displacement field and the topology are continuously varying unknown fields on the macroscale, and a regularization is included for well-posedness. We argue that prior solutions obtained from aligning trusses along the directions of principal stresses are included as a special case. The outlined approach results in heterogeneous truss architectures with a smoothly varying unit cell, enabling easy fabrication with a tunable length scale (the latter avoiding the ill-posedness stemming from classical nonconvex methods without an intrinsic length scale).


title: Atomic-scale Deformation Process of Glasses Unveiled by Stress-induced Structural Anisotropy

abstract: Experimentally resolving atomic-scale structural changes of a deformed glass remains challenging owing to the disordered nature of glass structure. Here, we show that the structural anisotropy emerges as a general hallmark for different types of glasses (metallic glasses, oxide glass, amorphous selenium, and polymer glass) after thermo-mechanical deformation, and it is highly correlates with local nonaffine atomic displacements detected by the high-energy X-ray diffraction technique. By analyzing the anisotropic pair density function, we unveil the atomic-level mechanism responsible for the plastic flow, which notably differs between metallic glasses and covalent glasses. The structural rearrangements in metallic glasses are mediated through cutting and formation of atomic bonds, which occurs in some localized inelastic regions embedded in elastic matrix, whereas that of covalent glasses is mediated through the rotation of atomic bonds or chains without bond length change, which occurs in a less localized manner.


title: Going Deeper into Permutation-Sensitive Graph Neural Networks

abstract: The invariance to permutations of the adjacency matrix, i.e., graph isomorphism, is an overarching requirement for Graph Neural Networks (GNNs). Conventionally, this prerequisite can be satisfied by the invariant operations over node permutations when aggregating messages. However, such an invariant manner may ignore the relationships among neighboring nodes, thereby hindering the expressivity of GNNs. In this work, we devise an efficient permutation-sensitive aggregation mechanism via permutation groups, capturing pairwise correlations between neighboring nodes. We prove that our approach is strictly more powerful than the 2-dimensional Weisfeiler-Lehman (2-WL) graph isomorphism test and not less powerful than the 3-WL test. Moreover, we prove that our approach achieves the linear sampling complexity. Comprehensive experiments on multiple synthetic and real-world datasets demonstrate the superiority of our model.


title: Study of Time Evolution of Thermal and Non-Thermal Emission from an M-Class Solar Flare

abstract: We conduct a wide-band X-ray spectral analysis in the energy range of 1.5-100 keV to study the time evolution of the M7.6 class flare of 2016 July 23, with the Miniature X-ray Solar Spectrometer (MinXSS) CubeSat and the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) spacecraft. With the combination of MinXSS for soft X-rays and RHESSI for hard X-rays, a non-thermal component and three-temperature multi-thermal component -- "cool" ($T \approx$ 3 MK), "hot" ($T \approx$ 15 MK), and "super-hot" ($T \approx$ 30 MK) -- were measured simultaneously. In addition, we successfully obtained the spectral evolution of the multi-thermal and non-thermal components with a 10 s cadence, which corresponds to the Alfv\'en time scale in the solar corona. We find that the emission measures of the cool and hot thermal components are drastically increasing more than hundreds of times and the super-hot thermal component is gradually appearing after the peak of the non-thermal emission. We also study the microwave spectra obtained by the Nobeyama Radio Polarimeters (NoRP), and we find that there is continuous gyro-synchrotron emission from mildly relativistic non-thermal electrons. In addition, we conducted a differential emission measure (DEM) analysis by using Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) and determine that the DEM of cool plasma increases within the flaring loop. We find that the cool and hot plasma components are associated with chromospheric evaporation. The super-hot plasma component could be explained by the thermalization of the non-thermal electrons trapped in the flaring loop.

categories: astro-ph.SR astro-ph.HE physics.plasm-ph physics.space-ph


title: SPH simulations of high-speed collisions between asteroids and comets

abstract: We studied impact processes by means of smoothed-particle hydrodynamics (SPH) simulations. The method was applied to modeling formation of main-belt families during the cometary bombardment (either early or late, ${\sim}\,3.85\,{\rm Gy}$ ago). If asteroids were bombarded by comets, as predicted by the Nice model, hundreds of asteroid families (catastrophic disruptions of diameter $D \ge 100\,{\rm km}$ bodies) should have been created, but the observed number is only 20. Therefore we computed a standard set of 125 simulations of collisions between representative $D = 100\,{\rm km}$ asteroids and high-speed icy projectiles (comets). According to our results, the largest remnant mass $M_{\rm lr}$ is similar as in low-speed collisions, due to appropriate scaling with the effective strength $Q_{\rm eff}$, but the largest fragment mass $M_{\rm lf}$ exhibits systematic differences - it is typically smaller for craterings and bigger for super-catastrophic events. This trend does not, however, explain the non-existence of old families. The respective parametric relations can be used in other statistical (Monte-Carlo) models to better understand collisions between asteroidal and cometary populations.


title: GLITCH: Automated Polyglot Security Smell Detection in Infrastructure as Code

abstract: Infrastructure as Code (IaC) is the process of managing IT infrastructure via programmable configuration files (also called IaC scripts). Like other software artifacts, IaC scripts may contain security smells, which are coding patterns that can result in security weaknesses. Automated analysis tools to detect security smells in IaC scripts exist, but they focus on specific technologies such as Puppet, Ansible, or Chef. This means that when the detection of a new smell is implemented in one of the tools, it is not immediately available for the technologies supported by the other tools -- the only option is to duplicate the effort. This paper presents an approach that enables consistent security smell detection across different IaC technologies. We conduct a large-scale empirical study that analyzes security smells on three large datasets containing 196,755 IaC scripts and 12,281,251 LOC. We show that all categories of security smells are identified across all datasets and we identify some smells that might affect many IaC projects. To conduct this study, we developed GLITCH, a new technology-agnostic framework that enables automated polyglot smell detection by transforming IaC scripts into an intermediate representation, on which different security smell detectors can be defined. GLITCH currently supports the detection of nine different security smells in scripts written in Ansible, Chef, or Puppet. We compare GLITCH with state-of-the-art security smell detectors. The results obtained not only show that GLITCH can reduce the effort of writing security smell analyses for multiple IaC technologies, but also that it has higher precision and recall than the current state-of-the-art tools.


title: On stationary inflection points in step responses

abstract: The step and impulse responses of a proper, rational transfer function are well-behaved analytic functions. We prove that such a response cannot have an inflection point such that the tangent at that point is parallel to the time-axis. Hence when a step or impulse response of a finite-dimensional LTI system crosses any given level, that crossing must be transversal, and can never be tangential.


title: Curvature-induced superconductivity enhancement for ultra-thin superconducting films

abstract: We derive the linearized Ginzburg-Landau (GL) equation for an ultra-thin superconducting film with curvature in a magnetic field. By introducing a novel transverse order parameter that varies slowly along the film, and applying the superconducting/vacuum boundary condition, we decouple the linearized GL equation into a transverse part and a surface part that includes the superconducting geometric potential (GP). The nucleation of the superconducting state in curved thin superconducting films can be equivalently described by the surface part equation. In the equivalent GL free energy of a curved superconducting film, the superconducting GP enables the film to remain in the superconducting state even when the superconducting parameter $ \alpha $ turns positive by further reducing the quadratic term of the order parameter. Furthermore, we numerically investigate the phase transition of a rectangle thin superconducting film bent around a cylindrical surface. Our numerical results show that the superconducting GP enhances the superconductivity of the curved film by weakening the effect of the magnetic field, and the increase of the critical temperature, in units of the bulk critical temperature, is equal to the product of the negative superconducting GP and the square of the zero-temperature coherence length, which agrees with our theoretical predictions.


title: Syntax-Guided Program Reduction for Understanding Neural Code Intelligence Models

abstract: Neural code intelligence (CI) models are opaque black-boxes and offer little insight on the features they use in making predictions. This opacity may lead to distrust in their prediction and hamper their wider adoption in safety-critical applications. Recently, input program reduction techniques have been proposed to identify key features in the input programs to improve the transparency of CI models. However, this approach is syntax-unaware and does not consider the grammar of the programming language. In this paper, we apply a syntax-guided program reduction technique that considers the grammar of the input programs during reduction. Our experiments on multiple models across different types of input programs show that the syntax-guided program reduction technique is faster and provides smaller sets of key tokens in reduced programs. We also show that the key tokens could be used in generating adversarial examples for up to 65% of the input programs.


title: WaveMix: A Resource-efficient Neural Network for Image Analysis

abstract: We propose a novel neural architecture for computer vision -- WaveMix -- that is resource-efficient and yet generalizable and scalable. While using fewer trainable parameters, GPU RAM, and computations, WaveMix networks achieve comparable or better accuracy than the state-of-the-art convolutional neural networks, vision transformers, and token mixers for several tasks. This efficiency can translate to savings in time, cost, and energy. To achieve these gains we used multi-level two-dimensional discrete wavelet transform (2D-DWT) in WaveMix blocks, which has the following advantages: (1) It reorganizes spatial information based on three strong image priors -- scale-invariance, shift-invariance, and sparseness of edges -- (2) in a lossless manner without adding parameters, (3) while also reducing the spatial sizes of feature maps, which reduces the memory and time required for forward and backward passes, and (4) expanding the receptive field faster than convolutions do. The whole architecture is a stack of self-similar and resolution-preserving WaveMix blocks, which allows architectural flexibility for various tasks and levels of resource availability. WaveMix establishes new benchmarks for segmentation on Cityscapes; and for classification on Galaxy 10 DECals, Places-365, five EMNIST datasets, and iNAT-mini and performs competitively on other benchmarks. Our code and trained models are publicly available.


title: One-Sided Repeated-Root Two-Dimensional Cyclic and Constacyclic Codes

abstract: In this paper, we study some repeated-root two-dimensional cyclic and constacyclic codes over a finite field $F=\mathbb{F}_q$. We obtain the generator matrices and generator polynomials of these codes and their duals. We also investigate when such codes are self-dual. Moreover, we prove that if there exists an asymptotically good family of one-sided repeated-root two-dimensional cyclic or constacyclic codes, then there exists an asymptotically good family of simple root two-dimensional cyclic or constacyclic codes with parameters at least as good as the first family. Furthermore, we show that several of the main results of the papers Rajabi and Khashyarmanesh (2018) and Sepasdar and Khashyarmanesh (2016) are not accurate and find other conditions needed for them to hold.


title: Enhancing Quality of Pose-varied Face Restoration with Local Weak Feature Sensing and GAN Prior

abstract: Facial semantic guidance (including facial landmarks, facial heatmaps, and facial parsing maps) and facial generative adversarial networks (GAN) prior have been widely used in blind face restoration (BFR) in recent years. Although existing BFR methods have achieved good performance in ordinary cases, these solutions have limited resilience when applied to face images with serious degradation and pose-varied (e.g., looking right, looking left, laughing, etc.) in real-world scenarios. In this work, we propose a well-designed blind face restoration network with generative facial prior. The proposed network is mainly comprised of an asymmetric codec and a StyleGAN2 prior network. In the asymmetric codec, we adopt a mixed multi-path residual block (MMRB) to gradually extract weak texture features of input images, which can better preserve the original facial features and avoid excessive fantasy. The MMRB can also be plug-and-play in other networks. Furthermore, thanks to the affluent and diverse facial priors of the StyleGAN2 model, we adopt it as the primary generator network in our proposed method and specially design a novel self-supervised training strategy to fit the distribution closer to the target and flexibly restore natural and realistic facial details. Extensive experiments on synthetic and real-world datasets demonstrate that our model performs superior to the prior art for face restoration and face super-resolution tasks.


title: Possible Systematic Rotation in the Mature Stellar Population of a $z=9.1$ Galaxy

abstract: We present new observations with the Atacama Large Millimeter/submillimeter Array for a gravitationally-lensed galaxy at $z=9.1$, MACS1149-JD1. [O III] 88-$\mu$m emission is detected at 10$\sigma$ with a spatial resolution of $\sim0.3$ kpc in the source plane, enabling the most distant morpho-kinematic study of a galaxy. The [O III] emission is distributed smoothly without any resolved clumps and shows a clear velocity gradient with $\Delta V_{\rm obs}/2\sigma_{\rm tot}=0.84\pm0.23$, where $\Delta V_{\rm obs}$ is the observed maximum velocity difference and $\sigma_{\rm tot}$ is the velocity dispersion measured in the spatially-integrated line profile, suggesting a rotating system. Assuming a geometrically thin self-gravitating rotation disk model, we obtain $V_{\rm rot}/\sigma_V=0.67_{-0.26}^{+0.73}$, where $V_{\rm rot}$ and $\sigma_V$ are the rotation velocity and velocity dispersion, respectively, still consistent with rotation. The resulting disk mass of $0.65_{-0.40}^{+1.37}\times10^{9}$ M$_\odot$ is consistent with being associated with the stellar mass identified with a 300 Myr-old stellar population independently indicated by a Balmer break in the spectral energy distribution. We conclude that the most of the dynamical mass is associated with the previously-identified mature stellar population that formed at $z\sim15$.


title: Signal strength and W-boson mass measurements as a probe of the electro-weak phase transition at colliders -- Snowmass White Paper

abstract: We consider an extension of the scalar sector of the Standard Model (SM) by an additional gauge singlet, which mixes with a part of the SM-like Higgs doublet. Within this model, parameter-space regions exist that can lead to a strong first-order electro-weak phase transition, a necessary condition for electro-weak baryogenesis. We discuss how such regions of the parameter space can be tested using the SM-like Higgs boson's signal strength measurements, as well as precision observables, such as e.g. the W-boson mass, at current and future colliders.


title: Deep Deconfounded Content-based Tag Recommendation for UGC with Causal Intervention

abstract: Traditional content-based tag recommender systems directly learn the association between user-generated content (UGC) and tags based on collected UGC-tag pairs. However, since a UGC uploader simultaneously creates the UGC and selects the corresponding tags, her personal preference inevitably biases the tag selections, which prevents these recommenders from learning the causal influence of UGCs' content features on tags. In this paper, we propose a deep deconfounded content-based tag recommender system, namely, DecTag, to address the above issues. We first establish a causal graph to represent the relations among uploader, UGC, and tag, where the uploaders are identified as confounders that spuriously correlate UGC and tag selections. Specifically, to eliminate the confounding bias, causal intervention is conducted on the UGC node in the graph via backdoor adjustment, where uploaders' influence on tags leaked through backdoor paths can be eliminated for causal effect estimation. Observing that adjusting the causal graph with do-calculus requires integrating the entire uploader space, which is infeasible, we design a novel Monte Carlo (MC)-based estimator with bootstrap, which can achieve asymptotic unbiasedness provided that uploaders for the collected UGCs are i.i.d. samples from the population. In addition, the MC estimator has the intuition of substituting the biased uploaders with a hypothetical random uploader from the population in the training phase, where deconfounding can be achieved in an interpretable manner. Finally, we establish a YT-8M-Causal dataset based on the widely used YouTube-8M dataset with causal intervention and propose an evaluation strategy accordingly to unbiasedly evaluate causal tag recommenders. Extensive experiments show that DecTag is more robust to confounding bias than state-of-the-art causal recommenders.


title: Finite element simulation for validation of multi-dipole line cusp magnetic field configuration for MPD

abstract: A Multi-dipole line cusp configured Plasma Device (MPD) having six electromagnets with embedded Vacoflux-50 as a core material has been operated with a capability to experimentally control the field-free region, the radial profile of magnetic field, and pole magnetic field by changing magnet current. For the validation of multi-dipole line cusp magnetic field (MMF) configuration in a 3-D geometry, a finite element simulation has been performed using COMSOL software. This paper presents 3-D magnetic field simulation results of multi-pole line cusp magnetic field configuration performed over the 1.2m length and 40cm diameter chamber in the vacuum condition. The simulation results show good agreement with the experimentally measured magnetic field profile. The performed magnetic field simulation results clearly capture that this configuration has full control over a null region (nearly field-free region) as well is capable to change the magnetic field values and radial variation of the magnetic field. Moreover, the magnetic field profiles over the end cross-section of the device have been discussed.


title: Study of Asymptotic Velocity in the Bondi-Hoyle Accretion Flows in the Domain of Kerr and 4-D Einstein-Gauss-Bonnet Gravities

abstract: Understanding the physical structures of the accreated matter very close to the black hole in quasars and active galactic nucleus (AGNs) is an important milestone to constrain the activities occurring in their centers. In this paper, we numerically investigate the effects of the asymptotic velocities on the physical structures of the accretion disk around the Kerr and Einstein-Gauss-Bonnet (EGB) rapidly rotating black holes. The Bondi-Hoyle accretion is considered with a falling gas towards the black hole in upstream region of the computational domain. The shock cones are naturally produced in the downstream part of the flow around both black holes. It is found that the structure of the cones and the amount of the accreated matter depend on asymptotic velocity $V_{\infty}$ (Mach number) and the types of the gravities (Kerr or EGB). Increasing the Mach number of the inflowing matter in the supersonic region causes the shock opening angle and accretion rates getting smaller because of the rapidly falling gas towards the black hole. The EGB gravity leads to an increase in the shock opening angle of the shock cones while the mass accretion rates $\dot{M}$ are decreasing in EGB gravity with a Gauss-Bonnet (GB) coupling constant $\alpha$. It is also confirmed that accretion rates and drag forces are significantly altered in the EGB gravity. Our numerical simulation results could be used to identify the accreation mechanism and physical properties of the accretion disk and black hole in the observed $X-$ rays such as NGC $1313$ $X-1$ and $1313$ $X-2$ and MAXI $J1803-298$.


title: CP violation in non-leptonic B decays as a portal to New Physics

abstract: CP violation offers a powerful probe for testing the Standard Model. In this endeavour, non-leptonic decays of B mesons play a key role. I will discuss benchmark decays and puzzling patterns emerging from theoretical analyses of the current data.


title: Hydrodynamic theory of two-dimensional incompressible polar active fluids with quenched and annealed disorder

abstract: We study the moving phase of two-dimensional (2D) incompressible polar active fluids in the presence of both quenched and annealed disorder. We show that long-range polar order persists even in this defect-ridden two-dimensional system. We obtain the large-distance, long-time scaling laws of the velocity fluctuations using three distinct dynamic renormalization group schemes. These are an uncontrolled one-loop calculation in exactly two dimensions, and two $d=(d_c-\epsilon)$-expansions to $O(\epsilon)$, obtained by two different analytic continuations of our 2D model to higher spatial dimensions: a ``hard" continuation which has $d_c={7\over 3}$, and a ``soft" continuation with $d_c={5\over 2}$. Surprisingly, the quenched and annealed parts of the velocity correlation function have the same anisotropy exponent and the relaxational and propagating parts of the dispersion relation have the same dynamic exponent in the nonlinear theory even though they are distinct in the linearized theory. This is due to anomalous hydrodynamics. Furthermore, all three renormalization schemes yield very similar values for the universal exponents, and, therefore, we expect the numerical values we predict for them to be highly accurate.


title: Isometric direct limits of bidual Banach spaces

abstract: Sequences of $n$-th order bidual Banach spaces, called tower systems and their direct and inverse limits are considered. Motivated by recent applications in uniform algebras, we introduce two functors: $\textrm{Dir}$ and $\textrm{Inv}$ assigning to Banach spaces (and to bounded linear operators) some new Banach spaces and operators. Of particular interest is the enormous "tower space" built over the space of continuous functions. We prove that the action of $\textrm{Dir}$ preserves direct sum decompositions. This functor preserves also spectra of operators, their Fredholmness and compactness properties. An application of these functors to the problem of location of supports of representing measures for function algebras is outlined in the last section.


title: Pressure-induced mixed states caused by spin-elastic interactions during first-order spin phase transition in spin crossover compounds

abstract: Recently, the possibility of exploiting the phenomenon of spin transition (ST) has been intensively investigated, therefore, it is particularly important to study the behavior of ST under various stimuli. Here, the shape and content of the intermediate phase of ST in Hoffmann-like compounds [Fe(Fpz)2M(CN)4](M = Pt, Pd) under external stimuli are studied. For this purpose, magnetic and Raman spectroscopy measurements were carried out. In pressure-induced spin transition (PIST), a mixture of high-spin and low-spin states appears, while in temperature-induced spin transition (TIST), a homogeneous state occurs. The first-order ST induced by pressure has a hysteresis, but is not abrupt. Whereas, the temperature-induced spin transition at ambient pressure is hysteretic and abrupt. To investigate this difference, we discuss using a thermodynamic model that considers elastic interactions, showing that the slope of the hysteresis loop is related to the appearance of internal pressure, which is related to the difference in sample compressibility under high spin and low spin states.


title: Regulating Matching Markets with Constraints: Data-driven Taxation

abstract: This paper develops a framework to conduct a counterfactual analysis to regulate matching markets with regional constraints that impose lower and upper bounds on the number of matches in each region. Our work is motivated by the Japan Residency Matching Program, in which the policymaker wants to guarantee the least number of doctors working in rural regions to achieve the minimum standard of service. Among the multiple possible policies that satisfy such constraints, a policymaker wants to choose the best. To this end, we develop a discrete choice model approach that estimates the utility functions of agents from observed data and predicts agents' behavior under different counterfactual policies. Our framework also allows the policymaker to design the welfare-maximizing tax scheme, which outperforms the policy currently used in practice. Furthermore, a numerical experiment illustrates how our method works.


title: Schauder regularity results in separable Hilbert spaces

abstract: We prove Schauder type estimates for solutions of stationary and evolution equations driven by weak generators of transition semigroups associated to a semilinear stochastic partial differential equations with values in a separable Hilbert space.


title: Quantum coherence on selectivity and transport of ion channels

abstract: Recently, it has been suggested that ion channel selectivity filter may exhibit quantum coherence, which may be appropriate to explain ion selection and conduction processes. Potassium channels play a vital role in many physiological processes. One of their main physiological functions is the efficient and highly selective transfer of K+ ions through the membranes into the cells. To do this, ion channels must be highly selective, allowing only certain ions to pass through the membrane, while preventing the others. The present research is an attempt to investigate the relationship between hopping rate and maintaining coherence in ion channels. Using the Lindblad equation to describe a three-level system, the results in different quantum regimes are examined. We studied the distillable coherence and the second order coherence function of the system. The oscillation of distillable coherence from zero, after the decoherence time, and also the behavior of the coherence function clearly show the point that the system is coherent in ion channels with high throughput rates.


title: Topological phase estimation method for reparameterized periodic functions

abstract: We consider a signal composed of several periods of a periodic function, of which we observe a noisy reparametrisation. The phase estimation problem consists of finding that reparametrisation, and, in particular, the number of observed periods. Existing methods are well-suited to the setting where the periodic function is known, or at least, simple. We consider the case when it is unknown and we propose an estimation method based on the shape of the signal. We use the persistent homology of sublevel sets of the signal to capture the temporal structure of its local extrema. We infer the number of periods in the signal by counting points in the persistence diagram and their multiplicities. Using the estimated number of periods, we construct an estimator of the reparametrisation. It is based on counting the number of sufficiently prominent local minima in the signal. This work is motivated by a vehicle positioning problem, on which we evaluated the proposed method.


title: Phonon and optical-roton branches of excitations of the Bose system

abstract: For a system of a large number of Bose particles, a chain of coupled equations for the averages of field operators is obtained. In the approximation where only the averages of one field operator and the averages of products of two operators at zero temperature are taken into account, there is derived a closed system of dynamic equations. Taking into account the finite range of the interaction potential between particles, the spectrum of elementary excitations of a many-particle Bose system is calculated, and it is shown that it has two branches: a sound branch and an optical branch with an energy gap at zero momentum. At high density, both branches are nonmonotonic and have the roton-like minima. The dispersion of the phonon part of the spectrum is considered. The performed calculations and analysis of experiments on neutron scattering allow to make a statement about the complex structure of the Landau dispersion curve in the superfluid He-4.


title: Einstein-Maxwell fields as solutions of higher-order theories

abstract: We study four-dimensional Einstein-Maxwell fields for which any higher-order corrections to the field equations effectively reduces to just a rescaling of the gravitational and the cosmological constant. These configurations are thus simultaneous solutions of (virtually) any modified theory of gravity coupled (possibly non-minimally) to any electrodynamics. In the case of non-null electromagnetic fields we provide a full characterization of such universal solutions, which correspond to a family of gravitational waves propagating in universes of the Levi-Civita--Bertotti--Robinson type. For null fields we first obtain a set of general necessary conditions, and then a full characterization for a special subfamily, which turns out to represent electromagnetic waves accompanied by pure radiation in the (anti-)Nariai background. The results are exemplified for the case of Born-Infeld, ModMax and Horndeski electrodynamics.


title: Relation-Specific Attentions over Entity Mentions for Enhanced Document-Level Relation Extraction

abstract: Compared with traditional sentence-level relation extraction, document-level relation extraction is a more challenging task where an entity in a document may be mentioned multiple times and associated with multiple relations. However, most methods of document-level relation extraction do not distinguish between mention-level features and entity-level features, and just apply simple pooling operation for aggregating mention-level features into entity-level features. As a result, the distinct semantics between the different mentions of an entity are overlooked. To address this problem, we propose RSMAN in this paper which performs selective attentions over different entity mentions with respect to candidate relations. In this manner, the flexible and relation-specific representations of entities are obtained which indeed benefit relation classification. Our extensive experiments upon two benchmark datasets show that our RSMAN can bring significant improvements for some backbone models to achieve state-of-the-art performance, especially when an entity have multiple mentions in the document.


title: Normality and associated primes of Closed neighborhood ideals and dominating ideals

abstract: In this paper, we first give some sufficient criteria for normality of monomial ideals. As applications, we show that closed neighborhood ideals of complete bipartite graphs are normal, and hence satisfy the (strong) persistence property. We also prove that dominating ideals of complete bipartite graphs are nearly normally torsion-free. In addition, we show that dominating ideals of $h$-wheel graphs, under certain condition, are normal.


title: Characterizing Tourist Daily Trip Chains Using Mobile Phone Big Data

abstract: Tourists tend to visit multiple destinations out of their variety-seeking motivations in their trips. Thus, it is critical to discover travel patterns involving multi-destinations in tourism research. Existing relevant research most relied on survey data or focused on citizens due to the lack of large-scale, fine-grained tourism datasets. Several scholars have mentioned the notion of trip chains, but few works have been done towards quantitatively identifying the structures of trip chains. In this paper, we propose a model for quantitatively characterizing tourist daily trip chains. After applying this model to tourist mobile phone big data, underlying tourist travel patterns are discovered. Through the framework, we find that: (1) Most "hybrid" (inter-city and intra-city) and "intra-city" (only intra-city) patterns can be captured by only 13 key trip chains relatively; (2) For two continuous days, almost all kinds of original chains have a rather high probability to transfer to either the first two transferred chains, or other infrequent chains in our study areas; (3) The principle of least efforts (PLE) affects tourists' structures of trip chains. We can use average degree and average travel distance to interpret tourist travel behavior (achieving tasks in PLE). This study not only demonstrate the complex daily travel trip chains from tourism big data, but also fill the gap in tourism literature on multi-destination trips by discovering significant and underlying patterns based on mobile datasets.


title: Deep Learning-based Spatially Explicit Emulation of an Agent-Based Simulator for Pandemic in a City

abstract: Agent-Based Models are very useful for simulation of physical or social processes, such as the spreading of a pandemic in a city. Such models proceed by specifying the behavior of individuals (agents) and their interactions, and parameterizing the process of infection based on such interactions based on the geography and demography of the city. However, such models are computationally very expensive, and the complexity is often linear in the total number of agents. This seriously limits the usage of such models for simulations, which often have to be run hundreds of times for policy planning and even model parameter estimation. An alternative is to develop an emulator, a surrogate model that can predict the Agent-Based Simulator's output based on its initial conditions and parameters. In this paper, we discuss a Deep Learning model based on Dilated Convolutional Neural Network that can emulate such an agent based model with high accuracy. We show that use of this model instead of the original Agent-Based Model provides us major gains in the speed of simulations, allowing much quicker calibration to observations, and more extensive scenario analysis. The models we consider are spatially explicit, as the locations of the infected individuals are simulated instead of the gross counts. Another aspect of our emulation framework is its divide-and-conquer approach that divides the city into several small overlapping blocks and carries out the emulation in them parallelly, after which these results are merged together. This ensures that the same emulator can work for a city of any size, and also provides significant improvement of time complexity of the emulator, compared to the original simulator.


title: On Characteristics of the Range of Integral Operators

abstract: We show that a positive operator between $L^p$-spaces is given by integration against a kernel function if and only if the image of each positive function has a lower semi-continuous representative with respect to a suitable topology. This is a consequence of a new characterization of kernel operators on general Banach lattices as those operators whose range can be represented over a fixed countable set of positive vectors. Similar results are shown to hold for operators that merely dominate a non-trivial kernel operator.


title: Deep neural networks overcome the curse of dimensionality in the numerical approximation of semilinear partial differential equations

abstract: We prove that deep neural networks are capable of approximating solutions of semilinear Kolmogorov PDE in the case of gradient-independent, Lipschitz-continuous nonlinearities, while the required number of parameters in the networks grow at most polynomially in both dimension $d \in \mathbb{N}$ and prescribed reciprocal accuracy $\varepsilon$. Previously, this has only been proven in the case of semilinear heat equations.


title: Incentive Mechanism Design for Emergency Frequency Control in Multi-Infeed Hybrid AC-DC System

abstract: In multi-infeed hybrid AC-DC (MIDC) systems, the emergency frequency control (EFC) with LCC-HVDC systems participating is of vital importance for system frequency stability. Nevertheless, when regional power systems are operated by different decision-makers, the LCC-HVDC systems and their connected AC systems might be unwilling to participate in the EFC due to the costs and losses. In this paper, to incentivize the LCC-HVDC systems and their connected adjacent AC systems to participate in the droop-based EFC, a novel control-parameter-based incentive mechanism is proposed, which can deal with various possible emergency frequency faults. Then, a non-cooperative-based incentive game model is formulated to implement the incentive mechanism in the MIDC system. An algorithm for seeking the Nash equilibrium is designed, and the uniqueness of Nash equilibrium is proven. Moreover, the individual rationality, incentive compatibility and social optimality of the proposed mechanism are analyzed and proven. The effectiveness of the proposed incentive mechanism is verified through a case study.


title: Agent-based Simulation of District-based Elections

abstract: In district-based elections, electors cast votes in their respective districts. In each district, the party with maximum votes wins the corresponding seat in the governing body. The election result is based on the number of seats won by different parties. In this system, locations of electors across the districts may severely affect the election result even if the total number of votes obtained by different parties remains unchanged. A less popular party may end up winning more seats if their supporters are suitably distributed spatially. This happens due to various regional and social influences on individual voters which modulate their voting choice. In this paper, we explore agent-based models for district-based elections, where we consider each elector as an agent, and try to represent their social and geographical attributes and political inclinations using probability distributions. This model can be used to simulate election results by Monte Carlo sampling. The models allow us to explore the full space of possible outcomes of an electoral setting, though they can also be calibrated to actual election results for suitable values of parameters. We use Approximate Bayesian Computation (ABC) framework to estimate model parameters. We show that our model can reproduce the results of elections held in India and USA, and can also produce counterfactual scenarios.