title: Charge Carrier Mediation and Ferromagnetism induced in MnBi6Te10 Magnetic Topological Insulators by antimony doping

abstract: A new kind of intrinsic magnetic topological insulators (MTI) MnBi2Te4 family have shed light on the observation of novel topological quantum effect such as quantum anomalous Hall effect (QAHE). However, the strong anti-ferromagnetic (AFM) coupling and high carrier concentration in the bulk hinder the practical applications. In closely related materials MnBi4Te7 and MnBi6Te10, the interlayer magnetic coupling is greatly suppressed by Bi2Te3 layer intercalation. However, AFM is still the ground state in these compounds. Here by magnetic and transport measurements, we demonstrate that Sb substitutional dopant plays a dual role in MnBi6Te10, which can not only adjust the charge carrier type and the concentration, but also induce the solid into a ferromagnetic (FM) ground state. AFM ground state region which is also close to the charge neutral point can be found in the phase diagram of Mn(SbxBi1-x)6Te10 when x ~ 0.25. An intrinsic FM-MTI candidate is thus demonstrated, and it may take a step further for the realization of high-quality and high-temperature QAHE and the related topological quantum effects in the future.


title: Optimism and Delays in Episodic Reinforcement Learning

abstract: There are many algorithms for regret minimisation in episodic reinforcement learning. This problem is well-understood from a theoretical perspective, providing that the sequences of states, actions and rewards associated with each episode are available to the algorithm updating the policy immediately after every interaction with the environment. However, feedback is almost always delayed in practice. In this paper, we study the impact of delayed feedback in episodic reinforcement learning from a theoretical perspective and propose two general-purpose approaches to handling the delays. The first involves updating as soon as new information becomes available, whereas the second waits before using newly observed information to update the policy. For the class of optimistic algorithms and either approach, we show that the regret increases by an additive term involving the number of states, actions, episode length, the expected delay and an algorithm-dependent constant. We empirically investigate the impact of various delay distributions on the regret of optimistic algorithms to validate our theoretical results.


title: An aggregation model of cockroaches with fast-or-slow motion dichotomy

abstract: We propose a mathematical model, namely a reaction-diffusion system, to describe social behaviour of cockroaches. An essential new aspect in our model is that the dispersion behaviour due to overcrowding effect is taken into account {as a counterpart to commonly studied aggregation}. This consideration leads to an intriguing new phenomenon which has not been observed in the literature. Namely, due to the competition between aggregation towards areas of higher concentration of pheromone and dispersion avoiding overcrowded areas, the cockroaches aggregate more at the transition area of pheromone. Moreover, we also consider the fast reaction limit where the switching rate between active and inactive subpopulations tends to infinity. By utilising improved duality and energy methods, together with the regularisation of heat operator, we prove that the weak solution of the reaction-diffusion system converges to that of a reaction-cross-diffusion system.


title: Axiomatic characterization of the $\chi^{2}$ dissimilarity measure

abstract: We axiomatically characterize the $\chi^{2}$ dissimilarity measure. To this end, we solve a new generalization of a functional equation discussed in Aczel (Lectures on functional equations and their applications, Academic Press, 1966).


title: Inexact proximal DC Newton-type method for nonconvex composite functions

abstract: We consider a class of difference-of-convex (DC) optimization problems where the objective function is the sum of a smooth function and a possible nonsmooth DC function. The application of proximal DC algorithms to address this problem class is well-known. In this paper, we combine a proximal DC algorithm with an inexact proximal Newton-type method to propose an inexact proximal DC Newton-type method. We demonstrate global convergence properties of the proposed method. In addition, we give a memoryless quasi-Newton matrix for scaled proximal mappings and consider a two-dimensional system of semi-smooth equations that arise in calculating scaled proximal mappings. To efficiently obtain the scaled proximal mappings, we adopt a semi-smooth Newton method to inexactly solve the system. Finally, we present some numerical experiments to investigate the efficiency of the proposed method, showing that the proposed method outperforms existing methods.


title: Microscopic derivation of a traffic flow model with a bifurcation

abstract: The goal of the paper is a rigorous derivation of a macroscopic traffic flow model with a bifurcation or a local perturbation from a microscopic one. The microscopic model is a simple follow-the-leader with random parameters. The random parameters are used as a statistical description of the road taken by a vehicle and its law of motion. The limit model is a deterministic and scalar Hamilton-Jacobi on a network with a flux limiter, the flux-limiter describing how much the bifurcation or the local perturbation slows down the vehicles. The proof of the existence of this flux limiter-the first one in the context of stochastic homogenization-relies on a concentration inequality and on a delicate derivation of a superadditive inequality.


title: Fingerprint Presentation Attack Detection by Channel-wise Feature Denoising

abstract: Due to the diversity of attack materials, fingerprint recognition systems (AFRSs) are vulnerable to malicious attacks. It is thus important to propose effective fingerprint presentation attack detection (PAD) methods for the safety and reliability of AFRSs. However, current PAD methods often exhibit poor robustness under new attack types settings. This paper thus proposes a novel channel-wise feature denoising fingerprint PAD (CFD-PAD) method by handling the redundant noise information ignored in previous studies. The proposed method learns important features of fingerprint images by weighing the importance of each channel and identifying discriminative channels and "noise" channels. Then, the propagation of "noise" channels is suppressed in the feature map to reduce interference. Specifically, a PA-Adaptation loss is designed to constrain the feature distribution to make the feature distribution of live fingerprints more aggregate and that of spoof fingerprints more disperse. Experimental results evaluated on the LivDet 2017 dataset showed that the proposed CFD-PAD can achieve a 2.53% average classification error (ACE) and a 93.83% true detection rate when the false detection rate equals 1.0% (TDR@FDR=1%). Also, the proposed method markedly outperforms the best single-model-based methods in terms of ACE (2.53% vs. 4.56%) and TDR@FDR=1%(93.83% vs. 73.32%), which demonstrates its effectiveness. Although we have achieved a comparable result with the state-of-the-art multiple-model-based methods, there still is an increase in TDR@FDR=1% from 91.19% to 93.83%. In addition, the proposed model is simpler, lighter and more efficient and has achieved a 74.76% reduction in computation time compared with the state-of-the-art multiple-model-based method. The source code is available at https://github.com/kongzhecn/cfd-pad.


title: Critical spin dynamics of Heisenberg ferromagnets revisited

abstract: We calculate the dynamic structure factor $S (\boldsymbol{k},\omega)$ in the paramagnetic regime of quantum Heisenberg ferromagnets for temperatures $T$ close to the critical temperature $T_c$ using our recently developed functional renormalization group approach to quantum spin systems. In $d=3$ dimensions we find that for small momenta $\boldsymbol{k}$ and frequencies $\omega$ the dynamic structure factor assumes the scaling form $S(\boldsymbol{k},\omega) = (\tau T G (\boldsymbol{k})/\pi)\Phi (k\xi, \omega\tau)$, where $ G (\boldsymbol{k})$ is the static spin-spin correlation function, $\xi$ is the correlation length, and the characteristic time-scale $\tau$ is proportional to $\xi^{5/2}$. We explicitly calculate the dynamic scaling function $\Phi (x,y)$ and find satisfactory agreement with neutron scattering experiments probing the critical spin dynamics in EuO and EuS. Precisely at the critical point where $\xi = \infty$ our result for the dynamic structure factor can be written as $S (\boldsymbol{k},\omega) = (\pi\omega_k)^{-1} T_c G (\boldsymbol{k}) \Psi_c (\omega/\omega_k)$, where $\omega_k \propto k^{5/2}$. We find that $\Psi_c(\nu)$ vanishes as $\nu^{-13/5}$ for large $\nu$, and as $\nu^{3/5}$ for small $\nu$. While the large-frequency behavior of $\Psi_c (\nu)$ is consistent with calculations based on mode-coupling theory and with perturbative renormalization group calculations to second order in $\epsilon = 6-d$, our result for small frequencies disagrees with previous calculations. We argue that up until now neither experiments nor numerical simulations are sufficiently accurate to determine the low-frequency behavior of $\Psi_c (\nu)$. We also calculate the low-temperature behavior of $S ( \boldsymbol{k},\omega)$ in one- and two dimensional ferromagnets and find that it satisfies dynamic scaling with exponent $z=2$ and exhibits a pseudogap for small frequencies.


title: MR-compatible loading device for assessment of heel pad internal tissue displacements under shearing load

abstract: In the last decade, the role of shearing loads has been increasingly suspected to play a determinant impact in the formation of deep pressure ulcers. In vivo observations of such deformations are complex to obtain. Previous studies only provide global measurements of such deformations without getting the quantitative values of the loads that generate these deformations. To study the role that shearing loads have in the aetiology of heel pressure ulcers, an MR-compatible device for the application of shearing and normal loads was designed. Magnetic resonance imaging is a key feature that allows to monitor deformations of soft tissues after loading in a non-invasive way. Measuring applied forces in an MR-environment is challenging due to the impossibility to use magnetic materials. In our device, forces are applied through the compression of springs made in polylactide. Shearing and normal loads were applied on the plantar skin of the human heel by means of an indenting plate while acquiring MR images. The device materials did not introduce any imaging artifact and allowed for high quality MR measurements permitting to identify the deformation of the internal components of the heel. The obtained subject-specific results are an original data set that can be used in validations for Finite Element analysis and therefore contribute to a better understanding of the factors involved in pressure ulcer development.


title: Covalent Functionalization of HiPco Single-Walled Carbon Nanotubes: Differences in the Oxidizing Action of H2SO4 and HNO3 during a Soft Oxidation Process

abstract: The results of a study on the evolution of HiPco single-walled carbon nanotubes during the oxidizing action of H2SO4 and HNO3 are presented. The process conditions used have been chosen so as to avoid any significant damage to the nanotube structure. The type and level of functionalization, the location of the grafted functions on the surface of the nanotube and the changes in morphological characteristics of the samples were examined by using a wide and complementary range of analytical techniques. We propose an explanation for the differences in the oxidizing action of sulfuric and nitric acids. The combined results allow us to suggest possible reaction mechanisms that occur on the surface of the nanotube.


