title: Triadic Exploration and Exploration with Multiple Experts

abstract: Formal Concept Analysis (FCA) provides a method called attribute exploration which helps a domain expert discover structural dependencies in knowledge domains that can be represented by a formal context (a cross table of objects and attributes). Triadic Concept Analysis is an extension of FCA that incorporates the notion of conditions. Many extensions and variants of attribute exploration have been studied but only few attempts at incorporating multiple experts have been made. In this paper we present triadic exploration based on Triadic Concept Analysis to explore conditional attribute implications in a triadic domain. We then adapt this approach to formulate attribute exploration with multiple experts that have different views on a domain.


title: Lattice Boltzmann Multicomponent Model for Direct-Writing Printing

abstract: We introduce a mesoscale approach for the simulation of multicomponent flows to model the direct-writing printing process, along with the early stage of ink deposition. As an application scenario, alginate solutions at different concentrations are numerically investigated alongside processing parameters, such as apparent viscosity, extrusion rate, and print head velocity. The present approach offers useful insights on the ink rheological effects upon printed products, susceptible to geometric accuracy and shear stress, by manufacturing processes such as the direct-writing printing for complex photonic circuitry, bio-scaffold fabrication, and tissue engineering.


title: The Incompressible Navier-Stokes-Fourier Limit from Boltzmann-Fermi-Dirac Equation

abstract: We study Boltzmann-Fermi-Dirac equation when quantum effects are taken into account in dilute gas dynamics. By employing new estimates on trilinear terms of collision kernels, we prove the global existence of the classical solution to Boltzmann-Fermi-Dirac equation near equilibrium. Furthermore, the limit from Boltzmann-Fermi-Dirac equation to incompressible Navier-Stokes-Fourier equations is justified rigorously. The corresponding formal analysis was given in the thesis of Zakrevskiy \cite{Zakrevskiy}


title: Remarks on holographic models of the Kerr-AdS$_{5}$ geometry

abstract: We study the low-temperature limit of scalar perturbations of the Kerr-AdS$_{5}$ black-hole for generic rotational parameters. We motivate the study by considering real-time holography of small black hole backgrounds. Using the isomonodromic technique, we show that corrections to the extremal limit can be encoded in the monodromy parameters of the Painlev\'e V transcendent, whose expansion is given in terms of irregular chiral conformal blocks. After discussing the contribution of the intermediate states to the quasi-normal modes, we perform a numerical analysis of the low-lying frequencies. We find that the fundamental mode is perturbatively stable at low temperatures for small black holes and that excited perturbations are superradiant, as expected from thermodynamical considerations. We close by considering the holographic interpretation of the unstable modes and the decaying process.


title: Measurement of the proton spin structure at long distances

abstract: Measuring the spin structure of protons and neutrons tests our understanding of how they arise from quarks and gluons, the fundamental building blocks of nuclear matter. At long distances the coupling constant of the strong interaction becomes large, requiring non-perturbative methods to calculate quantum chromodynamics processes, such as lattice gauge theory or effective field theories. Here we report proton spin structure measurements from scattering a polarized electron beam off polarized protons. The spin-dependent cross-sections were measured at large distances, corresponding to the region of low momentum transfer squared between 0.012 and 1.0 GeV$^2$. This kinematic range provides unique tests of chiral effective field theory predictions. Our results show that a complete description of the nucleon spin remains elusive, and call for further theoretical works, e.g. in lattice quantum chromodynamics. Finally, our data extrapolated to the photon point agree with the Gerasimov-Drell-Hearn sum rule, a fundamental prediction of quantum field theory that relates the anomalous magnetic moment of the proton to its integrated spin-dependent cross-sections.


title: An Evaluation of Cryptocurrency Payment Channel Networks and Their Privacy Implications

abstract: Cryptocurrencies redefined how money can be stored and transferred among users. However, independent of the amount being sent, public blockchain-based cryptocurrencies suffer from high transaction waiting times and fees. These drawbacks hinder the wide use of cryptocurrencies by masses. To address these challenges, payment channel network concept is touted as the most viable solution to be used for micro-payments. The idea is exchanging the ownership of money by keeping the state of the accounts locally. The users inform the blockchain rarely, which decreases the load on the blockchain. Specifically, payment channel networks can provide transaction approvals in seconds by charging a nominal fee proportional to the payment amount. Such attraction on payment channel networks inspired many recent studies which focus on how to design them and allocate channels such that the transactions will be secure and efficient. However, as payment channel networks are emerging and reaching large number of users, privacy issues are becoming more relevant that raise concerns about exposing not only individual habits but also businesses' revenues. In this paper, we first propose a categorization of the existing payment networks formed on top of blockchain-backed cryptocurrencies. After discussing several emerging attacks on user/business privacy in these payment channel networks, we qualitatively evaluate them based on a number of privacy metrics that relate to our case. Based on the discussions on the strengths and weaknesses of the approaches, we offer possible directions for research for the future of privacy based payment channel networks.


title: Exponential speedup of incoherent tunneling via dissipation

abstract: We study the escape rate of a particle in a metastable potential in presence of a dissipative bath coupled to the momentum of the particle. Using the semiclassical bounce technique, we find that this rate is exponentially enhanced. In particular, the influence of momentum dissipation depends on the slope of the barrier that the particle is tunneling through. We investigate also the influence of dissipative baths coupled to the position, and to the momentum of the particle, respectively. In this case the rate exhibits a non-monotonic behavior as a function of the dissipative coupling strengths. Remarkably, even in presence of position dissipation, momentum dissipation can enhance exponentially the escape rate in a large range of the parameter space. The influence of the momentum dissipation is also witnessed by the substantial increase of the average energy loss during inelastic (environment-assisted) tunneling.


title: Times of Arrival and Gauge Invariance

abstract: We revisit the arguments underlying two well-known arrival-time distributions in quantum mechanics, viz., the Aharonov-Bohm and Kijowski (ABK) distribution, applicable for freely moving particles, and the quantum flux (QF) distribution. An inconsistency in the original axiomatic derivation of Kijowski's result is pointed out, along with an inescapable consequence of the "negative arrival times" inherent to this proposal (and generalizations thereof). The ABK free-particle restriction is lifted in a discussion of an explicit arrival-time setup featuring a charged particle moving in a constant magnetic field. A natural generalization of the ABK distribution is in this case shown to be critically gauge-dependent. A direct comparison to the QF distribution, which does not exhibit this flaw, is drawn (its acknowledged drawback concerning the quantum backflow effect notwithstanding).


title: No-reference denoising of low-dose CT projections

abstract: Low-dose computed tomography (LDCT) became a clear trend in radiology with an aspiration to refrain from delivering excessive X-ray radiation to the patients. The reduction of the radiation dose decreases the risks to the patients but raises the noise level, affecting the quality of the images and their ultimate diagnostic value. One mitigation option is to consider pairs of low-dose and high-dose CT projections to train a denoising model using deep learning algorithms; however, such pairs are rarely available in practice. In this paper, we present a new self-supervised method for CT denoising. Unlike existing self-supervised approaches, the proposed method requires only noisy CT projections and exploits the connections between adjacent images. The experiments carried out on an LDCT dataset demonstrate that our method is almost as accurate as the supervised approach, while also outperforming the considered self-supervised denoising methods.


title: Asymptotics of some integrals involving modified Bessel and hyper-Bessel functions

abstract: We investigate the asymptotic expansion of integrals analogous to Ball's integral \[\int_0^\infty \left(\frac{\Gamma(1+\nu)|J_\nu(x)|}{(x/2)^\nu}\right)^{\!n}dx\] for large $n$ in which the Bessel function $J_\nu(x)$ is replaced by the modified Bessel functions $I_\nu(x)$ and $K_\nu(x)$ together with appropriate exponential factors $e^{\mp x}$, respectively. The above integral with $J_\nu(x)$ replaced by a hyper-Bessel function of the type recently discussed in Aktas {\it et al.} [The Ramanujan J., 2019] and taken over a finite interval determined by the first positive zero of the function is also considered for $n\to\infty$. We give the leading asymptotic behaviour of the hyper-Bessel function for $x\to+\infty$ in an appendix. Numerical examples are given to illustrate the accuracy of the various expansions obtained.


title: Digital twins based on bidirectional LSTM and GAN for modelling the COVID-19 pandemic

abstract: The outbreak of the coronavirus disease 2019 (COVID-19) has now spread throughout the globe infecting over 150 million people and causing the death of over 3.2 million people. Thus, there is an urgent need to study the dynamics of epidemiological models to gain a better understanding of how such diseases spread. While epidemiological models can be computationally expensive, recent advances in machine learning techniques have given rise to neural networks with the ability to learn and predict complex dynamics at reduced computational costs. Here we introduce two digital twins of a SEIRS model applied to an idealised town. The SEIRS model has been modified to take account of spatial variation and, where possible, the model parameters are based on official virus spreading data from the UK. We compare predictions from a data-corrected Bidirectional Long Short-Term Memory network and a predictive Generative Adversarial Network. The predictions given by these two frameworks are accurate when compared to the original SEIRS model data. Additionally, these frameworks are data-agnostic and could be applied to towns, idealised or real, in the UK or in other countries. Also, more compartments could be included in the SEIRS model, in order to study more realistic epidemiological behaviour.


title: Hybrid consistency and plausibility verification of product data according to FIC

abstract: The labelling of food products in the EU is regulated by the Food Information of Customers (FIC). Companies are required to provide the corresponding information regarding nutrients and allergens among others. With the rise of e-commerce more and more food products are sold online. There are often errors in the online product descriptions regarding the FIC-relevant information due to low data quality in the vendors' product data base. In this paper we propose a hybrid approach of both rule-based and machine learning to verify nutrient declaration and allergen labelling according to FIC requirements. Special focus is given to the problem of false negatives in allergen prediction since this poses a significant health risk to customers. Results show that a neural net trained on a subset of the ingredients of a product is capable of predicting the allergens contained with a high reliability.


title: The Wisdom of the Crowd and Higher-Order Beliefs

