title: Galaxy morphology rules out astrophysically relevant Hu-Sawicki $f(R)$ gravity

abstract: $f(R)$ is a paradigmatic modified gravity theory that typifies extensions to General Relativity with new light degrees of freedom and hence screened fifth forces between masses. These forces produce observable signatures in galaxy morphology, caused by a violation of the weak equivalence principle due to a differential impact of screening among galaxies' mass components. We compile statistical datasets of two morphological indicators -- offsets between stars and gas in galaxies and warping of stellar disks -- and use them to constrain the strength and range of a thin-shell-screened fifth force. This is achieved by applying a comprehensive set of upgrades to past work (Desmond et al 2018a,b): we construct a robust galaxy-by-galaxy Bayesian forward model for the morphological signals, including full propagation of uncertainties in the input quantities and marginalisation over an empirical model describing astrophysical noise. Employing more stringent data quality cuts than previously we find no evidence for a screened fifth force of any strength $\Delta G/G_\text{N}$ in the Compton wavelength range $0.3-8$ Mpc, setting a $1\sigma$ bound of $\Delta G/G_\text{N}<0.8$ at $\lambda_C=0.3$ Mpc that strengthens to $\Delta G/G_\text{N}<3\times10^{-5}$ at $\lambda_C=8$ Mpc. These are the tightest bounds to date beyond the Solar System by over an order of magnitude. For the Hu-Sawicki model of $f(R)$ with $n=1$ we require a background scalar field value $f_{R0} < 1.4 \times 10^{-8}$, forcing practically all astrophysical objects to be screened. We conclude that this model can have no relevance to astrophysics or cosmology.


title: A consistent variational formulation of Bishop nonlocal rods

abstract: Thick rods are employed in Nanotechnology to build modern electro mechanical systems. Design and optimization of such structures can be carried out by nonlocal continuum mechanics which is computationally convenient when compared to atomistic strategies. Bishop's kinematics is able to describe small-scale thick rods if a proper mathematical model of nonlocal elasticity is formulated to capture size effects. In all papers on the matter, nonlocal contributions are evaluated by replacing Eringen's integral convolution with the consequent (but not equivalent) differential equation governed by Helmholtz's differential operator. As notorious in integral equation theory, this replacement is possible for convolutions, defined in unbounded domains, governed by averaging kernels which are Green's functions of differential operators. Indeed, Eringen himself, in order to study nonlocal problems defined in unbounded domains, such as screw dislocations and wave propagations, suggested to replace integro-differential equations with differential conditions. A different scenario appears in Bishop rod mechanics where nonlocal integral convolutions are defined in bounded structural domains, so that Eringen's nonlocal differential equation has to be supplemented with additional boundary conditions. The objective is achieved by formulating the governing nonlocal equations by a proper variational statement. The new methodology provides an amendment of previous contributions in literature and is illustrated by investigating the elastostatic behavior of simple structural schemes. Exact solutions of Bishop rods are evaluated in terms of nonlocal parameter and cross-section gyration radius. Both hardening and softening structural responses are predictable with a suitable tuning of the parameters.


title: Direct observation of charge separation in an organic light harvesting system by femtosecond time-resolved XPS

abstract: The ultrafast dynamics of photon-to-charge conversion in an organic light harvesting system is studied by femtosecond time-resolved X-ray photoemission spectroscopy (TR-XPS) at the free-electron laser FLASH. This novel experimental technique provides site-specific information about charge separation and enables the monitoring of free charge carrier generation dynamics on their natural timescale, here applied to the model donor-acceptor system CuPc:C$_{60}$. A previously unobserved channel for exciton dissociation into mobile charge carriers is identified, providing the first direct, real-time characterization of the timescale and efficiency of charge generation from low-energy charge-transfer states in an organic heterojunction. The findings give strong support to the emerging realization that charge separation even from energetically disfavored excitonic states is contributing significantly, indicating new options for light harvesting in organic heterojunctions.


title: Moving object detection for visual odometry in a dynamic environment based on occlusion accumulation

abstract: Detection of moving objects is an essential capability in dealing with dynamic environments. Most moving object detection algorithms have been designed for color images without depth. For robotic navigation where real-time RGB-D data is often readily available, utilization of the depth information would be beneficial for obstacle recognition. Here, we propose a simple moving object detection algorithm that uses RGB-D images. The proposed algorithm does not require estimating a background model. Instead, it uses an occlusion model which enables us to estimate the camera pose on a background confused with moving objects that dominate the scene. The proposed algorithm allows to separate the moving object detection and visual odometry (VO) so that an arbitrary robust VO method can be employed in a dynamic situation with a combination of moving object detection, whereas other VO algorithms for a dynamic environment are inseparable. In this paper, we use dense visual odometry (DVO) as a VO method with a bi-square regression weight. Experimental results show the segmentation accuracy and the performance improvement of DVO in the situations. We validate our algorithm in public datasets and our dataset which also publicly accessible.


title: Poly-freeness of large even Artin groups

abstract: We prove that any large even Artin group is poly-free and that any even Artin group based on a triangle graph is also poly-free.


title: A polynomial size model with implicit SWAP gate counting for exact qubit reordering

abstract: Due to the physics behind quantum computing, quantum circuit designers must adhere to the constraints posed by the limited interaction distance of qubits. Existing circuits need therefore to be modified via the insertion of SWAP gates, which alter the qubit order by interchanging the location of two qubits' quantum states. We consider the Nearest Neighbor Compliance problem on a linear array, where the number of required SWAP gates is to be minimized. We introduce an Integer Linear Programming model of the problem of which the size scales polynomially in the number of qubits and gates. Furthermore, we solve $131$ benchmark instances to optimality using the commercial solver CPLEX. The benchmark instances are substantially larger in comparison to those evaluated with exact methods before. The largest circuits contain up to $18$ qubits or over $100$ quantum gates. This formulation also seems to be suitable for developing heuristic methods since (near) optimal solutions are discovered quickly in the search process.


title: Heavy Positive Ion Groups in Titan's Ionosphere from Cassini Plasma Spectrometer IBS Observations

abstract: Titan's ionosphere contains a plethora of hydrocarbons and nitrile cations and anions as measured by the Ion Neutral Mass Spectrometer and Cassini Plasma Spectrometer (CAPS) onboard the Cassini spacecraft. Data from the CAPS Ion Beam Spectrometer (IBS) sensor have been examined for five close encounters of Titan during 2009. The high relative velocity of Cassini with respect to the cold ions in Titan's ionosphere allows CAPS IBS to function as a mass spectrometer. Positive ion masses between 170 and 310 u/q are examined with ion mass groups identified between 170 and 275 u/q containing between 14 and 21 heavy (carbon/nitrogen/oxygen) atoms. These groups are the heaviest positive ion groups reported so far from the available in situ ion data at Titan. The ion group peaks are found to be consistent with masses associated with Polycyclic Aromatic Compounds (PAC), including Polycyclic Aromatic Hydrocarbon (PAH) and nitrogen-bearing polycyclic aromatic molecular ions. The ion group peak identifications are compared with previously proposed neutral PAHs and are found to be at similar masses, supporting a PAH interpretation. The spacing between the ion group peaks is also investigated, finding a spacing of 12 or 13 u/q indicating the addition of C or CH. Lastly, the occurrence of several ion groups is seen to vary across the five flybys studied, possibly relating to the varying solar radiation conditions observed across the flybys. These findings further the understanding between the low mass ions and the high mass negative ions, as well as with aerosol formation in Titan's atmosphere.


title: Dominant negative Poynting effect for soft tissues

abstract: We identify three distinct shearing modes for simple shear deformations of transversely isotropic soft tissue which allow for both positive and negative Poynting effects (that is, they require compressive and tensile lateral normal stresses, respectively, in order to maintain simple shear). The positive Poynting effect is that usually found for isotropic rubber. Here, specialisation of the general results to three strain-energy functions which are quadratic in the anisotropic invariants, linear in the isotropic strain invariants and consistent with the linear theory suggests that there are two Poynting effects which can accompany the shearing of soft tissue: a dominant negative effect in one mode of shear and a relatively small positive effect in the other two modes. We propose that the relative inextensibility of the fibres relative to the matrix is the primary mechanism behind this large negative Poynting effect.


title: Hardness and approximation of the Probabilistic p-Center problem under Pressure

abstract: The Probabilistic p-Center problem under Pressure (Min PpCP) is a variant of the usual p-Center problem we recently introduced in the context of wildfire management. The problem is to locate p shelters minimizing the maximum distance people will have to cover to reach the closest accessible shelter in case of fire. The landscape is divided into zones and is modeled as an edge-weighted graph with vertices corresponding to zones and edges corresponding to direct connections between two adjacent zones. The risk associated with fire outbreaks is modeled using a finite set of fire scenarios. Each scenario corresponds to a fire outbreak on a single zone (i.e., on a vertex) with the main consequence of modifying evacuation paths in two ways. First, an evacuation path cannot pass through the vertex on fire. Second, the fact that someone close to the fire may not take rational decisions when selecting a direction to escape is modeled using new kinds of evacuation paths. In this paper, for a given instance of Min PpCP defined by an edge-weighted graph G=(V,E,L) and an integer p, we characterize the set of feasible solutions of Min PpCP. We prove that Min PpCP cannot be approximated with a ratio less than 56/55 on subgrids (subgraphs of grids) of degree at most 3. Then, we propose some approximation results for Min PpCP. These results require approximation results for two variants of the (deterministic) Min p-Center problem called Min MAC p-Center and Min Partial p-Center.


title: Gent models for the inflation of spherical balloons