title: Attention Mechanisms in Computer Vision: A Survey

abstract: Humans can naturally and effectively find salient regions in complex scenes. Motivated by this observation, attention mechanisms were introduced into computer vision with the aim of imitating this aspect of the human visual system. Such an attention mechanism can be regarded as a dynamic weight adjustment process based on features of the input image. Attention mechanisms have achieved great success in many visual tasks, including image classification, object detection, semantic segmentation, video understanding, image generation, 3D vision, multi-modal tasks and self-supervised learning. In this survey, we provide a comprehensive review of various attention mechanisms in computer vision and categorize them according to approach, such as channel attention, spatial attention, temporal attention and branch attention; a related repository https://github.com/MenghaoGuo/Awesome-Vision-Attentions is dedicated to collecting related work. We also suggest future directions for attention mechanism research.


title: On the validation of pansharpening methods

abstract: Validation of the quality of pansharpening methods is a difficult task because the reference is not directly available. In the meantime, two main approaches have been established: validation in reduced resolution and original resolution. In the former approach it is still not clear how the data are to be processed to a lower resolution. Other open issues are related to the question which resolution and measures should be used. In the latter approach the main problem is how the appropriate measure should be selected. In the most comparison studies the results of both approaches do not correspond, that means in each case other methods are selected as the best ones. Thus, the developers of the new pansharpening methods still stand in the front of dilemma: how to perform a correct or appropriate comparison/evaluation/validation. It should be noted, that the third approach is possible, that is to perform the comparison of methods in a particular application with the usage of their ground truth. But this is not always possible, because usually developers are not working with applications. Moreover, it can be an additional computational load for a researcher in a particular application. In this paper some of the questions/problems raised above are approached/discussed. The following component substitution (CS) and high pass filtering (HPF) pansharpening methods with additive and multiplicative models and their enhancements such as haze correction, histogram matching, usage of spectral response functions (SRF), modulation transfer function (MTF) based lowpass filtering are investigated on remote sensing data of WorldView-2 and WorldView-4 sensors.


title: Low-Subpacketization Multi-Antenna Coded Caching for Dynamic Networks

abstract: Multi-antenna coded caching combines a global caching gain, proportional to the total cache size in the network, with an additional spatial multiplexing gain that stems from multiple transmitting antennas. However, classic centralized coded caching schemes are not suitable for dynamic networks as they require prior knowledge of the number of users to indicate what data should be cached at each user during the placement phase. On the other hand, fully decentralized schemes provide comparable gains to their centralized counterparts only when the number of users is very large. In this paper, we propose a novel multi-antenna coded caching scheme for dynamic networks, where instead of defining individual cache contents, we associate users with a limited set of predefined caching profiles. Then, during the delivery phase, we aim at achieving a combined caching and spatial multiplexing gain, comparable to a large extent with the ideal case of fully centralized schemes. The resulting scheme imposes small subpacketization and beamforming overheads, is robust under dynamic network conditions, and incurs small finite-SNR performance loss compared with centralized schemes.


title: Impact of Airways Geometry on Transport of Gases to Blood

abstract: Topological structure of bronchial trees affects transport of gases and aerosols in the respiratory system. We start by providing a quantitative assessment of the ability of the alternative tree representations to predict observable geometric and mechanistic characteristics, such as network resistance, dead space volume, and path length. Then we present a model of dynamic transport of oxygen and carbon dioxide along the airways, in the alveoli, across the alveolar membrane, and along the pulmonary blood capillaries. The model also accounts for the exchange of these two gases with blood in the capillaries, as well as for age, gender and other in-species characteristics. Our model's predictions are compared with corresponding observations, providing an additional venue to assess the validity of the existing representations of the lung's bronchial tree.


title: Recognizing Series-Parallel Matrices in Linear Time

abstract: A series-parallel matrix is a binary matrix that can be obtained from an empty matrix by successively adjoining rows or columns that are parallel to an existing row/column or have at most one 1-entry. Equivalently, series-parallel matrices are representation matrices of graphic matroids of series-parallel graphs, which can be recognized in linear time. We propose an algorithm that, for an m-by-n matrix A with k nonzeros, determines in expected $\mathcal{O}(m + n + k)$ time whether A is series-parallel, or returns a minimal non-series-parallel submatrix of A. We complement the developed algorithm by an efficient implementation and report about computational results.


title: Improved Decoding of Expander Codes

abstract: We study the classical expander codes, introduced by Sipser and Spielman \cite{SS96}. Given any constants $0< \alpha, \varepsilon < 1/2$, and an arbitrary bipartite graph with $N$ vertices on the left, $M < N$ vertices on the right, and left degree $D$ such that any left subset $S$ of size at most $\alpha N$ has at least $(1-\varepsilon)|S|D$ neighbors, we show that the corresponding linear code given by parity checks on the right has distance at least roughly $\frac{\alpha N}{2 \varepsilon }$. This is strictly better than the best known previous result of $2(1-\varepsilon ) \alpha N$ \cite{Sudan2000note, Viderman13b} whenever $\varepsilon < 1/2$, and improves the previous result significantly when $\varepsilon $ is small. Furthermore, we show that this distance is tight in general, thus providing a complete characterization of the distance of general expander codes. Next, we provide several efficient decoding algorithms, which vastly improve previous results in terms of the fraction of errors corrected, whenever $\varepsilon < \frac{1}{4}$. Finally, we also give a bound on the list-decoding radius of general expander codes, which beats the classical Johnson bound in certain situations (e.g., when the graph is almost regular and the code has a high rate). Our techniques exploit novel combinatorial properties of bipartite expander graphs. In particular, we establish a new size-expansion tradeoff, which may be of independent interests.


title: Quantitative stability of harmonic maps from $\mathbb{R}^2$ to $\mathbb{S}^2$ with higher degree

abstract: For degree $\pm 1$ harmonic maps from $\mathbb{R}^2$ (or $\mathbb{S}^2$) to $\mathbb{S}^2$, Bernand-Mantel, Muratov and Simon \cite{bernand2021quantitative} recently establish a uniformly quantitative stability estimate. Namely, for any map $u:\mathbb{R}^2\to \mathbb{S}^2$ with degree $\pm 1$, the discrepancy of its Dirichlet energy and $4\pi$ can linearly control the $\dot H^1$-difference of $u$ from the set of degree $\pm 1$ harmonic maps. Whether a similar estimate holds for harmonic maps with higher degree is unknown. In this paper, we prove that a similar quantitative stability result for higher degree is true only in local sense. Namely, given a harmonic map, a similar estimate holds if $u$ is already sufficiently near to it (modulo M\"{o}bius transform) and the bound in general depends on the given harmonic map. More importantly, we investigate an example of degree 2 case thoroughly, which shows that it fails to have a uniformly quantitative estimate like the degree $\pm 1$ case. This phenomenon show the striking difference of degree $\pm1$ ones and higher degree ones. Finally, we also conjecture a new uniformly quantitative stability estimate based on our computation.


title: AI in Human-computer Gaming: Techniques, Challenges and Opportunities

abstract: With breakthrough of the AlphaGo, human-computer gaming AI has ushered in a big explosion, attracting more and more researchers all around the world. As a recognized standard for testing artificial intelligence, various human-computer gaming AI systems (AIs) have been developed such as the Libratus, OpenAI Five and AlphaStar, beating professional human players. The rapid development of human-computer gaming AIs indicate a big step of decision making intelligence, and it seems that current techniques can handle very complex human-computer games. So, one natural question raises: what are the possible challenges of current techniques in human-computer gaming, and what are the future trends? To answer the above question, in this paper, we survey recent successful game AIs, covering board game AIs, card game AIs, first-person shooting game AIs and real time strategy game AIs. Through this survey, we 1) compare the main difficulties among different kinds of games and the corresponding techniques utilized for achieving professional human level AIs; 2) summarize the mainstream frameworks and techniques that can be properly relied on for developing AIs for complex human-computer gaming; 3) raise the challenges or drawbacks of current techniques in the successful AIs; and 4) try to point out future trends in human-computer gaming AIs. Finally, we hope this brief review can provide an introduction for beginners, and inspire insights for researchers in the field of AI in human-computer gaming.


title: CoReS: Compatible Representations via Stationarity

abstract: Compatible features enable the direct comparison of old and new learned features allowing to use them interchangeably over time. In visual search systems, this eliminates the need to extract new features from the gallery-set when the representation model is upgraded with novel data. This has a big value in real applications as re-indexing the gallery-set can be computationally expensive when the gallery-set is large, or even infeasible due to privacy or other concerns of the application. In this paper, we propose CoReS, a new training procedure to learn representations that are \textit{compatible} with those previously learned, grounding on the stationarity of the features as provided by fixed classifiers based on polytopes. With this solution, classes are maximally separated in the representation space and maintain their spatial configuration stationary as new classes are added, so that there is no need to learn any mappings between representations nor to impose pairwise training with the previously learned model. We demonstrate that our training procedure largely outperforms the current state of the art and is particularly effective in the case of multiple upgrades of the training-set, which is the typical case in real applications.


title: Dynamic Network Quantile Regression Model

abstract: We propose a dynamic network quantile regression model to investigate the quantile connectedness using a predetermined network information. We extend the existing network quantile autoregression model of Zhu et al. (2019b) by explicitly allowing the contemporaneous network effects and controlling for the common factors across quantiles. To cope with the endogeneity issue due to simultaneous network spillovers, we adopt the instrumental variable quantile regression (IVQR) estimation and derive the consistency and asymptotic normality of the IVQR estimator using the near epoch dependence property of the network process. Via Monte Carlo simulations, we confirm the satisfactory performance of the IVQR estimator across different quantiles under the different network structures. Finally, we demonstrate the usefulness of our proposed approach with an application to the dataset on the stocks traded in NYSE and NASDAQ in 2016.


title: Pseudo-domains in imaging data improve prediction of future disease status in multi-center studies

