abstract: In this paper we develop the Greedy Recombination Interpolation Method (GRIM) for finding sparse approximations of functions initially given as linear combinations of some (large) number of simpler functions. In a similar spirit to the CoSaMP algorithm, GRIM combines dynamic growth-based interpolation techniques and thinning-based reduction techniques. The dynamic growth-based aspect is a modification of the greedy growth utilised in the Generalised Empirical Interpolation Method (GEIM). A consequence of the modification is that our growth is not restricted to being one-per-step as it is in GEIM. The thinning-based aspect is carried out by recombination, which is the crucial component of the recent ground-breaking convex kernel quadrature method. GRIM provides the first use of recombination outside the setting of reducing the support of a measure. The sparsity of the approximation found by GRIM is controlled by the geometric concentration of the data in a sense that is related to a particular packing number of the data. We apply GRIM to a kernel quadrature task for the radial basis function kernel, and verify that its performance matches that of other contemporary kernel quadrature techniques.


title: Efficient Knowledge Compilation Beyond Weighted Model Counting

abstract: Quantitative extensions of logic programming often require the solution of so called second level inference tasks, i.e., problems that involve a third operation, such as maximization or normalization, on top of addition and multiplication, and thus go beyond the well-known weighted or algebraic model counting setting of probabilistic logic programming under the distribution semantics. We introduce Second Level Algebraic Model Counting (2AMC) as a generic framework for these kinds of problems. As 2AMC is to (algebraic) model counting what forall-exists-SAT is to propositional satisfiability, it is notoriously hard to solve. First level techniques based on Knowledge Compilation (KC) have been adapted for specific 2AMC instances by imposing variable order constraints on the resulting circuit. However, those constraints can severely increase the circuit size and thus decrease the efficiency of such approaches. We show that we can exploit the logical structure of a 2AMC problem to omit parts of these constraints, thus limiting the negative effect. Furthermore, we introduce and implement a strategy to generate a sufficient set of constraints statically, with a priori guarantees for the performance of KC. Our empirical evaluation on several benchmarks and tasks confirms that our theoretical results can translate into more efficient solving in practice. Under consideration for acceptance in TPLP.


title: Selecting a complete sample of blazars in sub-millimetre catalogues

abstract: The \textit{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS), that has covered about 642 sq. deg. in 5 bands from 100 to 500 $\mu\rm m$, allows a blind flux-limited selection of blazars at sub-mm wavelengths. However, blazars constitute a tiny fraction of H-ATLAS sources and therefore identifying them is not a trivial task. Using the data on known blazars detected by the H-ATLAS we have defined a locus for 500\,$\mu$m selected blazars and exploited it to select blazar candidates in the H-ATLAS fields. Candidates and known blazars in the H-ATLAS equatorial and South Galactic Pole fields were followed up with the Australia Telescope Compact Array (ATCA) or with the Karl G. Jansky Very Large Array (VLA), and matched with existing radio- and mm-catalogues to reconstruct the spectral behaviour over at least 6 orders of magnitude in frequency. We identified a selection approach that, combining the information in the sub-mm and radio domains, efficiently singles out genuine blazars. In this way, we identified a sample of 39 blazars brighter than $S_{500\mu\rm m} = 35\,$mJy in the H-ATLAS fields. Tests made cross-matching the H-ATLAS catalogues with large catalogues of blazar candidates indicate that the sample is complete. The derived counts are compared with model predictions finding good consistency with the C2Ex model and with estimates based on ALMA data.


title: On density of $Z_3$-flow-critical graphs

abstract: For an abelian group $\Gamma$, a graph $G$ is said to be $\Gamma$-flow-critical if $G$ does not admit a nowhere-zero $\Gamma$-flow, but for each edge $e\in E(G)$, the contraction $G/e$ has a nowhere-zero $\Gamma$-flow. A bound on the density of $Z_3$-flow-critical graphs drawn on a fixed surface is obtained, generalizing the planar case of the bound on the density of 4-critical graphs by Kostochka and Yancey.


title: Mitigating Hidden Confounding Effects for Causal Recommendation

abstract: Recommender systems suffer from confounding biases when there exist confounders affecting both item features and user feedback (e.g., like or not). Existing causal recommendation methods typically assume confounders are fully observed and measured, forgoing the possible existence of hidden confounders in real applications. For instance, product quality is a confounder since affecting both item prices and user ratings, but is hidden for the third-party e-commerce platform due to the difficulty of large-scale quality inspection; ignoring it could result in the bias effect of over-recommending high-price items. This work analyzes and addresses the problem from a causal perspective. The key lies in modeling the causal effect of item features on a user's feedback. To mitigate hidden confounding effects, it is compulsory but challenging to estimate the causal effect without measuring the confounder. Towards this goal, we propose a Hidden Confounder Removal (HCR) framework that leverages front-door adjustment to decompose the causal effect into two partial effects, according to the mediators between item features and user feedback. The partial effects are independent from the hidden confounder and identifiable. During training, HCR performs multi-task learning to infer the partial effects from historical interactions. We instantiate HCR for two scenarios and conduct experiments on three real-world datasets. Empirical results show that the HCR framework provides more accurate recommendations, especially for less-active users. We will release the code once accepted.


title: Computing Bend-Minimum Orthogonal Drawings of Plane Series-Parallel Graphs in Linear Time

abstract: A planar orthogonal drawing of a planar 4-graph G (i.e., a planar graph with vertex-degree at most four) is a crossing-free drawing that maps each vertex of G to a distinct point of the plane and each edge of $G$ to a sequence of horizontal and vertical segments between its end-points. A longstanding open question in Graph Drawing, dating back over 30 years, is whether there exists a linear-time algorithm to compute an orthogonal drawing of a plane 4-graph with the minimum number of bends. The term "plane" indicates that the input graph comes together with a planar embedding, which must be preserved by the drawing (i.e., the drawing must have the same set of faces as the input graph). In this paper, we positively answer the question above for the widely-studied class of series-parallel graphs. Our linear-time algorithm is based on a characterization of the planar series-parallel graphs that admit an orthogonal drawing without bends. This characterization is given in terms of the orthogonal spirality that each type of triconnected component of the graph can take; the orthogonal spirality of a component measures how much that component is "rolled-up" in an orthogonal drawing of the graph.


title: Factors influencing the energy gap in topological states of antiferromagnetic MnBi$_2$Te$_4$

abstract: The experimentally measured angle-resolved photoemission dispersion maps for MnBi$_{2}$Te$_{4}$ samples, which show different energy gaps at the Dirac point (DP), are compared with the results of theoretical calculations to find the conditions for the best agreement between theory and experiment. We have analyzed different factors which influence the Dirac gap width: (i) the surface van der Waals (SvdW) distance between the first and second septuple layers (SLs), (ii) the magnetic moment on Mn atoms, (iii) the spin-orbit coupling (SOC) strength for the surface Te and Bi atoms and related changes in the localization of the topological surface states (TSSs). It was shown that all these factors may change the gap width at the DP in a wide range from 5 to $\sim$90~meV. We show that the Dirac gap variation is mainly determined by the corresponding changes in the TSSs spatial distribution. The best agreement between the presented experimental data (with the Dirac gaps between $\sim$15 and 55~meV) and the calculations takes place for a slightly compressed SvdW interval (of about -3.5~\% compared to the bulk value) with modified SOC for surface atoms (that can occur in the presence of various defects in the near-surface region). We show that upon changing the values of the SvdW interval and surface SOC strength the TSSs spatial distribution shifts between the SLs with opposite magnetizations, which leads to a non-monotonic change in the Dirac gap size.


title: KGRGRL: A User's Permission Reasoning Method Based on Knowledge Graph Reward Guidance Reinforcement Learning

abstract: In general, multiple domain cyberspace security assessments can be implemented by reasoning user's permissions. However, while existing methods include some information from the physical and social domains, they do not provide a comprehensive representation of cyberspace. Existing reasoning methods are also based on expert-given rules, resulting in inefficiency and a low degree of intelligence. To address this challenge, we create a Knowledge Graph (KG) of multiple domain cyberspace in order to provide a standard semantic description of the multiple domain cyberspace. Following that, we proposed a user's permissions reasoning method based on reinforcement learning. All permissions in cyberspace are represented as nodes, and an agent is trained to find all permissions that user can have according to user's initial permissions and cyberspace KG. We set 10 reward setting rules based on the features of cyberspace KG in the reinforcement learning of reward information setting, so that the agent can better locate user's all permissions and avoid blindly finding user's permissions. The results of the experiments showed that the proposed method can successfully reason about user's permissions and increase the intelligence level of the user's permissions reasoning method. At the same time, the F1 value of the proposed method is 6% greater than that of the Translating Embedding (TransE) method.


title: Morse functions and contact convex surfaces

abstract: Let $f$ be a Morse function on a closed surface $\Sigma$ such that zero is a regular value and such that $f$ admits neither positive minima nor negative maxima. In this expository note, we show that $\Sigma\times \mathbb{R}$ admits an $\mathbb{R}$-invariant contact form $\alpha=fdt+\beta$ whose characteristic foliation along the zero section is (negative) weakly gradient-like with respect to $f$. The proof is self-contained and gives explicit constructions of any $\mathbb{R}$-invariant contact structure in $\Sigma \times \mathbb{R}$, up to isotopy. As an application, we give an alternative geometric proof of the homotopy classification of $\mathbb{R}$-invariant contact structures in terms of their dividing set.


title: Large eddy simulations of reacting and non-reacting transcritical fuel sprays using multiphase thermodynamics