abstract: We revisit an iconic deformation of non-linear elasticity: the inflation of a rubber spherical thin shell. We use the 3-parameter Mooney and Gent-Gent (GG) phenomenological models to explain the stretch-strain curve of a typical inflation, as these two models cover a wide spectrum of known models for rubber, including the Varga, Mooney-Rivlin, one-term Ogden, Gent-Thomas and Gent models. We find that the basic physics of inflation exclude the Varga, one-term Ogden and Gent-Thomas models. We find the link between the exact solution of non-linear elasticity and the membrane and Young-Laplace theories often used a priori in the literature. We compare the performance of both models on fitting the data for experiments on rubber balloons and animal bladder. We conclude that the GG model is the most accurate and versatile model on offer for the modelling of rubber balloon inflation.


title: DeltaGAN: Towards Diverse Few-shot Image Generation with Sample-Specific Delta

abstract: Learning to generate new images for a novel category based on only a few images, named as few-shot image generation, has attracted increasing research interest. Several state-of-the-art works have yielded impressive results, but the diversity is still limited. In this work, we propose a novel Delta Generative Adversarial Network (DeltaGAN), which consists of a reconstruction subnetwork and a generation subnetwork. The reconstruction subnetwork captures intra-category transformation, i.e., "delta", between same-category pairs. The generation subnetwork generates sample-specific "delta" for an input image, which is combined with this input image to generate a new image within the same category. Besides, an adversarial delta matching loss is designed to link the above two subnetworks together. Extensive experiments on five few-shot image datasets demonstrate the effectiveness of our proposed method.


title: A robust anisotropic hyperelastic formulation for the modelling of soft tissue

abstract: The Holzapfel-Gasser-Ogden (HGO) model for anisotropic hyperelastic behaviour of collagen fibre reinforced materials was initially developed to describe the elastic properties of arterial tissue, but is now used extensively for modelling a variety of soft biological tissues. Such materials can be regarded as incompressible, and when the incompressibility condition is adopted the strain energy \Psi of the HGO model is a function of one isotropic and two anisotropic deformation invariants. A compressible form (HGO-C model) is widely used in finite element simulations whereby the isotropic part of \Psi is decoupled into volumetric and isochoric parts and the anisotropic part of \Psi is expressed in terms of isochoric invariants. Here, by using three simple deformations (pure dilatation, pure shear and uniaxial stretch), we demonstrate that the compressible HGO-C formulation does not correctly model compressible anisotropic material behaviour, because the anisotropic component of the model is insensitive to volumetric deformation due to the use of isochoric anisotropic invariants. In order to correctly model compressible anisotropic behaviour we present a modified anisotropic (MA) model, whereby the full anisotropic invariants are used, so that a volumetric anisotropic contribution is represented. The MA model correctly predicts an anisotropic response to hydrostatic tensile loading, whereby a sphere deforms into an ellipsoid. It also computes the correct anisotropic stress state for pure shear and uniaxial deformation. To look at more practical applications, we developed a finite element user-defined material subroutine for the simulation of stent deployment in a slightly compressible artery. Significantly higher stress triaxiality and arterial compliance are computed when the full anisotropic invariants are used (MA model) instead of the isochoric form (HGO-C model).


title: Multichannel superconductivity of monolayer FeSe on SrTiO$_3$: Interplay of spin fluctuations and electron-phonon interaction

abstract: We investigate the effects of electron-phonon coupling, as well as of spin and charge fluctuations on the superconducting state in a single layer of FeSe on SrTiO$_3$ substrate. These three bosonic mediators of Cooper pairing are treated on equal footing in a multichannel, full-bandwidth, multiband, and anisotropic Eliashberg theory of the interacting state. Our self-consistent calculations show that an s-wave symmetry of the superconducting gap is compatible only with a complete absence of spin fluctuations. When spin fluctuations are present, the sign-changing nodeless d-wave pairing symmetry is always obtained, yet the essential ingredient for explaining the gap magnitude and critical temperature is still the interfacial electron-phonon interaction.


title: Completing $1/m_b^3$ corrections to non-leptonic bottom-to-up-quark decays

abstract: We compute the contributions of dimension six two-quark operators to the non-leptonic decay width of heavy hadrons due to the flavor changing bottom-to-up-quark transition in the heavy quark expansion. Analytical expressions for the Darwin term $\rho_D$ and the spin-orbit term $\rho_{\rm LS}$ are obtained with leading order accuracy.


title: Volume-Scaled Common Nearest Neighbor Clustering Algorithm with Free-Energy Hierarchy

abstract: The combination of Markov state modeling (MSM) and molecular dynamics (MD) simulations has been shown in recent years to be a valuable approach to unravel the slow processes of molecular systems with increasing complexity. While the algorithms for intermediate steps in the MSM workflow like featurization and dimensionality reduction have been specifically adapted for MD data sets, conventional clustering methods are generally applied for the discretization step. This work adds to recent efforts to develop specialized density-based clustering algorithms for the Boltzmann-weighted data from MD simulations. We introduce the volume-scaled common nearest neighbor (vs-CNN) clustering that is an adapted version of the common nearest neighbor (CNN) algorithm. A major advantage of the proposed algorithm is that the introduced density-based criterion directly links to a free-energy notion via Boltzmann inversion. Such a free-energy perspective allows for a straightforward hierarchical scheme to identify conformational clusters at different levels of a generally rugged free-energy landscape of complex molecular systems.

categories: physics.comp-ph cond-mat.soft physics.bio-ph physics.chem-ph


title: Revisiting Consensus-Based Energy-Management in Smart Grid with Transmission Losses and Directed Communication

abstract: We discovered a deficiency in Algorithm 1 and Theorem 3 of [1]. The algorithm called CEMA aims to solve an energy management problem distributively. However, by means of a counter example, we show that Theorem 2 and 3 of [1] contradict each other in the case of a valid scenario, proving that the suggested algorithm does not always find the optimum. Furthermore, we provide theoretic results, showing that Theorem 3 of [1] does not hold generally. At last, we provide a rectification by adjusting the algorithm and the corresponding proof of Theorem 3.


title: Nonlinear vibration localisation in a symmetric system of two coupled beams

abstract: We report nonlinear vibration localisation in a system of two symmetric weakly coupled nonlinear oscillators. A two degree-of-freedom model with piecewise linear stiffness shows bifurcations to localised solutions. An experimental investigation employing two weakly coupled beams touching against stoppers for large vibration amplitudes confirms the nonlinear localisation.


title: Differential molecule-cavity mode coupling in soft-cavities

abstract: The way molecules absorb, transfer, and emit light can be dramatically modified by coupling them to optical cavities. The extent of the modification is often defined by the cavity-molecule coupling strength. Evaluating this coupling strength for different types of modes supported by a cavity is crucial in designing cavities for molecule-cavity coupling. Here we probe a unique multimode cavity, a dielectric microsphere, also called a soft-cavity, which supports two distinct types of mode, dark-field scattering (DFS) modes and whispering gallery modes (WGM). Though seemingly similar, these modes show different characteristics such as spatial electric field profile, resonance line-width etc. We investigated coupling of a mono-layer of J-aggregated dye molecules and a dielectric plastic microsphere using two techniques, far-field excitation and evanescent excitation to generate DFS modes and WGMs respectively. We found that using WGMs we observe a clear signature of strong coupling, whereas with DFS modes we do not. We compared our experimental data to a simple coupled oscillator model and performed finite-element method based numerical simulations to provide a clearer understanding of our experimental findings.


title: Online Knowledge Base for Designing Shape-changing Interfaces using Modular Workshop Elements

abstract: Building and maintaining knowledge about specific interface technologies is a challenge. Current solutions include standard file-based document repositories, wikis, and other online tools. However, these solutions are often only available in intranets, become outdated and do not support the acquisition of knowledge in an efficient manner. The effort to gain an overview and detailed knowledge about novel interface technologies can be overwhelming and requires to research and read many technical reports and scientific publications. We propose to implement open source online knowledge bases that include building blocks for creating custom workshops to understand and apply the contained knowledge. We demonstrate this concept with a knowledge base for shape-changing interfaces (SCI-KB). The SCI-KB is hosted online at GitHub and fosters collaborative creation of knowledge elements accompanied by practical exercises and workshop elements that can be combined and adapted by individuals or groups of people new to the topic of shape-changing interfaces.


title: On the analyticity of solutions to non-linear elliptic partial differential systems

abstract: We give an easy proof of the fact that $C^\infty$ solutions to non-linear elliptic equations of second order $$ \phi(x, u, D u, D^2 u)=0 $$ are analytic. Following ideas of Kato, the proof uses an inductive estimate for suitable weighted derivatives. We then conclude the proof using Cauchy's method of majorants}.


title: Hearts for commutative noetherian rings: torsion pairs and derived equivalences

abstract: Over a commutative noetherian ring $R$, the prime spectrum controls, via the assignment of support, the structure of both $\mathsf{Mod}(R)$ and $\mathsf{D}(R)$. We show that, just like in $\mathsf{Mod}(R)$, the assignment of support classifies hereditary torsion pairs in the heart of any nondegenerate compactly generated $t$-structure of $\mathsf{D}(R)$. Moreover, we investigate whether these $t$-structures induce derived equivalences, obtaining a new source of Grothendieck categories which are derived equivalent to $\mathsf{Mod}(R)$.


title: Accelerating MPC by online detection of state space sets with common optimal feedback laws