abstract: In multi-center randomized clinical trials imaging data can be diverse due to acquisition technology or scanning protocols. Models predicting future outcome of patients are impaired by this data heterogeneity. Here, we propose a prediction method that can cope with a high number of different scanning sites and a low number of samples per site. We cluster sites into pseudo-domains based on visual appearance of scans, and train pseudo-domain specific models. Results show that they improve the prediction accuracy for steatosis after 48 weeks from imaging data acquired at an initial visit and 12-weeks follow-up in liver disease


title: Identifying urban air pollution hot-spots by dispersion modeling when data are scarce: application to diesel generators in Beirut, Lebanon

abstract: Diesel generators are emerging as community-initiated solutions to compensate for electricity shortage in cities marred by economical crisis and/or conflict. The resulting pollution distribution in dense urban environments is a major source of concern to the population. In the absence of periodic observations from properly distributed sensors, as is the case in Beirut, physically based computational modeling stand out as an effective tool for predicting the pollutant distribution in complex environments, and a cost-effective framework for investigating what-if scenarios and assessing mitigation strategies. Here, we present a Lagrangian transport model-based study of PM2.5 dispersion originating from a large number of diesel generators in Beirut. We explore large and small scale dispersion patterns in selected smalls domains and over the entire city. The scenarios considered investigate the impact of topography, atmospheric stability, presence of buildings, diesel generators distribution, and stacks elevations for representative meteorological conditions. Assessment of these scenarios is carried out in terms of small and large scale dispersion patterns and the mean concentration at street level and population exposure proxy indicators. We also report on the efficacy of elevating the stack height as a mitigation measure at different representative wind and atmospheric stability conditions.


title: Entanglement Classification via Operator Size

abstract: In this work, multipartite entanglement is classified by polynomials. I show that the operator size is closely related to the entanglement structure. Given a generic quantum state, I define a series of subspaces generated by operators of different sizes acting on it. The information about the entanglement is encoded into these subspaces. With the dimension of these subspaces as coefficients, I define a polynomial which I call the entanglement polynomial. The entanglement polynomial induces a homomorphism from quantum states to polynomials. It implies that we can characterize and find the building blocks of entanglement by polynomial factorization. Two states share the same entanglement polynomial if they are equivalent under the stochastic local operations and classical communication. To calculate the entanglement polynomial practically, I construct a series of states, called renormalized states, whose ranks are related to the coefficients of the entanglement polynomial.


title: Drone delivery: Reliable Cellular UAV Communication Using Multi-Operator Diversity

abstract: The market size of Unmanned Aerial Vehicles (UAVs, a.k.a drones) can reach up to 10% of the global market value. In particular, drone delivery is one of the most attractive applications. The growing number of drones requires appropriate traffic management systems that will rely on cellular networks. However, it has been shown in the literature that these networks cannot provide reliable communication due to low coverage probability and frequent handovers. This article presents a potential solution targeting these problems while requiring no modifications of the existing infrastructure. Namely, equipping the UAV with multiple cellular modems to connect to different providers' networks introduces network diversity resulting in 98% coverage probability at the flight altitude of 100 meters. In contrast, one network ensures only 80% coverage. At the same time, the size of the outage zones becomes up to ten times smaller and the frequency of harmful handovers is reduced to zero. The results are obtained with a physical-layer simulator utilizing a real urban 3D environment, cellular network parameters (e.g., site locations, antenna orientation and gains), and specific aerial channel models.


title: Operation of an archaeological lead PbWO$_4$ crystal to search for neutrinos from astrophysical sources with a Transition Edge Sensor

abstract: The experimental detection of the CE$\nu$NS allows the investigation of neutrinos and neutrino sources with all-flavor sensitivity. Given its large content in neutrons and stability, Pb is a very appealing choice as target element. The presence of the radioisotope $^{210}$Pb (T$_{1/2}\sim$22 yrs) makes natural Pb unsuitable for low-background, low-energy event searches. This limitation can be overcome employing Pb of archaeological origin, where several half-lives of $^{210}$Pb have gone by. We present results of a cryogenic measurement of a 15g PbWO$_4$ crystal, grown with archaeological Pb (older than $\sim$2000 yrs) that achieved a sub-keV nuclear recoil detection threshold. A ton-scale experiment employing such material, with a detection threshold for nuclear recoils of just 1 keV would probe the entire Milky Way for SuperNovae, with equal sensitivity for all neutrino flavors, allowing the study of the core of such exceptional events.


title: Second-virial theory for shape-persistent living polymers templated by discs

abstract: Living polymers composed of non-covalently bonded building blocks with weak backbone flexibility may self-assemble into thermoresponsive lyotropic liquid crystals. We demonstrate that the reversible polymer assembly and phase behavior can be controlled by the addition of (non-adsorbing) rigid colloidal discs which act as an entropic reorienting ``template" onto the supramolecular polymers. Using a particle-based second-virial theory that correlates the various entropies associated with the polymers and discs, we demonstrate that small fractions of discotic additives promote the formation of a polymer nematic phase. At larger disc concentrations, however, the phase is disrupted by collective disc alignment in favor of a discotic nematic fluid in which the polymers are dispersed anti-nematically. We show that the anti-nematic arrangement of the polymers generates a non-exponential molecular-weight distribution and stimulates the formation of oligomeric species. At sufficient concentrations the discs facilitate a liquid-liquid phase separation which can be brought into simultaneously coexistence with the two fractionated nematic phases, providing evidence for a four-fluid coexistence in reversible shape-dissimilar hard-core mixtures without cohesive interparticle forces. We stipulate the conditions under which such a phenomenon could be found in experiment.


title: AnimeCeleb: Large-Scale Animation CelebHeads Dataset for Head Reenactment

abstract: We present a novel Animation CelebHeads dataset (AnimeCeleb) to address an animation head reenactment. Different from previous animation head datasets, we utilize 3D animation models as the controllable image samplers, which can provide a large amount of head images with their corresponding detailed pose annotations. To facilitate a data creation process, we build a semi-automatic pipeline leveraging an open 3D computer graphics software with a developed annotation system. After training with the AnimeCeleb, recent head reenactment models produce high-quality animation head reenactment results, which are not achievable with existing datasets. Furthermore, motivated by metaverse application, we propose a novel pose mapping method and architecture to tackle a cross-domain head reenactment task. During inference, a user can easily transfer one's motion to an arbitrary animation head. Experiments demonstrate the usefulness of the AnimeCeleb to train animation head reenactment models, and the superiority of our cross-domain head reenactment model compared to state-of-the-art methods. Our dataset and code are available at https://github.com/kangyeolk/AnimeCeleb.


title: Systematic study of {\Delta}(1232) resonance excitations using single isobaric charge-exchange reactions induced by medium-mass projectiles of Sn

abstract: The fragment separator FRS has been for the first time used to measure the (n,p) and (p,n)-type isobaric charge-exchange cross sections of stable 112,124Sn isotopes accelerated at 1A GeV with an uncertainty of 3% and to separate quasi-elastic and inelastic components in the missing-energy spectra of the ejectiles. The inelastic contribution can be associated to the excitation of isobar {\Delta}(1232) resonances and to the pion emission in s-wave, both in the target and projectile nuclei, while the quasi-elastic contribution is associated to the nuclear spin-isospin response of nucleon-hole excitations. The data lead to interesting results where we observe a clear quenching of the quasi-elastic component and their comparisons to theoretical calculations demonstrate that the baryonic resonances can be excited in the target and projectile nuclei. To go further in this investigation, we propose to study the excitation of baryonic resonances taking advantage of the combination of high-resolving power magnetic spectrometers with the WASA calorimeter. These new measurements will allow us to determine the momenta of the ejectiles and pions emitted in coincidence after the single isobaric charge-exchange collisions, providing us unique opportunities to study the evolution of the baryonic resonance dynamics with the neutron-proton asymmetry through the use of exotic radioactive ion beams.


title: The Role of Self Interactions in the Cosmological Evolution of Warm Dark Matter

abstract: In this work we present a summary of recent studies on the effects of elastic self interactions in the evolution of Warm Dark Matter models (WDM), focusing on structure formation and the evolution of cosmological perturbations. We pay special attention to a particular class of sterile neutrino WDM known as $\nu$MSM and provide examples for the case of vector field self interactions. We calculate the effects of assuming self interacting dark matter in X-Ray astrophysical observations, in the formation of fermionic DM halos in (quasi) equilibrium states and in the evolution of DM perturbations in the early universe, assuming particle masses between $\mathcal{O}(1-100)$ keV. In the latter topic, we perform simulations using a modification to the public Boltzmann solver CLASS and compare our results with observations. We find self interactions to be an interesting addition to WDM models, which can alleviate tensions both present in standard CDM cosmology and regarding WDM itself, as well as provide an interesting avenue for DM halo formation.


title: Mean-field models of dynamics on networks via moment closure: an automated procedure

abstract: In the study of dynamics on networks, moment closure is a commonly used method to obtain low-dimensional evolution equations amenable to analysis. The variables in the evolution equations are mean counts of subgraph states and are referred to as moments. Due to interaction between neighbours, each moment equation is a function of higher-order moments, such that an infinite hierarchy of equations arises. Hence, the derivation requires truncation at a given order, and, an approximation of the highest-order moments in terms of lower-order ones,known as a closure formula. Recent systematic approximations have either restricted focus to closed moment equations for SIR epidemic spreading or to unclosed moment equations for arbitrary dynamics. In this paper, we develop a general procedure that automates both derivation and closure of arbitrary order moment equations for dynamics with nearest-neighbour interactions on undirected networks. Automation of the closure step was made possible by our generalised closure scheme,which systematically decomposes the largest subgraphs into their smaller components.We show that this decomposition is exact if these components form a tree, there is independence at distances beyond their graph diameter, and there is spatial homogeneity. Testing our method for SIS epidemic spreading on lattices and random networks confirms that biases are larger for networks with many short cycles in regimes with long-range dependence. A Mathematica package that automates the moment closure is available for download.


title: Irreducibility of a universal Prym-Brill-Noether locus

abstract: For genus $g = \frac{r(r+1)}{2}+1$, we prove that via the forgetful map, the universal Prym-Brill-Noether locus $\mathcal{R}^r_g$ has a unique irreducible component dominating the moduli space $\mathcal{R}_g$ of Prym curves.