abstract: Accurate simulations of high-pressure transcritical fuel sprays are essential for the design and optimization of next-generation gas turbines, internal combustion engines, and liquid propellant rocket engines. Most important and challenging is the accurate modelling of complex real-gas effects in high-pressure environments, especially the hybrid subcritical-to-supercritical mode of evaporation during the mixing of fuel and oxidizer. In this paper, we present a novel modeling framework for high-fidelity simulations of reacting and non-reacting transcritical fuel sprays. In this method, the high-pressure jet disintegration is modeled using a diffuse interface method with multiphase thermodynamics, which combines multi-component real-fluid kinetic and caloric state equations with vapor-liquid equilibrium calculations in order to compute thermodynamic properties of the mixture at transcritical pressures. All multiphase thermodynamic formulations are presented for general cubic state equations coupled with a rapid phase-equilibrium calculation method. The proposed method represents multiphase turbulent fluid flows at transcritical pressures without relying on any semi-empirical break-up and evaporation models. Combustion source terms are evaluated using a finite-rate chemistry model, including real-gas effects based on the fugacity of the species in the mixture. The adaptive local deconvolution method (ALDM) is used as a physically consistent turbulence model for large-eddy simulation (LES). LES results show a very good agreement with available experimental data for the reacting and non-reacting ECN Spray A at transcritical operating conditions.

categories: physics.flu-dyn nlin.AO physics.app-ph physics.comp-ph


title: Cyanopolyyne line survey towards high-mass star-forming regions with TMRT

abstract: We carried out a cyanopolyyne line survey towards a large sample of HMSFRs using the Shanghai Tian Ma 65m Radio Telescope (TMRT). Our sample consisted of 123 targets taken from the TMRT C band line survey. It included three kinds of sources, namely those with detection of the 6.7 GHz CH3OH maser alone, with detection of the radio recombination line (RRL) alone, and with detection of both (hereafter referred to as Maser-only, RRL-only, and Maser-RRL sources, respectively). We detected HC3N in 38 sources, HC5N in 11 sources, and HC7N in G24.790+0.084, with the highest detection rate being found for Maser-RRL sources and a very low detection rate found for RRL-only sources. Their column densities were derived using the rotational temperature measured from the NH3 lines. And we constructed and fitted the far-infrared (FIR) spectral energy distributions. Based on these, we derive their dust temperatures, H2 column densities, and abundances of cyanopolyynes relative to H2. The detection rate, the column density, and the relative abundance of HC3N increase from Maser-only to Maser-RRL sources and decrease from Maser-RRL to RRL-only sources. This trend is consistent with the proposed evolutionary trend of HC3N under the assumption that our Maser-only, Maser-RRL, and RRL-only sources correspond to massive young stellar objects, ultra-compact HII regions, and normal classical HII regions, respectively. Furthermore, a statistical analysis of the integrated line intensity and column density of HC3N and shock-tracing molecules (SiO, H2CO) enabled us to find positive correlations between them. This suggests that HC3N may be another tracer of shocks, and should therefore be the subject of further observations and corresponding chemical simulations. Our results indirectly support the idea that the neutral--neutral reaction between C2H2 and CN is the dominant formation pathway of HC3N.


title: Enhanced $2\pi$-periodic Aharonov-Bohm Effect as a Signature of Majorana Zero Modes Probed by Nonlocal Measurements

abstract: We propose the $2\pi$-periodic Aharonov-Bohm (AB) effect as a nonlocal probe of Majorana zero modes (MZMs) without the restriction of fermion parity. We demonstrate the enhancement of the AB effect, where the topological protection of MZMs yields amplified and robust Andreev reflection mediated by MZMs at multiple superconductor-normal metal junctions. We investigate the influence of trivial bound states and show that a nonlocal index enables a more explicit distinction between the trivial and topological bound states than local probes.


title: Packet Switching in Quantum Networks: A Path to Quantum Internet

abstract: Large-scale quantum networks with thousands of nodes require scalable network protocols and physical hardware to realize. In this work, we introduce packet switching as a new paradigm for quantum data transmission in both future and near-term quantum networks. We propose a classical-quantum data frame structure and explore methods of frame generation and processing. Further, we present conceptual designs for a quantum reconfigurable optical add-drop multiplexer to realize the proposed transmission scheme. Packet switching allows for a universal design for a next generation Internet where classical and quantum data share the same network protocols and infrastructure. In this new quantum networking paradigm, entanglement distribution, as with quantum key distribution, is an application built on top of the quantum network rather than as a network designed especially for those purposes. For analysis of the network model, we simulate the feasibility of quantum packet switching for some preliminary models of quantum key and entanglement distribution. Finally, we discuss how our model can be integrated with other network models toward a realization of a quantum Internet.


title: Graphene/SOI-based self-powered Schottky barrier photodiode array

abstract: We have fabricated 4-element Graphene/Silicon on Insulator (SOI) based Schottky barrier photodiode array ()PDA and investigated its optoelectronic device performance. In our device design, monolayer graphene is utilized as common electrode on lithographically defined linear array of n-type Si channels on SOI substrate. As revealed by wavelength resolved photocurrent spectroscopy measurements, each element in the PDA structure exhibited a maximum spectral responsivity comparable to that of commercially available photodiodes operating under self-powered operational mode. Time-dependent photocurrent spectroscopy measurements showed excellent photocurrent reversibility of the device. The study presented here is expected to offer exciting opportunities in terms of high value-added graphene/Si based PDA device applications such as multi-wavelength light measurement, level metering, high-speed photometry, position/motion detection, and more.


title: Classification of simple Harish-Chandra modules over the Ovsienko-Roger superalgebra

abstract: With the $\Omega$-operators for the Virasoro algebra \cite{BF} and the super Virasoro algebra in \cite{CL, CLL}, we get the $\Omega$-operators for the Ovsienko-Roger superalgebras in this paper and then use it to classify all simple cuspidal modules for the $\bZ$-graded and $\frac12\bZ$-graded Ovsienko-Roger superalgebras. By this result, we can easily classify all simple Harish-Chandra modules over some related Lie superalgebras, including the $N=1$ BMS$_3$ algebra, the super $W(2,2)$, etc.


title: Crowdsourced Hypothesis Generation and their Verification: A Case Study on Sleep Quality Improvement

abstract: A clinical study is often necessary for exploring important research questions; however, this approach is sometimes time and money consuming. Another extreme approach, which is to collect and aggregate opinions from crowds, provides a result drawn from the crowds' past experiences and knowledge. To explore a solution that takes advantage of both the rigid clinical approach and the crowds' opinion-based approach, we design a framework that exploits crowdsourcing as a part of the research process, whereby crowd workers serve as if they were a scientist conducting a "pseudo" prospective study. This study evaluates the feasibility of the proposed framework to generate hypotheses on a specified topic and verify them in the real world by employing many crowd workers. The framework comprises two phases of crowd-based workflow. In Phase 1 - the hypothesis generation and ranking phase - our system asks workers two types of questions to collect a number of hypotheses and rank them. In Phase 2 - the hypothesis verification phase - the system asks workers to verify the top-ranked hypotheses from Phase 1 by implementing one of them in real life. Through experiments, we explore the potential and limitations of the framework to generate and evaluate hypotheses about the factors that result in a good night's sleep. Our results on significant sleep quality improvement show the basic feasibility of our framework, suggesting that crowd-based research is compatible with experts' knowledge in a certain domain.


title: Tempestuous life beyond R500: X-ray view on the Coma cluster with SRG/eROSITA. II. Shock & Relic

abstract: This is the second paper in a series of studies of the Coma cluster using the SRG/eROSITA X-ray data obtained during the calibration and performance verification phase of the mission. Here, we focus on the region adjacent to the radio source 1253+275 (radio relic, RR, hereafter). We show that the X-ray surface brightness exhibits its steepest gradient at $\sim 79'$ ($\sim 2.2\,{\rm Mpc}\approx R_{200c}$), which is almost co-spatial to the outer edge of the RR. As in the case of several other relics, the Mach number of the shock derived from the X-ray surface brightness profile ($M_X\approx 1.9$) appears to be lower than needed to explain the slope of the integrated radio spectrum in the diffusive shock acceleration (DSA) model ($M_R\approx 3.5$) if the magnetic field is uniform and the radiative losses are fast. However, the shock geometry is plausibly much more complicated than a spherical wedge centered on the cluster, given the non-trivial correlation between radio, X-ray, and SZ images. While the complicated shock geometry alone might cause a negative bias in $M_X$, we speculate on a few other possibilities that may affect the $M_X$-$M_R$ relation, including the shock substructure that might be modified by the presence of non-thermal filaments stretching across the shock and the propagation of relativistic electrons along the non-thermal filaments with a strong magnetic field. We also discuss the "history" of the radio galaxy NGC4789, which is located ahead of the relic in the context of the Coma-NGC4839 merger scenario.


title: Concavity property of minimal $L^{2}$ integrals with Lebesgue measurable gain VII -- Negligible weights

abstract: In this article, we present characterizations of the concavity property of minimal $L^2$ integrals with negligible weights degenerating to linearity on the fibrations over open Riemann surfaces and the fibrations over products of open Riemann surfaces. As applications, we obtain characterizations of the holding of equality in optimal jets $L^2$ extension problem with negligible weights on the fibrations over open Riemann surfaces and the fibrations over products of open Riemann surfaces.


title: Optimization of Sparse Sensor Placement for Estimation of Wind Direction and Surface Pressure Distribution Using Time-Averaged Pressure-Sensitive Paint Data on Automobile Model

abstract: This study proposes a method for predicting the wind direction against the simple automobile model (Ahmed model) and the surface pressure distributions on it by using data-driven optimized sparse pressure sensors. Positions of sparse pressure sensor pairs on the Ahmed model were selected for estimation of the yaw angle and reconstruction of pressure distributions based on the time-averaged surface pressure distributions database of various yaw angles, whereas the symmetric sensors in the left and right sides of the model were assumed. The surface pressure distributions were obtained by pressure-sensitive paint measurements. Three algorithms for sparse sensor selection based on the greedy algorithm were applied, and the sensor positions were optimized. The sensor positions and estimation accuracy of yaw angle and pressure distributions of three algorithms were compared and evaluated. The results show that a few optimized sensors can accurately predict the yaw angle and the pressure distributions.