abstract: Model predictive control (MPC) samples a generally unknown and complicated feedback law point by point. The solution for the current state $x$ contains, however, more information than only the optimal signal $u$ for this particular state. In fact, it provides an optimal affine feedback law $x\rightarrow u(x)$ on a polytope $\Pi \subset \mathbb{R}^n$, i.e., on a full-dimensional state space set. It is an obvious idea to reuse this affine feedback law as long as possible. Reusing it on its polytope $\Pi$ is too conservative, however, because any $\Pi$ is a state space set with a common affine law $x\rightarrow (u_0^\prime (x), \dots, u_{N-1}^\prime (x))^\prime\in\mathbb{R}^{Nm}$ for the entire horizon $N$. We show a simple criterion exists for identifying the polytopes that have a common $x\rightarrow u_0(x)$, but may differ with respect to $u_1(x), \dots, u_{N-1}(x)$. Because this criterion is too computationally expensive for an online use, we introduce a simple heuristics for the fast construction of a subset of the polytopes of interest. Computational examples show (i) a considerable fraction of QPs can be avoided (10% to 40%) and (ii) the heuristics results in a reduction very close to the maximum one that could be achieved if the explicit solution was available. We stress the proposed approach is intended for use in online MPC and it does not require the explicit solution.


title: On the Capacity Enlargement of Gaussian Broadcast Channels with Passive Noisy Feedback

abstract: It is well known that the capacity region of an average transmit power constrained Gaussian Broadcast Channel (GBC) with independent noise realizations at the receivers is enlarged by the presence of causal noiseless feedback. Capacity region enlargement is also known to be possible by using only passive noisy feedback, when the GBC has identical noise variances at the receivers. The last fact remains true even when the feedback noise variance is very high, and available only from one of the receivers. While such capacity enlargements are feasible for several other feedback models in the Gaussian BC setting, it is also known that feedback does not change the capacity region for physically degraded broadcast channels. In this paper, we consider a two user GBC with independent noise realizations at the receivers, where the feedback links from the receivers are corrupted by independent additive Gaussian noise processes. We investigate the set of four noise variances, two forward and two feedback, for which no capacity enlargement is possible. A sharp characterization of this region is derived, i.e., any quadruple outside the presented region will lead to a capacity enlargement, whereas quadruples inside will leave the capacity region unchanged. Our results lead to the conclusion that when the forward noise variances are different, too noisy a feedback from one of the receivers alone is not always beneficial for enlarging the capacity region, be it from the stronger user or the weaker one, in sharp contrast to the case of equal forward noise variances.


title: Real-space multiple scattering theory for superconductors with impurities

abstract: We implement the Bogoliubov-de Gennes (BdG) equation in real-space using the screened Korringa-Kohn-Rostoker (KKR) method. This allows us to solve, self-consistently, the superconducting state for 3d crystals including substitutional impurities with a full normal-state DFT band structure. We apply the theoretical framework to bulk Nb with impurities. Without impurities, Nb has an anisotropic gap structure with two distinct peaks around the Fermi level. In the presence of non-magnetic impurities those peaks are broadened due to the scattering between the two bulk superconducting gaps, however the peaks remain separated. As a second example of self-consistent real-space solutions of the BdG equations we examine superconducting clusters embedded within a non-superconducting bulk metallic host. This allows us to estimate the coherence length of the superconductor and we show that, within our framework, the coherence length of the superconductor is related to the inverse of the gap size, just as in bulk BCS theory.


title: Seifert fibrations of lens spaces over non-orientable bases