title: Active microrheology in corrugated channels: comparison of thermal and colloidal baths

abstract: The dynamics of colloidal suspension confined within porous materials strongly differs from that in the bulk. In particular, within porous materials, the presence of boundaries with complex shapes entangles the longitudinal and transverse degrees of freedom inducing a coupling between the transport of the suspension and the density inhomogeneities induced by the walls. Colloidal suspension confined within model porous media are characterized by means of active microrheology where a net force is applied on a single colloid (tracer particle) whose transport properties are then studied. The trajectories provided by active microrheology are exploited to determine the local transport coefficients. In order to asses the role of the colloid-colloid interactions we compare the case of a tracer embedded in a colloidal suspension to the case of a tracer suspended in an ideal bath. Our results show that the friction coefficient increases and the passage time distribution widens upon increasing the corrugation of the channel. These features are obtained for a tracer suspended in a (thermalized) colloidal bath as well as for the case of an ideal thermal bath. These results highlight the relevance of the confinement on the transport and show a mild dependence on the colloidal/thermal bath. Finally, we rationalize our numerical results with a semi-analytical model. Interestingly, the predictions of the model are quantitatively reliable for mild external forces, hence providing a reliable tool for predicting the transport across porous materials.


title: Multimodal Generalized Zero Shot Learning for Gleason Grading using Self-Supervised Learning

abstract: Gleason grading from histopathology images is essential for accurate prostate cancer (PCa) diagnosis. Since such images are obtained after invasive tissue resection quick diagnosis is challenging under the existing paradigm. We propose a method to predict Gleason grades from magnetic resonance (MR) images which are non-interventional and easily acquired. We solve the problem in a generalized zero-shot learning (GZSL) setting since we may not access training images of every disease grade. Synthetic MRI feature vectors of unseen grades (classes) are generated by exploiting Gleason grades' ordered nature through a conditional variational autoencoder (CVAE) incorporating self-supervised learning. Corresponding histopathology features are generated using cycle GANs, and combined with MR features to predict Gleason grades of test images. Experimental results show our method outperforms competing feature generating approaches for GZSL, and comes close to performance of fully supervised methods.


title: Enumerating Minimal Separators in Ranked Order

abstract: Let $G$ be an $n$-vertex graph, and $s,t$ vertices of $G$. We present an efficient algorithm which enumerates the set of minimal $st$-separators of $G$ in ascending order of cardinality, with a delay of $O(n^{3.5})$ per separator. In particular, we present an algorithm that lists, in ascending order of cardinality, all minimal separators with at most $k$ vertices. In that case, we show that the delay of the enumeration algorithm is $O(kn^{2.5})$ per separator. Our process is based on a new method that can decide, in polynomial time, whether the set of minimal separators under certain inclusion, exclusion, and cardinality constraints is empty.


title: The Possibilistic Horn Non-Clausal Knowledge Bases

abstract: Posibilistic logic is the most extended approach to handle uncertain and partially inconsistent information. Regarding normal forms, advances in possibilistic reasoning are mostly focused on clausal form. Yet, the encoding of real-world problems usually results in a non-clausal (NC) formula and NC-to-clausal translators produce severe drawbacks that heavily limit the practical performance of clausal reasoning. Thus, by computing formulas in its original NC form, we propose several contributions showing that notable advances are also possible in possibilistic non-clausal reasoning. {\em Firstly,} we define the class of {\em Possibilistic Horn Non-Clausal Knowledge Bases,} or $\mathcal{\overline{H}}_\Sigma$, which subsumes the classes: possibilistic Horn and propositional Horn-NC. $\mathcal{\overline{H}}_\Sigma $ is shown to be a kind of NC analogous of the standard Horn class. {\em Secondly}, we define {\em Possibilistic Non-Clausal Unit-Resolution,} or $ \mathcal{UR}_\Sigma $, and prove that $ \mathcal{UR}_\Sigma $ correctly computes the inconsistency degree of $\mathcal{\overline{H}}_\Sigma $members. $\mathcal{UR}_\Sigma $ had not been proposed before and is formulated in a clausal-like manner, which eases its understanding, formal proofs and future extension towards non-clausal resolution. {\em Thirdly}, we prove that computing the inconsistency degree of $\mathcal{\overline{H}}_\Sigma $ members takes polynomial time. Although there already exist tractable classes in possibilistic logic, all of them are clausal, and thus, $\mathcal{\overline{H}}_\Sigma $ turns out to be the first characterized polynomial non-clausal class within possibilistic reasoning.


title: Towards a Classification of Multi-Faced Independence: A Representation-Theoretic Approach

abstract: We attack the classification problem of multi-faced independences, the first non-trivial example being Voiculescu's bi-freeness. While the present paper does not achieve a complete classification, it formalizes the idea of lifting an operator on a pre-Hilbert space in a "universal" way to a larger product space, which is key for the construction of (old and new) examples. It will be shown how universal lifts can be used to construct very well-behaved (multi-faced) independences in general. Furthermore, we entirely classify universal lifts to the tensor product and to the free product of pre-Hilbert spaces. Our work brings to light surprising new examples of 2-faced independences. Most noteworthy, for many known 2-faced independences, we find that they admit continuous deformations within the class of 2-faced independences, showing in particular that, in contrast with the single faced case, this class is infinite (and even uncountable).


title: Joint FCLT for Sample Quantile and Measures of Dispersion for Functionals of Mixing Processes

abstract: In this paper, we establish a joint (bivariate) functional central limit theorem of the sample quantile and the $r$-th absolute centred sample moment for functionals of mixing processes. More precisely, we consider $L_2$-near epoch dependent processes that are functionals of either $\phi$-mixing or absolutely regular processes. The general results we obtain can be used for two classes of popular and important processes in applications: The class of augmented GARCH($p$,$q$) processes with independent and identically distributed innovations (including many GARCH variations used in practice) and the class of ARMA($p$,$q$) processes with mixing innovations (including, e.g., ARMA-GARCH processes). For selected examples, we provide exact conditions on the moments and parameters of the process for the joint asymptotics to hold.


title: On maximal solvable extensions of a pure non-characteristically nilpotent Lie algebra

abstract: In this paper we introduce the notion of pure non-characteristically nilpotent Lie algebra and under a condition we prove that a complex maximal extension of a finite-dimensional pure non-characteristically nilpotent Lie algebra is isomorphic to a semidirect sum of the nilradical and its maximal torus. We also prove that such solvable Lie algebras are complete and we specify a subclass of the maximal solvable extensions of pure non-characteristically nilpotent Lie algebras that have trivial cohomology group. Some comparisons with the results obtained earlier are given.


title: Comparison of heavy ion transport simulations: Ag+Ag collisions at Elab = 1.58 AGeV

abstract: We compare the microscopic transport models UrQMD, PHSD, PHQMD, and SMASH to make predictions for the upcoming Ag+Ag data at $E_\mathrm{lab}=1.58$~$A$GeV ($\sqrt{s_\mathrm{NN}}=2.55$~GeV) by the HADES collaboration. We study multiplicities, spectra and effective source temperatures of protons, $\pi^{\pm,0}$, $K^\pm$, the $\eta$, $\Lambda+\Sigma^0$ and the $\Xi^-$ within these models. Despite variations in the detailed implementation of the dynamics in the different models, the employed transport approaches all show consistent multiplicities of the bulk of investigated hadrons. The main differences are in the $\Xi^-$ production, which is treated differently between UrQMD/SMASH on one side employing high mass resonance states with explicit decays to $\mathrm{Resonance}\rightarrow \Xi+K+K$ in contrast to PHSD/PHQMD which account only non-resonant $\Xi$ production channels. A comparison of the spectra, summarized by effective source temperatures, shows that all models provide similar source temperatures around $T_\mathrm{source}=80-95$~MeV, and show substantial radial flow on the order of $\langle v_T\rangle=0.22c-0.3c$ even for such a small system.


title: Is $\chi_{c1}(3872)$ generated from string breaking?

abstract: We show, from a diabatic analysis of Lattice results for string breaking, that mixing of $Q\bar{Q}$ with open-flavor meson-meson configurations may be expressed through a mixing potential which is order $1 / m_{Q}$. A relation between the minimum string breaking energy gap and the string tension comes out naturally. Using this relation, and matching the energy gap for $b\bar{b}$ with Lattice QCD data, we study the mixing in the $c\bar{c}$ case without any additional parameter. A $1^{++}$ bound state very close below the ${D}^{0}\bar{D}^{\ast 0}$ threshold, in perfect correspondence with $\chi_{c1}(3872)$, is predicted.


title: On Russell typicality in Set Theory

abstract: By Tzouvaras, a set is nontypical in the Russell sense, if it belongs to a countable ordinal definable set. The class HNT of all hereditarily nontypical sets satisfies all axioms of ZF and the double inclusion HOD$\subseteq$HNT$\subseteq$V holds. Several questions about the nature of such sets, recently proposed by Tzouvaras, are solved in this paper. In particular, a model of ZFC is presented in which HOD$\subsetneqq$HNT$\subsetneqq$V, and another model of ZFC, in which HNT does not satisfy the axiom of choice.


title: Lack of other molecules in CO-rich debris discs: is it primordial or secondary gas?

abstract: The nature of the gas in CO-rich debris discs remains poorly understood, as it could either be a remnant from the earlier Class II phase or of secondary origin, driven by the destruction of icy planetesimals. The aim of this paper was to elucidate the origin of the gas content in the debris discs via various simple molecules that are often detected in the less-evolved Class II discs. We present millimetre molecular line observations of nine circumstellar discs around A-type stars: four CO-rich debris discs (HD 21997, HD 121617, HD 131488, HD 131835) and five old Herbig Ae protoplanetary discs (HD 139614, HD 141569, HD 142666, HD 145718, HD 100453). The sources were observed with the Atacama Large Millimeter/submillimeter Array (ALMA) in Bands 5 and 6 with 1-2'' resolution. The Herbig Ae discs are detected in the CO isotopologues, CN, HCN, HCO$^+$, C2H, and CS lines. In contrast, only CO isotopologues are detected in the debris discs, showing a similar amount of CO to that found in the Herbig Ae protoplanetary discs. Using chemical and radiative transfer modelling, we show that the abundances of molecules other than CO in debris discs are expected to be very low. We consider multiple sets of initial elemental abundances with various degrees of H2 depletion. We find that the HCO+ lines should be the second brightest after the CO lines and that their intensities strongly depend on the overall CO/H$_2$ ratio of the gas. However, even in the ISM-like scenario, the simulated HCO+ emission remains weak as required by our non-detections.