title: Residual Local Feature Network for Efficient Super-Resolution

abstract: Deep learning based approaches has achieved great performance in single image super-resolution (SISR). However, recent advances in efficient super-resolution focus on reducing the number of parameters and FLOPs, and they aggregate more powerful features by improving feature utilization through complex layer connection strategies. These structures may not be necessary to achieve higher running speed, which makes them difficult to be deployed to resource-constrained devices. In this work, we propose a novel Residual Local Feature Network (RLFN). The main idea is using three convolutional layers for residual local feature learning to simplify feature aggregation, which achieves a good trade-off between model performance and inference time. Moreover, we revisit the popular contrastive loss and observe that the selection of intermediate features of its feature extractor has great influence on the performance. Besides, we propose a novel multi-stage warm-start training strategy. In each stage, the pre-trained weights from previous stages are utilized to improve the model performance. Combined with the improved contrastive loss and training strategy, the proposed RLFN outperforms all the state-of-the-art efficient image SR models in terms of runtime while maintaining both PSNR and SSIM for SR. In addition, we won the first place in the runtime track of the NTIRE 2022 efficient super-resolution challenge. Code will be available at https://github.com/fyan111/RLFN.


title: There is no universal acceleration scale in galaxies

abstract: Recently, many studies seem to reveal the existence of some correlations between dark matter and baryonic matter. In particular, the unexpected tight Radial Acceleration Relation (RAR) discovered in rotating galaxies has caught much attention. The RAR suggests the existence of a universal and fundamental acceleration scale in galaxies, which seems to challenge the $\Lambda$CDM model and favour some modified gravity theories. A large debate about whether RAR is compatible with the $\Lambda$CDM model has arisen. Here, by analysing the high quality velocity dispersion profiles of 13 E0-type elliptical galaxies in the SDSS-IV MaNGA sample and assuming a power-law function of radius $r$ for the 3-dimensional velocity dispersion in each galaxy, we report the RAR for E0-type elliptical galaxies and we show that the resultant RAR has more than $5\sigma$ deviations from the RAR in late-type galaxies. This new RAR provides an independent probe to falsify the existence of any universal acceleration scale in galaxies. Our result significantly challenges those modified gravity theories which suggest the existence of any universal acceleration scale.


title: The use of deep learning in interventional radiotherapy (brachytherapy): a review with a focus on open source and open data

abstract: Deep learning advanced to one of the most important technologies in almost all medical fields. Especially in areas, related to medical imaging it plays a big role. However, in interventional radiotherapy (brachytherapy) deep learning is still in an early phase. In this review, first, we investigated and scrutinised the role of deep learning in all processes of interventional radiotherapy and directly related fields. Additionally we summarised the most recent developments. To reproduce results of deep learning algorithms both source code and training data must be available. Therefore, a second focus of this work was on the analysis of the availability of open source, open data and open models. In our analysis, we were able to show that deep learning plays already a major role in some areas of interventional radiotherapy, but is still hardly presented in others. Nevertheless, its impact is increasing with the years, partly self-propelled but also influenced by closely related fields. Open source, data and models are growing in number but are still scarce and unevenly distributed among different research groups. The reluctance in publishing code, data and models limits reproducibility and restricts evaluation to mono-institutional datasets. Summarised, deep learning will change positively the workflow of interventional radiotherapy but there is room for improvement when it comes to reproducible results and standardised evaluation methods.


title: Fairness in Participatory Budgeting via Equality of Resources

abstract: We introduce a family of normative principles to assess fairness in the context of participatory budgeting. These principles are based on the fundamental idea that budget allocations should be fair in terms of the resources invested into meeting the wishes of individual voters. This is in contrast to earlier proposals that are based on specific assumptions regarding the satisfaction of voters with a given budget allocation. We analyse these new principles in axiomatic, algorithmic, and experimental terms.


title: Learning-Based Orchestration for Dynamic Functional Split and Resource Allocation in vRANs

abstract: One of the key benefits of virtualized radio access networks (vRANs) is network management flexibility. However, this versatility raises previously-unseen network management challenges. In this paper, a learning-based zero-touch vRAN orchestration framework (LOFV) is proposed to jointly select the functional splits and allocate the virtualized resources to minimize the long-term management cost. First, testbed measurements of the behaviour between the users' demand and the virtualized resource utilization are collected using a centralized RAN system. The collected data reveals that there are non-linear and non-monotonic relationships between demand and resource utilization. Then, a comprehensive cost model is proposed that takes resource overprovisioning, declined demand, instantiation and reconfiguration into account. Moreover, the proposed cost model also captures different routing and computing costs for each split. Motivated by our measurement insights and cost model, LOFV is developed using a model-free reinforcement learning paradigm. The proposed solution is constructed from a combination of deep Q-learning and a regression-based neural network that maps the network state and users' demand into split and resource control decisions. Our numerical evaluations show that LOFV can offer cost savings by up to 69\% of the optimal static policy and 45\% of the optimal fully dynamic policy.


title: Fair Shares: Feasibility, Domination and Incentives

abstract: We consider fair allocation of a set $M$ of indivisible goods to $n$ equally-entitled agents, with no monetary transfers. Every agent $i$ has a valuation $v_i$ from some given class of valuation functions. A share $s$ is a function that maps a pair $(v_i,n)$ to a value, with the interpretation that if an allocation of $M$ to $n$ agents fails to give agent $i$ a bundle of value at least equal to $s(v_i,n)$, this serves as evidence that the allocation is not fair towards $i$. For such an interpretation to make sense, we would like the share to be feasible, meaning that for any valuations in the class, there is an allocation that gives every agent at least her share. The maximin share was a natural candidate for a feasible share for additive valuations. However, Kurokawa, Procaccia and Wang [2018] show that it is not feasible. We initiate a systematic study of the family of feasible shares. We say that a share is \emph{self maximizing} if truth-telling maximizes the implied guarantee. We show that every feasible share is dominated by some self-maximizing and feasible share. We seek to identify those self-maximizing feasible shares that are polynomial time computable, and offer the highest share values. We show that a SM-dominating feasible share -- one that dominates every self-maximizing (SM) feasible share -- does not exist for additive valuations (and beyond). Consequently, we relax the domination property to that of domination up to a multiplicative factor of $\rho$ (called $\rho$-dominating). For additive valuations we present shares that are feasible, self-maximizing and polynomial-time computable. For $n$ agents we present such a share that is $\frac{2n}{3n-1}$-dominating. For two agents we present such a share that is $(1 - \epsilon)$-dominating. Moreover, for these shares we present poly-time algorithms that compute allocations that give every agent at least her share.


title: Probe the color screening in proton-nucleus collisions with complex potentials

abstract: Color screening and parton inelastic scattering modify the heavy-quark antiquark potential in the medium that consists of particles from quantum chromodynamics (QCD), leading to suppression of quarkonium production in relativistic heavy-ion collisions. Due to small charm/anti-charm ($c\bar{c}$) pair production number in proton-nucleus (pA) collisions, the correlation between different $c\bar{c}$ pairs is negligible, which makes the Schr\"odinger equation viable for tracking the evolution of only one $c\bar{c}$ pair. We employ the time-dependent Schr\"odinger equation with in-medium $c\bar{c}$ potential to study the evolution of charmonium wave functions in the hydrodynamic like QCD medium produced in pA collisions. We explore different parametrizations of real and imaginary parts of $c\bar{c}$ potential and calculate the nuclear modification factors ($R_{\rm pA}$) of $J/\psi$ and $\psi(2S)$ in $\sqrt{s_{NN}}=5.02$ TeV energy p-Pb collisions at Large Hadron Collider (LHC). Comparing a strong and a weak screening scenario with experimental data in this approach, we arrive at the conclusion that the color screening is weak at temperature close to deconfined phase transition. Moreover, the imaginary part of the potential is crucial to describe the experimental data which is consistent with widely studied semi-classical approaches where the dissociation rates are essential.


title: A scalable deep learning approach for solving high-dimensional dynamic optimal transport

abstract: The dynamic formulation of optimal transport has attracted growing interests in scientific computing and machine learning, and its computation requires to solve a PDE-constrained optimization problem. The classical Eulerian discretization based approaches suffer from the curse of dimensionality, which arises from the approximation of high-dimensional velocity field. In this work, we propose a deep learning based method to solve the dynamic optimal transport in high dimensional space. Our method contains three main ingredients: a carefully designed representation of the velocity field, the discretization of the PDE constraint along the characteristics, and the computation of high dimensional integral by Monte Carlo method in each time step. Specifically, in the representation of the velocity field, we apply the classical nodal basis function in time and the deep neural networks in space domain with the H1-norm regularization. This technique promotes the regularity of the velocity field in both time and space such that the discretization along the characteristic remains to be stable during the training process. Extensive numerical examples have been conducted to test the proposed method. Compared to other solvers of optimal transport, our method could give more accurate results in high dimensional cases and has very good scalability with respect to dimension. Finally, we extend our method to more complicated cases such as crowd motion problem.


title: MOA-2019-BLG-008Lb: a new microlensing detection of an object at the planet/brown dwarf boundary

abstract: We report on the observations, analysis and interpretation of the microlensing event MOA-2019- BLG-008. The observed anomaly in the photometric light curve is best described through a binary lens model. In this model, the source did not cross caustics and no finite source effects were observed. Therefore the angular Einstein ring radius cannot be measured from the light curve alone. However, the large event duration, t E about 80 days, allows a precise measurement of the microlensing parallax. In addition to the constraints on the angular radius and the apparent brightness I s of the source, we employ the Besancon and GalMod galactic models to estimate the physical properties of the lens. We find excellent agreement between the predictions of the two Galactic models: the companion is likely a resident of the brown dwarf desert with a mass Mp about 30 MJup and the host is a main sequence dwarf star. The lens lies along the line of sight to the Galactic Bulge, at a distance of less then4 kpc. We estimate that in about 10 years, the lens and source will be separated by 55 mas, and it will be possible to confirm the exact nature of the lensing system by using high-resolution imaging from ground or space-based observatories.