abstract: We classify the Seifert fibrations of lens spaces where the base orbifold is non-orientable. This is an addendum to our earlier paper `Seifert fibrations of lens spaces'. We correct Lemma 4.1 of that paper and fill the gap in the classification that resulted from the erroneous lemma.


title: Improving Linux-Kernel Tests for LockDoc with Feedback-driven Fuzzing

abstract: LockDoc is an approach to extract locking rules for kernel data structures from a dynamic execution trace recorded while the system is under a benchmark load. These locking rules can e.g. be used to locate synchronization bugs. For high rule precision and thorough bug finding, the approach heavily depends on the choice of benchmarks: They must trigger the execution of as much code as possible in the kernel subsystem relevant for the targeted data structures. However, existing test suites such as those provided by the Linux Test Project (LTP) only achieve -- in the case of LTP -- about 35 percent basic-block coverage for the VFS subsystem, which is the relevant subsystem when extracting locking rules for filesystem-related data structures. In this article, we discuss how to complement the LTP suites to improve the code coverage for our LockDoc scenario. We repurpose syzkaller -- a coverage-guided fuzzer with the goal to validate the robustness of kernel APIs -- to 1) not aim for kernel crashes, and to 2) maximize code coverage for a specific kernel subsystem. Thereby, we generate new benchmark programs that can be run in addition to the LTP, and increase VFS basic-block coverage by 26.1 percent.


title: Typestates to Automata and back: a tool

abstract: Development of software is an iterative process. Graphical tools to represent the relevant entities and processes can be helpful. In particular, automata capture well the intended execution flow of applications, and are thus behind many formal approaches, namely behavioral types. Typestate-oriented programming allow us to model and validate the intended protocol of applications, not only providing a top-down approach to the development of software, but also coping well with compositional development. Moreover, it provides important static guarantees like protocol fidelity and some forms of progress. Mungo is a front-end tool for Java that associates a typestate describing the valid orders of method calls to each class, and statically checks that the code of all classes follows the prescribed order of method calls. To assist programming with Mungo, as typestates are textual descriptions that are terms of an elaborate grammar, we developed a tool that bidirectionally converts typestates into an adequate form of automata, providing on one direction a visualization of the underlying protocol specified by the typestate, and on the reverse direction a way to get a syntactically correct typestate from the more intuitive automata representation.


title: Probably Approximately Correct Explanations of Machine Learning Models via Syntax-Guided Synthesis

abstract: We propose a novel approach to understanding the decision making of complex machine learning models (e.g., deep neural networks) using a combination of probably approximately correct learning (PAC) and a logic inference methodology called syntax-guided synthesis (SyGuS). We prove that our framework produces explanations that with a high probability make only few errors and show empirically that it is effective in generating small, human-interpretable explanations.


title: Edge dynamics of an Integer Quantum Hall system

abstract: In this master thesis work the linear and non linear edge dynamics of a non-interacting system of fermions in a Integer Quantum Hall state is theoretically and numerically studied.


title: The Boon and Bane of Cross-Signing: Shedding Light on a Common Practice in Public Key Infrastructures

abstract: Public Key Infrastructures (PKIs) with their trusted Certificate Authorities (CAs) provide the trust backbone for the Internet: CAs sign certificates which prove the identity of servers, applications, or users. To be trusted by operating systems and browsers, a CA has to undergo lengthy and costly validation processes. Alternatively, trusted CAs can cross-sign other CAs to extend their trust to them. In this paper, we systematically analyze the present and past state of cross-signing in the Web PKI. Our dataset (derived from passive TLS monitors and public CT logs) encompasses more than 7 years and 225 million certificates with 9.3 billion trust paths. We show benefits and risks of cross-signing. We discuss the difficulty of revoking trusted CA certificates where, worrisome, cross-signing can result in valid trust paths to remain after revocation; a problem for non-browser software that often blindly trusts all CA certificates and ignores revocations. However, cross-signing also enables fast bootstrapping of new CAs, e.g., Let's Encrypt, and achieves a non-disruptive user experience by providing backward compatibility. In this paper, we propose new rules and guidance for cross-signing to preserve its positive potential while mitigating its risks.


title: Oumuamuas passing through molecular clouds

abstract: The detections of 1I/Oumuamua and 2I/Borisov within just two years demonstrate impressively that interstellar objects (ISOs) must be common in the Milky Way. Once released from their parent system, these ISOs travel for Gyr through interstellar space. While often imagined as empty, interstellar space contains gas and dust most prominent in the form of molecular clouds. Performing numerical simulations, we test how often ISOs cross such molecular clouds. We find that the ISOs pass amazingly often through molecular clouds. In the solar neighbourhood, ISOs typically spend 0.1-0.2% of their journey inside molecular clouds, for relative slow ISOs ($<$ 5 km/s) this can increase to 1-2%, equivalent to 10 - 20 Myr per Gyr. Thus the dynamically youngest ISOs spend the longest time in molecular clouds. In other words, molecular clouds must mainly contain relatively young ISOs ($<$ 1-2 Gyr). Thus the half-life of the seeding process by ISOs is substantially shorter than a stellar lifetime. The actual amount of time spent in MCs decreases with distance to the Galactic Centre. We find that ISOs pass so often through MCs that backtracing their path to find their parent star beyond 250 Myr seems beyond the point. Besides, we give a first estimate of the ISO density depending on the galactic centre distance based on the stellar distribution.


title: The forbidden region for random zeros: appearance of quadrature domains

abstract: Our main discovery is a surprising interplay between quadrature domains on the one hand, and the zero process of the Gaussian Entire Function (GEF) on the other. Specifically, consider the GEF conditioned on the rare hole event that there are no zeros in a given large Jordan domain. We show that in the natural scaling limit, a quadrature domain enclosing the hole emerges as a forbidden region, where the zero density vanishes. Moreover, we give a description of those holes for which the forbidden region is a disk. The connecting link between random zeros and potential theory is supplied by a constrained extremal problem for the Zeitouni-Zelditch functional. To solve this problem, we recast it in terms of a seemingly novel obstacle problem, where the solution is forced to be harmonic inside the hole.


title: Document-level Neural Machine Translation with Document Embeddings

abstract: Standard neural machine translation (NMT) is on the assumption of document-level context independent. Most existing document-level NMT methods are satisfied with a smattering sense of brief document-level information, while this work focuses on exploiting detailed document-level context in terms of multiple forms of document embeddings, which is capable of sufficiently modeling deeper and richer document-level context. The proposed document-aware NMT is implemented to enhance the Transformer baseline by introducing both global and local document-level clues on the source end. Experiments show that the proposed method significantly improves the translation performance over strong baselines and other related studies.


title: Dealing with Incompatibilities among Procedural Goals under Uncertainty

abstract: By considering rational agents, we focus on the problem of selecting goals out of a set of incompatible ones. We consider three forms of incompatibility introduced by Castelfranchi and Paglieri, namely the terminal, the instrumental (or based on resources), and the superfluity. We represent the agent's plans by means of structured arguments whose premises are pervaded with uncertainty. We measure the strength of these arguments in order to determine the set of compatible goals. We propose two novel ways for calculating the strength of these arguments, depending on the kind of incompatibility that exists between them. The first one is the logical strength value, it is denoted by a three-dimensional vector, which is calculated from a probabilistic interval associated with each argument. The vector represents the precision of the interval, the location of it, and the combination of precision and location. This type of representation and treatment of the strength of a structured argument has not been defined before by the state of the art. The second way for calculating the strength of the argument is based on the cost of the plans (regarding the necessary resources) and the preference of the goals associated with the plans. Considering our novel approach for measuring the strength of structured arguments, we propose a semantics for the selection of plans and goals that is based on Dung's abstract argumentation theory. Finally, we make a theoretical evaluation of our proposal.


title: Numerical Methods to Compute Stresses and Displacements from Cellular Forces: Application to the Contraction of Tissue

abstract: We consider a mathematical model for wound contraction, which is based on solving a momentum balance under the assumptions of isotropy, homogeneity, Hooke's Law, infinitesimal strain theory and point forces exerted by cells. However, point forces, described by Dirac Delta distributions lead to a singular solution, which in many cases may cause trouble to finite element methods due to a low degree of regularity. Hence, we consider several alternatives to address point forces, that is, whether to treat the region covered by the cells that exert forces as part of the computational domain or as 'holes' in the computational domain. The formalisms develop into the immersed boundary approach and the 'hole approach', respectively. Consistency between these approaches is verified in a theoretical setting, but also confirmed computationally. However, the 'hole approach' is much more expensive and complicated for its need of mesh adaptation in the case of migrating cells while it increases the numerical accuracy, which makes it hard to adapt to the multi-cell model. Therefore, for multiple cells, we consider the polygon that is used to approximate the boundary of cells that exert contractile forces. It is found that a low degree of polygons, in particular triangular or square shaped cell boundaries, already give acceptable results in engineering precision, so that it is suitable for the situation with a large amount of cells in the computational domain.


title: A study of irreducible components of Springer fibers using quiver varieties

abstract: It is a remarkable theorem by Maffei--Nakajima that the Slodowy variety, which is a subvariety of the resolution of the nilpotent cone, can be realized as a Nakajima quiver variety of type A. However, the isomorphism is rather implicit as it takes to solve a system of equations in which the variables are linear maps. In this paper, we construct solutions to this system under certain assumptions. This establishes an explicit and efficient way to compute the image of a complete flag contained in the Slodowy variety under the Maffei--Nakajima isomorphism and describe these flags in terms of quiver representations. As Slodowy varieties contain Springer fibers naturally, we can use these results to provide an explicit description of the irreducible components of two-row Springer fibers in terms of a family of kernel relations via quiver representations.


title: Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

abstract: Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events.

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


title: Analysis of the Convergence Speed of the Arimoto-Blahut Algorithm by the Second Order Recurrence Formula

abstract: In this paper, we investigate the convergence speed of the Arimoto-Blahut algorithm. For many channel matrices the convergence is exponential, but for some channel matrices it is slower than exponential. By analyzing the Taylor expansion of the defining function of the Arimoto-Blahut algorithm, we will make the conditions clear for the exponential or slower convergence. The analysis of the slow convergence is new in this paper. Based on the analysis, we will compare the convergence speed of the Arimoto-Blahut algorithm numerically with the values obtained in our theorems for several channel matrices. The purpose of this paper is a complete understanding of the convergence speed of the Arimoto-Blahut algorithm.


title: Effects of neutrino magnetic moment and charge radius constraints and medium modifications of the nucleon form factors on the neutrino mean free path in dense matter

abstract: Effects of neutrino charge radius and magnetic moment constraints obtained from the astrophysical observations and reactor experiments, as well as in-medium modifications of the weak and electromagnetic nucleon form factors of the matter on the neutrino electroweak interaction with dense matter, are estimated. We use a relativistic mean-field model for the in-medium effective nucleon mass and quark-meson coupling model for nucleon form factors. We analyze the neutrino scattering cross section and mean free path in cold nuclear matter. We find that the increase of the cross section relative to that without neutrino form factors results in the decrease of the neutrino mean free path when neutrino form factors and the in-medium modifications of the nucleon weak and electromagnetic form factors are simultaneously considered. The quenching of the neutrino mean free path is evaluated to be about 12-58\% for the values of $\mu_\nu = 3 \times 10^{-12} \mu_B$ and $R_\nu = 3.5 \times 10^{-5}~\textrm{MeV}^{-1}$ compared with that obtained for the $\mu_\nu =0$ and $R_\nu =0$. The decrease of the neutrino mean free path is expected to decelerate the cooling of neutron stars. Each contribution of the neutrino form factors to the neutrino mean free path is discussed.


title: Randomized Reduced Forward Models for Efficient Metropolis--Hastings MCMC, with Application to Subsurface Fluid Flow and Capacitance Tomography

abstract: Bayesian modelling and computational inference by Markov chain Monte Carlo (MCMC) is a principled framework for large-scale uncertainty quantification, though is limited in practice by computational cost when implemented in the simplest form that requires simulating an accurate computer model at each iteration of the MCMC. The delayed acceptance Metropolis--Hastings MCMC leverages a reduced model for the forward map to lower the compute cost per iteration, though necessarily reduces statistical efficiency that can, without care, lead to no reduction in the computational cost of computing estimates to a desired accuracy. Randomizing the reduced model for the forward map can dramatically improve computational efficiency, by maintaining the low cost per iteration but also avoiding appreciable loss of statistical efficiency. Randomized maps are constructed by a posteriori adaptive tuning of a randomized and locally-corrected deterministic reduced model. Equivalently, the approximated posterior distribution may be viewed as induced by a modified likelihood function for use with the reduced map, with parameters tuned to optimize the quality of the approximation to the correct posterior distribution. Conditions for adaptive MCMC algorithms allow practical approximations and algorithms that have guaranteed ergodicity for the target distribution. Good statistical and computational efficiencies are demonstrated in examples of calibration of large-scale numerical models of geothermal reservoirs and electrical capacitance tomography.


title: Blowing up solutions for supercritical Yamabe problems on manifolds with non umbilic boundary

abstract: We build blowing-up solutions for a supercritical perturbation of the Yamabe problem on manifolds with boundary, provided the dimension of the manifold is n>6 and the trace-free part of the second fundamental form is non-zero everywhere on the boundary.


title: Modelling of limitations of bulk heterojunction architecture in organic solar cells II: 3d model

abstract: Polymer solar cells are considered as very promising candidates for development of photovoltaics of the future. They are cheap and easy to fabricate. However, up to now, they possess fundamental drawback: low effectiveness. In the most popular BHJ (bulk heterojunction) architecture the actual long-standing top efficiency is about 12\% (recent achievements about 15\%). One ask the question how fundamental this limitation is, as certain theoretical considerations suggest that it should be about two times higher. In our paper we analyze the `geometric factor' as one of possible explanation of relatively low efficiency of BHJ architecture. More precisly, we calculate the effective area of the donor-acceptor border in the random mixture of donor and acceptor nanocrystals and further compare it with an ideal 'brush' architecture. In our previous calculation for the two dimensional model, we have found that the maximal value of geometric factor was about 40\%. In the actual three dimensional model, it turned out that both architectures give very close value of the effective area. So the geometric factor seems to be not significant as a factor limiting efficiency. Implications of this fact are discussed: we list two other factors (mentioned but not thoroughly discussed in literature) which can be responsible for limitations of efficiency of BHJ architecture. We estimate their scale, and suggest that these limitations are inevitable, or at least very hard to overcome. We suggest that return to layer architecture could radically improve efficiency limitations -- however, to make breakthrough, materials with large exciton diffusion length have to be invented.


title: Hybrid Digital-Analog Beamforming and MIMO Radar with OTFS Modulation

abstract: Motivated by future automotive applications, we study some joint radar target detection and parameter estimation problems where the transmitter, equipped with a mono-static MIMO radar, wishes to detect multiple targets and then estimate their respective parameters, while simultaneously communicating information data using orthogonal time frequency space (OTFS) modulation. Assuming that the number of radio frequency chains is smaller than the number of antennas over the mmWave frequency band, we design hybrid digital-analog beamforming at the radar transmitter adapted to different operating phases. The first scenario considers a wide angular beam in order to perform the target detection and parameter estimation, while multicasting a common message to all possible active users. The second scenario considers narrow angular beams to send information streams individually to the already detected users and simultaneously keep tracking of their respective parameters. Under this setup, we propose an efficient maximum likelihood scheme combined with hybrid beamforming to jointly perform target detection and parameter estimation. Our numerical results demonstrate that the proposed algorithm is able to reliably detect multiple targets with a sufficient number of antennas and achieves the Cram\'er-Rao lower bound for radar parameter estimation such as delay, Doppler and angle-of-arrival (AoA).


title: Automated Extended Aperture Photometry for K2 RR Lyrae stars

abstract: Light curves for RR Lyrae stars can be difficult to obtain properly in the K2 mission due to the similarities between the timescales of the observed physical phenomena and the instrumental signals appearing in the data. We developed a new photometric method called Extended Aperture Photometry (EAP), a key element of which is to extend the aperture to an optimal size to compensate for the motion of the telescope and to collect all available flux from the star before applying further corrections. We determined the extended apertures for individual stars by hand so far. Now we managed to automate the pipeline that we intend to use for the nearly four thousand RR Lyrae targets observed in the K2 mission. We present the outline of our pipeline and make some comparisons to other photometric solutions.