abstract: The classic wisdom-of-the-crowd problem asks how a principal can "aggregate" information about the unknown state of the world from agents without understanding the information structure among them. We propose a new simple procedure called Population-Mean-Based Aggregation to achieve this goal. The procedure only requires eliciting agents' beliefs about the state, and also eliciting some agents' expectations of the average belief in the population. We show that this procedure fully aggregates information: in an infinite population, it always infers the true state of the world. The procedure can accommodate correlations in agents' information, misspecified beliefs, any finite number of possible states of the world, and only requires very weak assumptions on the information structure.


title: Lyapunov exponent ISCO and Kolmogorov Senai entropy for Kerr Kiselev black hole

abstract: Geodesic motion has significant characteristics of space-time. We calculate the principle Lyapunov exponent (LE), which is the inverse of the instability timescale associated with this geodesics and Kolmogorov-Senai (KS) entropy for our rotating Kerr-Kiselev (KK) black hole. We have investigate the existence of stable/unstable equatorial circular orbits via LE and KS entropy for time-like and null circular geodesics. We have shown that both LE and KS entropy can be written in terms of the radial equation of innermost stable circular orbit (ISCO) for time-like circular orbit. Also, we computed the equation marginally bound circular orbit, which gives the radius (smallest real root) of marginally bound circular orbit (MBCO). We found that the null circular geodesics has larger angular frequency than time-like circular geodesics ($Q_o > Q_{\sigma}$). Thus, null-circular geodesics provides the fastest way to circulate KK black holes. Further, it is also to be noted that null circular geodesics has shortest orbital period $(T_{photon}< T_{ISCO})$ among the all possible circular geodesics. Even null circular geodesics traverses fastest than any stable time-like circular geodesics other than the ISCO.


title: Disease Prediction with a Maximum Entropy Method

abstract: In this paper, we propose a maximum entropy method for predicting disease risks. It is based on a patient's medical history with diseases coded in ICD-10 which can be used in various cases. The complete algorithm with strict mathematical derivation is given. We also present experimental results on a medical dataset, demonstrating that our method performs well in predicting future disease risks and achieves an accuracy rate twice that of the traditional method. We also perform a comorbidity analysis to reveal the intrinsic relation of diseases.


title: OmiEmbed: a unified multi-task deep learning framework for multi-omics data

abstract: High-dimensional omics data contains intrinsic biomedical information that is crucial for personalised medicine. Nevertheless, it is challenging to capture them from the genome-wide data due to the large number of molecular features and small number of available samples, which is also called 'the curse of dimensionality' in machine learning. To tackle this problem and pave the way for machine learning aided precision medicine, we proposed a unified multi-task deep learning framework named OmiEmbed to capture biomedical information from high-dimensional omics data with the deep embedding and downstream task modules. The deep embedding module learnt an omics embedding that mapped multiple omics data types into a latent space with lower dimensionality. Based on the new representation of multi-omics data, different downstream task modules were trained simultaneously and efficiently with the multi-task strategy to predict the comprehensive phenotype profile of each sample. OmiEmbed support multiple tasks for omics data including dimensionality reduction, tumour type classification, multi-omics integration, demographic and clinical feature reconstruction, and survival prediction. The framework outperformed other methods on all three types of downstream tasks and achieved better performance with the multi-task strategy comparing to training them individually. OmiEmbed is a powerful and unified framework that can be widely adapted to various application of high-dimensional omics data and has a great potential to facilitate more accurate and personalised clinical decision making.


title: Multimodal-Aware Weakly Supervised Metric Learning with Self-weighting Triplet Loss

abstract: In recent years, we have witnessed a surge of interests in learning a suitable distance metric from weakly supervised data. Most existing methods aim to pull all the similar samples closer while push the dissimilar ones as far as possible. However, when some classes of the dataset exhibit multimodal distribution, these goals conflict and thus can hardly be concurrently satisfied. Additionally, to ensure a valid metric, many methods require a repeated eigenvalue decomposition process, which is expensive and numerically unstable. Therefore, how to learn an appropriate distance metric from weakly supervised data remains an open but challenging problem. To address this issue, in this paper, we propose a novel weakly supervised metric learning algorithm, named MultimoDal Aware weakly supervised Metric Learning (MDaML). MDaML partitions the data space into several clusters and allocates the local cluster centers and weight for each sample. Then, combining it with the weighted triplet loss can further enhance the local separability, which encourages the local dissimilar samples to keep a large distance from the local similar samples. Meanwhile, MDaML casts the metric learning problem into an unconstrained optimization on the SPD manifold, which can be efficiently solved by Riemannian Conjugate Gradient Descent (RCGD). Extensive experiments conducted on 13 datasets validate the superiority of the proposed MDaML.


title: Directive Explanations for Actionable Explainability in Machine Learning Applications

abstract: This paper investigates the prospects of using directive explanations to assist people in achieving recourse of machine learning decisions. Directive explanations list which specific actions an individual needs to take to achieve their desired outcome. If a machine learning model makes a decision that is detrimental to an individual (e.g. denying a loan application), then it needs to both explain why it made that decision and also explain how the individual could obtain their desired outcome (if possible). At present, this is often done using counterfactual explanations, but such explanations generally do not tell individuals how to act. We assert that counterfactual explanations can be improved by explicitly providing people with actions they could use to achieve their desired goal. This paper makes two contributions. First, we present the results of an online study investigating people's perception of directive explanations. Second, we propose a conceptual model to generate such explanations. Our online study showed a significant preference for directive explanations ($p<0.001$). However, the participants' preferred explanation type was affected by multiple factors, such as individual preferences, social factors, and the feasibility of the directives. Our findings highlight the need for a human-centred and context-specific approach for creating directive explanations.


title: DNN Based Beam Selection in mmW Heterogeneous Networks

abstract: We consider a heterogeneous cellular network wherein multiple small cell millimeter wave (mmW) base stations (BSs) coexist with legacy sub-6GHz macro BSs. In the mmW band, small cells use multiple narrow beams to ensure sufficient coverage and User Equipments (UEs) have to select the best small cell and the best beam in order to access the network. This process usually based on exhaustive search may introduce unacceptable latency. In order to address this issue, we rely on the sub-6GHz macro BS support and propose a deep neural network (DNN) architecture that utilizes basic components from the Channel State Information (CSI) of sub-6GHz network as input features. The output of the DNN is the mmW BS and beam selection that can provide the best communication performance. In the set of features, we avoid using the UE location, which may not be readily available for every device. We formulate a mmW BS selection and beam selection problem as a classification and regression problem respectively and propose a joint solution using a branched neural network. The numerical comparison with the conventional exhaustive search results shows that the proposed design demonstrate better performance than exhaustive search in terms of latency with at least 85\% accuracy.


title: Continuous Random Variable Estimation is not Optimal for the Witsenhausen Counterexample

abstract: Optimal design of distributed decision policies can be a difficult task, illustrated by the famous Witsenhausen counterexample. In this paper we characterize the optimal control designs for the vector-valued setting assuming that it results in an internal state that can be described by a continuous random variable which has a probability density function. More specifically, we provide a genie-aided outer bound that relies on our previous results for empirical coordination problems. This solution turns out to be not optimal in general, since it consists of a time-sharing strategy between two linear schemes of specific power. It follows that the optimal decision strategy for the original scalar Witsenhausen problem must lead to an internal state that cannot be described by a continuous random variable which has a probability density function.


title: Spectral extrema of graphs with fixed size: cycles and complete bipartite graphs

abstract: Nikiforov [Some inequalities for the largest eigenvalue of a graph, Combin. Probab. Comput. 179--189] showed that if $G$ is $K_{r+1}$-free then the spectral radius $\rho(G)\leq\sqrt{2m(1-1/r)}$, which implies that $G$ contains $C_3$ if $\rho(G)>\sqrt{m}$. In this paper, we follow this direction on determining which subgraphs will be contained in $G$ if $\rho(G)> f(m)$, where $f(m)\sim\sqrt{m}$ as $m\rightarrow \infty$. We first show that if $\rho(G)\geq \sqrt{m}$, then $G$ contains $K_{2,r+1}$ unless $G$ is a star; and $G$ contains either $C_3^+$ or $C_4^+$ unless $G$ is a complete bipartite graph, where $C_t^+$ denotes the graph obtained from $C_t$ and $C_3$ by identifying an edge. Secondly, we prove that if $\rho(G)\geq{\frac12+\sqrt{m-\frac34}}$, then $G$ contains pentagon and hexagon unless $G$ is a book; and if $\rho(G)>{\frac12(k-\frac12)+\sqrt{m+\frac14(k-\frac12)^2}}$, then $G$ contains $C_t$ for every $t\leq 2k+2$. In the end, some related conjectures are provided for further research.


title: Observability of the relative motion from inertial data in kinematic chains

abstract: Real-time motion tracking of kinematic chains is a key prerequisite in the control of, e.g., robotic actuators and autonomous vehicles and also has numerous biomechanical applications. In recent years, it has been shown that, by placing inertial sensors on segments that are connected by rotational joints, the motion of that kinematic chain can be tracked accurately. These methods specifically avoid using magnetometer measurements, which are known to be unreliable since the magnetic field at the different sensor locations is typically different. They rely on the assumption that the motion of the kinematic chain is sufficiently rich to assure observability of the relative pose. However, a formal investigation of this crucial requirement has not yet been presented, and no specific conditions for observability have so far been given. In this work, we present an observability analysis and show that the relative pose of the body segments is indeed observable under a very mild condition on the motion. We support our results by simulation studies, in which we employ a state estimator that neither uses magnetometer measurements nor additional sensors and does not impose assumptions on the accelerometer to measure only the direction of gravity, nor on the range of motion or degrees of freedom of the joints. We investigate the effect of the amount of excitation and of stationary periods in the data on the accuracy of the estimates. We then use experimental data from two mechanical joints as well as from a human gait experiment to validate the observability criterion in practice and to show that small excitation levels are sufficient for obtaining accurate estimates even in the presence of time periods during which the motion is not observable.


title: Combinatorial and number-theoretic properties of generic reals