title: Adaptive VEM: Stabilization-Free A Posteriori Error Analysis and Contraction Property

abstract: In the present paper we initiate the challenging task of building a mathematically sound theory for Adaptive Virtual Element Methods (AVEMs). Among the realm of polygonal meshes, we restrict our analysis to triangular meshes with hanging nodes in 2d -- the simplest meshes with a systematic refinement procedure that preserves shape regularity and optimal complexity. A major challenge in the a posteriori error analysis of AVEMs is the presence of the stabilization term, which is of the same order as the residual-type error estimator but prevents the equivalence of the latter with the energy error. Under the assumption that any chain of recursively created hanging nodes has uniformly bounded length, we show that the stabilization term can be made arbitrarily small relative to the error estimator provided the stabilization parameter of the scheme is sufficiently large. This quantitative estimate leads to stabilization-free upper and lower a posteriori bounds for the energy error. This novel and crucial property of VEMs hinges on the largest subspace of continuous piecewise linear functions and the delicate interplay between its coarser scales and the finer ones of the VEM space. An important consequence for piecewise constant data is a contraction property between consecutive loops of AVEMs, which we also prove. Our results apply to $H^1$-conforming (lowest order) VEMs of any kind, including the classical and enhanced VEMs.


title: Symbolic Music Loop Generation with VQ-VAE

abstract: Music is a repetition of patterns and rhythms. It can be composed by repeating a certain number of bars in a structured way. In this paper, the objective is to generate a loop of 8 bars that can be used as a building block of music. Even considering musical diversity, we assume that music patterns familiar to humans can be defined in a finite set. With explicit rules to extract loops from music, we found that discrete representations are sufficient to model symbolic music sequences. Among VAE family, musical properties from VQ-VAE are better observed rather than other models. Further, to emphasize musical structure, we have manipulated discrete latent features to be repetitive so that the properties are more strengthened. Quantitative and qualitative experiments are extensively conducted to verify our assumptions.


title: Calculating Question Similarity is Enough: A New Method for KBQA Tasks

abstract: Knowledge Base Question Answering (KBQA) aims to answer natural language questions with the help of an external knowledge base. The core idea is to find the link between the internal knowledge behind questions and known triples of the knowledge base. Traditional KBQA task pipelines contain several steps, including entity recognition, entity linking, answering selection, etc. In this kind of pipeline methods, errors in any procedure will inevitably propagate to the final prediction. To address this challenge, this paper proposes a Corpus Generation - Retrieve Method (CGRM) with Pre-training Language Model (PLM) for the KBQA task. The major novelty lies in the design of the new method, wherein our approach, the knowledge enhanced T5 (kT5) model aims to generate natural language QA pairs based on Knowledge Graph triples and directly solve the QA by retrieving the synthetic dataset. The new method can extract more information about the entities from PLM to improve accuracy and simplify the processes. We test our method on NLPCC-ICCPOL 2016 KBQA dataset, and the results show that our method improves the performance of KBQA and the out straight-forward method is competitive with the state-of-the-art.


title: Role of Time-Reversal Symmetry in the Dynamical Response of "One-Way" Nonlinear Devices

abstract: We study the role of time-reversal symmetry on the dynamical response of nonlinear optical systems that behave as unidirectional ("one-way") devices. It is shown that lossless nonlinear materials, despite being nonreciprocal, are typically time-reversal invariant. This property raises an apparent paradox because time-reversal invariant systems are forcibly bi-directional. Here, we present a solution for this conundrum, and theoretically explain why the "one-way" behavior can indeed be compatible with the time-reversal invariance. It is found that in the time-reversed problem the incident waves have a variation in time that is generally incompatible with the adiabatic approximation. Due to this reason the adiabatic approximation fails to predict the bi-directional nature of nonlinear system. We discuss the implications of this finding in the performance of practical nonlinear "one-way" devices.


title: Long-time behavior of several point particles in a 1D viscous compressible fluid

abstract: We study the long-time behavior of \textit{several} point particles in a 1D viscous compressible fluid. It is shown that the velocities of the point particles all obey the power law $t^{-3/2}$. This result extends author's previous works on the long-time behavior of a \textit{single} point particle. New difficulties arise in the derivation of pointwise estimates of Green's functions due to infinite reflections of waves in-between the point particles. In particular, the differential equation technique used in previous works alone does not suffice. We overcome this by carefully analyzing the structure of Green's functions in the Laplace variable, especially their asymptotic and analyticity properties.


title: Colossal Brownian yet non-Gaussian diffusion in a periodic potential: impact of nonequilibrium noise amplitude statistics

abstract: Last year in [Phys. Rev. E 102, 042121 (2020)] the authors studied an overdamped dynamics of nonequilibrium noise driven Brownian particle dwelling in a spatially periodic potential and discovered a novel class of Brownian, yet non-Gaussian diffusion. The mean square displacement of the particle grows linearly with time and the probability density for the particle position is Gaussian, however, the corresponding distribution for the increments is non-Gaussian. The latter property induces the colossal enhancement of diffusion, significantly exceeding the well known effect of giant diffusion. Here we considerably extend the above predictions by investigating the influence of nonequilibrium noise amplitude statistics on the colossal Brownian, yet non-Gaussian diffusion. The tail of amplitude distribution crucially impacts both the magnitude of diffusion amplification as well as Gaussianity of the position and increments statistics. Our results carry profound consequences for diffusive behaviour in nonequilibrium settings such as living cells in which diffusion is a central transport mechanism.


title: How does the background atmosphere affect the onset of the runaway greenhouse ?

abstract: As the insolation of an Earth-like (exo)planet with a large amount of water increases, its surface and atmospheric temperatures also increase, eventually leading to a catastrophic runaway greenhouse transition. While some studies have shown that the onset of the runaway greenhouse may be delayed due to an overshoot of the outgoing longwave radiation (OLR) -- compared to the Simpson-Nakajima threshold -- by radiatively inactive gases, there is still no consensus on whether this is occurring and why. Here, we used a suite of 1D radiative-convective models to study the runaway greenhouse transition, with particular emphasis on taking into account the radical change in the amount of water vapour (from trace gas to dominant gas). The aim of this work is twofold: first, to determine the most important physical processes and parametrisations affecting the OLR; and second, to propose reference OLR curves for N$_2$+H$_2$O atmospheres. Through multiple sensitivity tests, we list and select the main important physical processes and parametrisations that need to be accounted for in 1D radiative-convective models to compute an accurate estimate of the OLR for N$_2$+H$_2$O atmospheres. The reference OLR curve is computed with a 1D model built according to the sensitivity tests. These tests also allow us to interpret the diversity of results already published in the literature. Moreover, we provide a correlated-k table able to reproduce line-by-line calculations with high accuracy. We find that the transition between an N$_2$-dominated atmosphere and an H$_2$O-dominated atmosphere induces an overshoot of the OLR compared to the (pure H$_2$O) Simpson-Nakajima asymptotic limit. This overshoot is first due to a transition between foreign and self-broadening of the water absorption lines, and second to a transition between dry and moist adiabatic lapse rates.


title: Shifts of prepotentials (with an appendix by Michele Vergne)

abstract: We study the dynamics of supersymmetric theories in five dimensions obtained by compactifications of M-theory on a Calabi-Yau threefold X. For a compact X, this is determined by the geometry of X, in particular the Kahler class dependence of the volume of X determines the effective couplings of vector multiplets. Rigid supersymmetry emerges in the limit of divergent volume, prompting the study of the structure of Duistermaat-Heckman formula and its generalizations for non-compact toric Kahler manifolds. Our main tool is the set of finite-difference equations obeyed by equivariant volumes and their quantum versions. We also discuss a physical application of these equations in the context of seven-dimensional gauge theories, extending and clarifying our previous results. The appendix by M. Vergne provides an alternative local proof of the shift equation.


title: Improved Gibbs samplers for Cosmic Microwave Background power spectrum estimation

abstract: We study different variants of the Gibbs sampler algorithm from the perspective of their applicability to the estimation of power spectra of the cosmic microwave background (CMB) anisotropies. These include approaches studied earlier in the CMB literature as well as new ones which are proposed in this work. We demonstrate all these variants on full and cut sky simulations and compare their performance, assessing both their computational and statistical efficiency. For this we employ a consistent comparison metric, an effective sample size (ESS) per second, commonly used in this context in the statistical literature. We show that one of the proposed approaches, referred to as Centered overrelax, which capitalizes on additional, auxiliary variables to minimize computational time needed per sample, and uses overrelaxation to decorrelate subsequent samples, performs better than the standard Gibbs sampler by a factor between one and two orders of magnitude in the nearly full-sky, satellite-like cases. It therefore potentially provides an interesting alternative to the currently favored approaches.


title: Instant tissue field and magnetic susceptibility mapping from MR raw phase using Laplacian enabled deep neural networks

abstract: Quantitative susceptibility mapping (QSM) is a valuable MRI post-processing technique that quantifies the magnetic susceptibility of body tissue from phase data. However, the traditional QSM reconstruction pipeline involves multiple non-trivial steps, including phase unwrapping, background field removal, and dipole inversion. These intermediate steps not only increase the reconstruction time but amplify noise and errors. This study develops a large-stencil Laplacian preprocessed deep learning-based neural network for near instant quantitative field and susceptibility mapping (i.e., iQFM and iQSM) from raw MR phase data. The proposed iQFM and iQSM methods were compared with established reconstruction pipelines on simulated and in vivo datasets. In addition, experiments on patients with intracranial hemorrhage and multiple sclerosis were also performed to test the generalization of the novel neural networks. The proposed iQFM and iQSM methods yielded comparable results to multi-step methods in healthy subjects while dramatically improving reconstruction accuracies on intracranial hemorrhages with large susceptibilities. The reconstruction time was also substantially shortened from minutes using multi-step methods to only 30 milliseconds using the trained iQFM and iQSM neural networks.