title: Determining Number of Kneser Graphs: Exact Values and Improved Bounds

abstract: The determining number of a graph $G = (V,E)$ is the minimum cardinality of a set $S\subseteq V$ such that pointwise stabilizer of $S$ under the action of $Aut(G)$ is trivial. In this paper, we provide some improved upper and lower bounds on the determining number of Kneser graphs. Moreover, we provide the exact value of the determining number for some subfamilies of Kneser graphs.


title: Extremely precise age and metallicity of the open cluster NGC 2506 using detached eclipsing binaries

abstract: Accurate stellar parameters of stars in open clusters can help constrain models of stellar structure and evolution. Here we wish to determine the age and metallicity content of the open cluster NGC 2506. To this end we investigated three detached eclipsing binaries (DEBs; V2032, V4, and V5) for which we determined their masses and radii, as well as four red giant branch stars for which we determined their effective temperatures, surface gravities, and metallicities. Three of the stars in the DEBs have masses close to the cluster turn-off mass, allowing for extremely precise age determination. Comparing the values for the masses and radii of the binaries to BaSTI isochrones we estimated a cluster age of $2.01 \pm 0.10$ Gyr. This does depend on the models used in the comparison, where we have found that the inclusion of convective core-overshooting is necessary to properly model the cluster. From red giant branch stars we determined values for the effective temperatures, the surface gravities, and the metallicities. From these we find a cluster metallicity of $-0.36 \pm 0.10$ dex. Using this value and the values for the effective temperatures we determine the reddening to be E$(b - y) = 0.057 \pm 0.004$ mag. Furthermore, we derived the distance to the cluster from Gaia parallaxes and found $3.101 \pm 0.017$ kpc, and we have performed a radial velocity membership determination for stars in the field of the cluster. Finally, we report on the detection of oscillation signals in $\gamma$ Dor and $\delta$ Scuti members in data from the TESS mission, including the possible detection of solar-like oscillations in two of the red giants.


title: Additive Models for Symmetric Positive-Definite Matrices, Riemannian Manifolds and Lie groups

abstract: In this paper an additive regression model for a symmetric positive-definite matrix valued response and multiple scalar predictors is proposed. The model exploits the abelian group structure inherited from either the Log-Cholesky metric or the Log-Euclidean framework that turns the space of symmetric positive-definite matrices into a Riemannian manifold and further a bi-invariant Lie group. The additive model for responses in the space of symmetric positive-definite matrices with either of these metrics is shown to connect to an additive model on a tangent space. This connection not only entails an efficient algorithm to estimate the component functions but also allows to generalize the proposed additive model to general Riemannian manifolds that might not have a Lie group structure. Optimal asymptotic convergence rates and normality of the estimated component functions are also established. Numerical studies show that the proposed model enjoys superior numerical performance, especially when there are multiple predictors. The practical merits of the proposed model are demonstrated by analyzing diffusion tensor brain imaging data.


title: Cough Against COVID: Evidence of COVID-19 Signature in Cough Sounds

abstract: Testing capacity for COVID-19 remains a challenge globally due to the lack of adequate supplies, trained personnel, and sample-processing equipment. These problems are even more acute in rural and underdeveloped regions. We demonstrate that solicited-cough sounds collected over a phone, when analysed by our AI model, have statistically significant signal indicative of COVID-19 status (AUC 0.72, t-test,p <0.01,95% CI 0.61-0.83). This holds true for asymptomatic patients as well. Towards this, we collect the largest known(to date) dataset of microbiologically confirmed COVID-19 cough sounds from 3,621 individuals. When used in a triaging step within an overall testing protocol, by enabling risk-stratification of individuals before confirmatory tests, our tool can increase the testing capacity of a healthcare system by 43% at disease prevalence of 5%, without additional supplies, trained personnel, or physical infrastructure


title: Multi-source Data Mining for e-Learning

abstract: Data mining is the task of discovering interesting, unexpected or valuable structures in large datasets and transforming them into an understandable structure for further use . Different approaches in the domain of data mining have been proposed, among which pattern mining is the most important one. Pattern mining mining involves extracting interesting frequent patterns from data. Pattern mining has grown to be a topic of high interest where it is used for different purposes, for example, recommendations. Some of the most common challenges in this domain include reducing the complexity of the process and avoiding the redundancy within the patterns. So far, pattern mining has mainly focused on the mining of a single data source. However, with the increase in the amount of data, in terms of volume, diversity of sources and nature of data, mining multi-source and heterogeneous data has become an emerging challenge in this domain. This challenge is the main focus of our work where we propose to mine multi-source data in order to extract interesting frequent patterns.


title: Commands 4 Autonomous Vehicles (C4AV) Workshop Summary