abstract: We discuss some properties of Cohen and random reals. We show that they belong to any definable partition regular family, and hence they satisfy most "largeness" properties studied in Ramsey theory. We determine their position in the Mahler's classification of the reals and using it, we get some information about Liouville numbers. We also show that they are wild in the sense of o-minimality, i.e., they define the set of integers.


title: BATTPOWER Application: Large-Scale Integration of EVs in an Active Distribution Grid -- A Norwegian Case Study

abstract: In this paper, we introduce an application of a high-performance MultiPeriod AC Optimal Power Flow (MPOPF) solver, called "BATTPOWER", to simulate active distribution grids for a near-future scenario. A large-scale Norwegian distribution grid along with a large population of Electric Vehicles (EV) are here taken as the case-study. We suggest and analyse three operational strategies (in terms of control of charge scheduling fleet of EV) for the Distribution System Operator (DSO): (a) uncoordinated/dumb charge scheduling, (b) coordinated charge scheduling with the objective of energy cost-minimisation without operational constraints of the grid, and (c) coordinated charge scheduling with the objective of energy cost-minimisation along with the operational constraints of the grid. The results demonstrate that the uncoordinated charging would lead to: 1) overloading of lines and transformers when the share of EVs is above 20%, and 2) higher operational costs than the proposed control strategies of (b) and (c). In strategy (b) operational line/transformer limits are violated when the populations of EVs are growing above 36%. This implies that current market design must be altered to allow active control of a large proportion of DERs within grid operational limits to achieve cost minimization at system level. To our knowledge, the work presented in this paper is the first ever attempt to do a comprehensive analysis of the impact of EV charging demand on a real Norwegian distribution grid. Moreover, the inference of the analysis says that the Norwegian distribution networks are more prone to congestion problems than the voltage problems for the EV demand which includes a smart charging scheme accounting for grid conditions.


title: Equivalence between angular spectrum-based and multipole expansion-based formulas of the acoustic radiation force and torque

abstract: Two main methods have been proposed to derive the acoustical radiation force and torque applied by an arbitrary acoustic field on a particle: The first one relies on the plane wave angular spectrum decomposition of the incident field (see [Sapozhnikov and Bailey, J. Acoust. Soc. Am. 133, 661 (2013)] for the force and [Gong and Baudoin, J. Acoust. Soc. Am. 148, 3131 (2020)] for the torque), while the second one relies on the decomposition of the incident field into a sum of spherical waves, the so-called multipole expansion (see [Silva, J. Acoust. Soc. Am. 130, 3541 (2011)] and [Baresh et al., J. Acoust. Soc. Am. 133, 25 (2013)] for the force, and [Silva et al., EPL 97, 54003 (2012)] and [Gong et al., Phys. Rev. Applied 11, 064022 (2019)] for the torque). In this paper, we formally establish the equivalence between the expressions obtained with these two methods for both the force and torque.


title: Certifying Differential Equation Solutions from Computer Algebra Systems in Isabelle/HOL

abstract: The Isabelle/HOL proof assistant has a powerful library for continuous analysis, which provides the foundation for verification of hybrid systems. However, Isabelle lacks automated proof support for continuous artifacts, which means that verification is often manual. In contrast, Computer Algebra Systems (CAS), such as Mathematica and SageMath, contain a wealth of efficient algorithms for matrices, differential equations, and other related artifacts. Nevertheless, these algorithms are not verified, and thus their outputs cannot, of themselves, be trusted for use in a safety critical system. In this paper we integrate two CAS systems into Isabelle, with the aim of certifying symbolic solutions to ordinary differential equations. This supports a verification technique that is both automated and trustworthy.


title: Hierarchical Multi-head Attentive Network for Evidence-aware Fake News Detection