title: Prompting to Distill: Boosting Data-Free Knowledge Distillation via Reinforced Prompt

abstract: Data-free knowledge distillation (DFKD) conducts knowledge distillation via eliminating the dependence of original training data, and has recently achieved impressive results in accelerating pre-trained language models. At the heart of DFKD is to reconstruct a synthetic dataset by inverting the parameters of the uncompressed model. Prior DFKD approaches, however, have largely relied on hand-crafted priors of the target data distribution for the reconstruction, which can be inevitably biased and often incompetent to capture the intrinsic distributions. To address this problem, we propose a prompt-based method, termed as PromptDFD, that allows us to take advantage of learned language priors, which effectively harmonizes the synthetic sentences to be semantically and grammatically correct. Specifically, PromptDFD leverages a pre-trained generative model to provide language priors and introduces a reinforced topic prompter to control data synthesis, making the generated samples thematically relevant and semantically plausible, and thus friendly to downstream tasks. As shown in our experiments, the proposed method substantially improves the synthesis quality and achieves considerable improvements on distillation performance. In some cases, PromptDFD even gives rise to results on par with those from the data-driven knowledge distillation with access to the original training data.


title: A Two-Phase Method for Production Planning and Machine Speed Optimization Problem

abstract: Textile industry is becoming a highly competitive area with the increase in demand for textile products. Since expanding the production capacity is not always feasible, optimizing the existing system is more practical. In particular, we consider a felt production system of a textile factory operating in Turkey in this study. We aim to minimize the production costs by optimizing machine operating speeds as well as building an efficient production lot sizing plan within the planning horizon. In this direction, we propose the Lot Sizing and Machine Speed (LSMS) nonlinear model to determine the optimal unit processing times and production quantities while minimizing the work-in-process and end item inventories by changing the machine operating speeds dynamically according to the demands. Since LSMS nonlinear optimization problem is NP-hard, we design a Two-Phase heuristic which iteratively processes a linear programming model by utilizing a commercial solver at each phase. We intensively test our Two-Phase heuristic via randomly generated demand, planning horizon and machine-hour capacity scenarios. Our computational experiments show that the introduced Two-Phase heuristic can find the local-optimal results in acceptable amount of time.


title: A model aggregation approach for high-dimensional large-scale optimization

abstract: Bayesian optimization (BO) has been widely used in machine learning and simulation optimization. With the increase in computational resources and storage capacities in these fields, high-dimensional and large-scale problems are becoming increasingly common. In this study, we propose a model aggregation method in the Bayesian optimization (MamBO) algorithm for efficiently solving high-dimensional large-scale optimization problems. MamBO uses a combination of subsampling and subspace embeddings to collectively address high dimensionality and large-scale issues; in addition, a model aggregation method is employed to address the surrogate model uncertainty issue that arises when embedding is applied. This surrogate model uncertainty issue is largely ignored in the embedding literature and practice, and it is exacerbated when the problem is high-dimensional and data are limited. Our proposed model aggregation method reduces these lower-dimensional surrogate model risks and improves the robustness of the BO algorithm. We derive an asymptotic bound for the proposed aggregated surrogate model and prove the convergence of MamBO. Benchmark numerical experiments indicate that our algorithm achieves superior or comparable performance to other commonly used high-dimensional BO algorithms. Moreover, we apply MamBO to a cascade classifier of a machine learning algorithm for face detection, and the results reveal that MamBO finds settings that achieve higher classification accuracy than the benchmark settings and is computationally faster than other high-dimensional BO algorithms.


title: The invariant action for solitonic 5-branes

abstract: We construct the full effective action including DBI and WZ terms for solitonic 5-branes covariant under T-duality. The result is a completion of results known in the literature to a full T-duality covariant expression. The covariant WZ action includes previously omitted R-R terms. The obtained full covariant effective action reproduces the one obtained by S-duality from the D5-brane upon the correct choice of the covariant charge.


title: Strong $\mathbb A^1$-invariance of $\mathbb A^1$-connected components of reductive algebraic groups

abstract: We show that the sheaf of $\mathbb A^1$-connected components of a reductive algebraic group over a perfect field is strongly $\mathbb A^1$-invariant. As a consequence, torsors under such groups give rise to $\mathbb A^1$-fiber sequences. We also show that sections of $\mathbb A^1$-connected components of anisotropic, semisimple, simply connected algebraic groups over an arbitrary field agree with their $R$-equivalence classes, thereby removing the perfectness assumption in the previously known results about the characterization of isotropy in terms of affine homotopy invariance of Nisnevich locally trivial torsors.


title: The Smallest Hard Trees

abstract: We find an orientation of a tree with 20 vertices such that the corresponding fixed-template constraint satisfaction problem (CSP) is NP-complete, and prove that for every orientation of a tree with fewer vertices the corresponding CSP can be solved in polynomial time. We also compute the smallest tree that is NL-hard (assuming L is not NL), the smallest tree that cannot be solved by arc consistency, and the smallest tree that cannot be solved by Datalog. Our experimental results also support a conjecture of Bulin concerning a question of Hell, Nesetril and Zhu, namely that "easy trees lack the ability to count". Most proofs are computer-based and make use of the most recent universal-algebraic theory about the complexity of finite-domain CSPs. However, further ideas are required because of the huge number of orientations of trees. In particular, we use the well-known fact that it suffices to study orientations of trees that are cores and show how to efficiently decide whether a given orientation of a tree is a core using the arc-consistency procedure. Moreover, we present a method to generate orientations of trees that are cores that works well in practice. In this way we found interesting examples for the open research problem to classify finite-domain CSPs in NL.


title: Specification is Law: Safe Creation and Upgrade of Ethereum Smart Contracts

abstract: Smart contracts are the building blocks of the "code is law" paradigm: the smart contract's code indisputably describes how its assets are to be managed - once it is created, its code is typically immutable. Faulty smart contracts present the most significant evidence against the practicality of this paradigm; they are well-documented and resulted in assets worth vast sums of money being compromised. To address this issue, the Ethereum community proposed (i) tools and processes to audit/analyse smart contracts, and (ii) design patterns implementing a mechanism to make contract code mutable. Individually, (i) and (ii) only partially address the challenges raised by the "code is law" paradigm. In this paper, we combine elements from (i) and (ii) to create a systematic framework that moves away from "code is law" and gives rise to a new "specification is law" paradigm. It allows contracts to be created and upgraded but only if they meet a corresponding formal specification. The framework is centered around \emph{a trusted deployer}: an off-chain service that formally verifies and enforces this notion of conformance. We have prototyped this framework, and investigated its applicability to contracts implementing two widely used Ethereum standards: the ERC20 Token Standard and ERC1155 Multi Token Standard, with promising results.


title: Quasinormal modes of Kerr-like black bounce spacetime

abstract: The quasinormal modes of the Kerr-like black bounce spacetime under the scalar field perturbations are studied. We derive the effective potential of scalar perturbation in Kerr-like black bounce spacetime, and study the influence of Kerr-like black bounce spacetime parameters on quasinormal modes by using WKB method and P\"{o}schl-Teller potential approximation. We find that Kerr-like black bounce spacetime has the analogous double-peaked potential as Schwarzschild-like black bounce spacetime, and mass of the scalar particle has a non-negligible effect on quasinormal modes of Kerr-like black bounce spacetime.


title: Wasserstein t-SNE

abstract: Scientific datasets often have hierarchical structure: for example, in surveys, individual participants (samples) might be grouped at a higher level (units) such as their geographical region. In these settings, the interest is often in exploring the structure on the unit level rather than on the sample level. Units can be compared based on the distance between their means, however this ignores the within-unit distribution of samples. Here we develop an approach for exploratory analysis of hierarchical datasets using the Wasserstein distance metric that takes into account the shapes of within-unit distributions. We use t-SNE to construct 2D embeddings of the units, based on the matrix of pairwise Wasserstein distances between them. The distance matrix can be efficiently computed by approximating each unit with a Gaussian distribution, but we also provide a scalable method to compute exact Wasserstein distances. We use synthetic data to demonstrate the effectiveness of our Wasserstein t-SNE, and apply it to data from the 2017 German parliamentary election, considering polling stations as samples and voting districts as units. The resulting embedding uncovers meaningful structure in the data.


title: Quantitative Discourse Cohesion Analysis of Scientific Scholarly Texts using Multilayer Networks

abstract: Discourse cohesion facilitates text comprehension and helps the reader form a coherent narrative. In this study, we aim to computationally analyze the discourse cohesion in scientific scholarly texts using multilayer network representation and quantify the writing quality of the document. Exploiting the hierarchical structure of scientific scholarly texts, we design section-level and document-level metrics to assess the extent of lexical cohesion in text. We use a publicly available dataset along with a curated set of contrasting examples to validate the proposed metrics by comparing them against select indices computed using existing cohesion analysis tools. We observe that the proposed metrics correlate as expected with the existing cohesion indices. We also present an analytical framework, CHIAA (CHeck It Again, Author), to provide pointers to the author for potential improvements in the manuscript with the help of the section-level and document-level metrics. The proposed CHIAA framework furnishes a clear and precise prescription to the author for improving writing by localizing regions in text with cohesion gaps. We demonstrate the efficacy of CHIAA framework using succinct examples from cohesion-deficient text excerpts in the experimental dataset.


title: Flocking dynamics mediated by weighted social networks

abstract: We study the effects of animal social networks with a weighted pattern of interactions on the flocking transition exhibited by models of self-organized collective motion. Considering a model representing dynamics on a one-dimensional substrate, application of a heterogeneous mean-field theory provides a phase diagram as function of the heterogeneity of the network connections and the correlations between weights and degree. In this diagram we observe two phases, one corresponding to the presence of a transition and other to a transition suppressed in an always ordered system, already observed in the non-weighted case. Interestingly, a third phase, with no transition in an always disordered state, is also obtained. These predictions, numerically recovered in computer simulations, are also fulfilled for the more realistic Vicsek model, with movement in a two-dimensional space. Additionally, we observe at finite network sizes the presence of a maximum threshold for particular weight configurations, indicating that it is possible to tune weights to achieve a maximum resilience to noise effects. Simulations in real weighted animal social networks show that, in general, the presence of weights diminishes the value of the flocking threshold, thus increasing the fragility of the flocking state. The shift in the threshold is observed to depend on the heterogeneity of the weight pattern.