abstract: The task of visual grounding requires locating the most relevant region or object in an image, given a natural language query. So far, progress on this task was mostly measured on curated datasets, which are not always representative of human spoken language. In this work, we deviate from recent, popular task settings and consider the problem under an autonomous vehicle scenario. In particular, we consider a situation where passengers can give free-form natural language commands to a vehicle which can be associated with an object in the street scene. To stimulate research on this topic, we have organized the \emph{Commands for Autonomous Vehicles} (C4AV) challenge based on the recent \emph{Talk2Car} dataset (URL: https://www.aicrowd.com/challenges/eccv-2020-commands-4-autonomous-vehicles). This paper presents the results of the challenge. First, we compare the used benchmark against existing datasets for visual grounding. Second, we identify the aspects that render top-performing models successful, and relate them to existing state-of-the-art models for visual grounding, in addition to detecting potential failure cases by evaluating on carefully selected subsets. Finally, we discuss several possibilities for future work.


title: GuessTheMusic: Song Identification from Electroencephalography response

abstract: The music signal comprises of different features like rhythm, timbre, melody, harmony. Its impact on the human brain has been an exciting research topic for the past several decades. Electroencephalography (EEG) signal enables non-invasive measurement of brain activity. Leveraging the recent advancements in deep learning, we proposed a novel approach for song identification using a Convolution Neural network given the electroencephalography (EEG) responses. We recorded the EEG signals from a group of 20 participants while listening to a set of 12 song clips, each of approximately 2 minutes, that were presented in random order. The repeating nature of Music is captured by a data slicing approach considering brain signals of 1 second duration as representative of each song clip. More specifically, we predict the song corresponding to one second of EEG data when given as input rather than a complete two-minute response. We have also discussed pre-processing steps to handle large dimensions of a dataset and various CNN architectures. For all the experiments, we have considered each participant's EEG response for each song in both train and test data. We have obtained 84.96\% accuracy for the same. The performance observed gives appropriate implication towards the notion that listening to a song creates specific patterns in the brain, and these patterns vary from person to person.


title: Simplicial persistence of financial markets: filtering, generative processes and portfolio risk

abstract: We introduce simplicial persistence, a measure of time evolution of network motifs in subsequent temporal layers. We observe long memory in the evolution of structures from correlation filtering, with a two regime power law decay in the number of persistent simplicial complexes. Null models of the underlying time series are tested to investigate properties of the generative process and its evolutional constraints. Networks are generated with both TMFG filtering technique and thresholding showing that embedding-based filtering methods (TMFG) are able to identify higher order structures throughout the market sample, where thresholding methods fail. The decay exponents of these long memory processes are used to characterise financial markets based on their stage of development and liquidity. We find that more liquid markets tend to have a slower persistence decay. This is in contrast with the common understanding that developed markets are more random. We find that they are indeed less predictable for what concerns the dynamics of each single variable but they are more predictable for what concerns the collective evolution of the variables. This could imply higher fragility to systemic shocks.

categories: q-fin.ST physics.soc-ph q-fin.PM q-fin.RM


title: Point Forces in Elasticity Equation and Their Alternatives in Multi Dimensions

abstract: We consider several mathematical issues regarding models that simulate forces exerted by cells. Since the size of cells is much smaller than the size of the domain of computation, one often considers point forces, modelled by Dirac Delta distributions on boundary segments of cells. In the current paper, we treat forces that are directed normal to the cell boundary and that are directed toward the cell centre. Since it can be shown that there exists no smooth solution, at least not in $H^1$ for solutions to the governing momentum balance equation, we analyse the convergence and quality of the approximation. Furthermore, the expected finite element problems that we get necessitate scrutinizing alternative model formulations, such as the use of smoothed Dirac distributions, or the so-called smoothed particle approach as well as the so-called 'hole' approach where cellular forces are modelled through the use of (natural) boundary conditions. In this paper, we investigate and attempt to quantify the conditions for consistency between the various approaches. This has resulted in error analyses in the $H^1$-norm of the numerical solution based on Galerkin principles that entail Lagrangian basis functions. The paper also addresses well-posedness in terms of existence and uniqueness. The current analysis has been performed for the linear steady-state (hence neglecting inertia and damping) momentum equations under the assumption of Hooke's law.


title: $\sigma^2$R Loss: a Weighted Loss by Multiplicative Factors using Sigmoidal Functions

abstract: In neural networks, the loss function represents the core of the learning process that leads the optimizer to an approximation of the optimal convergence error. Convolutional neural networks (CNN) use the loss function as a supervisory signal to train a deep model and contribute significantly to achieving the state of the art in some fields of artificial vision. Cross-entropy and Center loss functions are commonly used to increase the discriminating power of learned functions and increase the generalization performance of the model. Center loss minimizes the class intra-class variance and at the same time penalizes the long distance between the deep features inside each class. However, the total error of the center loss will be heavily influenced by the majority of the instances and can lead to a freezing state in terms of intra-class variance. To address this, we introduce a new loss function called sigma squared reduction loss ($\sigma^2$R loss), which is regulated by a sigmoid function to inflate/deflate the error per instance and then continue to reduce the intra-class variance. Our loss has clear intuition and geometric interpretation, furthermore, we demonstrate by experiments the effectiveness of our proposal on several benchmark datasets showing the intra-class variance reduction and overcoming the results obtained with center loss and soft nearest neighbour functions.


title: A Finite Elements Strategy for Spread Contract Valuation Via Associated PIDE

abstract: We study an efficient strategy based on finite elements to value spread options on commodities whose underlying assets follow a dynamic described by a certain class of two-dimensional Levy models by solving their associated partial integro-differential equation (PIDE). To this end we consider a Galerkin approximation in space along with an implicit scheme for time evolution. Diffusion and drift in the associated operator are discretized using an exact Gaussian quadrature, while the integral part corresponding to jumps is approximated using the symbol method recently introduced in the literature. A system with blocked Toeplitz with Toeplitz blocks (BTTB) matrix is efficiently solved via biconjugate stabilized gradient method (BICSTAB) with a circulant pre-conditioner at each time step. The technique is applied to the pricing of \textit{crack} spread options between the prices of futures RBOB gasoline (reformulated blendstock for oxygenate blending) and West Texas Intermediate(WTI) oil in NYMEX.


title: Designing Compact Features for Remote Stroke Rehabilitation Monitoring using Wearable Accelerometers

abstract: Stroke is known as a major global health problem, and for stroke survivors it is key to monitor the recovery levels. However, traditional stroke rehabilitation assessment methods (such as the popular clinical assessment) can be subjective and expensive, and it is also less convenient for patients to visit clinics in a high frequency. To address this issue, in this work based on wearable sensing and machine learning techniques, we develop an automated system that can predict the assessment score in an objective manner. With wrist-worn sensors, accelerometer data is collected from 59 stroke survivors in free-living environments for a duration of 8 weeks, and we map the week-wise accelerometer data(3 days per week) to the assessment score by developing signal processing and predictive model pipeline. To achieve this, we propose two types of new features, which can encode the rehabilitation information from both paralysed and non-paralysed sides while suppressing the high level noises such as irrelevant daily activities. Based on the proposed features, we further develop the longitudinal mixed-effects model with Gaussian process prior (LMGP), which can model the random effects caused by different subjects and time slots (during the 8 weeks). Comprehensive experiments are conducted to evaluate our system on both acute and chronic patients, and the promising results suggest its effectiveness.


title: Formation of the cosmic-ray halo: Galactic spectrum of primary cosmic rays

abstract: A self-consistent model of a one-dimensional cosmic-ray (CR) halo around the Galactic disk is formulated with the restriction to a minimum number of free parameters. It is demonstrated that the turbulent cascade of MHD waves does not necessarily play an essential role in the halo formation. Instead, an increase of the Alfven velocity with distance to the disk leads to an efficient generic mechanism of the turbulent redshift, enhancing CR scattering by the self-generated MHD waves. As a result, the calculated size of the CR halo at lower energies is determined by the halo sheath, an energy-dependent region around the disk beyond which the CR escape becomes purely advective. At sufficiently high energies, the halo size is set by the characteristic thickness of the ionized gas distribution. The calculated Galactic spectrum of protons shows a remarkable agreement with observations, reproducing the position of spectral break at ~ 0.6 TeV and the spectral shape up to ~ 10 TeV.


title: Rainbow Spacetime from a Nonlocal Gravitational Uncertainty Principle

abstract: Occurrence of spacetime singularities is one of the peculiar features of Einstein gravity, signalling limitation on probing short distances in spacetime. This alludes to the existence of a fundamental length scale in nature. On contrary, Heisenberg quantum uncertainty relation seems to allow for probing arbitrarily small length scales. To reconcile these two conflicting ideas in line with a well known framework of quantum gravity, several modifications of Heisenberg algebra have been proposed. However, it has been extensively argued that such a minimum length would introduce nonlocality in theories of quantum gravity. In this Letter, we analyze a previously proposed deformation of the Heisenberg algebra (i.e. $p \rightarrow p (1 + \lambda p^{-1})$) for a particle confined in a box subjected to a gravitational field. For the problem in hand, such deformation seems to yield an energy-dependent behavior of spacetime in a way consistent with gravity's rainbow, hence demonstrating a connection between non-locality and gravity's rainbow.


title: SciBERT-based Semantification of Bioassays in the Open Research Knowledge Graph

abstract: As a novel contribution to the problem of semantifying biological assays, in this paper, we propose a neural-network-based approach to automatically semantify, thereby structure, unstructured bioassay text descriptions. Experimental evaluations, to this end, show promise as the neural-based semantification significantly outperforms a naive frequency-based baseline approach. Specifically, the neural method attains 72% F1 versus 47% F1 from the frequency-based method.


title: A computationally efficient physiologically comprehensive 3D-0D closed-loop model of the heart and circulation

abstract: Computer models of cardiac electro-mechanics (EM) show promise as an effective means for quantitative analysis of clinical data and, potentially, for predicting therapeutic responses.realize such advanced applications methodological key challenges must be addressed. Enhanced computational efficiency and robustness is crucial to facilitate, within tractable time frames, model personalization, the simulation of prolonged observation periods under a broad range of conditions, and physiological completeness encompassing therapy-relevant mechanisms is needed to endow models with predictive capabilities beyond the mere replication of observations. Here, we introduce a universal feature-complete cardiac EM modeling framework that builds on a flexible method for coupling a 3D model of bi-ventricular EM to the physiologically comprehensive 0D CircAdapt model representing atrial mechanics and closed-loop circulation. A detailed mathematical description is given and efficiency, robustness, and accuracy of numerical scheme and solver implementation are evaluated. After parameterization and stabilization of the coupled 3D-0D model to a limit cycle under baseline conditions, the model's ability to replicate physiological behaviors is demonstrated, by simulating the transient response to alterations in loading conditions and contractility, as induced by experimental protocols used for assessing systolic and diastolic ventricular properties. Mechanistic completeness and computational efficiency of this novel model render advanced applications geared towards predicting acute outcomes of EM therapies feasible.


title: Differentiation of the Wright functions with respect to parameters and other results

abstract: In this survey we discuss derivatives of the Wright functions (of the first and the second kind) with respect to parameters. Differentiation of these functions leads to infinite power series with coefficient being quotients of the digamma (psi) and gamma functions. Only in few cases it is possible to obtain the sums of these series in a closed form. Functional form of the power series resembles those derived for the Mittag-Leffler functions. If the Wright functions are treated as the generalized Bessel functions, differentiation operations can be expressed in terms of the Bessel functions and their derivatives with respect to the order. It is demonstrated that in many cases it is possible to derive the explicit form of the Mittag-Leffler functions by performing simple operations with the Laplace transforms of the Wright functions. The Laplace transform pairs of the both kinds of the Wright functions are discussed for particular values of the parameters. Some transform pairs serve to obtain functional limits by applying the shifted Dirac delta function.


title: Improving the spatial resolution of a BOTDA sensor using deconvolution algorithm

abstract: Spatial resolution improvement from an acquired measurement using long pulse is developed for Brillouin optical time domain analysis (BOTDA) systems based on the total variation deconvolution algorithm. The frequency dependency of Brillouin gain temporal envelope is investigated by simulation, and its impact on the recovered results of deconvolution algorithm is thoroughly analyzed. To implement a reliable deconvolution process, differential pulse-width pair (DPP) technique is utilized to effectively eliminate the systematic BFS distortion stemming from the frequency dependency of temporal envelope. The width of the pulse pairs should be larger than 40 ns as is analyzed theoretically and verified experimentally. It has been demonstrated that the proposed method can realize flexible adjustment of spatial resolution with enhanced signal-to-noise ratio (SNR) from an established measurement with long pump pulse. In the experiment, the spatial resolution is increased to 0.5 m and 1 m with high measurement accuracy by using the deconvolution algorithm from the measurement of 60/40 ns DPP signals. Compared with the raw DPP results with the same spatial resolution, 9.2 dB and 8.4 dB SNR improvements are obtained for 0.5 m and 1 m spatial resolution respectively, thanks to the denoising capability of the total variation deconvolution algorithm. The impact of sampling rate on the recovery results is also studied. The proposed sensing system allows for distortion-free Brillouin distributed sensing with higher spatial resolution and enhanced SNR from the conventional DPP setup with long pulse pairs.


title: HDGlab: An open-source implementation of the hybridisable discontinuous Galerkin method in MATLAB

abstract: This paper presents HDGlab, an open source MATLAB implementation of the hybridisable discontinuous Galerkin (HDG) method. The main goal is to provide a detailed description of both the HDG method for elliptic problems and its implementation available in HDGlab. Ultimately, this is expected to make this relatively new advanced discretisation method more accessible to the computational engineering community. HDGlab presents some features not available in other implementations of the HDG method that can be found in the free domain. First, it implements high-order polynomial shape functions up to degree nine, with both equally-spaced and Fekete nodal distributions. Second, it supports curved isoparametric simplicial elements in two and three dimensions. Third, it supports non-uniform degree polynomial approximations and it provides a flexible structure to devise degree adaptivity strategies. Finally, an interface with the open-source high-order mesh generator Gmsh is provided to facilitate its application to practical engineering problems.


title: Blocking total dominating sets via edge contractions

abstract: In this paper, we study the problem of deciding whether the total domination number of a given graph $G$ can be reduced using exactly one edge contraction (called 1-Edge Contraction($\gamma_t$)). We focus on several graph classes and determine the computational complexity of this problem. By putting together these results, we manage to obtain a complete dichotomy for $H$-free graphs.


title: Monte Carlo Tree Search Based Tactical Maneuvering

abstract: In this paper we explore the application of simultaneous move Monte Carlo Tree Search (MCTS) based online framework for tactical maneuvering between two unmanned aircrafts. Compared to other techniques, MCTS enables efficient search over long horizons and uses self-play to select best maneuver in the current state while accounting for the opponent aircraft tactics. We explore different algorithmic choices in MCTS and demonstrate the framework numerically in a simulated 2D tactical maneuvering application.


title: Acoustic radiation force and torque on spheroidal particles in an ideal cylindrical chamber

abstract: We theoretically investigate how the acoustic radiation force and torque arise on a small spheroidal particle immersed in a nonviscous fluid inside an ideal cylindrical chamber. The ideal chamber comprises a hard top and bottom (rigid boundary condition), and a soft or hard lateral wall. By assuming the particle is much smaller than the acoustic wavelength, we present analytical expressions of the radiation force and torque caused by an acoustic wave of arbitrary shape. Unlike previous results, these expressions are given relative to a fixed laboratory frame. Our model is showcased for analyzing the behavior of an elongated metallic microspheroid (with a 10 : 1 aspect ratio) in a half-wavelength acoustofluidic chamber with a few millimeters diameter. The results show the radiation torque aligns the microspheroid along the nodal plane, and the radiation force causes a translational motion with a speed of up to one body length per second. At last, we discuss the implications of this study to propelled nanorods by ultrasound.


title: Deep learning denoising for EOG artifacts removal from EEG signals

abstract: There are many sources of interference encountered in the electroencephalogram (EEG) recordings, specifically ocular, muscular, and cardiac artifacts. Rejection of EEG artifacts is an essential process in EEG analysis since such artifacts cause many problems in EEG signals analysis. One of the most challenging issues in EEG denoising processes is removing the ocular artifacts where Electrooculographic (EOG), and EEG signals have an overlap in both frequency and time domains. In this paper, we build and train a deep learning model to deal with this challenge and remove the ocular artifacts effectively. In the proposed scheme, we convert each EEG signal to an image to be fed to a U-NET model, which is a deep learning model usually used in image segmentation tasks. We proposed three different schemes and made our U-NET based models learn to purify contaminated EEG signals similar to the process used in the image segmentation process. The results confirm that one of our schemes can achieve a reliable and promising accuracy to reduce the Mean square error between the target signal (Pure EEGs) and the predicted signal (Purified EEGs).


title: Automatic Differentiation to Simultaneously Identify Nonlinear Dynamics and Extract Noise Probability Distributions from Data

abstract: The sparse identification of nonlinear dynamics (SINDy) is a regression framework for the discovery of parsimonious dynamic models and governing equations from time-series data. As with all system identification methods, noisy measurements compromise the accuracy and robustness of the model discovery procedure. In this work, we develop a variant of the SINDy algorithm that integrates automatic differentiation and recent time-stepping constrained motivated by Rudy et al. for simultaneously (i) denoising the data, (ii) learning and parametrizing the noise probability distribution, and (iii) identifying the underlying parsimonious dynamical system responsible for generating the time-series data. Thus within an integrated optimization framework, noise can be separated from signal, resulting in an architecture that is approximately twice as robust to noise as state-of-the-art methods, handling as much as 40% noise on a given time-series signal and explicitly parametrizing the noise probability distribution. We demonstrate this approach on several numerical examples, from Lotka-Volterra models to the spatio-temporal Lorenz 96 model. Further, we show the method can identify a diversity of probability distributions including Gaussian, uniform, Gamma, and Rayleigh.


title: Disordered complex networks: energy optimal lattices and persistent homology

abstract: Disordered complex networks are of fundamental interest as stochastic models for information transmission over wireless networks. Well-known networks based on the Poisson point process model have limitations vis-a-vis network efficiency, whereas strongly correlated alternatives, such as those based on random matrix spectra (RMT), have tractability and robustness issues. In this work, we demonstrate that network models based on random perturbations of Euclidean lattices interpolate between Poisson and rigidly structured networks, and allow us to achieve the best of both worlds : significantly improve upon the Poisson model in terms of network efficacy measured by the Signal to Interference plus Noise Ratio (abbrv. SINR) and the related concept of coverage probabilities, at the same time retaining a considerable measure of mathematical and computational simplicity and robustness to erasure and noise. We investigate the optimal choice of the base lattice in this model, connecting it to the celebrated problem optimality of Euclidean lattices with respect to the Epstein Zeta function, which is in turn related to notions of lattice energy. This leads us to the choice of the triangular lattice in 2D and face centered cubic lattice in 3D. We demonstrate that the coverage probability decreases with increasing strength of perturbation, eventually converging to that of the Poisson network. In the regime of low disorder, we approximately characterize the statistical law of the coverage function. In 2D, we determine the disorder strength at which the PTL and the RMT networks are the closest measured by comparing their network topologies via a comparison of their Persistence Diagrams . We demonstrate that the PTL network at this disorder strength can be taken to be an effective substitute for the RMT network model, while at the same time offering the advantages of greater tractability.

categories: eess.SP cond-mat.dis-nn cs.IT cs.LG math.IT math.PR


title: Symbol functions for symmetric frameworks

abstract: We prove a variant of the well-known result that intertwiners for the bilateral shift on `$\ell^2(Z)$ are unitarily equivalent to multiplication operators on $L^2(T)$. This enables us to unify and extend fundamental aspects of rigidity theory for bar-joint frameworks with an abelian symmetry group. In particular, we formulate the symbol function for a wide class of frameworks and show how to construct generalised rigid unit modes in a variety of new contexts.


title: Resolution Enhancement in Protein NMR Spectra by Deconvolution with Compressed Sensing Reconstruction

abstract: Multidimensional NMR spectroscopy is one of the basic tools for determining the structure of biomolecules. Unfortunately, the resolution of the spectra is often limited by inter-nuclear couplings. This limitation cannot be overcome by common ways of increasing resolution, i.e. non-uniform sampling (NUS) followed by compressed sensing (CS) reconstruction. In this paper, we show how to enrich CS processing with virtual decoupling leading to an increase in resolution, sensitivity, and overall quality of NUS reconstruction. A mathematical description of the decoupling by deconvolution approach explains the effects of noise, modulation of the sampling schedule, and reveals relation with the underlying assumption of the CS. The gain in resolution and sensitivity is demonstrated for the basic experiment used for protein backbone assignment 3D HNCA applied to two large protein systems: intrinsically disordered 441-residue Tau and a 509-residue globular bacteriophytochrome fragment.


title: Short dated smile under Rough Volatility: asymptotics and numerics

abstract: In [Precise Asymptotics for Robust Stochastic Volatility Models; Ann. Appl. Probab. 2021] we introduce a new methodology to analyze large classes of (classical and rough) stochastic volatility models, with special regard to short-time and small noise formulae for option prices, using the framework [Bayer et al; A regularity structure for rough volatility; Math. Fin. 2020]. We investigate here the fine structure of this expansion in large deviations and moderate deviations regimes, together with consequences for implied volatility. We discuss computational aspects relevant for the practical application of these formulas. We specialize such expansions to prototypical rough volatility examples and discuss numerical evidence.


title: Surface Josephson plasma waves in a high-temperature superconductor

abstract: Electron density oscillations with acoustic dispersions and sustained at boundaries between different media provide information about surface and interface properties of hetero-structures. In ultra-thin metallic films these plasmonic excitations are heavily damped. Superconductivity is predicted to reduce dissipation allowing detection of these resonances. Emerging low-loss interface Cooper-pair waves have been studied before, however, the observation of surface-confined Josephson plasmons has remained elusive. Here, we report on generation and coupling to these excitations in an ultrathin single-crystal film of high-temperature superconductor La1.85Sr0.15CuO4. The film becomes brighter than Au below the critical temperature when probed with sub-gap THz photons. We show that the enhanced signal in the superconducting state, which can be visualized with a spatial resolution better than {\lambda}/3,000, originates from near-field coupling of light to surface Josephson plasmons. Our results open a path towards non-invasive investigation of enhanced superconductivity in artificial multilayers, buried interface states in topological hetero-structures, and non-linear phenomena in Josephson devices.


title: Improved Coding over Sets for DNA-Based Data Storage

abstract: Error-correcting codes over sets, with applications to DNA storage, are studied. The DNA-storage channel receives a set of sequences, and produces a corrupted version of the set, including sequence loss, symbol substitution, symbol insertion/deletion, and limited-magnitude errors in symbols. Various parameter regimes are studied. New bounds on code parameters are provided, which improve upon known bounds. New codes are constructed, at times matching the bounds up to lower-or der terms or small constant factors.


title: Calmness of the solution mapping of Navier-Stokes problems modeled by hemivariational inequalities

abstract: The main purpose of this paper is to find conditions for Holder calmness of the solution mapping, viewed as a function of the boundary data, of a hemivariational inequality governed by the Navier-Stokes operator. To this end, a more abstract model is studied first: a class of parametric equilibrium problems defined by trifunctions. The presence of trifunctions allows the extension of the monotonicity notions and of the duality principle in the theory of equilibrium problems.


title: Resolution of a bent-crystal spectrometer for X-ray free electron laser pulses: diamond vs. silicon

abstract: The resolution function of a spectrometer based on a strongly bent single crystal (bending radius of 10 cm or less) is evaluated. It is shown that the resolution is controlled by two parameters, (i) the ratio of the lattice spacing of the chosen reflection to the crystal thickness and (ii) a single parameter comprising crystal thickness, its bending radius, and anisotropic elastic constants of the chosen crystal. Diamond, due to its unique elastic properties, can provide notably higher resolution than silicon. The results allow to optimize the parameters of bent crystal spectrometers for the hard X-ray free electron laser sources.


title: Real-Time Optimization Meets Bayesian Optimization and Derivative-Free Optimization: A Tale of Modifier Adaptation

abstract: This paper investigates a new class of modifier-adaptation schemes to overcome plant-model mismatch in real-time optimization of uncertain processes. The main contribution lies in the integration of concepts from the areas of Bayesian optimization and derivative-free optimization. The proposed schemes embed a physical model and rely on trust-region ideas to minimize risk during the exploration, while employing Gaussian process regression to capture the plant-model mismatch in a non-parametric way and drive the exploration by means of acquisition functions. The benefits of using an acquisition function, knowing the process noise level, or specifying a nominal process model are illustrated on numerical case studies, including a semi-batch photobioreactor optimization problem.


title: FarsTail: A Persian Natural Language Inference Dataset

abstract: Natural language inference (NLI) is known as one of the central tasks in natural language processing (NLP) which encapsulates many fundamental aspects of language understanding. With the considerable achievements of data-hungry deep learning methods in NLP tasks, a great amount of effort has been devoted to develop more diverse datasets for different languages. In this paper, we present a new dataset for the NLI task in the Persian language, also known as Farsi, which is one of the dominant languages in the Middle East. This dataset, named FarsTail, includes 10,367 samples which are provided in both the Persian language as well as the indexed format to be useful for non-Persian researchers. The samples are generated from 3,539 multiple-choice questions with the least amount of annotator interventions in a way similar to the SciTail dataset. A carefully designed multi-step process is adopted to ensure the quality of the dataset. We also present the results of traditional and state-of-the-art methods on FarsTail including different embedding methods such as word2vec, fastText, ELMo, BERT, and LASER, as well as different modeling approaches such as DecompAtt, ESIM, HBMP, and ULMFiT to provide a solid baseline for the future research. The best obtained test accuracy is 83.38% which shows that there is a big room for improving the current methods to be useful for real-world NLP applications in different languages. We also investigate the extent to which the models exploit superficial clues, also known as dataset biases, in FarsTail, and partition the test set into easy and hard subsets according to the success of biased models. The dataset is available at https://github.com/dml-qom/FarsTail


title: A bounded operator approach to technical indicators without lag

abstract: In the framework of technical analysis for algorithmic trading we use a linear algebra approach in order to define classical technical indicators as bounded operators of the space $l^\infty(\mathbb{N})$. This more abstract view enables us to define in a very simple way the no-lag versions of these tools. Then we apply our results to a basic trading system in order to compare the classical Elder's impulse system with its no-lag version and the so-called Nyquist-Elder's impulse system.


title: PSZSPT: a joint Planck and SPT-SZ cluster catalog

abstract: We present the first cluster catalog extracted from combined space-based (Planck) and ground-based (South Pole Telescope; SPT-SZ) millimeter data. We developed and applied a matched multi-filter (MMF) capable of dealing with the different transfer functions and resolutions of the two datasets. We verified that it produces results consistent with publications from Planck and SPT collaborations when applied on the datasets individually. We also verified that Planck and SPT-SZ cluster fluxes are consistent with each other. When applied blindly to the combined dataset, the MMF generated a catalog of 419 detections ($S/N>5$), of which 323 are already part of the SPT-SZ or PSZ2 catalogs; 54 are new SZ detections, which have been identified in other catalogs or surveys; and 42 are new unidentified candidates. The MMF takes advantage of the complementarity of the two datasets, Planck being particularly useful for detecting clusters at a low redshift ($z<0.3$), while SPT is efficient at finding higher redshift ($z>0.3$) sources. This work represents a proof of concept that blind cluster extraction can be performed on combined, inhomogeneous millimeter datasets acquired from space and ground. This result is of prime importance for planned ground-based cosmic microwave background (CMB) experiments (e.g., Simons Observatory, CMB-S4) and envisaged CMB space missions (e.g., PICO, Backlight) that will detect hundreds of thousands of clusters in the low mass regime ($M_{500} \leqslant 10^{14} M_\odot$), for which the various sources of intra-cluster emission (gas, dust, synchrotron) will be of the same order of magnitude and hence require broad ground and space frequency coverage with a comparable spatial resolution for adequate separation.


title: Equivalence of three classical algorithms with quantum side information: Privacy amplification, error correction, and data compression

abstract: Privacy amplification (PA) is an indispensable component in classical and quantum cryptography. Error correction (EC) and data compression (DC) algorithms are also indispensable in classical and quantum information theory. We here study these three algorithms (PA, EC, and DC) in the presence of quantum side information, and show that they all become equivalent in the one-shot scenario. As an application of this equivalence, we take previously known security bounds of PA, and translate them into coding theorems for EC and DC which have not been obtained previously. Further, we apply these results to simplify and improve our previous result that the two prevalent approaches to the security proof of quantum key distribution (QKD) are equivalent. We also propose a new method to simplify the security proof of QKD.


title: Pedestrian Motion Tracking by Using Inertial Sensors on the Smartphone

abstract: Inertial Measurement Unit (IMU) has long been a dream for stable and reliable motion estimation, especially in indoor environments where GPS strength limits. In this paper, we propose a novel method for position and orientation estimation of a moving object only from a sequence of IMU signals collected from the phone. Our main observation is that human motion is monotonous and periodic. We adopt the Extended Kalman Filter and use the learning-based method to dynamically update the measurement noise of the filter. Our pedestrian motion tracking system intends to accurately estimate planar position, velocity, heading direction without restricting the phone's daily use. The method is not only tested on the self-collected signals, but also provides accurate position and velocity estimations on the public RIDI dataset, i.e., the absolute transmit error is 1.28m for a 59-second sequence.


title: Densely Guided Knowledge Distillation using Multiple Teacher Assistants

abstract: With the success of deep neural networks, knowledge distillation which guides the learning of a small student network from a large teacher network is being actively studied for model compression and transfer learning. However, few studies have been performed to resolve the poor learning issue of the student network when the student and teacher model sizes significantly differ. In this paper, we propose a densely guided knowledge distillation using multiple teacher assistants that gradually decreases the model size to efficiently bridge the large gap between the teacher and student networks. To stimulate more efficient learning of the student network, we guide each teacher assistant to every other smaller teacher assistants iteratively. Specifically, when teaching a smaller teacher assistant at the next step, the existing larger teacher assistants from the previous step are used as well as the teacher network. Moreover, we design stochastic teaching where, for each mini-batch, a teacher or teacher assistants are randomly dropped. This acts as a regularizer to improve the efficiency of teaching of the student network. Thus, the student can always learn salient distilled knowledge from the multiple sources. We verified the effectiveness of the proposed method for a classification task using CIFAR-10, CIFAR-100, and ImageNet. We also achieved significant performance improvements with various backbone architectures such as ResNet, WideResNet, and VGG.


title: Generalized distance to a simplex and a new geometrical method for portfolio optimization

abstract: Risk aversion plays a significant and central role in investors' decisions in the process of developing a portfolio. In this framework of portfolio optimization we determine the portfolio that possesses the minimal risk by using a new geometrical method. For this purpose, we elaborate an algorithm that enables us to compute any generalized Euclidean distance to a standard simplex. With this new approach, we are able to treat the case of portfolio optimization without short-selling in its entirety, and we also recover in geometrical terms the well-known results on portfolio optimization with allowed short-selling. Then, we apply our results in order to determine which convex combination of the CAC 40 stocks possesses the lowest risk: not only we get a very low risk compared to the index, but we also get a return rate that is almost three times better than the one of the index.


title: Steady State Of Random Dynamical Systems

abstract: Random dynamical systems (RDS) evolve by a dynamical rule chosen independently with a certain probability, from a given set of deterministic rules. These dynamical systems in an interval reach a steady state with a unique well-defined probability density only under certain conditions, namely Pelikan's criterion. We investigate and characterize the steady state of a bounded RDS when Pelikan's criterion breaks down. In this regime, the system is attracted to a common fixed point (CFP) of all the maps, which is attractive for at least one of the constituent mapping functions. If there are many such fixed points, the initial density is shared among the CFPs; we provide a mapping of this problem with the well known hitting problem of random walks and find the relative weights at different CFPs. The weights depend upon the initial distribution.


title: Radon daughter removal from PTFE surfaces and its application in liquid xenon detectors

abstract: Long-lived radon daughters are a critical background source in experiments searching for low-energy rare events. Originating from radon in ambient air, radioactive polonium, bismuth and lead isotopes plate-out on materials that are later employed in the experiment. In this paper, we examine cleaning procedures for their capability to remove radon daughters from PTFE surfaces, a material often used in liquid xenon TPCs. We found a large difference between the removal efficiency obtained for the decay chains of $^{222}$Rn and $^{220}$Rn, respectively. This indicates that the plate-out mechanism has an effect on the cleaning success. While the long-lived $^{222}$Rn daughters could be reduced by a factor of ~2, the removal of $^{220}$Rn daughters was up to 10 times more efficient depending on the treatment. Furthermore, the impact of a nitric acid based PTFE cleaning on the liquid xenon purity is investigated in a small-scale liquid xenon TPC.


title: Residual Spatial Attention Network for Retinal Vessel Segmentation

abstract: Reliable segmentation of retinal vessels can be employed as a way of monitoring and diagnosing certain diseases, such as diabetes and hypertension, as they affect the retinal vascular structure. In this work, we propose the Residual Spatial Attention Network (RSAN) for retinal vessel segmentation. RSAN employs a modified residual block structure that integrates DropBlock, which can not only be utilized to construct deep networks to extract more complex vascular features, but can also effectively alleviate the overfitting. Moreover, in order to further improve the representation capability of the network, based on this modified residual block, we introduce the spatial attention (SA) and propose the Residual Spatial Attention Block (RSAB) to build RSAN. We adopt the public DRIVE and CHASE DB1 color fundus image datasets to evaluate the proposed RSAN. Experiments show that the modified residual structure and the spatial attention are effective in this work, and our proposed RSAN achieves the state-of-the-art performance.