title: Enhancement of the efficiency of acousto-optic Bragg diffraction due to optical activity. A case of Pb$_5$Ge$_3$O$_{11}$ crystals

abstract: We show that the existence of optical activity in an optical material can lead to essential enhancement of acousto-optic (AO) figure of merit for this material. The reason is that the ellipticity of interacting optical eigenwaves approaches unity near the optic axis and so additional elasto-optic (EO) tensor components with relatively large values become involved into the effective EO coefficient. We demonstrate on the example of lead germanate crystals, Pb$_5$Ge$_3$O$_{11}$, that the increase in the efficiency of AO diffraction manifests itself for all the types of isotropic and anisotropic interactions, whenever the incident optical wave propagates close to the optic axis. We find that, in the particular case of diffraction in the interaction plane XZ of Pb$_5$Ge$_3$O$_{11}$ crystals, the maximal enhancement of the AO figure of merit takes place under conditions of the types V and VI of isotropic diffraction, with the AO figure of merit increasing from zero up to 13.3x10$^{-15}$ s$^3$/kg, and the type IX of anisotropic diffraction when the AO figure of merit increases more than twice (from 12.5x10$^{-15}$ up to 26.5x10$^{-15}$ s$^3$/kg). The maximal AO efficiency in the XZ interaction plane is reached at the types I and II of isotropic AO interactions. In these cases the AO figure of merit increases from 6.8x10$^{-15}$ up to 37.9x10$^{-15}$ s$^3$/kg and from 31.1x10$^{-15}$ to 37.9x10$^{-15}$ s$^3$/kg, respectively.


title: Versatile Inverse Reinforcement Learning via Cumulative Rewards

abstract: Inverse Reinforcement Learning infers a reward function from expert demonstrations, aiming to encode the behavior and intentions of the expert. Current approaches usually do this with generative and uni-modal models, meaning that they encode a single behavior. In the common setting, where there are various solutions to a problem and the experts show versatile behavior this severely limits the generalization capabilities of these methods. We propose a novel method for Inverse Reinforcement Learning that overcomes these problems by formulating the recovered reward as a sum of iteratively trained discriminators. We show on simulated tasks that our approach is able to recover general, high-quality reward functions and produces policies of the same quality as behavioral cloning approaches designed for versatile behavior.


title: Fast Axiomatic Attribution for Neural Networks

abstract: Mitigating the dependence on spurious correlations present in the training dataset is a quickly emerging and important topic of deep learning. Recent approaches include priors on the feature attribution of a deep neural network (DNN) into the training process to reduce the dependence on unwanted features. However, until now one needed to trade off high-quality attributions, satisfying desirable axioms, against the time required to compute them. This in turn either led to long training times or ineffective attribution priors. In this work, we break this trade-off by considering a special class of efficiently axiomatically attributable DNNs for which an axiomatic feature attribution can be computed with only a single forward/backward pass. We formally prove that nonnegatively homogeneous DNNs, here termed $\mathcal{X}$-DNNs, are efficiently axiomatically attributable and show that they can be effortlessly constructed from a wide range of regular DNNs by simply removing the bias term of each layer. Various experiments demonstrate the advantages of $\mathcal{X}$-DNNs, beating state-of-the-art generic attribution methods on regular DNNs for training with attribution priors.


title: Local minimality of $\mathbb{R}^N$-valued and $\mathbb{S}^N$-valued Ginzburg-Landau vortex solutions in the unit ball $B^N$

abstract: We study the existence, uniqueness and minimality of critical points of the form $m_{\varepsilon,\eta}(x) = (f_{\varepsilon,\eta}(|x|)\frac{x}{|x|}, g_{\varepsilon,\eta}(|x|))$ of the functional \[ E_{\varepsilon,\eta}[m] = \int_{B^N} \Big[\frac{1}{2} |\nabla m|^2 + \frac{1}{2\varepsilon^2} (1 - |m|^2)^2 + \frac{1}{2\eta^2} m_{N+1}^2\Big]\,dx \] for $m=(m_1, \dots, m_N, m_{N+1}) \in H^1(B^N,\mathbb{R}^{N+1})$ with $m(x) = (x,0)$ on $\partial B^N$. We establish a necessary and sufficient condition on the dimension $N$ and the parameters $\varepsilon$ and $\eta$ for the existence of an escaping vortex solution $(f_{\varepsilon,\eta}, g_{\varepsilon,\eta})$ with $g_{\varepsilon,\eta}> 0$. We also establish its uniqueness and local minimality. In the limiting case $\eta = 0$, we prove the local minimality of the degree-one vortex solution for the Ginzburg-Landau (GL) energy for every $\varepsilon > 0$ and $N \geq 2$. Similarly, when $\varepsilon = 0$, we prove the local minimality of the degree-one escaping vortex solution to an $\mathbb{S}^N$-valued GL model arising in micromagnetics for every $\eta > 0$ and $2 \leq N \leq 6$.


title: Where to Drill Next? A Dual-Weighted Approach to Adaptive Optimal Design of Groundwater Surveys

abstract: We present a novel approach to adaptive optimal design of groundwater surveys - a methodology for choosing the location of the next monitoring well. Our dual-weighted approach borrows ideas from Bayesian Optimisation and goal-oriented error estimation to propose the next monitoring well, given that some data is already available from existing wells. Our method is distinct from other optimal design strategies in that it does not rely on Fisher Information and it instead directly exploits the posterior uncertainty and the expected solution to a dual (or adjoint) problem to construct an acquisition function that optimally reduces the uncertainty in the model as a whole and some engineering quantity of interest in particular. We demonstrate our approach in the context of 2D groundwater flow example and show that the dual-weighted approach outperforms the baseline approach with respect to reducing the error in the posterior estimate of the quantity of interest.


title: NeuralPDE: Modelling Dynamical Systems from Data

abstract: Many physical processes such as weather phenomena or fluid mechanics are governed by partial differential equations (PDEs). Modelling such dynamical systems using Neural Networks is an active research field. However, current methods are still very limited, as they do not exploit the knowledge about the dynamical nature of the system, require extensive prior knowledge about the governing equations or are limited to linear or first-order equations. In this work we make the observation that the Method of Lines used to solve PDEs can be represented using convolutions which makes convolutional neural networks (CNNs) the natural choice to parametrize arbitrary PDE dynamics. We combine this parametrization with differentiable ODE solvers to form the NeuralPDE Model, which explicitly takes into account the fact that the data is governed by differential equations. We show in several experiments on toy and real-world data that our model consistently outperforms state-of-the-art models used to learn dynamical systems.


title: A Data Quarantine Model to Secure Data in Edge Computing

abstract: Edge computing provides an agile data processing platform for latency-sensitive and communication-intensive applications through a decentralized cloud and geographically distributed edge nodes. Gaining centralized control over the edge nodes can be challenging due to security issues and threats. Among several security issues, data integrity attacks can lead to inconsistent data and intrude edge data analytics. Further intensification of the attack makes it challenging to mitigate and identify the root cause. Therefore, this paper proposes a new concept of data quarantine model to mitigate data integrity attacks by quarantining intruders. The efficient security solutions in cloud, ad-hoc networks, and computer systems using quarantine have motivated adopting it in edge computing. The data acquisition edge nodes identify the intruders and quarantine all the suspected devices through dimensionality reduction. During quarantine, the proposed concept builds the reputation scores to determine the falsely identified legitimate devices and sanitize their affected data to regain data integrity. As a preliminary investigation, this work identifies an appropriate machine learning method, Linear Discriminant Analysis (LDA), for dimensionality reduction. The LDA results in 72.83% quarantine accuracy and 0.9 seconds training time, which is efficient than other state-of-the-art methods. In future, this would be implemented and validated with ground truth data.


title: Improving needle visibility in LED-based photoacoustic imaging using deep learning with semi-synthetic datasets

abstract: Photoacoustic imaging has shown great potential for guiding minimally invasive procedures by accurate identification of critical tissue targets and invasive medical devices (such as metallic needles). The use of light emitting diodes (LEDs) as the excitation light sources accelerates its clinical translation owing to its high affordability and portability. However, needle visibility in LED-based photoacoustic imaging is compromised primarily due to its low optical fluence. In this work, we propose a deep learning framework based on U-Net to improve the visibility of clinical metallic needles with a LED-based photoacoustic and ultrasound imaging system. To address the complexity of capturing ground truth for real data and the poor realism of purely simulated data, this framework included the generation of semi-synthetic training datasets combining both simulated data to represent features from the needles and in vivo measurements for tissue background. Evaluation of the trained neural network was performed with needle insertions into blood-vessel-mimicking phantoms, pork joint tissue ex vivo and measurements on human volunteers. This deep learning-based framework substantially improved the needle visibility in photoacoustic imaging in vivo compared to conventional reconstruction by suppressing background noise and image artefacts, achieving 5.8 and 4.5 times improvements in signal-to-noise ratio (SNR) and the modified Hausdorff distance (MHD) respectively. Thus, the proposed framework could be helpful for reducing complications during percutaneous needle insertions by accurate identification of clinical needles in photoacoustic imaging.


title: Reliability Assessment for Distribution Systems with Embedded Microgrids