title: The role of dineutrino modes in the search for new physics

abstract: Dineutrino modes offer promising searches for new physics. The potential aspects of these modes are reviewed in detail. Performing a proper combination of them, novel tests of the SM symmetries are derived. Different phenomenological applications are worked out, including charm, beauty and kaons, which result in novel tests of lepton universality and charged lepton flavour conservation with flavour-summed dineutrino observables, in addition to improved bounds on $\tau\,\ell$ couplings with $\ell=e,\mu,\tau$.


title: Code Smells Detection via Modern Code Review: A Study of the OpenStack and Qt Communities

abstract: Code review that detects and locates defects and other quality issues plays an important role in software quality control. One type of issue that may impact the quality of software is code smells. Yet, little is known about the extent to which code smells are identified during modern code review. To investigate the concept behind code smells identified in modern code review and what actions reviewers suggest and developers take in response to the identified smells, we conducted a study of code smells in code reviews by analyzing reviews from four large open source projects from the OpenStack (Nova and Neutron) and Qt (Qt Base and Qt Creator) communities. We manually checked a total of 25,415 code review comments obtained by keywords search and random selection, and identified 1,539 smell-related reviews. Our analysis found that 1) code smells were not commonly identified in code reviews, 2) smells were usually caused by violation of coding conventions, 3) reviewers usually provided constructive feedback, including fixing (refactoring) recommendations to help developers remove smells, 4) developers generally followed those recommendations and actioned the changes, 5) once identified by reviewers, it usually takes developers less than one week to fix the smells, and 6) the main reason why developers chose to ignore the identified smells is not worth fixing the smell. Our results suggest that: 1) developers should closely follow coding conventions in their projects to avoid introducing code smells, 2) review-based detection of code smells is perceived to be a trustworthy approach by developers, mainly because reviews are context-sensitive (as reviewers are more aware of the context of the code given that they are part of the project's development team), and 3) program context needs to be fully considered in order to make a decision of whether to fix the identified code smell immediately.


title: Reachability Constrained Reinforcement Learning

abstract: Constrained reinforcement learning (CRL) has gained significant interest recently, since safety constraints satisfaction is critical for real-world problems. However, existing CRL methods constraining discounted cumulative costs generally lack rigorous definition and guarantee of safety. In contrast, in the safe control research, safety is defined as persistently satisfying certain state constraints. Such persistent safety is possible only on a subset of the state space, called feasible set, where an optimal largest feasible set exists for a given environment. Recent studies incorporate feasible sets into CRL with energy-based methods such as control barrier function (CBF), safety index (SI), and leverage prior conservative estimations of feasible sets, which harms the performance of the learned policy. To deal with this problem, this paper proposes the reachability CRL (RCRL) method by using reachability analysis to establish the novel self-consistency condition and characterize the feasible sets. The feasible sets are represented by the safety value function, which is used as the constraint in CRL. We use the multi-time scale stochastic approximation theory to prove that the proposed algorithm converges to a local optimum, where the largest feasible set can be guaranteed. Empirical results on different benchmarks validate the learned feasible set, the policy performance, and constraint satisfaction of RCRL, compared to CRL and safe control baselines.


title: Problem Decomposition and Multi-shot ASP Solving for Job-shop Scheduling

abstract: The Job-shop Scheduling Problem (JSP) is a well-known and challenging combinatorial optimization problem in which tasks sharing a machine are to be arranged in a sequence such that encompassing jobs can be completed as early as possible. In this paper, we propose problem decomposition into time windows whose operations can be successively scheduled and optimized by means of multi-shot Answer Set Programming (ASP) solving. Decomposition aims to split highly complex scheduling tasks into better manageable sub-problems with a balanced number of operations so that good quality or even optimal partial solutions can be reliably found in a small fraction of runtime. Problem decomposition must respect the precedence of operations within their jobs and partial schedules optimized by time windows should yield better global solutions than obtainable in similar runtime on the entire instance. We devise and investigate a variety of decomposition strategies in terms of the number and size of time windows as well as heuristics for choosing their operations. Moreover, we incorporate time window overlapping and compression techniques into the iterative scheduling process to counteract window-wise optimization limitations restricted to partial schedules. Our experiments on JSP benchmark sets of several sizes show that successive optimization by multi-shot ASP solving leads to substantially better schedules within the runtime limit than global optimization on the full problem, where the gap increases with the number of operations to schedule. While the obtained solution quality still remains behind a state-of-the-art Constraint Programming system, our multi-shot solving approach comes closer the larger the instance size, demonstrating good scalability by problem decomposition.


title: Dynamical spontaneous scalarization in Einstein-Maxwell-scalar theory

abstract: We study the linear instability and the nonlinear dynamical evolution of the Reissner-Nordstr\"om (RN) black hole in the Einstein-Maxwell-scalar theory in asymptotic flat spacetime. We focus on the coupling function $f(\phi)=e^{-b\phi^2}$ which allows both the scalar-free RN solution and scalarized black hole solution. We first present the evolution of system parameters during dynamic scalarization. For parameter regions where spontaneous scalarization occurs, we find that the evolution of the scalar field at the horizon is dominated by the fundamental unstable mode from linear analysis at early times. At late times, the nonlinear evolution can be viewed as the perturbation of scalarized black holes.


title: Filtrations on the globalization of twisted D-modules over Dedekind schemes

abstract: Fabian Januszewski and the author established the theory of twisted D-modules over general base schemes. In this short note, we construct a $K$-invariant positive exhaustive filtration on the globalization of the twisted D-module on a smooth quasi-compact $K$-scheme over a Dedekind scheme $S$ obtained by the direct image of a $K$-equivariant twisted integrable connection along a $K$-equivariant closed immersion from a smooth proper $K$-scheme $Y$ with $K$ a smooth $S$-affine group scheme, whose $p$th associated graded $\mathcal{O}_S$-module is locally free of finite rank for every integer $p$. In particular, the $k$-module of its global sections is projective if $S$ is affine with coordinate ring $k$.


title: The AI Teacher Test: Measuring the Pedagogical Ability of Blender and GPT-3 in Educational Dialogues

abstract: How can we test whether state-of-the-art generative models, such as Blender and GPT-3, are good AI teachers, capable of replying to a student in an educational dialogue? Designing an AI teacher test is challenging: although evaluation methods are much-needed, there is no off-the-shelf solution to measuring pedagogical ability. This paper reports on a first attempt at an AI teacher test. We built a solution around the insight that you can run conversational agents in parallel to human teachers in real-world dialogues, simulate how different agents would respond to a student, and compare these counterpart responses in terms of three abilities: speak like a teacher, understand a student, help a student. Our method builds on the reliability of comparative judgments in education and uses a probabilistic model and Bayesian sampling to infer estimates of pedagogical ability. We find that, even though conversational agents (Blender in particular) perform well on conversational uptake, they are quantifiably worse than real teachers on several pedagogical dimensions, especially with regard to helpfulness (Blender: {\Delta} ability = -0.75; GPT-3: {\Delta} ability = -0.93).


title: Ferrimagnetism in stable non-metal covalent organic framework

abstract: We synthesized a pure organic non-metal crystalline covalent organic framework TAPA-BTD-COF by bottom-up Schiff base chemical reaction. And this imine-based COF is stable in aerobic condition and room-temperature. We discovered that this TAPA-BTD-COF exhibited strong magneticity in 300 K generating magnetic hysteresis loop in M-H characterization and giant chimol up to 0.028. And we further conducted zero-field cooling and field-cooling measurement of M-T curves. The as-synthesized materials showed a large chi/mol up to 0.028 in 300 K and increasing to 0.037 in 4.0 K with 200 Oe measurement field. The TAPA-BTD-COF 1/chimol~T curve supported its ferrimagnetism, with an intrinsic delta temperature as -33.03 K by extrapolating the 1/chimol~T curve. From the continuously increasing slope of 1/chimol~T, we consider that this TAPA-BTD-COF belongs to ferrimagnetic other than antiferromagnetic materials. And the large chimol value 0.028 at 300 K and 0.037 at 4.0 K also supported this, since common antiferromagnetic materials possess chimol in the range of 10-5 to 10-3 as weak magnetics other than strong magnetic materials such as ferrimagnetics and ferromagnetics. Since this material is purely non-metal organic polymer, the possibility of d-block and f-block metal with unpaired-electron induced magnetism can be excluded. Besides, since the COF does not involve free-radical monomer in the processes of synthesis, we can also exclude the origin of free-radical induced magnetism. According to recent emerging flat-band strong correlated exotic electron property, this unconventional phenomenon may relate to n-type doping on the flat-band locating in the CBM, thus generating highly-localized electron with infinite effective mass and exhibiting strong correlation, which accounts for this non-trivial strong and stable ferrimagneticity at room-temperature and aerobic atmospheric conditions.


title: Web of resonances and possible path of evolution of the small Uranian satellites

abstract: Satellite systems around giant planets are immersed in a region of complex resonant configurations. Understanding the role of satellite resonances contributes to comprehending the dynamical processes in planetary formation and posterior evolution. Our main goal is to analyse the resonant structure of small moons around Uranus and propose different scenarios able to describe the current configuration of these satellites. We focus our study on the external members of the regular satellites interior to Miranda, namely Rosalind, Cupid, Belinda, Perdita, Puck, and Mab, respectively. We use N-body integrations to perform dynamical maps to analyse their dynamics and proximity to two-body and three-body mean-motion resonances (MMR). We found a complicated web of low-order resonances amongst them. Employing analytical prescriptions, we analysed the evolution by gas drag and type-I migration in a circumplanetary disc (CPD) to explain different possible histories for these moons. We also model the tidal evolution of these satellites using some crude approximations and found possible paths that could lead to MMRs crossing between pairs of moons. Finally, our simulations show that each mechanism can generate significant satellite radial drift leading to possible resonant capture, depending on the distances and sizes.


title: Effect of Sediment Form and Form Distribution on Porosity: A Simulation Study Based on the Discrete Element Method

abstract: Porosity is one of the key properties of dense particle packings like sediment deposits and is influenced by a multitude of grain characteristics such as their size distribution and shape. In the present work, we focus on the form, a specific aspect of the overall shape, of sedimentary grains in order to investigate and quantify its effect on porosity, ultimately deriving novel porosity-prediction models. To this end, we develop a robust and accurate simulation tool based on the discrete element method which we validate against laboratory experiments. Utilizing digital representations of actual sediment from the Rhine river, we first study packings that are composed of particles with a single form. There, the porosity is found to be mainly determined by the inverse equancy, i.e., the ratio of the longest to the smallest form-defining axis. Only for small ratios, additional shape-related properties become relevant, as revealed by a direct comparison to packings of form-equivalent ellipsoids. Since sediment naturally features form mixtures, we extend our simulation tool to study sediment packings with normally-distributed forms. In agreement with our single form studies, the porosity depends primarily on the inverse of the mean equancy. By supplying additional information about a second form factor and the standard deviations, we derive an accurate model for porosity prediction. Due to its simplicity, it can be readily applied to sediment packings for which some measurements of flatness and elongation, the two most common form factors, are available.


title: Stopping Rules for Gradient Methods for Non-Convex Problems with Additive Noise in Gradient

abstract: We study the gradient method under the assumption that an additively inexact gradient is available for, generally speaking, non-convex problems. The non-convexity of the objective function, as well as the use of an inexactness specified gradient at iterations, can lead to various problems. For example, the trajectory of the gradient method may be far enough away from the starting point. On the other hand, the unbounded removal of the trajectory of the gradient method in the presence of noise can lead to the removal of the trajectory of the method from the desired exact solution. The results of investigating the behavior of the trajectory of the gradient method are obtained under the assumption of the inexactness of the gradient and the condition of gradient dominance. It is well known that such a condition is valid for many important non-convex problems. Moreover, it leads to good complexity guarantees for the gradient method. A rule of early stopping of the gradient method is proposed. Firstly, it guarantees achieving an acceptable quality of the exit point of the method in terms of the function. Secondly, the stopping rule ensures a fairly moderate distance of this point from the chosen initial position. In addition to the gradient method with a constant step, its variant with adaptive step size is also investigated in detail, which makes it possible to apply the developed technique in the case of an unknown Lipschitz constant for the gradient. Some computational experiments have been carried out which demonstrate effectiveness of the proposed stopping rule for the investigated gradient methods.


title: Heri-Graphs: A Workflow of Creating Datasets for Multi-modal Machine Learning on Graphs of Heritage Values and Attributes with Social Media

abstract: Values (why to conserve) and Attributes (what to conserve) are essential concepts of cultural heritage. Recent studies have been using social media to map values and attributes conveyed by public to cultural heritage. However, it is rare to connect heterogeneous modalities of images, texts, geo-locations, timestamps, and social network structures to mine the semantic and structural characteristics therein. This study presents a methodological workflow for constructing such multi-modal datasets using posts and images on Flickr for graph-based machine learning (ML) tasks concerning heritage values and attributes. After data pre-processing using state-of-the-art ML models, the multi-modal information of visual contents and textual semantics are modelled as node features and labels, while their social relationships and spatiotemporal contexts are modelled as links in Multi-Graphs. The workflow is tested in three cities containing UNESCO World Heritage properties - Amsterdam, Suzhou, and Venice, which yielded datasets with high consistency for semi-supervised learning tasks. The entire process is formally described with mathematical notations, ready to be applied in provisional tasks both as ML problems with technical relevance and as urban/heritage study questions with societal interests. This study could also benefit the understanding and mapping of heritage values and attributes for future research in global cases, aiming at inclusive heritage management practices.


title: Finite-horizon Equilibria for Neuro-symbolic Concurrent Stochastic Games

abstract: We present novel techniques for neuro-symbolic concurrent stochastic games, a recently proposed modelling formalism to represent a set of probabilistic agents operating in a continuous-space environment using a combination of neural network based perception mechanisms and traditional symbolic methods. To date, only zero-sum variants of the model were studied, which is too restrictive when agents have distinct objectives. We formalise notions of equilibria for these models and present algorithms to synthesise them. Focusing on the finite-horizon setting, and (global) social welfare subgame-perfect optimality, we consider two distinct types: Nash equilibria and correlated equilibria. We first show that an exact solution based on backward induction may yield arbitrarily bad equilibria. We then propose an approximation algorithm called frozen subgame improvement, which proceeds through iterative solution of nonlinear programs. We develop a prototype implementation and demonstrate the benefits of our approach on two case studies: an automated car-parking system and an aircraft collision avoidance system.


title: SQ-VAE: Variational Bayes on Discrete Representation with Self-annealed Stochastic Quantization

abstract: One noted issue of vector-quantized variational autoencoder (VQ-VAE) is that the learned discrete representation uses only a fraction of the full capacity of the codebook, also known as codebook collapse. We hypothesize that the training scheme of VQ-VAE, which involves some carefully designed heuristics, underlies this issue. In this paper, we propose a new training scheme that extends the standard VAE via novel stochastic dequantization and quantization, called stochastically quantized variational autoencoder (SQ-VAE). In SQ-VAE, we observe a trend that the quantization is stochastic at the initial stage of the training but gradually converges toward a deterministic quantization, which we call self-annealing. Our experiments show that SQ-VAE improves codebook utilization without using common heuristics. Furthermore, we empirically show that SQ-VAE is superior to VAE and VQ-VAE in vision- and speech-related tasks.


title: A Neuro-Symbolic ASP Pipeline for Visual Question Answering

abstract: We present a neuro-symbolic visual question answering (VQA) pipeline for CLEVR, which is a well-known dataset that consists of pictures showing scenes with objects and questions related to them. Our pipeline covers (i) training neural networks for object classification and bounding-box prediction of the CLEVR scenes, (ii) statistical analysis on the distribution of prediction values of the neural networks to determine a threshold for high-confidence predictions, and (iii) a translation of CLEVR questions and network predictions that pass confidence thresholds into logic programs so that we can compute the answers using an ASP solver. By exploiting choice rules, we consider deterministic and non-deterministic scene encodings. Our experiments show that the non-deterministic scene encoding achieves good results even if the neural networks are trained rather poorly in comparison with the deterministic approach. This is important for building robust VQA systems if network predictions are less-than perfect. Furthermore, we show that restricting non-determinism to reasonable choices allows for more efficient implementations in comparison with related neuro-symbolic approaches without loosing much accuracy. This work is under consideration for acceptance in TPLP.


title: Magnetic neutron scattering from spherical nanoparticles with Neel surface anisotropy: Analytical treatment

abstract: The magnetization profile and the related magnetic small-angle neutron scattering cross section of a single spherical nanoparticle with Neel surface anisotropy is analytically investigated. We employ a Hamiltonian that comprises the isotropic exchange interaction, an external magnetic field, a uniaxial magnetocrystalline anisotropy in the core of the particle, and the Neel anisotropy at the surface. Using a perturbation approach, the determination of the magnetization profile can be reduced to a Helmholtz equation with Neumann boundary condition, whose solution is represented by an infinite series in terms of spherical harmonics and spherical Bessel functions. From the resulting infinite series expansion, we analytically calculate the Fourier transform, which is algebraically related to the magnetic small-angle neutron scattering cross section. The approximate analytical solution is compared to the numerical solution using the Landau-Lifshitz equation, which accounts for the full nonlinearity of the problem.


title: Stable Matching with Multilayer Approval Preferences: Approvals can be Harder than Strict Preferences

abstract: We study stable matching problems where agents have multilayer preferences: There are $\ell$ layers each consisting of one preference relation for each agent. Recently, Chen et al. [EC '18] studied such problems with strict preferences, establishing four multilayer adaptions of classical notions of stability. We follow up on their work by analyzing the computational complexity of stable matching problems with multilayer approval preferences. We consider eleven stability notions derived from three well-established stability notions for stable matchings with ties and the four adaptions proposed by Chen et al. For each stability notion, we show that the problem of finding a stable matching is either polynomial-time solvable or NP-hard. Furthermore, we examine the influence of the number of layers and the desired "degree of stability" on the problems' complexity. Somewhat surprisingly, we discover that assuming approval preferences instead of strict preferences does not considerably simplify the situation (and sometimes even makes polynomial-time solvable problems NP-hard).


title: Dynamical emergence of a Kosterlitz-Thouless transition in a disordered Bose gas following a quench

abstract: We study the dynamical evolution of a two-dimensional Bose gas after a disorder potential quench. Depending on the initial conditions, the system evolves either to a thermal or a superfluid state. Using extensive quasi-exact numerical simulations, we show that the two phases are separated by a Kosterlitz-Thouless transition. The thermalization time is shown to be longer in the superfluid phase, but no critical slowing down is observed at the transition. The long-time phase diagram is well reproduced by a simple theoretical model. The spontaneous emergence of Kosterlitz-Thouless transitions following a quench is a generic phenomenon that should arise both in the context of non-equilibrium quantum gases and nonlinear, classical wave systems.


title: Magnetic neutron scattering from spherical nanoparticles with Neel surface anisotropy: Atomistic simulations

abstract: We consider a dilute ensemble of randomly-oriented noninteracting spherical nanomagnets and investigate its magnetization structure and ensuing neutron-scattering response by numerically solving the Landau-Lifshitz equation. Taking into account the isotropic exchange interaction, an external magnetic field, a uniaxial magnetic anisotropy for the particle core, and in particular the Neel surface anisotropy, we compute the magnetic small-angle neutron scattering cross section and pair-distance distribution function from the obtained equilibrium spin structures. The numerical results are compared to the well-known analytical expressions for uniformly magnetized particles and provide guidance to the experimentalist. Moreover, the effect of a particle-size distribution function is modeled.


title: The interactions of SARS-CoV-2 with co-circulating pathogens: Epidemiological implications and current knowledge gaps

abstract: Despite the availability of effective vaccines, the persistence of SARS-CoV-2 suggests that co-circulation with other pathogens and resulting multi-epidemics may become increasingly frequent. To better forecast and control the risk of such multi-epidemics, it is essential to elucidate the potential interactions of SARS-CoV-2 with other pathogens; these interactions, however, remain poorly defined. Here, we aimed to review the current body of evidence about SARS-CoV-2 interactions. To study pathogen interactions in a systematic way, we first developed a general framework to capture their major components: sign, strength, symmetry, duration, and mechanism. We then reviewed the experimental evidence from animal models about SARS-CoV-2 interactions. Of the 14 studies identified, 11 focused on the outcomes of co-infection with non-attenuated influenza A viruses and generally demonstrated that co-infection increased disease severity compared with either mono-infection. By contrast, the effect of co-infection on the viral load of either virus was variable and inconsistent across studies. Next, we reviewed the epidemiological evidence about SARS-CoV-2 interactions in human populations. Although numerous studies were identified, only few were specifically designed to infer interaction and many were prone to multiple biases, including confounding. Nevertheless, their results suggested that influenza and pneumococcal conjugate vaccinations were associated with reduced riskof SARS-CoV-2 infection. Finally, we formulated simple transmission models of SARS-CoV-2 co-circulation with a viral or a bacterial pathogen, showing how they can naturally incorporate the proposed framework. More generally, we argue that such models, when designed with an integrative and multidisciplinary perspective, will be invaluable tools to resolve the substantial uncertainties that remain about SARS-CoV-2 interactions.


title: Towards on-sky adaptive optics control using reinforcement learning

abstract: The direct imaging of potentially habitable Exoplanets is one prime science case for the next generation of high contrast imaging instruments on ground-based extremely large telescopes. To reach this demanding science goal, the instruments are equipped with eXtreme Adaptive Optics (XAO) systems which will control thousands of actuators at a framerate of kilohertz to several kilohertz. Most of the habitable exoplanets are located at small angular separations from their host stars, where the current XAO systems' control laws leave strong residuals.Current AO control strategies like static matrix-based wavefront reconstruction and integrator control suffer from temporal delay error and are sensitive to mis-registration, i.e., to dynamic variations of the control system geometry. We aim to produce control methods that cope with these limitations, provide a significantly improved AO correction and, therefore, reduce the residual flux in the coronagraphic point spread function. We extend previous work in Reinforcement Learning for AO. The improved method, called PO4AO, learns a dynamics model and optimizes a control neural network, called a policy. We introduce the method and study it through numerical simulations of XAO with Pyramid wavefront sensing for the 8-m and 40-m telescope aperture cases. We further implemented PO4AO and carried out experiments in a laboratory environment using MagAO-X at the Steward laboratory. PO4AO provides the desired performance by improving the coronagraphic contrast in numerical simulations by factors 3-5 within the control region of DM and Pyramid WFS, in simulation and in the laboratory. The presented method is also quick to train, i.e., on timescales of typically 5-10 seconds, and the inference time is sufficiently small (< ms) to be used in real-time control for XAO with currently available hardware even for extremely large telescopes.


title: Peridynamic modeling for impact failure of wet concrete considering the influence of saturation

abstract: In this paper, a modified intermediately homogenized peridynamic (IH-PD) model for analyzing impact failure of wet concrete has been presented under the configuration of ordinary state-based peridynamic theory. The meso-structural properties of concrete are linked to the macroscopic mechanical behavior in the IH-PD model, where the heterogeneity of concrete is taken into account, and the calculation cost does not increase. Simultaneously, the porosity of concrete is considered, which is implemented by deleting the bond between two material points, as well as the influence of porosity on the mechanical properties of concrete. Moreover, the effective bulk and shear modulus of cement mortar in wet concrete (saturated and unsaturated concrete) are calculated respectively. The dynamic model for wet concrete is described from three aspects: strength, dynamic increase factor, and equation of state. Validation of the proposed model is established through analyzing some benchmark tests and comparing with the corresponding experiment and other available numerical results.


title: An Effective Transformer-based Solution for RSNA Intracranial Hemorrhage Detection Competition

abstract: We present an effective method for Intracranial Hemorrhage Detection (IHD) which exceeds the performance of the winner solution in RSNA-IHD competition (2019). Meanwhile, our model only takes quarter parameters and ten percent FLOPs compared to the winner's solution. The IHD task needs to predict the hemorrhage category of each slice for the input brain CT. We review the top-5 solutions for the IHD competition held by the Radiological Society of North America(RSNA) in 2019. Nearly all the top solutions rely on 2D convolutional networks and sequential models (Bidirectional GRU or LSTM) to extract intra-slice and inter-slice features, respectively. All the top solutions enhance the performance by leveraging the model ensemble, and the model number varies from 7 to 31. In the past years, since much progress has been made in the computer vision regime especially Transformer-based models, we introduce the Transformer-based techniques to extract the features in both intra-slice and inter-slice views for IHD tasks. Additionally, a semi-supervised method is embedded into our workflow to further improve the performance. The code is available in the manuscript.


title: Heroes, Villains, and Victims, and GPT-3: Automated Extraction of Character Roles Without Training Data

abstract: This paper shows how to use large-scale pre-trained language models to extract character roles from narrative texts without training data. Queried with a zero-shot question-answering prompt, GPT-3 can identify the hero, villain, and victim in diverse domains: newspaper articles, movie plot summaries, and political speeches.


title: Evolution of Braneworld Kerr-Newman Naked Singularities

abstract: We study evolution of the braneworld Kerr--Newman (K-N) naked singularities, namely their mass $M$ , spin $a$, and tidal charge $b$ characterizing the role of the bulk space, due to matter in-falling from Keplerian accretion disk. We construct the evolution in two limiting cases applied to the tidal charge. In the first case we assume $b$ = const during the evolution, in the second one we assume that the dimensionless tidal charge $\beta \equiv b/M^2$ = const. For positive values of the tidal charge the evolution is equivalent to the case of the standard K-N naked singularity under accretion of electrically neutral matter. We demonstrate that counter-rotating accretion always converts a K-N naked singularity into an extreme K-N black hole and that the corotating accretion leads to variety of outcomes. The conversion to an extreme K-N black hole is possible for naked singularity with dimensionless tidal charge $\beta < 0.25$, and $\beta \in (0.25, 1)$ with sufficiently low spin. In other cases the accretion ends in a transcendental state. For $0.25 < \beta < 1$ this is a mining unstable K-N naked singularity enabling formally unlimited energy extraction from the naked singularity. In the case of $\beta > 1$, the corotating accretion creates unlimited torodial structure of mater orbiting the naked singularity. Both non-standard outcomes of the corotating accretion imply a transcendence of such naked singularity due to nonlinear gravitational effects.


title: Appropriate reduction of the posterior distribution in fully Bayesian inversions

abstract: Bayesian inversion generates a posterior distribution of model parameters from an observation equation and prior information both weighted by hyperparameters. The prior is also introduced for the hyperparameters in fully Bayesian inversions and enables us to evaluate both the model parameters and hyperparameters probabilistically by the joint posterior. However, even in a linear inverse problem, it is unsolved how we should extract useful information on the model parameters from the joint posterior. This study presents a theoretical exploration into the appropriate dimensionality reduction of the joint posterior in the fully Bayesian inversion. We classify the ways of probability reduction into the following three categories focused on the marginalisation of the joint posterior: (1) using the joint posterior without marginalisation, (2) using the marginal posterior of the model parameters and (3) using the marginal posterior of the hyperparameters. First, we derive several analytical results that characterise these categories. One is a suite of semianalytic representations of the probability maximisation estimators for respective categories in the linear inverse problem. The mode estimators of categories (1) and (2) are found asymptotically identical for a large number of data and model parameters. We also prove the asymptotic distributions of categories (2) and (3) delta-functionally concentrate on their probability peaks, which predicts two distinct optimal estimates of the model parameters. Second, we conduct a synthetic test and find an appropriate reduction is realised by category (3), typified by Akaike's Bayesian information criterion (ABIC). The other reduction categories are shown inappropriate for the case of many model parameters, where the probability concentration of the marginal posterior of the model parameters is found no longer to mean the central limit theorem...


title: Ensemble Kalman Inversion method for an inverse problem in soil-structure interaction

abstract: The interaction between the foundation structures and the soil has been developed for many engineering applications. For the determination of the stress in foundation structure it is needed to determine the influence of the stiffness of soil with respect to the displacement w of the deformable plate (direct problem), and viceversa, how the stiffness of the foundation structure affects the resulting subsidence (inverse problem). In this paper, we deal with the Winkler mathematical model and propose to use an efficient Ensemble Kalman Inversion scheme (EKI) that regularizes iteratively the ill-posedness of the inverse problem. It is a regularizing optimizer used in Bayesian inverse problems that samples particles in pseudo-time introducing a motion due to the movement of these particles. The EKI algorithm converges to the solution of an optimization problem that minimizes the objective function. In this context we show how to reconstruct the Winkler subgrade reaction coefficient of a rectangular thin plate loaded with an existing building by using the EKI methodology combined by the finite difference method (FDM) to discretize the biharmonic operator of the governing equations.


title: State of Health Estimation of Lithium-Ion Batteries in Vehicle-to-Grid Applications Using Recurrent Neural Networks for Learning the Impact of Degradation Stress Factors

abstract: This work presents an effective state of health indicator to indicate lithium-ion battery degradation based on a long short-term memory (LSTM) recurrent neural network (RNN) coupled with a sliding-window. The developed LSTM RNN is able to capture the underlying long-term dependencies of degraded cell capacity on battery degradation stress factors. The learning performance was robust when there was sufficient training data, with an error of < 5% if more than 1.15 years worth of data was supplied for training.


title: Autonomous Open-Ended Learning of Tasks with Non-Stationary Interdependencies

abstract: Autonomous open-ended learning is a relevant approach in machine learning and robotics, allowing the design of artificial agents able to acquire goals and motor skills without the necessity of user assigned tasks. A crucial issue for this approach is to develop strategies to ensure that agents can maximise their competence on as many tasks as possible in the shortest possible time. Intrinsic motivations have proven to generate a task-agnostic signal to properly allocate the training time amongst goals. While the majority of works in the field of intrinsically motivated open-ended learning focus on scenarios where goals are independent from each other, only few of them studied the autonomous acquisition of interdependent tasks, and even fewer tackled scenarios where goals involve non-stationary interdependencies. Building on previous works, we tackle these crucial issues at the level of decision making (i.e., building strategies to properly select between goals), and we propose a hierarchical architecture that treating sub-tasks selection as a Markov Decision Process is able to properly learn interdependent skills on the basis of intrinsically generated motivations. In particular, we first deepen the analysis of a previous system, showing the importance of incorporating information about the relationships between tasks at a higher level of the architecture (that of goal selection). Then we introduce H-GRAIL, a new system that extends the previous one by adding a new learning layer to store the autonomously acquired sequences of tasks to be able to modify them in case the interdependencies are non-stationary. All systems are tested in a real robotic scenario, with a Baxter robot performing multiple interdependent reaching tasks.


title: Resemblance of the power-law scaling behavior of a non-Markovian and nonlinear point processes

abstract: We analyze the statistical properties of a temporal point process driven by a confined fractional Brownian motion. The event count distribution and power spectral density of this non--Markovian point process exhibit power--law scaling. We show that a nonlinear Markovian point process can reproduce the same scaling behavior. This result indicates a possible link between nonlinearity and apparent non--Markovian behavior.


title: El Logaritmo Integral: N\'umeros primos y algo m\'as

abstract: We show the relevance of the logarithmic integral function in the development of mathematics in the first half of the 19th century. Its importance involved first level mathematicians such as Euler, Gauss, Bessel, Riemann. Our perspective is the result of a detailed study of the original sources. We manage to establish the timeline of how the advances took place. In particular, our study vindicates the contributions of Bessel (in collaboration with Gauss). The logarithmic integral in Gauss's mind played a fundamental role in several fields such as complex analysis (Cauchy's theorem), numerical methods (Gaussian quadrature) as well as its best known relation with the distribution of prime numbers. We study this last aspect in detail starting from the works by Legendre and Tchebycheff up to Riemann.


title: A Framework to Map VMAF with the Probability of Just Noticeable Difference between Video Encoding Recipes

abstract: Just Noticeable Difference (JND) model developed based on Human Vision System (HVS) through subjective studies is valuable for many multimedia use cases. In the streaming industries, it is commonly applied to reach a good balance between compression efficiency and perceptual quality when selecting video encoding recipes. Nevertheless, recent state-of-the-art deep learning based JND prediction model relies on large-scale JND ground truth that is expensive and time consuming to collect. Most of the existing JND datasets contain limited number of contents and are limited to a certain codec (e.g., H264). As a result, JND prediction models that were trained on such datasets are normally not agnostic to the codecs. To this end, in order to decouple encoding recipes and JND estimation, we propose a novel framework to map the difference of objective Video Quality Assessment (VQA) scores, i.e., VMAF, between two given videos encoded with different encoding recipes from the same content to the probability of having just noticeable difference between them. The proposed probability mapping model learns from DCR test data, which is significantly cheaper compared to standard JND subjective test. As we utilize objective VQA metric (e.g., VMAF that trained with contents encoded with different codecs) as proxy to estimate JND, our model is agnostic to codecs and computationally efficient. Throughout extensive experiments, it is demonstrated that the proposed model is able to estimate JND values efficiently.


title: Light transfer transitions beyond higher-order exceptional points in parity-time and anti-parity-time symmetric waveguide arrays

abstract: We propose two non-Hermitian arrays consisting of $N=2l+1$ waveguides and exhibiting parity-time ($\mathcal{PT}$) or anti-$\mathcal{PT}$ symmetry for investigating light transfer dynamics based on $N$th-order exceptional points (EPs). The $\mathcal{PT}$-symmetric array supports two $N$th-order EPs separating an unbroken and a broken phase with real and imaginary eignvalues, respectively. Light transfer dynamics in this array exhibits radically different behaviors, i.e. a unidirectional oscillation behavior in the unbroken phase, an edge-towards localization behavior in the broken phase, and a center-towards localization behavior just at $N$th-order EPs. The anti-$\mathcal{PT}$-symmetric array supports also two $N$th-order EPs separating an unbroken and a broken phase, which refer however to imaginary and real eigenvalues, respectively. Accordingly, light transfer dynamics in this array exhibits a center-towards localization behavior in the unbroken phase and an origin-centered oscillation behavior in the broken phase. These nontrivial light transfer behaviors and their controlled transitions are not viable for otherwise split lower-order EPs and depend on the underlying $SU(2)$ symmetry of spin-$l$ matrices.


title: DMRF-UNet: A Two-Stage Deep Learning Scheme for GPR Data Inversion under Heterogeneous Soil Conditions

abstract: Traditional ground-penetrating radar (GPR) data inversion leverages iterative algorithms which suffer from high computation costs and low accuracy when applied to complex subsurface scenarios. Existing deep learning-based methods focus on the ideal homogeneous subsurface environments and ignore the interference due to clutters and noise in real-world heterogeneous environments. To address these issues, a two-stage deep neural network (DNN), called DMRF-UNet, is proposed to reconstruct the permittivity distributions of subsurface objects from GPR B-scans under heterogeneous soil conditions. In the first stage, a U-shape DNN with multi-receptive-field convolutions (MRF-UNet1) is built to remove the clutters due to inhomogeneity of the heterogeneous soil. Then the denoised B-scan from the MRF-UNet1 is combined with the noisy B-scan to be inputted to the DNN in the second stage (MRF-UNet2). The MRF-UNet2 learns the inverse mapping relationship and reconstructs the permittivity distribution of subsurface objects. To avoid information loss, an end-to-end training method combining the loss functions of two stages is introduced. A wide range of subsurface heterogeneous scenarios and B-scans are generated to evaluate the inversion performance. The test results in the numerical experiment and the real measurement show that the proposed network reconstructs the permittivities, shapes, sizes, and locations of subsurface objects with high accuracy. The comparison with existing methods demonstrates the superiority of the proposed methodology for the inversion under heterogeneous soil conditions.


title: Weakly-supervised Biomechanically-constrained CT/MRI Registration of the Spine

abstract: CT and MRI are two of the most informative modalities in spinal diagnostics and treatment planning. CT is useful when analysing bony structures, while MRI gives information about the soft tissue. Thus, fusing the information of both modalities can be very beneficial. Registration is the first step for this fusion. While the soft tissues around the vertebra are deformable, each vertebral body is constrained to move rigidly. We propose a weakly-supervised deep learning framework that preserves the rigidity and the volume of each vertebra while maximizing the accuracy of the registration. To achieve this goal, we introduce anatomy-aware losses for training the network. We specifically design these losses to depend only on the CT label maps since automatic vertebra segmentation in CT gives more accurate results contrary to MRI. We evaluate our method on an in-house dataset of 167 patients. Our results show that adding the anatomy-aware losses increases the plausibility of the inferred transformation while keeping the accuracy untouched.


title: Limit of solutions for semilinear Hamilton-Jacobi equations with degenerate viscosity

abstract: In the paper we prove the convergence of viscosity solutions $u_{\lambda}$ as $\lambda\rightarrow0_+$ for the parametrized degenerate viscous Hamilton-Jacobi equation \[ H(x,d_x u, \lambda u)=\alpha(x)\Delta u,\quad \alpha(x)\geq 0,\quad x\in \mathbb T^n \] under suitable convex and monotonic conditions on $H: T^*M\times\mathbb R\rightarrow\mathbb R$. Such a limit can be characterized in terms of stochastic Mather measures associated with the critical equation \[ H(x,d_x u,0)=\alpha(x)\Delta u. \]


title: Weighted approximation in higher-dimensional missing digit sets

abstract: In this note, we use the mass transference principle for rectangles, recently obtained by Wang and Wu (Math. Ann., 2021), to study the Hausdorff dimension of sets of "weighted $\Psi$-well-approximable" points in certain self-similar sets in $\mathbb{R}^{d}$. Specifically, we investigate weighted $\Psi$-well-approximable points in "missing digit" sets in $\mathbb{R}^{d}$. The sets we consider are natural generalisations of Cantor-type sets in $\mathbb{R}$ to higher dimensions and include, for example, four corner Cantor sets (or Cantor dust) in the plane with contraction ratio $\frac{1}{n}$ with $n \in \mathbb{N}$.


title: On some generalized inverses and partial orders in $\ast$-rings

abstract: Let $\mathcal{R}$ be a unital ring with involution. The notions of 1MP-inverse and MP1-inverse are extended from $M_{m,n}(\mathbb{C)}$, the set of all $m\times n $ matrices over $\mathbb{C}$, to the set $\mathcal{R}% ^{\dagger}$ of all Moore-Penrose invertible elements in $\mathcal{R}$. We study partial orders on $\mathcal{R}^{\dagger}$ that are induced by 1MP-inverses and MP1-inverses. We also extend to the setting of Rickart $\ast $-rings the concept of another partial order, called the plus order, which has been recently introduced on the set of all bounded linear operators between Hilbert spaces. Properties of these relations are investigated and some known results are thus generalized.