abstract: The widespread of fake news and misinformation in various domains ranging from politics, economics to public health has posed an urgent need to automatically fact-check information. A recent trend in fake news detection is to utilize evidence from external sources. However, existing evidence-aware fake news detection methods focused on either only word-level attention or evidence-level attention, which may result in suboptimal performance. In this paper, we propose a Hierarchical Multi-head Attentive Network to fact-check textual claims. Our model jointly combines multi-head word-level attention and multi-head document-level attention, which aid explanation in both word-level and evidence-level. Experiments on two real-word datasets show that our model outperforms seven state-of-the-art baselines. Improvements over baselines are from 6\% to 18\%. Our source code and datasets are released at \texttt{\url{https://github.com/nguyenvo09/EACL2021}}.


title: The connection between warm carbon chain chemistry and interstellar irradiation of star-forming cores

abstract: Some observations of warm carbon chain chemistry (WCCC) cores indicate that they are often located near the edges of molecular clouds. This finding may suggest that WCCC is promoted in star-forming cores exposed to radiation from the interstellar medium. We aim to investigate the chemistry of carbon chains in such a core. A chemical simulation of a gas parcel in a low-mass star-forming core with a full level of irradiation by interstellar photons and cosmic rays was compared to a simulation of a core receiving only one-tenth of such irradiation. In the full irradiation model, the abundances of carbon chains were found to be higher by a factor of few to few hundred, compared to the model with low irradiation. Higher carbon-chain abundances in the prestellar stage and, presumably, in the extended circumstellar envelope, arise because of irradiation of gas and dust by interstellar photons and cosmic rays. A full standard rate of cosmic-ray induced ionization is essential for a high carbon-chain abundance peak to occur in the circumstellar envelope, which is heated by the protostar (the "true" WCCC phenomenon). The full irradiation model has lower abundances of complex organic molecules than the low-irradiation model. We conclude that WCCC can be caused by exposure of a star-forming core to interstellar radiation, or even just to cosmic rays. The Appendix describes an updated accurate approach for calculating the rate of cosmic-ray induced desorption.


title: Neural Network Potential Energy Surface for the low temperature Ring Polymer Molecular Dynamics of the H2CO + OH reaction

abstract: A new potential energy surface (PES) and dynamical study are presented of the reactive process between H2CO + OH towards the formation of HCO + H2O and HCOOH + H. In this work a source of spurious long range interactions in symmetry adapted neural network (NN) schemes is identified, what prevents their direct application for low temperature dynamical studies. For this reason, a partition of the PES into a diabatic matrix plus a NN many body term has been used fitted with a novel artificial neural networks scheme that prevents spurious asymptotic interactions. Quasi-classical trajectory and ring polymer molecular dynamics (RPMD) studies have been carried on this PES to evaluate the rate constant temperature dependence for the different reactive processes, showing a good agreement with the available experimental data. Of special interest is the analysis of the previously identified trapping mechanism in the RPMD study, which can be attributed to spurious resonances associated to excitations of the normal modes of the ring polymer.


title: A Novel Graphene-Based Circulator with Multi-layer Triangular Post for THz Region

abstract: This article proposes a novel three-port circulator with a triangular graphene-based post for the THz region. This new circulator is formed by three 120 symmetrical metal-based waveguides with a multi-layer triangular graphene-based post. The anisotropic feature for circulation is provided by magnetically-biased graphene in the triangular post. The DC magnetic bias is applied in the z-direction. The magnetized, triangular graphene sheet supports hybrid TM-TE plasmons. To realize the proposed circulator, the structure has been simulated in COMSOL software. In our simulation results, isolation of -40 dB with a transmission loss of -3.5 dB is obtained at the central frequency of 5 THz for a specific design. The bandwidth of the simulated circulator is reported 7.25% with respect to the isolation level of -15 dB. It has been shown that the scattering parameters of the proposed circulator can be changed by altering the chemical potential of the graphene and DC magnetic bias, which makes this circulator a tunable device to be utilized in various plasmonic systems.


title: Force-Directed Layout of Order Diagrams using Dimensional Reduction

abstract: Order diagrams allow human analysts to understand and analyze structural properties of ordered data. While an experienced expert can create easily readable order diagrams, the automatic generation of those remains a hard task. In this work, we adapt force-directed approaches, which are known to generate aesthetically-pleasing drawings of graphs, to the realm of order diagrams. Our algorithm ReDraw thereby embeds the order in a high dimension and then iteratively reduces the dimension until a two-dimensional drawing is achieved. To improve aesthetics, this reduction is equipped with two force-directed steps where one optimizes on distances of nodes and the other on distances of lines in order to satisfy a set of a priori fixed conditions. By respecting an invariant about the vertical position of the elements in each step of our algorithm we ensure that the resulting drawings satisfy all necessary properties of order diagrams. Finally, we present the results of a user study to demonstrate that our algorithm outperforms comparable approaches on drawings of lattices with a high degree of distributivity.


title: Impossibility of Partial Recovery in the Graph Alignment Problem

abstract: Random graph alignment refers to recovering the underlying vertex correspondence between two random graphs with correlated edges. This can be viewed as an average-case and noisy version of the well-known graph isomorphism problem. For the correlated Erd\"os-R\'enyi model, we prove an impossibility result for partial recovery in the sparse regime, with constant average degree and correlation, as well as a general bound on the maximal reachable overlap. Our bound is tight in the noiseless case (the graph isomorphism problem) and we conjecture that it is still tight with noise. Our proof technique relies on a careful application of the probabilistic method to build automorphisms between tree components of a subcritical Erd\"os-R\'enyi graph.


title: Variation of Stratifications From Toric GIT

abstract: When a reductive group acts on an algebraic variety, a linearized ample line bundle induces a stratification on the variety where the strata are ordered by the degrees of instability. In this paper, we study variation of stratifications coming from the group actions in the GIT quotient construction for projective toric varieties. Cox showed that each projective toric variety is a GIT quotient of an affine space by a diagonalizable group with respect to linearizations that come from ample divisors on the toric variety. We provide a sufficient conditions for two ample divisors to induce the same stratification and formulate two types of walls in the ample cone that completely capture two kinds of variations. We also prove that the variation is intrinsic to the primitive collections and the relations among ray generators of the fans.


title: Regular integral models for Shimura varieties of orthogonal type

abstract: We consider Shimura varieties for orthogonal or spin groups acting on hermitian symmetric domains of type IV. We give regular p-adic integral models for these varieties over odd primes p at which the level subgroup is the connected stabilizer of a vertex lattice in the orthogonal space. Our construction is obtained by combining results of Kisin and the first author with an explicit presentation and resolution of a corresponding local model.


title: A Simulation Study of Electric Field Engineering with Multi-Level Pinned Photodiodes for Fast and Complete Charge Transfer

abstract: In a CMOS image sensor pixel, fast and complete charge transfer from pinned photodiode (PPD) is desired and necessary in some applications. In special cases such as time-of-flight imaging or large pinned photodiodes, the PPD potential well shape highly affects the charge transfer performance and should be engineered carefully. In the present work, a PPD structure named multi-level PPD is introduced and examined through simulation study. Moreover, a fast and effective way to analyze the pinning process for a lag-free design is introduced. It is concluded that the proposed PPD achieves fast and complete charge transfer without additional implementation masks or process steps. The proposed PPD is compared with a similar conventional rectangular pixel and 31% reduction in the charge transfer time is observed.


title: Local well-posedness and parabolic smoothing of solutions of fully nonlinear third-order equations on the torus

abstract: We study the initial value problem of fully nonlinear third-order equations on the torus. Under some conditions on the nonlinearity and the data, we prove that the equation behaves like a parabolic one: there exists a unique local solution in one direction of time that is infinitely smooth and the problem in not well-posed in the other direction. Under other conditions on the nonlinearity and the data, we prove that the equation behaves like a dispersive one: there exists a unique local solution in both directions of time.


title: TricycleGAN: Unsupervised Image Synthesis and Segmentation Based on Shape Priors

abstract: Medical image segmentation is routinely performed to isolate regions of interest, such as organs and lesions. Currently, deep learning is the state of the art for automatic segmentation, but is usually limited by the need for supervised training with large datasets that have been manually segmented by trained clinicians. The goal of semi-superised and unsupervised image segmentation is to greatly reduce, or even eliminate, the need for training data and therefore to minimze the burden on clinicians when training segmentation models. To this end we introduce a novel network architecture for capable of unsupervised and semi-supervised image segmentation called TricycleGAN. This approach uses three generative models to learn translations between medical images and segmentation maps using edge maps as an intermediate step. Distinct from other approaches based on generative networks, TricycleGAN relies on shape priors rather than colour and texture priors. As such, it is particularly well-suited for several domains of medical imaging, such as ultrasound imaging, where commonly used visual cues may be absent. We present experiments with TricycleGAN on a clinical dataset of kidney ultrasound images and the benchmark ISIC 2018 skin lesion dataset.


title: HMC, an Algorithms in Data Mining, the Functional Analysis approach

abstract: The main purpose of this paper is to facilitate the communication between the Analytic, Probabilistic and Algorithmic communities. We present a proof of convergence of the Hamiltonian (Hybrid) Monte Carlo algorithm from the point of view of the Dynamical Systems, where the evolving objects are densities of probability distributions and the tool are derived from the Functional Analysis.


title: Mobility spaces and geodesics for the n-sphere

abstract: We introduce an algebraic system which can be used as a model for spaces with geodesic paths between any two of their points. This new algebraic structure is based on the notion of mobility algebra which has recently been introduced as a model for the unit interval of real numbers. We show that there is a strong connection between modules over a ring and affine mobility spaces over a mobility algebra. However, geodesics in general fail to be affine thus giving rise to the new algebraic structure of mobility space. We show that the so called formula for spherical linear interpolation, which gives geodesics on the n-sphere, is an example of a mobility space over the unit interval mobility algebra.


title: Skorohod and Stratonovich integrals for controlled processes

abstract: Given a continuous Gaussian process $x$ which gives rise to a $p$-geometric rough path for $p\in (2,3)$, and a general continuous process $y$ controlled by $x$, under proper conditions we establish the relationship between the Skorohod integral $\int_0^t y_s {\mathrm{d}}^\diamond x_s$ and the Stratonovich integral $\int_0^t y_s {\mathrm{d}} {\mathbf x}_s$. Our strategy is to employ the tools from rough paths theory and Malliavin calculus to analyze discrete sums of the integrals.


title: Invertible DenseNets with Concatenated LipSwish

abstract: We introduce Invertible Dense Networks (i-DenseNets), a more parameter efficient extension of Residual Flows. The method relies on an analysis of the Lipschitz continuity of the concatenation in DenseNets, where we enforce invertibility of the network by satisfying the Lipschitz constant. Furthermore, we propose a learnable weighted concatenation, which not only improves the model performance but also indicates the importance of the concatenated weighted representation. Additionally, we introduce the Concatenated LipSwish as activation function, for which we show how to enforce the Lipschitz condition and which boosts performance. The new architecture, i-DenseNet, out-performs Residual Flow and other flow-based models on density estimation evaluated in bits per dimension, where we utilize an equal parameter budget. Moreover, we show that the proposed model out-performs Residual Flows when trained as a hybrid model where the model is both a generative and a discriminative model.


title: The Infrared Emission and Vigorous Star Formation of Low-redshift Quasars

abstract: The star formation activity of the host galaxies of active galactic nuclei (AGNs) provides valuable insights into the complex interconnections between black hole growth and galaxy evolution. A major obstacle arises from the difficulty of estimating accurate star formation rates in the presence of a strong AGN. Analyzing the $1-500\, \mu m$ spectral energy distributions and high-resolution mid-infrared spectra of low-redshift ($z < 0.5$) Palomar-Green quasars with bolometric luminosity $\sim 10^{44.5}-10^{47.5}\rm\,erg\,s^{-1}$, we find, from comparison with an independent star formation rate indicator based on [Ne II] 12.81$\, \mu m$ and [Ne III] 15.56$\, \mu m$, that the torus-subtracted, total infrared ($8-1000\, \mu m$) emission yields robust star formation rates in the range $\sim 1-250\,M_\odot\,{\rm yr^{-1}}$. Combined with available stellar mass estimates, the vast majority ($\sim 75\%-90\%$) of the quasars lie on or above the main sequence of local star-forming galaxies, including a significant fraction ($\sim 50\%-70\%$) that would qualify as starburst systems. This is further supported by the high star formation efficiencies derived from the gas content inferred from the dust masses. Inspection of high-resolution Hubble Space Telescope images reveals a wide diversity of morphological types, including a number of starbursting hosts that have not experienced significant recent dynamical perturbations. The origin of the high star formation efficiency is unknown.


title: Active Boundary Loss for Semantic Segmentation

abstract: This paper proposes a novel active boundary loss for semantic segmentation. It can progressively encourage the alignment between predicted boundaries and ground-truth boundaries during end-to-end training, which is not explicitly enforced in commonly used cross-entropy loss. Based on the predicted boundaries detected from the segmentation results using current network parameters, we formulate the boundary alignment problem as a differentiable direction vector prediction problem to guide the movement of predicted boundaries in each iteration. Our loss is model-agnostic and can be plugged in to the training of segmentation networks to improve the boundary details. Experimental results show that training with the active boundary loss can effectively improve the boundary F-score and mean Intersection-over-Union on challenging image and video object segmentation datasets.


title: Covid-19 risk factors: Statistical learning from German healthcare claims data

abstract: We analyse prior risk factors for severe, critical or fatal courses of Covid-19 based on a retrospective cohort using claims data of the AOK Bayern. As our main methodological contribution, we avoid prior grouping and pre-selection of candidate risk factors. Instead, fine-grained hierarchical information from medical classification systems for diagnoses, pharmaceuticals and procedures are used, resulting in more than 33,000 covariates. Our approach has better predictive ability than well-specified morbidity groups but does not need prior subject-matter knowledge. The methodology and estimated coefficients are made available to decision makers to prioritize protective measures towards vulnerable subpopulations and to researchers who like to adjust for a large set of confounders in studies of individual risk factors.


title: On the Eptihermal Neutron Energy Limit for Accelerator-Based Boron Neutron Capture Therapy (AB-BNCT): Study and Impact of New Energy Limits

abstract: Background and purpose: Accelerator-Based Boron Neutron Capture Therapy is a radiotherapy based on compact accelerator neutron sources requiring an epithermal neutron field for tumour irradiations. Neutrons of 10 keV are considered as the maximum optimised energy to treat deep-seated tumours. We investigated, by means of Monte Carlo simulations, the epithermal range from 10 eV to 10 keV in order to optimise the maximum epithermal neutron energy as a function of the tumour depth. Methods: A Snyder head phantom was simulated and mono-energetic neutrons with 4 different incident energies were used: 10 eV, 100 eV, 1 keV and 10 keV. $^{10}$B capture rates and absorbed dose composition on every tissue were calculated to describe and compare the effects of lowering the maximum epithermal energy. The Therapeutic Gain (TG) was estimated considering the whole brain volume. Results: For tumours seated at 4 cm depth, 10 eV, 100 eV and 1 keV neutrons provided respectively 54 $\%$, 36 $\%$ and 18 $\%$ increase on the TG compared to 10 keV neutrons. Neutrons with energies between 10 eV and 1 keV provided higher TG than 10 keV neutrons for tumours seated up to 6.4 cm depth inside the head. The size of the tumour does not change these results. Conclusions: Using lower epithermal energy neutrons for AB-BNCT tumour irradiation could improve treatment efficacy, delivering more therapeutic dose while reducing the dose in healthy tissues. This could lead to new Beam Shape Assembly designs in order to optimise the BNCT irradiation.


title: Geometric quantization via cotangent models

abstract: In this article we give a universal model for geometric quantization associated to a real polarization given by an integrable system with non-degenerate singularities. This universal model goes one step further than the previous cotangent models by both considering singular orbits and adding to the cotangent models a model for the prequantum line bundle. These singularities are generic in the sense that are given by Morse-type functions and include elliptic, hyperbolic and focus-focus singularities. Examples of systems admitting such singularities are toric, semitoric and almost toric manifolds, as well as physical systems such as the coupling of harmonic oscillators, the spherical pendulum or the reduction of the Euler's equations of the rigid body on $T^*(SO(3))$ to a sphere. Our geometric quantization formulation coincides with the previous models away from the singularities and corrects former models for hyperbolic and focus-focus singularities cancelling out the infinite dimensional contributions obtained by former approaches. The geometric quantization models provided here match the classical physical methods for mechanical systems such as the spherical pendulum. Our cotangent models obey a local-to-global principle and can be glued to determine the geometric quantization of the global systems even if the global symplectic classification of the systems is not known in general.


title: Additive Average Schwarz Method for Elliptic Mortar Finite Element Problems with Highly Heterogeneous Coefficients

abstract: In this paper, we extend the additive average Schwarz method to solve second order elliptic boundary value problems with heterogeneous coefficients inside the subdomains and across their interfaces by the mortar technique, where the mortar finite element discretization is on nonmatching meshes. In this two-level method, we enrich the coarse space in two different ways, i.e., by adding eigenfunctions of two variants of the generalized eigenvalue problems. We prove that the condition numbers of the systems of algebraic equations resulting from the extended additive average Schwarz method, corresponding to both coarse spaces, are of the order O(H/h) and independent of jumps in the coefficients, where H and h are the mesh parameters.


title: Conformal Methods in General Relativity with application to Conformal Cyclic Cosmology: A minicourse given at the IXth IMLG Warsaw 2018

abstract: In these lectures my aim is to review enough of conformal differential geometry in four dimensions to give an account of Penrose's conformal cyclic geometry.


title: Biases in Moments of Dirichlet Coefficients of Elliptic Curve Families

abstract: Elliptic curves arise in many important areas of modern number theory. One way to study them is take local data, the number of solutions modulo $p$, and create an $L$-function. The behavior of this global object is related to two of the seven Clay Millennial Problems: the Birch and Swinnerton-Dyer Conjecture and the Generalized Riemann Hypothesis. We study one-parameter families over $\mathbb{Q}(T)$, which are of the form $y^2=x^3+A(T)x+B(T)$, with non-constant $j$-invariant. We define the r-th moment of an elliptic curve to be $A_{r,E}(p) := \frac1{p} \sum_{t \bmod p} a_t(p)^r$, where $a_t(p)$ is $p$ minus the number of solutions to $y^2 = x^3 + A(t)x + B(t) \bmod p$. Rosen and Silverman showed biases in the first moment equal the rank of the Mordell-Weil group of rational solutions. Michel proved that $pA_{2,E}(p)=p^2+O(p^{3/2})$. Based on several special families where computations can be done in closed form, Miller in his thesis conjectured that the largest lower-order term in the second moment that does not average to $0$ is on average negative. He further showed that such a negative bias has implications in the distribution of zeros of the elliptic curve $L$-function near the central point. To date, evidence for this conjecture is limited to special families. In addition to studying some additional families where the calculations can be done in closed form, we also systematically investigate families of various rank. These are the first general tests of the conjecture; while we cannot in general obtain closed form solutions, we discuss computations which support or contradict the conjecture. We then generalize to higher moments, and see evidence that the bias continues in the even moments.


title: Separable Joint Blind Deconvolution and Demixing

abstract: Blind deconvolution and demixing is the problem of reconstructing convolved signals and kernels from the sum of their convolutions. This problem arises in many applications, such as blind MIMO. This work presents a separable approach to blind deconvolution and demixing via convex optimization. Unlike previous works, our formulation allows separation into smaller optimization problems, which significantly improves complexity. We develop recovery guarantees, which comply with those of the original non-separable problem, and demonstrate the method performance under several normalization constraints.


title: Continuity of Monge-Amp{\`e}re Potentials in Big Cohomology Classes

abstract: Extending DiNezza-Lu's approach to the setting of big cohomology classes, we prove that solutions of degenerate complex Monge-Amp{\`e}re equations on compact K{\"a}hler manifolds are continuous on a Zariski open set. This allows us to show that singular K{\"a}hler-Einstein metrics on log canonical varieties of general type have continuous potentials on the ample locus outside of the non-klt part.


title: EFloat: Entropy-coded Floating Point Format for Compressing Vector Embedding Models

abstract: In a large class of deep learning models, including vector embedding models such as word and database embeddings, we observe that floating point exponent values cluster around a few unique values, permitting entropy based data compression. Entropy coding compresses fixed-length values with variable-length codes, encoding most probable values with fewer bits. We propose the EFloat compressed floating point number format that uses a variable field boundary between the exponent and significand fields. EFloat uses entropy coding on exponent values and signs to minimize the average width of the exponent and sign fields, while preserving the original FP32 exponent range unchanged. Saved bits become part of the significand field increasing the EFloat numeric precision by 4.3 bits on average compared to other reduced-precision floating point formats. EFloat makes 8-bit and even smaller floats practical without sacrificing the exponent range of a 32-bit floating point representation. We currently use the EFloat format for saving memory capacity and bandwidth consumption of large vector embedding models such as those used for database embeddings. Using the RMS error as metric, we demonstrate that EFloat provides higher accuracy than other floating point formats with equal bit budget. The EF12 format with 12-bit budget has less end-to-end application error than the 16-bit BFloat16. EF16 with 16-bit budget has an RMS-error 17 to 35 times less than BF16 RMS-error for a diverse set of embedding models. When making similarity and dissimilarity queries, using the NDCG ranking metric, EFloat matches the result quality of prior floating point representations with larger bit budgets.


title: ProxyFAUG: Proximity-based Fingerprint Augmentation

abstract: The proliferation of data-demanding machine learning methods has brought to light the necessity for methodologies which can enlarge the size of training datasets, with simple, rule-based methods. In-line with this concept, the fingerprint augmentation scheme proposed in this work aims to augment fingerprint datasets which are used to train positioning models. The proposed method utilizes fingerprints which are recorded in spacial proximity, in order to perform fingerprint augmentation, creating new fingerprints which combine the features of the original ones. The proposed method of composing the new, augmented fingerprints is inspired by the crossover and mutation operators of genetic algorithms. The ProxyFAUG method aims to improve the achievable positioning accuracy of fingerprint datasets, by introducing a rule-based, stochastic, proximity-based method of fingerprint augmentation. The performance of ProxyFAUG is evaluated in an outdoor Sigfox setting using a public dataset. The best performing published positioning method on this dataset is improved by 40% in terms of median error and 6% in terms of mean error, with the use of the augmented dataset. The analysis of the results indicate a systematic and significant performance improvement at the lower error quartiles, as indicated by the impressive improvement of the median error.


title: Diving inside a hairy black hole

abstract: We investigate the interior of the Einstein-Gauss-Bonnet charged black-hole with scalar hair. We find a variety of dynamical epochs, with the particular important feature that the Cauchy horizon is not present. This makes the violation of the no-hair theorem a possible tool to understand how might the strong cosmic censorship conjecture work.


title: Fractionally Log-Concave and Sector-Stable Polynomials: Counting Planar Matchings and More

abstract: We show fully polynomial time randomized approximation schemes (FPRAS) for counting matchings of a given size, or more generally sampling/counting monomer-dimer systems in planar, not-necessarily-bipartite, graphs. While perfect matchings on planar graphs can be counted exactly in polynomial time, counting non-perfect matchings was shown by [Jer87] to be #P-hard, who also raised the question of whether efficient approximate counting is possible. We answer this affirmatively by showing that the multi-site Glauber dynamics on the set of monomers in a monomer-dimer system always mixes rapidly, and that this dynamics can be implemented efficiently on downward-closed families of graphs where counting perfect matchings is tractable. As further applications of our results, we show how to sample efficiently using multi-site Glauber dynamics from partition-constrained strongly Rayleigh distributions, and nonsymmetric determinantal point processes. In order to analyze mixing properties of the multi-site Glauber dynamics, we establish two notions for generating polynomials of discrete set-valued distributions: sector-stability and fractional log-concavity. These notions generalize well-studied properties like real-stability and log-concavity, but unlike them robustly degrade under useful transformations applied to the distribution. We relate these notions to pairwise correlations in the underlying distribution and the notion of spectral independence introduced by [ALO20], providing a new tool for establishing spectral independence based on geometry of polynomials. As a byproduct of our techniques, we show that polynomials avoiding roots in a sector of the complex plane must satisfy what we call fractional log-concavity; this extends a classic result established by [Gar59] who showed homogeneous polynomials that have no roots in a half-plane must be log-concave over the positive orthant.


title: Semi-device-independent certification of entanglement in superdense coding

abstract: Superdense coding is a paradigmatic protocol in quantum information science, employing a quantum communication channel to send classical information more efficiently. As we show here, it can be understood as a particular case of a prepare and measure experiment, a scenario that has attracted growing attention for its fundamental and practical applications. Formulating superdense coding as a prepare and measure scenario allows us to provide a semi-device-independent witness of entanglement that significantly improves over previous tests. Furthermore, we also show how to adapt our results into self-testing of maximally entangled states and also provide a semidefinite program formulation allowing to efficiently optimize, for any shared quantum state, the probability of success in the superdense coding protocol.


title: Matching Impatient and Heterogeneous Demand and Supply

abstract: Service platforms must determine rules for matching heterogeneous demand (customers) and supply (workers) that arrive randomly over time and may be lost if forced to wait too long for a match. Our objective is to maximize the cumulative value of matches, minus costs incurred when demand and supply wait. We develop a fluid model, that approximates the evolution of the stochastic model, and captures explicitly the nonlinear dependence between the amount of demand and supply waiting and the distribution of their patience times, also known as reneging or abandonment times in the literature. The fluid model invariant states approximate the steady-state mean queue-lengths in the stochastic system, and, therefore, can be used to develop an optimization problem whose optimal solution provides matching rates between demand and supply types that are asymptotically optimal (on fluid scale, as demand and supply rates grow large). We propose a discrete review matching policy that asymptotically achieves the optimal matching rates. We further show that when the aforementioned matching optimization problem has an optimal extreme point solution, which occurs when the patience time distributions have increasing hazard rate functions, a state-independent priority policy, that ranks the edges on the bipartite graph connecting demand and supply, is asymptotically optimal. A key insight from this analysis is that the ranking critically depends on the patience time distributions, and may be different for different distributions even if they have the same mean, demonstrating that models assuming, e.g., exponential patience times for tractability, may lack robustness. Finally, we observe that when holding costs are zero, a discrete review policy, that does not require knowledge of inter-arrival and patience time distributions, is asymptotically optimal.


title: Fine-tuning deep learning model parameters for improved super-resolution of dynamic MRI with prior-knowledge

abstract: Dynamic imaging is a beneficial tool for interventions to assess physiological changes. Nonetheless during dynamic MRI, while achieving a high temporal resolution, the spatial resolution is compromised. To overcome this spatio-temporal trade-off, this research presents a super-resolution (SR) MRI reconstruction with prior knowledge based fine-tuning to maximise spatial information while reducing the required scan-time for dynamic MRIs. An U-Net based network with perceptual loss is trained on a benchmark dataset and fine-tuned using one subject-specific static high resolution MRI as prior knowledge to obtain high resolution dynamic images during the inference stage. 3D dynamic data for three subjects were acquired with different parameters to test the generalisation capabilities of the network. The method was tested for different levels of in-plane undersampling for dynamic MRI. The reconstructed dynamic SR results after fine-tuning showed higher similarity with the high resolution ground-truth, while quantitatively achieving statistically significant improvement. The average SSIM of the lowest resolution experimented during this research (6.25~\% of the k-space) before and after fine-tuning were 0.939 $\pm$ 0.008 and 0.957 $\pm$ 0.006 respectively. This could theoretically result in an acceleration factor of 16, which can potentially be acquired in less than half a second. The proposed approach shows that the super-resolution MRI reconstruction with prior-information can alleviate the spatio-temporal trade-off in dynamic MRI, even for high acceleration factors.


title: Tachyonic Preheating in Palatini $R^2$ Inflation

abstract: We study preheating in the Palatini formalism with a quadratic inflaton potential and an added $\alpha R^2$ term. In such models, the oscillating inflaton field repeatedly returns to the plateau of the Einstein frame potential, on which the tachyonic instability fragments the inflaton condensate within less than an e-fold. We find that tachyonic preheating takes place when $\alpha \gtrsim 10^{13}$ and that the energy density of the fragmented field grows with the rate $\Gamma/H \approx 0.011 \times \alpha^{0.31}$. The model extends the family of plateau models with similar preheating behaviour. Although it contains non-canonical quartic kinetic terms in the Einstein frame, we show that, in the first approximation, these can be neglected during both preheating and inflation.


title: Adiabatic waveforms for extreme mass-ratio inspirals via multivoice decomposition in time and frequency

abstract: We compute adiabatic waveforms for extreme mass-ratio inspirals (EMRIs) by "stitching" together a long inspiral waveform from a sequence of waveform snapshots, each of which corresponds to a particular geodesic orbit. We show that the complicated total waveform can be regarded as a sum of "voices." Each voice evolves in a simple way on long timescales, a property which can be exploited to efficiently produce waveform models that faithfully encode the properties of EMRI systems. We look at examples for a range of different orbital geometries: spherical orbits, equatorial eccentric orbits, and one example of generic (inclined and eccentric) orbits. To our knowledge, this is the first calculation of a generic EMRI waveform that uses strong-field radiation reaction. We examine waveforms in both the time and frequency domains. Although EMRIs evolve slowly enough that the stationary phase approximation (SPA) to the Fourier transform is valid, the SPA calculation must be done to higher order for some voices, since their instantaneous frequency can change from chirping forward ($\dot f > 0$) to chirping backward ($\dot f < 0$). The approach we develop can eventually be extended to more complete EMRI waveform models, for example to include effects neglected by the adiabatic approximation such as the conservative self force and spin-curvature coupling.


title: Magnon-spinon dichotomy in the Kitaev hyperhoneycomb $\beta$-Li$_2$IrO$_3$

abstract: The family of edge-sharing tri-coordinated iridates and ruthenates has emerged in recent years as a major platform for Kitaev spin liquid physics, where spins fractionalize into emergent magnetic fluxes and Majorana fermions with Dirac-like dispersions. While such exotic states are usually pre-empted by long-range magnetic order at low temperatures, signatures of Majorana fermions with long coherent times have been predicted to manifest at intermediate and higher energy scales, similar to the observation of spinons in quasi-1D spin chains. Here we present a Resonant Inelastic X-ray Scattering study of the magnetic excitations of the hyperhoneycomb iridate $\beta$-Li$_2$IrO$_3$ under a magnetic field with a record-high-resolution spectrometer. At low-temperatures, dispersing spin waves can be resolved around the predicted intertwined incommensurate spiral and field-induced zigzag orders, whose excitation energy reaches a maximum of 16meV. A 2T magnetic field softens the dispersion around ${\bf Q}=0$. The behavior of the spin waves under magnetic field is consistent with our semiclassical calculations for the ground state and the dynamical spin structure factor, which further predicts that the ensued intertwined uniform states remain robust up to very high fields (100 T). Most saliently, the low-energy magnon-like mode is superimposed by a broad continuum of excitations, centered around 35meV and extending up to 100meV. This high-energy continuum survives up to at least 300K -- well above the ordering temperature of 38K -- and gives evidence for pairs of long-lived Majorana fermions of the proximate Kitaev spin liquid.


title: Coherent ray-wave structured light based on (helical) Ince-Gaussian modes

abstract: The topological evolution of classic eigenmodes including Hermite-Laguerre-Gaussian and (helical) InceGaussian modes is exploited to construct coherent state modes, which unifies the representations of travelingwave (TW) and standing-wave (SW) ray-wave structured light for the first time and realizes the TW-SW unified ray-wave geometric beam with topology of raytrajectories splitting effect, breaking the boundary of TW and SW structured light. We experimentally generate these new modes with high purity and dynamic control by digital holography method, revealing potential applications in optical manipulation and communication.


title: Resonance suppression of the r-mode instability in superfluid neutron stars: Accounting for muons and entrainment

abstract: We calculate the finite-temperature r-mode spectrum of a superfluid neutron star accounting for both muons in the core and the entrainment between neutrons and protons. We show that the standard perturbation scheme, considering the rotation rate as an expansion parameter, breaks down in this case. We develop an original perturbation scheme which circumvents this problem by treating both the perturbations due to rotation and (weak) entrainment simultaneously. Applying this scheme, we propose a simple method for calculating the superfluid r-mode eigenfrequency in the limit of vanishing rotation rate. We also calculate the r-mode spectrum at finite rotation rate for realistic microphysics input (adopting, however, the Newtonian framework and Cowling approximation when considering perturbed oscillation equations) and show that the normal r-mode exhibits resonances with superfluid r-modes at certain values of temperatures and rotation frequencies in the parameter range relevant to neutron stars in low-mass X-ray binaries (LMXBs). This turns the recently suggested phenomenological model of resonance r-mode stabilization into a quantitative theory, capable of explaining observations. A strong dependence of resonance rotation rates and temperatures on the neutron superfluidity model allows us to constrain the latter by confronting our calculations with the observations of neutron stars in LMXBs.


title: RoI Tanh-polar Transformer Network for Face Parsing in the Wild

abstract: Face parsing aims to predict pixel-wise labels for facial components of a target face in an image. Existing approaches usually crop the target face from the input image with respect to a bounding box calculated during pre-processing, and thus can only parse inner facial Regions of Interest~(RoIs). Peripheral regions like hair are ignored and nearby faces that are partially included in the bounding box can cause distractions. Moreover, these methods are only trained and evaluated on near-frontal portrait images and thus their performance for in-the-wild cases has been unexplored. To address these issues, this paper makes three contributions. First, we introduce iBugMask dataset for face parsing in the wild, which consists of 21,866 training images and 1,000 testing images. The training images are obtained by augmenting an existing dataset with large face poses. The testing images are manually annotated with $11$ facial regions and there are large variations in sizes, poses, expressions and background. Second, we propose RoI Tanh-polar transform that warps the whole image to a Tanh-polar representation with a fixed ratio between the face area and the context, guided by the target bounding box. The new representation contains all information in the original image, and allows for rotation equivariance in the convolutional neural networks~(CNNs). Third, we propose a hybrid residual representation learning block, coined HybridBlock, that contains convolutional layers in both the Tanh-polar space and the Tanh-Cartesian space, allowing for receptive fields of different shapes in CNNs. Through extensive experiments, we show that the proposed method improves the state-of-the-art for face parsing in the wild and does not require facial landmarks for alignment.


title: On the stability of the martingale optimal transport problem: A set-valued map approach

abstract: Continuity of the value of the martingale optimal transport problem on the real line w.r.t. its marginals was recently established in Backhoff-Veraguas and Pammer [2] and Wiesel [21]. We present a new perspective of this result using the theory of set-valued maps. In particular, using results from Beiglb\"ock, Jourdain, Margheriti, and Pammer [5], we show that the set of martingale measures with fixed marginals is continuous, i.e., lower- and upper hemicontinuous, w.r.t. its marginals. Moreover, we establish compactness of the set of optimizers as well as upper hemicontinuity of the optimizers w.r.t. the marginals.


title: Dynamical phases and quantum correlations in an emitter-waveguide system with feedback

abstract: We investigate the creation and control of emergent collective behavior and quantum correlations using feedback in an emitter-waveguide system using a minimal model. Employing homodyne detection of photons emitted from a laser-driven emitter ensemble into the modes of a waveguide allows to generate intricate dynamical phases. In particular, we show the emergence of a time-crystal phase, the transition to which is controlled by the feedback strength. Feedback enables furthermore the control of many-body quantum correlations, which become manifest in spin squeezing in the emitter ensemble. Developing a theory for the dynamics of fluctuation operators we discuss how the feedback strength controls the squeezing and investigate its temporal dynamics and dependence on system size. The largely analytical results allow to quantify spin squeezing and fluctuations in the limit of large number of emitters, revealing critical scaling of the squeezing close to the transition to the time-crystal. Our study corroborates the potential of integrated emitter-waveguide systems -- which feature highly controllable photon emission channels -- for the exploration of collective quantum phenomena and the generation of resources, such as squeezed states, for quantum enhanced metrology.


title: Edge Computing: A Systematic Mapping Study

abstract: Edge computing is a novel computing paradigm that extends cloud resources at the edge of the network to tackle the problem of communication latency in latency-sensitive applications. For the last decades, there have been many efforts dedicated to this field. The dramatic growth in the publications, and the great attention of the research community, have made it necessary to conduct a Systematic Mapping Study (SMS) to enable researchers get a better view of the field. A three-tier search method is considered in this work. In this method, we defined some quality criteria to extract appropriate search spaces and studies. Using this methodology, we select 112 search spaces out of 805 ones, and by searching in these search spaces we select 1440 studies out of 8725. In our SMS, 8 research questions have been designed and answered to identify the main topics, architectures, techniques, etc. in the field of edge computing.


title: An electromagnetic physics constructor for low energy polarised X-/gamma ray transport in Geant4

abstract: The production, application, and/or measurement of polarised X-/gamma rays are key to the fields of synchrotron science and X-/gamma-ray astronomy. The design, development and optimisation of experimental equipment utilised in these fields typically relies on the use of Monte Carlo radiation transport modelling toolkits such as Geant4. In this work the Geant4 "G4LowEPPhysics" electromagnetic physics constructor has been reconfigured to offer a "best set" of electromagnetic physics models for studies exploring the transport of low energy polarised X-/gamma rays. An overview of the physics models implemented in "G4LowEPPhysics", and it's experimental validation against Compton X-ray polarimetry measurements of the BL38B1 beamline at the SPring-8 synchrotron (Sayo, Japan) is reported. "G4LowEPPhysics" is shown to be able to reproduce the experimental results obtained at the BL38B1 beamline (SPring-8) to within a level of accuracy on the same order as Geant4's X-/gamma ray interaction cross-sectional data uncertainty (approximately $\pm$ 5 \%).


title: $\mathcal{N}=1$ supersymmetric three-dimensional QED in the large-$N_f$ limit and applications to super-graphene

abstract: We study $\mathcal{N}=1$ supersymmetric three-dimensional Quantum Electrodynamics with $N_f$ two-component fermions. Due to the infra-red (IR) softening of the photon, $\ep$-scalar and photino propagators, the theory flows to an interacting fixed point deep in the IR, $p_E \ll e^2 N_f/8$, where $p_E$ is the euclidean momentum and $e$ the electric charge. At next-to-leading order in the $1/N_f$-expansion, we find that the flow of the dimensionless effective coupling constant $\overline{\al}$ is such that: $\overline{\al} \ra 8/\big(N_f \,(1+C/N_f)\big) \approx (8/N_f)(1-0.4317/N_f)$ where $C= 2\,(12-\pi^2)/\pi^2$. Hence, the non-trivial IR fixed point is stable with respect to quantum corrections. Various properties of the theory are explored and related via a mapping to the ones of a $\mathcal{N}=1$ model of super-graphene. In particular, we derive the interaction correction coefficient to the optical conductivity of super-graphene, $C_{\rm sg} = (12-\pi^2)/(2\pi) = 0.3391$, which is six times larger than in the non-supersymmetric case, $C_{\rm g} = (92-9\pi^2)/(18\pi) = 0.0561$.


title: Data-to-text Generation with Macro Planning

abstract: Recent approaches to data-to-text generation have adopted the very successful encoder-decoder architecture or variants thereof. These models generate text which is fluent (but often imprecise) and perform quite poorly at selecting appropriate content and ordering it coherently. To overcome some of these issues, we propose a neural model with a macro planning stage followed by a generation stage reminiscent of traditional methods which embrace separate modules for planning and surface realization. Macro plans represent high level organization of important content such as entities, events and their interactions; they are learnt from data and given as input to the generator. Extensive experiments on two data-to-text benchmarks (RotoWire and MLB) show that our approach outperforms competitive baselines in terms of automatic and human evaluation.


title: Central extensions of some linear cycle sets

abstract: For each member $\mathcal{A}$ of a family of linear cycle sets whose underlying abelian group is cyclic of order a power of a prime number, we compute all the central extensions of $\mathcal{A}$ by an arbitrary abelian group.


title: Water reactions on reconstructed rutile TiO$_2$: a DFT / DFTB approach

abstract: Far from being conclusively understood, the reactive interaction of water with rutile does still present a challenge to atomistic modelling techniques rooted on quantum mechanics. We show that static geometries of stoichiometric TiO$_2$/water interfaces can be described well by Density Functional Tight Binding (DFTB). However, this method needs further improvements to reproduce the low dissociation propensity of H$_2$O after adsorption predicted by Density Functional Theory (DFT). A reliable description of the surface reactivity of water is fundamental to investigate the non-stoichiometric reconstruction of the (001) facet rich in Ti interstitials. Calculations based on (DFT) predict the transition temperature for the onset of reconstruction in remarkable agreement with experiments and suggest that this surface, in contact with liquid water, can promote spontaneous H$_2$O splitting and formation of H$_2$ molecules.


title: Normalizing flows for microscopic many-body calculations: an application to the nuclear equation of state

abstract: Normalizing flows are a class of machine learning models used to construct a complex distribution through a bijective mapping of a simple base distribution. We demonstrate that normalizing flows are particularly well suited as a Monte Carlo integration framework for quantum many-body calculations that require the repeated evaluation of high-dimensional integrals across smoothly varying integrands and integration regions. As an example, we consider the finite-temperature nuclear equation of state. An important advantage of normalizing flows is the ability to build highly expressive models of the target integrand, which we demonstrate enables precise evaluations of the nuclear free energy and its derivatives. Furthermore, we show that a normalizing flow model trained on one target integrand can be used to efficiently calculate related integrals when the temperature, density, or nuclear force is varied. This work will support future efforts to build microscopic equations of state for numerical simulations of supernovae and neutron star mergers that employ state-of-the-art nuclear forces and many-body methods.


title: Multiple Criss-Cross Insertion and Deletion Correcting Codes

abstract: This paper investigates the problem of correcting multiple criss-cross insertions and deletions in arrays. More precisely, we study the unique recovery of $n \times n$ arrays affected by $t$-criss-cross deletions defined as any combination of $t_r$ row and $t_c$ column deletions such that $t_r + t_c = t$ for a given $t$. We show an equivalence between correcting $t$-criss-cross deletions and $t$-criss-cross insertions and show that a code correcting $t$-criss-cross insertions/deletions has redundancy at least $tn + t \log n - \log(t!)$. Then, we present an existential construction of $t$-criss-cross insertion/deletion correcting code with redundancy bounded from above by $tn + \mathcal{O}(t^2 \log^2 n)$. The main ingredients of the presented code construction are systematic binary $t$-deletion correcting codes and Gabidulin codes. The first ingredient helps locating the indices of the inserted/deleted rows and columns, thus transforming the insertion/deletion-correction problem into a row/column erasure-correction problem which is then solved using the second ingredient.


title: Minimum-Complexity Failure Correction in Linear Arrays via Compressive Processing

abstract: Given an array with defective elements, failure correction (FC) aims at finding a new set of weights for the working elements so that the properties of the original pattern can be recovered. Unlike several FC techniques available in the literature, which update all the working excitations, the Minimum-Complexity Failure Correction (MCFC) problem is addressed in this paper. By properly reformulating the FC problem, the minimum number of corrections of the whole excitations of the array is determined by means of an innovative Compressive Processing (CP) technique in order to afford a pattern as close as possible to the original one (i.e., the array without failures). Selected examples, from a wide set of numerical test cases, are discussed to assess the effectiveness of the proposed approach as well as to compare its performance with other competitive state-of-the-art techniques in terms of both pattern features and number of corrections.


title: Adversarial Attacks and Defenses in Physiological Computing: A Systematic Review

abstract: Physiological computing uses human physiological data as system inputs in real time. It includes, or significantly overlaps with, brain-computer interfaces, affective computing, adaptive automation, health informatics, and physiological signal based biometrics. Physiological computing increases the communication bandwidth from the user to the computer, but is also subject to various types of adversarial attacks, in which the attacker deliberately manipulates the training and/or test examples to hijack the machine learning algorithm output, leading to possible user confusion, frustration, injury, or even death. However, the vulnerability of physiological computing systems has not been paid enough attention to, and there does not exist a comprehensive review on adversarial attacks to them. This paper fills this gap, by providing a systematic review on the main research areas of physiological computing, different types of adversarial attacks and their applications to physiological computing, and the corresponding defense strategies. We hope this review will attract more research interests on the vulnerability of physiological computing systems, and more importantly, defense strategies to make them more secure.


title: Feedback Capacity of Parallel ACGN Channels and Kalman Filter: Power Allocation with Feedback

abstract: In this paper, we relate the feedback capacity of parallel additive colored Gaussian noise (ACGN) channels to a variant of the Kalman filter. By doing so, we obtain lower bounds on the feedback capacity of such channels, as well as the corresponding feedback (recursive) coding schemes, which are essentially power allocation policies with feedback, to achieve the bounds. The results are seen to reduce to existing lower bounds in the case of a single ACGN feedback channel, whereas when it comes to parallel additive white Gaussian noise (AWGN) channels with feedback, the recursive coding scheme reduces to a feedback "water-filling" power allocation policy.

categories: cs.IT cs.LG cs.SY eess.SP eess.SY math.IT math.OC


title: Identifying regions of inhomogeneities in spatial processes via an M-RA and mixture priors

abstract: Soils have been heralded as a hidden resource that can be leveraged to mitigate and address some of the major global environmental challenges. Specifically, the organic carbon stored in soils, called Soil Organic Carbon (SOC), can, through proper soil management, help offset fuel emissions, increase food productivity, and improve water quality. As collecting data on SOC is costly and time consuming, not much data on SOC is available, although understanding the spatial variability in SOC is of fundamental importance for effective soil management. In this manuscript, we propose a modeling framework that can be used to gain a better understanding of the dependence structure of a spatial process by identifying regions within a spatial domain where the process displays the same spatial correlation range. To achieve this goal, we propose a generalization of the Multi-Resolution Approximation (M-RA) modeling framework of Katzfuss (2017) originally introduced as a strategy to reduce the computational burden encountered when analyzing massive spatial datasets. To allow for the possibility that the correlation of a spatial process might be characterized by a different range in different subregions of a spatial domain, we provide the M-RA basis functions weights with a two-component mixture prior with one of the mixture components a shrinking prior. We call our approach the mixture M-RA. Application of the mixture M-RA model to both stationary and non-stationary data shows that the mixture M-RA model can handle both types of data, can correctly establish the type of spatial dependence structure in the data (e.g. stationary vs not), and can identify regions of local stationarity.


title: Computational identification of significant actors in paintings through symbols and attributes

abstract: The automatic analysis of fine art paintings presents a number of novel technical challenges to artificial intelligence, computer vision, machine learning, and knowledge representation quite distinct from those arising in the analysis of traditional photographs. The most important difference is that many realist paintings depict stories or episodes in order to convey a lesson, moral, or meaning. One early step in automatic interpretation and extraction of meaning in artworks is the identifications of figures (actors). In Christian art, specifically, one must identify the actors in order to identify the Biblical episode or story depicted, an important step in understanding the artwork. We designed an automatic system based on deep convolutional neural networks and simple knowledge database to identify saints throughout six centuries of Christian art based in large part upon saints symbols or attributes. Our work represents initial steps in the broad task of automatic semantic interpretation of messages and meaning in fine art.


title: Vesicle shape transformations driven by confined active filaments

abstract: In active matter systems, deformable boundaries provide a mechanism to organize internal active stresses and perform work on the external environment. To study a minimal model of such a system, we perform particle-based simulations of an elastic vesicle containing a collection of polar active filaments. The interplay between the active stress organization due to interparticle interactions and that due to the deformability of the confinement leads to a variety of filament spatiotemporal organizations that have not been observed in bulk systems or under rigid confinement, including highly-aligned rings and caps. In turn, these filament assemblies drive dramatic and tunable transformations of the vesicle shape and its dynamics. We present simple scaling models that reveal the mechanisms underlying these emergent behaviors and yield design principles for engineering active materials with targeted shape dynamics.


title: Deep learning-based synthetic-CT generation in radiotherapy and PET: a review

abstract: Recently, deep learning (DL)-based methods for the generation of synthetic computed tomography (sCT) have received significant research attention as an alternative to classical ones. We present here a systematic review of these methods by grouping them into three categories, according to their clinical applications: I) To replace CT in magnetic resonance (MR)-based treatment planning. II) Facilitate cone-beam computed tomography (CBCT)-based image-guided adaptive radiotherapy. III) Derive attenuation maps for the correction of positron emission tomography (PET). Appropriate database searching was performed on journal articles published between January 2014 and December 2020. The DL methods' key characteristics were extracted from each eligible study, and a comprehensive comparison among network architectures and metrics was reported. A detailed review of each category was given, highlighting essential contributions, identifying specific challenges, and summarising the achievements. Lastly, the statistics of all the cited works from various aspects were analysed, revealing the popularity and future trends, and the potential of DL-based sCT generation. The current status of DL-based sCT generation was evaluated, assessing the clinical readiness of the presented methods.


title: Electronic and dynamical properties of CeRh$_{2}$As$_{2}$: Role of Rh$_{2}$As$_{2}$ layers and expected hidden orbital order

abstract: Recently discovered heavy fermion CeRh$_{2}$As$_{2}$ compound crystallizes in the nonsymmorphic $P4/nmm$ symmetry, which enables the occurrence of topological protection. Experimental results show that this material exhibits unusual behavior, which is manifested by the appearance of two superconducting phases. In this work, we uncover and discuss a role of Rh$_{2}$As$_{2}$ layers and their impact on the electronic and dynamical properties of the system. The location of Ce atoms between two non-equivalent layers allows for the realization of hidden orbital order. We point out that the electronic band structure around the Fermi level is associated mostly with Ce $4f$ and Rh $4d$ orbitals and suggest the occurrence of the Lifshitz transition induced by the external magnetic field. We discuss also the role played by the $f$--$d$ orbital hybridization in the electronic band structure.


title: A Pointwise Inequality for Derivatives of Solutions of the Heat Equation in Bounded Domains

abstract: Let $u(t,x)$ be a solution of the heat equation in $\mathbb{R}^n$. Then, each $k-$th derivative also solves the heat equation and satisfies a maximum principle, the largest $k-$th derivative of $u(t,x)$ cannot be larger than the largest $k-$th derivative of $u(0,x)$. We prove an analogous statement for the solution of the heat equation on bounded domains $\Omega \subset \mathbb{R}^n$ with Dirichlet boundary conditions. As an application, we give a new and fairly elementary proof of the sharp growth of the second derivatives of Laplacian eigenfunction $-\Delta \phi_k = \lambda_k \phi_k$ with Dirichlet conditions on smooth domains $\Omega \subset \mathbb{R}^n$.


title: Insight of the Green's function as a defect state in a boundary value problem

abstract: A new perspective of the Green's function in a boundary value problem as the only eigenstate in an auxiliary formulation is introduced. In this treatment, the Green's function can be perceived as a defect state in the presence of a $\delta$-function potential, the height of which depends on the Green's function itself. This approach is illustrated in one-dimensional and two-dimensional Helmholtz equation problems, with an emphasis on systems that are open and have a non-Hermitian potential. We then draw an analogy between the Green's function obtained this way and a chiral edge state circumventing a defect in a topological lattice, which shines light on the local minimum of the Green's function at the source position.


title: SIMUS: an open-source simulator for medical ultrasound imaging. Part I: theory & examples

abstract: Background and Objective: Computational ultrasound imaging has become a well-established methodology in the ultrasound community. Simulations of ultrasound sequences and images allow the study of innovative techniques in terms of emission strategy, beamforming, and probe design. There is a wide spectrum of software dedicated to ultrasound imaging, each having its specificities in its applications and the numerical method. Methods: We describe in this two-part paper a new ultrasound simulator (SIMUS) for MATLAB, which belongs to the MATLAB UltraSound Toolbox (MUST). The SIMUS software simulates acoustic pressure fields and radiofrequency RF signals for uniform linear or convex probes. SIMUS is an open-source software whose features are 1) ease of use, 2) time-harmonic analysis, 3) pedagogy. The main goal was to offer a comprehensive turnkey tool, along with a detailed theory for pedagogical and research purposes. Results: This article describes in detail the underlying linear theory of SIMUS and provides examples of simulated acoustic fields and ultrasound images. The accompanying article (part II) is devoted to the comparison of SIMUS with several software: Field II, k-Wave, FOCUS, and the Verasonics simulator. The MATLAB open codes for the simulator SIMUS are distributed under the terms of the GNU Lesser General Public License, and can be downloaded from https://www.biomecardio.com/MUST. Conclusions: The simulations described in the accompanying paper (Part II) show that SIMUS can be used for realistic simulations in medical ultrasound imaging.


title: On the number of fixed points of the map $\gamma$

abstract: We recursively define a sequence $\{F_{n,k}\}_{n,k\in\mathbb N }$ and we prove that such sequence contains only the symbols $\{0,1\}.$ We investigate some number-theoretic properties of such sequence and of the way it can be generated. The number $F_n$ can be interpreted as the number of fixed points of semilength $n$ of the map $\gamma$ introduced by Barnabei et al. Our results partially answer conjectures posed by Cori et al.


title: Spin-orbit-entangled electronic phases in 4$d$ and 5$d$ transition-metal compounds

abstract: Complex oxides with $4d$ and $5d$ transition-metal ions recently emerged as a new paradigm in correlated electron physics, due to the interplay between spin-orbit coupling and electron interactions. For $4d$ and $5d$ ions, the spin-orbit coupling, $\zeta$, can be as large as 0.2-0.4 eV, which is comparable with and often exceeds other relevant parameters such as Hund's coupling $J_{\rm H}$, noncubic crystal field splitting $\Delta$, and the electron hopping amplitude $t$. This gives rise to a variety of spin-orbit-entangled degrees of freedom and, crucially, non-trivial interactions between them that depend on the $d$-electron configuration, the chemical bonding, and the lattice geometry. Exotic electronic phases often emerge, including spin-orbit assisted Mott insulators, quantum spin liquids, excitonic magnetism, multipolar orderings and correlated topological semimetals. This paper provides a selective overview of some of the most interesting spin-orbit-entangled phases that arise in $4d$ and $5d$ transition-metal compounds.


title: Hawkes Processes on Graphons

abstract: We propose a novel framework for modeling multiple multivariate point processes, each with heterogeneous event types that share an underlying space and obey the same generative mechanism. Focusing on Hawkes processes and their variants that are associated with Granger causality graphs, our model leverages an uncountable event type space and samples the graphs with different sizes from a nonparametric model called {\it graphon}. Given those graphs, we can generate the corresponding Hawkes processes and simulate event sequences. Learning this graphon-based Hawkes process model helps to 1) infer the underlying relations shared by different Hawkes processes; and 2) simulate event sequences with different event types but similar dynamics. We learn the proposed model by minimizing the hierarchical optimal transport distance between the generated event sequences and the observed ones, leading to a novel reward-augmented maximum likelihood estimation method. We analyze the properties of our model in-depth and demonstrate its rationality and effectiveness in both theory and experiments.


title: Analytic characterization of high dimension weighted special atom spaces

abstract: Special atom spaces have been around for quite awhile since the introduction of atoms by R. Coifman in his seminal paper who led to another proof that the dual of the Hardy space $H^1$ is in fact the space of functions of bounded means oscillations (BMO). Special atom spaces enjoy quite a few attributes of their own, among which the fact that they have an analytic extension to the unit disc. Recently, an extension of special atom spaces to higher dimensions was proposed, making ripe the possible exploration of the above extension in higher dimensions. In this paper we propose an analytic characterization of special atom spaces in higher dimensions.