abstract: The electrical distribution system is moving toward a more decentralized, complex, and dynamic system. The system is experiencing a higher penetration of distributed energy resources, flexible resources, and active end-users, leading to an active distribution system. If these components are used actively, they can have a positive effect on the distribution system's reliability. This paper aims to investigate how a microgrid with renewable energy sources and energy storage might influence the reliability of electricity supply in a radially operated distribution system. The reliability is investigated from both the distribution system and the microgrid perspective with the application of different scenarios. The study includes an extensive sensitivity analysis. In this study, we use our developed reliability assessment tool for distribution networks. The tool is an open available software where the example network and datasets are embedded. The reliability assessment tool is based on Monte Carlo simulations where load flow calculations are included to capture the behavior of the system and system components. In this paper, the reliability tool is extended to include microgrids and microgrid controller opportunities that handle interaction strategies with the distribution network. The model is demonstrated on the IEEE 33-bus network. The result is confirmed through statistical testing showing the statistical significance of providing support from the microgrid on the distribution system's reliability.


title: Emergent Weyl Fermions in an Orbital Multipolar Ordering Phase

abstract: Multipolar orderings in degenerate orbital systems offer unique opportunities for emergent topological phases. The phase diagram of interacting spinless fermions in a $p$-band diamond lattice at unit filling is first studied to elucidate the essential role of orbital multipolar orderings in the evolution of multifold degenerate band nodes. The free band structure around the Brillouin zone center is described by two quadratic band nodes each with a threefold degeneracy, which are spanned by the bonding and anti-bonding $p$-orbital multiplets, respectively. Upon switching on interactions, the triply degenerate band node is split into a pair of Weyl fermions with opposite chirality due to the onset of orbital multipolar orderings. Further raising interactions ultimately drives the system into an insulating phase with the orbital quadrupolar ordering. Our study is then generalized to spin-$1/2$ fermions, which has direct relevance with solid-state materials. The system develops full spin polarization through a ferromagnetic transition at tiny interactions, leaving the remaining orbital sector activated. The ensuing transitions take place in the orbital sector as a natural consequence, qualitatively recovering the phase diagram of spinless fermions. Our findings shed new light on the realization of emergent novel fermions with a prospect being a frontier at the confluence of topology, orbital physics and strong correlation.


title: Electron recoil effect in electrically tunable MoSe2 monolayers

abstract: Radiative recombination of excitons dressed by the interactions with free charge carriers often occurs under simultaneous excitation of either electrons or holes to unbound states. This phenomenon, known as the electron recoil effect, manifests itself in pronounced, asymmetric spectral lineshapes of the resulting emission. We study the electron recoil effect experimentally in electrically-tunable monolayer semiconductors and derive it theoretically using both trion and Fermi-polaron pictures. Time-resolved analysis of the recoil lineshapes is employed to access transient, non-equilibrium states of the exciton-carrier complexes. We demonstrate cooling of the initially overheated populations on the picosecond timescales and reveal the impact of lattice temperature and free carrier density. Both thermally activated phonons and the presence of free charges are shown to accelerate equilibration. Finally, we find strong correlations between relaxation times from recoil analysis and luminescence rise times, providing a consistent interpretation for the initial dynamics of trion/Fermi-polaron states.


title: FastFlow: Unsupervised Anomaly Detection and Localization via 2D Normalizing Flows

abstract: Unsupervised anomaly detection and localization is crucial to the practical application when collecting and labeling sufficient anomaly data is infeasible. Most existing representation-based approaches extract normal image features with a deep convolutional neural network and characterize the corresponding distribution through non-parametric distribution estimation methods. The anomaly score is calculated by measuring the distance between the feature of the test image and the estimated distribution. However, current methods can not effectively map image features to a tractable base distribution and ignore the relationship between local and global features which are important to identify anomalies. To this end, we propose FastFlow implemented with 2D normalizing flows and use it as the probability distribution estimator. Our FastFlow can be used as a plug-in module with arbitrary deep feature extractors such as ResNet and vision transformer for unsupervised anomaly detection and localization. In training phase, FastFlow learns to transform the input visual feature into a tractable distribution and obtains the likelihood to recognize anomalies in inference phase. Extensive experimental results on the MVTec AD dataset show that FastFlow surpasses previous state-of-the-art methods in terms of accuracy and inference efficiency with various backbone networks. Our approach achieves 99.4% AUC in anomaly detection with high inference efficiency.


title: Enhanced production of $^{8}$Be nuclei in relativistic nuclei fragmentation

abstract: This article is dedicated to the results of early measuring fragmentation of relativistic $^{16}$O, $^{22}$Ne, $^{28}$Si, and $^{197}$Au nuclei in a nuclear track emulsion. It has been found that there is a contribution of the unstable $^8$Be nucleus decays to $\alpha$-particle multiplicities. These renewed measurements of nuclear track emulsion exposed to $^{84}$Kr nuclei at the energy of 950 MeV per nucleon have been analyzed and given below.


title: Contrastive Representation Learning with Trainable Augmentation Channel

abstract: In contrastive representation learning, data representation is trained so that it can classify the image instances even when the images are altered by augmentations. However, depending on the datasets, some augmentations can damage the information of the images beyond recognition, and such augmentations can result in collapsed representations. We present a partial solution to this problem by formalizing a stochastic encoding process in which there exist a tug-of-war between the data corruption introduced by the augmentations and the information preserved by the encoder. We show that, with the infoMax objective based on this framework, we can learn a data-dependent distribution of augmentations to avoid the collapse of the representation.


title: Quadratic speedup of global search using a biased crossover of two good solutions

abstract: The minimisation of cost functions is crucial in various optimisation fields. However, identifying their global minimum remains challenging owing to the huge computational cost incurred. This work analytically expresses the computational cost to identify an approximate global minimum for a class of cost functions defined under a high-dimensional discrete state space. Then, we derive an optimal global search scheme that minimises the computational cost. Mathematical analyses demonstrate that a combination of the gradient descent algorithm and the selection and crossover algorithm--with a biased crossover weight--maximises the search efficiency. Remarkably, its computational cost is of the square root order in contrast to that of the conventional gradient descent algorithms, indicating a quadratic speedup of global search. We corroborate this proposition using numerical analyses of the travelling salesman problem. The simple computational architecture and minimal computational cost of the proposed scheme are highly desirable for biological organisms and neuromorphic hardware.


title: Finite-time dynamical phase transition in non-equilibrium relaxation

abstract: We uncover a finite-time dynamical phase transition in the thermal relaxation of a mean-field magnetic model. The phase transition manifests itself as a cusp singularity in the probability distribution of the magnetisation that forms at a critical time. The transition is due to a sudden switch in the dynamics, characterised by a dynamical order parameter. We derive a dynamical Landau theory for the transition that applies to a range of systems with scalar, parity-invariant order parameters. Close to criticalilty, our theory reveals an exact mapping between the dynamical and equilibrium phase transitions of the magnetic model, and implies critical exponents of mean-field type. We argue that interactions between nearby saddle points, neglected at the mean-field level, may lead to critical, spatiotemporal fluctuations of the order parameter, and thus give rise to novel, dynamical critical phenomena.


title: Geometrical scaling for light flavor hadrons

abstract: As it is well known by now, the pre-partonic phase in hadron collisions is successfully described by the Color Glass Condensate (CGC) approach. Previous studies, based on experimental data obtained on a wide range of energies at the Relativistic Heavy Ion Collider (RHIC) and at the Large Hadron Collider (LHC) for $\pi^+$, $K^+$ and $p$, evidenced that observables related to the dynamics of the collision, i.e. the mean transverse momentum ($\langle p_T \rangle$), the slope of the $\langle p_T \rangle$ dependence on the mass of the hadrons and the average transverse flow velocity obtained from the simultaneous fits of the $p_T$ spectra of the particles with the Boltzmann-Gibbs Blast Wave (BGBW) expression, scale rather well as a function of the square root of the ratio of the particle density over unit of rapidity to the overlapping area of the colliding nuclei ($\sqrt{(dN/dy)/S_\perp}$), the relevant scale in the gluon saturation picture. Results of a similar study extended to strange and multi-strange hadrons, for both proton-proton (\textit{pp}) and heavy-ion (\textit{A-A}) collision systems are presented in the present paper. The similarities and differences in the behaviour of strange hadrons relative to non-strange hadrons are discussed.


title: Reachability analysis of neural networks using mixed monotonicity

abstract: This paper presents a new reachability analysis approach to compute interval over-approximations of the output set of feedforward neural networks with input uncertainty. We adapt to neural networks an existing mixed-monotonicity method for the reachability analysis of dynamical systems and apply it to each partial network within the main network. This ensures that the intersection of the obtained results is the tightest interval over-approximation of the output of each layer that can be obtained using mixed-monotonicity on any partial network decomposition. Unlike other tools in the literature focusing on small classes of piecewise-affine or monotone activation functions, the main strength of our approach is its generality: it can handle neural networks with any Lipschitz-continuous activation function. In addition, the simplicity of our framework allows users to very easily add unimplemented activation functions, by simply providing the function, its derivative and the global argmin and argmax of the derivative. Our algorithm is compared to five other interval-based tools (Interval Bound Propagation, ReluVal, Neurify, VeriNet, CROWN) on both existing benchmarks and two sets of small and large randomly generated networks for four activation functions (ReLU, TanH, ELU, SiLU).


title: AutoGMap: Learning to Map Large-scale Sparse Graphs on Memristive Crossbars

abstract: The sparse representation of graphs has shown great potential for accelerating the computation of graph applications (e.g., Social Networks, Knowledge Graphs) on traditional computing architectures (CPU, GPU, or TPU). But the exploration of large-scale sparse graph computing on processing-in-memory (PIM) platforms (typically with memristive crossbars) is still in its infancy. To implement the computation or storage of large-scale or batch graphs on memristive crossbars, a natural assumption is that a large-scale crossbar is demanded, but with low utilization. Some recent works question this assumption, to avoid the waste of storage and computational resource, the fixed-size or progressively scheduled ''block partition'' schemes are proposed. However, these methods are coarse-grained or static, and are not effectively sparsity-aware. This work proposes the dynamic sparsity-aware mapping scheme generating method that models the problem with a sequential decision-making model, and optimizes it by reinforcement learning (RL) algorithm (REINFORCE). Our generating model (LSTM, combined with the dynamic-fill scheme) generates remarkable mapping performance on a small-scale graph/matrix data (complete mapping costs 43% area of the original matrix) and two large-scale matrix data (costing 22.5% area on qh882 and 17.1% area on qh1484). Our method may be extended to sparse graph computing on other PIM architectures, not limited to the memristive device-based platforms.


title: Analyzing the Behavior and Financial Status of Soccer Fans from a Mobile Phone Network Perspective: Euro 2016, a Case Study

abstract: In this study, Call Detail Records (CDRs), covering Budapest, for the month of June in 2016 has been analyzed. During this observation period, the 2016 UEFA European Football Championship took place, which affected significantly the habit of the residents, despite the fact that not a single match was played in the city. We evaluated the fans' behavior in Budapest, during and after the Hungarian matches, and found that the mobile phone network activity reflects the football fans' behavior, demonstrating the potential of mobile phone network data within a social sensing system. The Call Detail Records are enriched with mobile phone properties to analyze the subscribers' devices. Applying the device information (Type Allocation Code) from the activity records, the Subscriber Identity Modules, that do not operate in cell phones are omitted from mobility analyses, allowing to focus on people. The mobile phone price is proposed and evaluated as a socioeconomic indicator, and correlation between the phone price and the mobility customs have been found. We also found that, beside the cell phone price, the subscriber age and the subscription type also have an effect on the mobility. On the other hand, these do not seem to affect the interest in football.


title: A unified strategy to compute some special functions of number-theoretic interest

abstract: We introduce an algorithm to compute the functions belonging to a suitable set ${\mathscr F}$ defined as follows: $f\in {\mathscr F}$ means that $f(s,x)$, $s\in A\subset {\mathbb R}$ being fixed and $x>0$, has a power series expansion centred at $x_0=1$ with convergence radius greater or equal than $1$; moreover, it satisfies a functional equation of step $1$ and the Euler-Maclaurin summation formula can be applied to $f$. Denoting the Euler gamma-function as $\Gamma$, we will show that, for $x>0$, $\log \Gamma(x)$, the digamma function $\psi(x)$, the polygamma functions $\psi^{(w)}(x)$, $w\in {\mathbb N}$, $w\ge1$, and, for $s>1$ being fixed, the Hurwitz $\zeta(s,x)$-function and its first partial derivative $\frac{\partial\zeta}{\partial s}(s,x)$ are in ${\mathscr F}$. In all these cases the coefficients of the involved power series will depend on the values of $\zeta(u)$, $u>1$, where $\zeta$ is the Riemann zeta-function. As a by-product, we will also show how to compute the Dirichlet $L$-functions $L(s,\chi)$ and $L^\prime(s,\chi)$, $s> 1$, $\chi$ being a primitive Dirichlet character, by inserting the reflection formulae of $\zeta(s,x)$ and $\frac{\partial\zeta}{\partial s}(s,x)$ into the first step of the Fast Fourier Transform algorithm. Moreover, we will obtain some new formulae and algorithms for the Dirichlet $\beta$-function and for the Catalan constant $G$. Finally, we will study the case of the Bateman $G$-function and of the alternating Hurwitz zeta-function, also known as the $\eta$-function; we will show that, even if they are not in ${\mathscr F}$, our approach can be adapted to handle them too. In the last section we will also describe some tests that show a performance gain with respect to a standard multiprecision implementation of $\zeta(s,x)$ and $\frac{\partial\zeta}{\partial s}(s,x)$, $s>1$, $x>0$.


title: The ROSMA-PTV concept and the impact of rotational setup uncertainties on the tumor control probability in canine brain tumors

abstract: Purpose: In this modelling study we pursued two main goals. The first was to establish a new CTV-to-PTV expansion which considers the closest and most critical organ at risk (OAR). The second goal was to investigate the impact of the planning target volume (PTV) margin size on the tumor control probability (TCP) dependent on the geometrical setup uncertainties. The aim was to achieve a smaller margin expansion close to the OAR while allowing a moderately larger expansion in less critical areas further away from the OAR and most importantly while maintaining the TCP. Methods and Materials: Imaging data of radiation therapy (RT) plans from pet dogs which had undergone RT for brain tumor were used to estimate the clinic specific rotational setup uncertainties. A Monte-Carlo methodology used to quantify the implications of those uncertainties on the TCP. A computational CTV-to-PTV expansion method based on probability density was established. All required software modules were developed and integrated into a software package which directly interacts with the Varian Eclipse treatment planning system. Results: The rotational setup uncertainties were obtained from 44 dog patients. Their frequency distributions were analysed, and the mean values and standard deviations were determined. To investigate the impact on the TCP, several uniform and non-isotropic PTV targets were created. Standardized RT plans with equal optimization constraints were defined and automatically applied and calculated on these targets. The resulting TCPs were then evaluated and compared. Conclusion: The non-isotropic margins were found to result in larger TCPs with smaller margin excess volume. Further we have presented an additional application of the newly established CTV-to-PTV expansion method for radiation therapy of the spinal axis of human patients.


title: SDSS J1451+2709 a normal blue quasar but mis-classified as a HII galaxy in the BPT diagram by flux ratios of narrow emission lines

abstract: In the manuscript, we discuss properties of SDSS J1451+2709, a normal blue quasar but mis-classified as a HII galaxy in the BPT diagram (called as a mis-classified quasar). The emission lines around H$\alpha$ and H$\beta$ are well measured by different model functions with broad Balmer lines being described by Gaussian or Lorentz functions, in the SDSS spectrum in 2007 and in the KPNO spectrum in 1990. After considering variations of broad emission lines, different model functions lead to different fluxes of narrow emission lines, but the different narrow emission line flux ratios lead the SDSS J1451+2709 to be classified as a HII galaxy in the BPT diagram. In order to explain the unique properties of the mis-classified quasar SDSS J1451+2709 in the BPT diagram, two methods are proposed, the starforming contributions and compressed NLRs with high electron densities near to critical densities of forbidden emission lines. Unfortunately, the two methods cannot be preferred in the SDSS J1451+2709, further efforts are necessary to find the physical origin of the unique properties of the mis-classified quasar SDSS J1451+2709 in the BPT diagram. Meanwhile, there are not quite different long-term variabilities of SDSS J1451+2709 from the normal quasars. The mis-classified quasar SDSS J1451+2709, an extremely unique case or a special case among the normal quasars, could provide further clues on the applications of BPT diagrams to the normal broad line AGN and to narrow emission line objects, indicating part of narrow emission line HII galaxies actually harbouring central AGN activities.


title: Levitating frogs, machine learning and elliptic integrals

abstract: We present the calculation of the stability region of a perfect diamagnet levitated in a magnetic field created by a circular current loop. We make use of the machine learning technique of automatic differentiation to illustrate the calculations with the elliptic integrals involved.


title: Determination of spin chirality using x-ray magnetic circular dichroism

abstract: A 3-fold symmetric kagome lattice that has negative spin chirality can give a non-zero x-ray magnetic circular dichroism (XMCD) signal, despite that the total spin moment amounts to zero. This is explained by a hitherto unnoticed rule for the rotational symmetry invariance of the XMCD signal. A necessary condition is the existence of an anisotropic XMCD signal for the local magnetic atom, which can arise from a spin anisotropy either in the ground state or the final state. The angular dependence of the XMCD as a function of the beam direction has an unusual behavior. The maximum dichroism is not aligned along the spin direction, but depends on the relative orientation of the spin with respect to the atomic orientation. Therefore, different geometries can result in the same angular dependence, and the spin direction can only be determined if the atomic orientation is known. The consequences for the x-ray magneto-optical sum rules are given. The integrated XMCD signals are proportional to the anisotropy in the orbital moment and the magnetic dipole term, where the isotropic spin moment drops out.


title: Theoretical Guarantees for the Statistical Finite Element Method

abstract: The statistical finite element method (StatFEM) is an emerging probabilistic method that allows observations of a physical system to be synthesised with the numerical solution of a PDE intended to describe it in a coherent statistical framework, to compensate for model error. This work presents a new theoretical analysis of the statistical finite element method demonstrating that it has similar convergence properties to the finite element method on which it is based. Our results constitute a bound on the Wasserstein-2 distance between the ideal prior and posterior and the StatFEM approximation thereof, and show that this distance converges at the same mesh-dependent rate as finite element solutions converge to the true solution. Several numerical examples are presented to demonstrate our theory, including an example which test the robustness of StatFEM when extended to nonlinear quantities of interest.


title: Fast adjoint differentiation of chaos via computing unstable perturbations of transfer operators

abstract: We devise the fast adjoint response algorithm for the gradient of physical measures (long-time-average statistics) of discrete-time hyperbolic chaos with respect to many system parameters. Its cost is independent of the number of parameters. The algorithm transforms our new theoretical tools, the adjoint shadowing lemma and the equivariant divergence formula, into the form of progressively computing $u$ many bounded vectors on one orbit. Here $u$ is the unstable dimension. We demonstrate our algorithm on an example difficult for previous methods, a system with random noise, and a system of a discontinuous map. We also give a short formal proof of the equivariant divergence formula. Compared to the better-known finite-element method, our algorithm is not cursed by dimensionality of the phase space (typical real-life systems have very high dimensions), since it samples by one orbit. Compared to the ensemble/stochastic method, our algorithm is not cursed by the butterfly effect, since the recursive relations in our algorithm is bounded.


title: A massless boundary component mode synthesis method for elastodynamic contact problems

abstract: We propose to combine the ideas of mass redistribution and component mode synthesis. More specifically, we employ the MacNeal method, which readily leads to a singular mass matrix, and an accordingly modified version of the Craig-Bampton method. Besides obtaining a massless boundary, we achieve a drastic reduction of the mathematical model order in this way compared to the parent finite element model. Contact is modeled using set-valued laws and time stepping is carried out with a semi-explicit scheme. We assess the method's computational performance by a series of benchmarks, including both frictionless and frictional contact. The results indicate that the proposed method achieves excellent energy conservation properties and superior convergence behavior. It reduces the spurious oscillations and decreases the computational effort by about 1-2 orders of magnitude compared to the current state of the art (mass-carrying component mode synthesis method). We believe that the computational performance and favorable energy conservation properties will be valuable for the prediction of vibro-impact processes and physical damping.