title: Fast computation of p-values for the permutation test based on Pearson's correlation coefficient and other statistical tests

abstract: Permutation tests are among the simplest and most widely used statistical tools. Their p-values can be computed by a straightforward sampling of permutations. However, this way of computing p-values is often so slow that it is replaced by an approximation, which is accurate only for part of the interesting range of parameters. Moreover, the accuracy of the approximation can usually not be improved by increasing the computation time. We introduce a new sampling-based algorithm which uses the fast Fourier transform to compute p-values for the permutation test based on Pearson's correlation coefficient. The algorithm is practically and asymptotically faster than straightforward sampling. Typically, its complexity is logarithmic in the input size, while the complexity of straightforward sampling is linear. The idea behind the algorithm can also be used to accelerate the computation of p-values for many other common statistical tests. The algorithm is easy to implement, but its analysis involves results from the representation theory of the symmetric group.


title: Omega Theorems for The Twisted Divisor Function

abstract: For a fixed $\theta\neq 0$, we define the twisted divisor function $$ \tau(n, \theta):=\sum_{d\mid n}d^{i\theta}\ .$$ In this article we consider the error term $\Delta(x)$ in the following asymptotic formula $$ \sum_{n\leq x}^*|\tau(n, \theta)|^2=\omega_1(\theta)x\log x + \omega_2(\theta)x\cos(\theta\log x) +\omega_3(\theta)x + \Delta(x),$$ where $\omega_i(\theta)$ for $i=1, 2, 3$ are constants depending only on $\theta$. We obtain $$\Delta(T)=\Omega\left(T^{\alpha(T)}\right) \text{ where } \alpha(T) =\frac{3}{8}-\frac{c}{(\log T)^{1/8}} \text{ and } c>0,$$ along with an $\Omega$-bound for the Lebesgue measure of the set of points where the above estimate holds.


title: Fitting a recurrent dynamical neural network to neural spiking data: Tackling with the sigmoidal gain function issues

abstract: This is a continuation of a recent study on the modeling of the information coding in sensory system in the brain. The data from a sensory neurons are available as discrete spike timings with no amplitude information. In the simulations, these are generated from a data generator model which has certain differences from the model being estimated. The model under consideration is simpler than the one used as a data generator as it has no sigmoidal gain function parameters. This choice is based on a suggestion from a recent study which states that inclusion of gain functions to the estimation algorithm increases issues such as parameter confounding which leads to performance degradation issues like increased or unpredictable variance changes with different stimuli configurations. To resolve this issue we consider a more generic model that has no sigmoidal gain functions to be estimated. Like that of the first research, we applied a Fourier series stimulus to both data generator and the estimated network. In order to test the performance of the proposed scheme, we have repeated the simulations with different sample sizes, stimulus component counts, amplitude and base frequencies. Mean values of estimation are presented as tables and the statistical analysis results are presented in graphical forms. In addition, since we do not have any true parameter data (for the estimated model) we compare the firing rate responses of both data generator and estimated models to different stimuli. It appeared that the responses of both models are almost the same.


title: Varying Driver Velocity Fields in Photospheric MHD Wave Simulations

abstract: Torsional motions are ubiquitous in the solar atmosphere. In this work, we perform 3D numerical simulations which mimic a vortex-type photospheric driver with a Gaussian spatial profile. This driver is implemented to excite MHD waves in an axially symmetric, 3D magnetic flux tube embedded in a realistic solar atmosphere. The Gaussian width of the driver is varied and the resulting perturbations are compared. Velocity vectors were decomposed into parallel, perpendicular and azimuthal components with respect to pre-defined magnetic flux surfaces. These components correspond broadly to the fast, slow and Alfven modes, respectively. From these velocities the corresponding wave energy fluxes are calculated, allowing us to estimate the contribution of each mode to the energy flux. For the narrowest driver ($0.15$ Mm) the parallel component accounts for $\sim 55 - 65\%$ of the flux. This contribution increases smoothly with driver width up to nearly $90\%$ for the widest driver ($0.35$ Mm). The relative importance of the perpendicular and azimuthal components decrease at similar rates. The azimuthal energy flux varied between $\sim 35\%$ for the narrowest driver and $< 10\%$ for the widest one. Similarly, the perpendicular flux was $\sim 25 - 10\%$. We also demonstrate that the fast mode corresponds to the sausage wave in our simulations. Our results therefore show that the fast sausage wave is easily excited by this driver and that it carries the majority of the energy transported. For this vortex-type driver the Alfven wave does not contribute a significant amount of energy.


title: Low-cost error mitigation by symmetry verification

abstract: We investigate the performance of error mitigation via measurement of conserved symmetries on near-term devices. We present two protocols to measure conserved symmetries during the bulk of an experiment, and develop a zero-cost post-processing protocol which is equivalent to a variant of the quantum subspace expansion. We develop methods for inserting global and local symetries into quantum algorithms, and for adjusting natural symmetries of the problem to boost their mitigation against different error channels. We demonstrate these techniques on two- and four-qubit simulations of the hydrogen molecule (using a classical density-matrix simulator), finding up to an order of magnitude reduction of the error in obtaining the ground state dissociation curve.


title: PCach: The Case for Pre-Caching your Mobile Data

abstract: We present PCach, a smartphone-based approach for relieving the congestion in cellular network resulting from the exponential growth in mobile data traffic. The basic idea underlying PCach is simple: use WiFi to proactively cache content on the smartphone's memory, which otherwise would have been delivered through the cellular network. However, it leads to several challenging questions, including how much mobile data actually flows through cellular networks, how much data can be pre-cached, and when and what to pre-cache. We address these questions progressively using a thorough analysis of user data collected from our purpose-built crowdsensing Android application, actively utilized by 45 users for periods dating back to July 2014. Our analysis shows that the median smartphone user transfers 15% of their data via the cellular network and that 80\% of it can be pre-cached via WiFi. To capitalize on these results, we draw on a careful analysis of the measurement data to introduce an algorithm that can run stand-alone on off-the-shelf smartphones and predict with good accuracy when and what to pre-cache.


title: Quantum optimal control of the dissipative production of a maximally entangled state

abstract: Entanglement generation can be robust against noise in approaches that deliberately incorporate dissipation into the system dynamics. The presence of additional dissipation channels may, however, limit fidelity and speed of the process. Here we show how quantum optimal control techniques can be used to both speed up the entanglement generation and increase the fidelity in a realistic setup, whilst respecting typical experimental limitations. For the example of entangling two trapped ion qubits [Lin et al., Nature 504, 415 (2013)], we find an improved fidelity by simply optimizing the polarization of the laser beams utilized in the experiment. More significantly, an alternate combination of transitions between internal states of the ions, when combined with optimized polarization, enables faster entanglement and decreases the error by an order of magnitude.


title: Prescribed mean curvature surfaces in the product spaces $\mathbb{M}^2(\kappa)\times\mathbb{R}$; Height estimates and classification results for properly embedded surfaces

abstract: The aim of this paper is to extend classic results of the theory of CMC surfaces in the product spaces to the class of immersed surfaces in $\mathbb{M}^2(\kappa)\times\mathbb{R}$ whose mean curvature is given as a $C^1$ function depending on their angle function. We cover topics such as the existence of a priori curvature estimates for graphs, height estimates for horizontal and vertical compact graphs and a structure-type result, which classifies all the simply connected, properly embedded surfaces with finite topology.


title: Optimizing the tight-binding parametrization of the quasi-one-dimensional superconductor K2Cr3As3

abstract: We study the tight-binding dispersion of the recently discovered superconductor K2Cr3As3, obtained from Wannier projection of Density Functional Theory (DFT) results. In order to establish quantitatively the actual degree of quasi-one dimensionality of this compound, we analyze the electronic band structure for two reduced sets of hopping parameters: one restricted to the Cr-As tubes and another one retaining a minimal number of in-plane hoppings. The corresponding total and local density of states of the compound are also computed with the aim of assessing the tight-binding results with respect to the DFT ones. We find a quite good agreement with the DFT results for the more extended set of hopping parameters, especially for what concerns the orbitals that dominate at the Fermi level. Accordingly, we conclude that one cannot avoid taking into account in-plane hoppings up to the next-nearest-neighbors cells even only to describe correctly the Fermi surface cuts and the populations along the kz direction. Such a choice of a minimal number of hopping parameters directly reflects in the possibility of correctly describing correlations and magnetic interactions.


title: Judging the Judges: A General Framework for Evaluating the Performance of International Sports Judges

abstract: The monitoring of judges and referees in sports has become an important topic due to the increasing media exposure of international sporting events and the large monetary sums involved. In this article, we present a method to assess the accuracy of sports judges and estimate their bias. Our method is broadly applicable to all sports where panels of judges evaluate athletic performances on a finite scale. We analyze judging scores from eight different sports with comparable judging systems: diving, dressage, figure skating, freestyle skiing (aerials), freestyle snowboard (halfpipe, slopestyle), gymnastics, ski jumping and synchronized swimming. With the notable exception of dressage, we identify, for each aforementioned sport, a general and accurate pattern of the intrinsic judging error as a function of the performance level of the athlete. This intrinsic judging inaccuracy is heteroscedastic and can be approximated by a quadratic curve, indicating increased consensus among judges towards the best athletes. Using this observation, the framework developed to assess the performance of international gymnastics judges is applicable to all these sports: we can evaluate the performance of judges compared to their peers and distinguish cheating from unintentional misjudging. Our analysis also leads to valuable insights about the judging practices of the sports under consideration. In particular, it reveals a systemic judging problem in dressage, where judges disagree on what constitutes a good performance.


title: FINJ: A Fault Injection Tool for HPC Systems

abstract: We present FINJ, a high-level fault injection tool for High-Performance Computing (HPC) systems, with a focus on the management of complex experiments. FINJ provides support for custom workloads and allows generation of anomalous conditions through the use of fault-triggering executable programs. FINJ can also be integrated seamlessly with most other lower-level fault injection tools, allowing users to create and monitor a variety of highly-complex and diverse fault conditions in HPC systems that would be difficult to recreate in practice. FINJ is suitable for experiments involving many, potentially interacting nodes, making it a very versatile design and evaluation tool.


title: Fluctuations of random Motzkin paths

abstract: It is known that after scaling a random Motzkin path converges to a Brownian excursion. We prove that the fluctuations of the counting processes of the ascent steps, the descent steps and the level steps converge jointly to linear combinations of two independent processes: a Brownian motion and a Brownian excursion. The proofs rely on the Laplace transforms and an integral representation based on an identity connecting non-crossing pair partitions and joint moments of an explicit non-homogeneous Markov process.


title: Fast cosine transform for FCC lattices

abstract: Voxel representation and processing is an important issue in a broad spectrum of applications. E.g., 3D imaging in biomedical engineering applications, video game development and volumetric displays are often based on data representation by voxels. By replacing the standard sampling lattice with a face-centered lattice one can obtain the same sampling density with less sampling points and reduce aliasing error, as well. We introduce an analog of the discrete cosine transform for the facecentered lattice relying on multivariate Chebyshev polynomials. A fast algorithm for this transform is deduced based on algebraic signal processing theory and the rich geometry of the special unitary Lie group of degree four.


title: From Monolith to Microservices: A Classification of Refactoring Approaches

abstract: While the recently emerged Microservices architectural style is widely discussed in literature, it is difficult to find clear guidance on the process of refactoring legacy applications. The importance of the topic is underpinned by high costs and effort of a refactoring process which has several other implications, e.g. overall processes (DevOps) and team structure. Software architects facing this challenge are in need of selecting an appropriate strategy and refactoring technique. One of the most discussed aspects in this context is finding the right service granularity to fully leverage the advantages of a Microservices architecture. This study first discusses the notion of architectural refactoring and subsequently compares 10 existing refactoring approaches recently proposed in academic literature. The approaches are classified by the underlying decomposition technique and visually presented in the form of a decision guide for quick reference. The review yielded a variety of strategies to break down a monolithic application into independent services. With one exception, most approaches are only applicable under certain conditions. Further concerns are the significant amount of input data some approaches require as well as limited or prototypical tool support.


title: A lower Wegner estimate and bounds on the spectral shift function for continuum random Schr\"odinger operators

abstract: We prove a strictly positive, locally uniform lower bound on the density of states (DOS) of continuum random Schr\"odinger operators on the entire spectrum, i.e. we show that the DOS does not have a zero within the spectrum. This follows from a lower Wegner estimate for finite-volume continuum random Schr\"odinger operators. We assume throughout iid random variables and the single-site distribution having a Lebesgue density bounded from below on its support. The main mathematical novelty in this paper are pointwise-in-energy bounds on the expectation of the spectral shift function at all energies for these operators where we mainly focus on perturbations corresponding to a change from Dirichlet to Neumann boundary conditions along the boundary of a cube. We show that the bound scales with the area of the hypersurface where the boundary conditions are changed. We also prove bounds on the averaged spectral shift function for perturbations by bounded and compactly supported multiplication operators.


title: Origin of the hemispheric asymmetry of solar activity

abstract: The frequency spectrum of the hemispheric asymmetry of solar activity shows enhanced power for the period ranges around 8.5 years and between 30 and 50 years. This can be understood as the sum and beat periods of the superposition of two dynamo modes: a dipolar mode with a (magnetic) period of about 22 years and aquadrupolar mode with a period between 13 and 15 years. An updated Babcock-Leighton-type dynamo model with weak driving as indicated by stellar observations shows an excited dipole mode and a damped quadrupole mode in the correct range of periods. Random excitation of the quadrupole by stochastic fluctuations of the source term for the poloidal field leads to a time evolution of activity and asymmetry that is consistent with the observational results.


title: A note on Ramsey numbers for Berge-G hyper graphs

abstract: For a graph G=(V,E), a hypergraph H is called Berge-G if there is a bijection f from E(G) to E(H) such that for each e in E(G), e is a subset of f(e). The set of all Berge-G hypergraphs is denoted B(G). For integers k>1, r>1, and a graph G, let the Ramsey number R_r(B(G), k) be the smallest integer n such that no matter how the edges of a complete r-uniform n-vertex hypergraph are colored with k colors, there is a copy of a monochromatic Berge-G subhypergraph. Furthermore, let R(B(G),k) be the smallest integer n such that no matter how all subsets an n-element set are colored with k colors, there is a monochromatic copy of a Berge-G hypergraph. We give an upper bound for R_r(B(G),k) in terms of graph Ramsey numbers. In particular, we prove that when G becomes acyclic after removing some vertex, R_r(B(G),k)\le 4k|V(G)|+r-2, in contrast with classical multicolor Ramsey numbers. When G is a triangle or a K_4, we find sharper bounds and some exact results and determine some `small' Ramsey numbers: k/2 - o(k) < R_3(B(K_3)), k) < 3k/4+ o(k), For any odd integer t\neq 3, R(B(K_3),2^t-1)=t+2, 2^{ck} < R_3(B(K_4),k) < e(1+o(1))(k-1)k!, R_3(B(K_3),2)=R_3(B(K_3),3)=5, R_3(B(K_3),4)=6, R_3(B(K_3),5)=7, R_3(B(K_3),6)=8, R_3(B(K_3,8)=9, R_3(B(K_4),2)=6.


title: Differential of metric valued Sobolev maps

abstract: We introduce a notion of differential of a Sobolev map between metric spaces. The differential is given in the framework of tangent and cotangent modules of metric measure spaces, developed by the first author. We prove that our notion is consistent with Kirchheim's metric differential when the source is a Euclidean space, and with the abstract differential provided by the first author when the target is $\mathbb{R}$.


title: Renormalon-free definition of the gluon condensate within the large-$\beta_0$ approximation

abstract: We propose a clear definition of the gluon condensate within the large-$\beta_0$ approximation as an attempt toward a systematic argument on the gluon condensate. We define the gluon condensate such that it is free from a renormalon uncertainty, consistent with the renormalization scale independence of each term of the operator product expansion (OPE), and an identical object irrespective of observables. The renormalon uncertainty of $\mathcal{O}(\Lambda^4)$, which renders the gluon condensate ambiguous, is separated from a perturbative calculation by using a recently suggested analytic formulation. The renormalon uncertainty is absorbed into the gluon condensate in the OPE, which makes the gluon condensate free from the renormalon uncertainty. As a result, we can define the OPE in a renormalon-free way. Based on this renormalon-free OPE formula, we discuss numerical extraction of the gluon condensate using the lattice data of the energy density operator defined by the Yang--Mills gradient flow.


title: Differential geometry of line of curvature on parametric hypersurfaces in $\mathbb{E}^4$

abstract: The purpose of this paper is, first, to give an algorithm that enables to obtain the lines of curvature on parametric hypersurfaces in Euclidean 4-space, and then, to obtain the curvatures of such lines by using the extended Darboux frame along the curve.


title: A Better Baseline for AVA

abstract: We introduce a simple baseline for action localization on the AVA dataset. The model builds upon the Faster R-CNN bounding box detection framework, adapted to operate on pure spatiotemporal features - in our case produced exclusively by an I3D model pretrained on Kinetics. This model obtains 21.9% average AP on the validation set of AVA v2.1, up from 14.5% for the best RGB spatiotemporal model used in the original AVA paper (which was pretrained on Kinetics and ImageNet), and up from 11.3 of the publicly available baseline using a ResNet101 image feature extractor, that was pretrained on ImageNet. Our final model obtains 22.8%/21.9% mAP on the val/test sets and outperforms all submissions to the AVA challenge at CVPR 2018.


title: Exactness of Bohr-Sommerfeld quantisation for two non-central potentials

abstract: In this paper we demonstrate the integrability of the Hamilton-Jacobi equation for two non-central potentials in spherical polar coordinates, and present complete solutions for the classically bound orbits. We then show that the semiclassical method of Bohr-Sommerfeld quantisation exactly reproduces the bound state spectra of the corresponding quantum mechanical Schr\"odinger equations. One of these potentials has previously been analysed in parabolic coordinates; the results for the other are, to the authors' best knowledge, original.


title: A Linear Constrained Optimization Benchmark For Probabilistic Search Algorithms: The Rotated Klee-Minty Problem

abstract: The development, assessment, and comparison of randomized search algorithms heavily rely on benchmarking. Regarding the domain of constrained optimization, the number of currently available benchmark environments bears no relation to the number of distinct problem features. The present paper advances a proposal of a scalable linear constrained optimization problem that is suitable for benchmarking Evolutionary Algorithms. By comparing two recent EA variants, the linear benchmarking environment is demonstrated.


title: Third and fourth order well-balanced schemes for the shallow water equations based on the CWENO reconstruction

abstract: High order finite volume schemes for conservation laws are very useful in applications, due to their ability to compute accurate solutions on quite coarse meshes and with very few restrictions on the kind of cells employed in the discretization. For balance laws, the ability to approximate up to machine precision relevant steady states allows the scheme to compute accurately, also on coarse meshes, small perturbations of such states. In this paper we propose third and fourth order accurate finite volume schemes for the shallow water equations. The schemes have the well-balanced property thanks to a path-conservative approach applied to an appropriate non-conservative reformulation of the equations. High order accuracy is achieved by designing truly two-dimensional reconstruction procedures of the CWENO type. The novel schemes are tested for accuracy, well-balancing and shown to maintain posivity of the water height on wet/dry transitions. Finally they are applied to simulate the Tohoku 2011 tsunami event.


title: Construction of a quotient ring of $\mathbb{Z}_2\mathcal{F}$ in which a binomial $1 + w$ is invertible using small cancellation methods

abstract: We apply small cancellation methods originating from group theory to investigate the structure of a quotient ring $\mathbb{Z}_2\mathcal{F} / \mathcal{I}$, where $\mathbb{Z}_2\mathcal{F}$ is the group algebra of the free group $\mathcal{F}$ over the field $\mathbb{Z}_2$, and the ideal $\mathcal{I}$ is generated by a single trinomial $1 + v + vw$, where $v$ is a complicated word depending on $w$. In $\mathbb{Z}_2\mathcal{F} / \mathcal{I}$ we have $(1 + w)^{-1} = v$, so $1 + w$ becomes invertible. We construct an explicit linear basis of $\mathbb{Z}_2\mathcal{F} / \mathcal{I}$ (thus showing that $\mathbb{Z}_2\mathcal{F} / \mathcal{I}\neq 0$). This is the first step in constructing rings with exotic properties.


title: Tied links and invariants for singular links

abstract: Tied links and the tied braid monoid were introduced recently by the authors and used to define new invariants for classical links. Here, we give a version purely algebraic-combinatoric of tied links. With this new version we prove that the tied braid monoid has a decomposition like a semi--direct group product. By using this decomposition we reprove the Alexander and Markov theorem for tied links; also, we introduce the tied singular knots, the tied singular braid monoid and certain families of Homflypt type invariants for tied singular links; these invariants are five-variables polynomials. Finally, we study the behavior of these invariants; in particular, we show that our invariants distinguish non isotopic singular links indistinguishable by the Paris-Rabenda invariant.


title: Pion-kaon femtoscopy in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV measured with ALICE

abstract: Femtoscopic correlations between charged pions and kaons for different charge combinations are measured in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV with ALICE at the LHC. The three-dimensional pion-kaon ($\pi-{\rm K}$) correlation functions and double ratios in the out-side-long pair rest frame are studied in different centrality bins. The $\pi-{\rm K}$ femtoscopic source size parameter ($R_{\rm out}$) and emission asymmetry ($\mu_{\rm out}$) are extracted. It is observed that the average source size of the system and the emission asymmetry between pions and kaons increase from peripheral to central events.


title: Loosely-Coupled Semi-Direct Monocular SLAM

abstract: We propose a novel semi-direct approach for monocular simultaneous localization and mapping (SLAM) that combines the complementary strengths of direct and feature-based methods. The proposed pipeline loosely couples direct odometry and feature-based SLAM to perform three levels of parallel optimizations: (1) photometric bundle adjustment (BA) that jointly optimizes the local structure and motion, (2) geometric BA that refines keyframe poses and associated feature map points, and (3) pose graph optimization to achieve global map consistency in the presence of loop closures. This is achieved in real-time by limiting the feature-based operations to marginalized keyframes from the direct odometry module. Exhaustive evaluation on two benchmark datasets demonstrates that our system outperforms the state-of-the-art monocular odometry and SLAM systems in terms of overall accuracy and robustness.


title: On the multiwavelength properties of several {\gamma}-ray detected narrow-line Seyfert 1 galaxies

abstract: The {\gamma}-ray detection from several radio-loud (RL) narrow-line Seyfert 1 (NLS1) galaxies has enabled us to study powerful relativistic jets in active galactic nuclei (AGNs) with smaller black hole masses and higher accretion rates than classical blazars. However, the sample of those {\gamma}-ray detected NLS1s available is still not large enough for a comprehensive and statistical study. We provide a summary of our detections and follow-up studies of three {\gamma}-ray-emitting NLS1s: SDSS J211852.96-073227.5 with flaring {\gamma}-ray radiation (Yang et al. 2018, Paliya et al. 2018) and SDSS J122222.55+041315.7 with the highest redshift by far (z~1) (Yao et al. 2015a), along with the prototype 1H 0323+342 (Zhou et al. 2007, Yao et al. 2015b). And we will discuss their multiwavelength properties and variability properties, including implications from high-energy observations in {\gamma}-rays and X-rays, infrared and radio properties, and correlated variability between several wavebands.


title: Pseudo-scalar {\boldm $ q\bar q $} bound states at finite temperatures within a Dyson-Schwinger--Bethe-Salpeter approach

abstract: The combined Dyson-Schwinger--Bethe-Salpeter equations are employed at non-zero temperature. The truncations refer to a rainbow-ladder approximation augmented with an interaction kernel which facilitates a special temperature dependence. At low temperatures, $T \to 0$, we recover a quark propagator from the Dyson-Schwinger (gap) equation which delivers, e.g. mass functions $B$, quark renormalization wave function $A$, and two-quark condensate $\la q \bar q \ra$ smoothly interpolating to the $T = 0$ results, despite the broken O(4) symmetry in the heat bath and discrete Matsubara frequencies. Besides the Matsubara frequency difference entering the interaction kernel, often a Debye screening mass term is introduced when extending the $T = 0$ kernel to non-zero temperatures. At larger temperatures, however, we are forced to drop this Debye mass in the infra-red part of the longitudinal interaction kernel to keep the melting of the two-quark condensate in a range consistent with lattice QCD results. Utilizing that quark propagator for the first few hundred fermion Matsubara frequencies we evaluate the Bethe-Salpeter vertex function in the pseudo-scalar $ q \bar q$ channel for the lowest boson Matsubara frequencies and find a competition of $ q \bar q$ bound states and quasi-free two-quark states at $T = {\cal O}$ (100 MeV). This indication of pseudo-scalar meson dissociation below the anticipated QCD deconfinement temperature calls for an improvement of the approach, which is based on an interaction adjusted to the meson spectrum at $T = 0$.


title: Concurrent Learning of Semantic Relations

abstract: Discovering whether words are semantically related and identifying the specific semantic relation that holds between them is of crucial importance for NLP as it is essential for tasks like query expansion in IR. Within this context, different methodologies have been proposed that either exclusively focus on a single lexical relation (e.g. hypernymy vs. random) or learn specific classifiers capable of identifying multiple semantic relations (e.g. hypernymy vs. synonymy vs. random). In this paper, we propose another way to look at the problem that relies on the multi-task learning paradigm. In particular, we want to study whether the learning process of a given semantic relation (e.g. hypernymy) can be improved by the concurrent learning of another semantic relation (e.g. co-hyponymy). Within this context, we particularly examine the benefits of semi-supervised learning where the training of a prediction function is performed over few labeled data jointly with many unlabeled ones. Preliminary results based on simple learning strategies and state-of-the-art distributional feature representations show that concurrent learning can lead to improvements in a vast majority of tested situations.


title: High-temperature Dust Condensation around an AGB Star: Evidence from a Highly Pristine Presolar Corundum

abstract: Corundum ($\alpha$-Al$_{2}$O$_{3}$) and amorphous or metastable Al$_{2}$O$_{3}$ are common components of circumstellar dust observed around O-rich asymptotic giant branch (AGB) stars and found in primitive meteorites. We report a detailed isotopic and microstructural investigation of a unique presolar corundum grain, QUE060, identified in an acid residue of the Queen Alexandra Range 97008 (LL3.05) meteorite. Based on its O and Mg isotopic compositions, this 1.4 $\mu$m diameter grain formed in a low- or intermediate-mass AGB star. It has four developed rhombohedral $\{$011$\}$ faces of corundum and a rough, rounded face with cavities. High Mg contents (Mg/Al $>$ 0.004) are due to the decay of radioactive $^{26}$Al. No spinel (MgAl$_{2}$O$_{4}$) inclusions that might have exsolved from the corundum are observed, but there are several high-Mg domains with modulated structures. The subhedral shape of grain QUE060 is the first clear evidence that corundum condenses and grows to micrometer sizes in the extended atmospheres around AGB stars. The flat faces indicate that grain QUE060 experienced little modification by gas-grain and grain-grain collisions in the interstellar medium (ISM) and solar nebula. The Mg distribution in its structure indicates that grain QUE060 has not experienced any severe heating events since the exhaustion of $^{26}$Al. However, it underwent at least one very transient heating event to form the high-Mg domains. A possible mechanism for producing this transient event, as well as the one rough surface and cavity, is a single grain-grain collision in the ISM. These results indicate that grain QUE060 is the most pristine circumstellar corundum studied to date.


title: Topological crystalline insulators from stacked graphene layers

abstract: In principle the stacking of different two-dimensional (2D) materials allows the construction of 3D systems with entirely new electronic properties. Here we propose to realize topological crystalline insulators (TCI) protected by mirror symmetry in heterostructures consisting of graphene monolayers separated by two-dimensional polar spacers. The polar spacers are arranged such that they can induce an alternating doping and/or spin-orbit coupling in the adjacent graphene sheets. When spin-orbit coupling dominates, the non-trivial phase arises due to the fact that each graphene sheet enters a quantum spin-Hall phase. Instead, when the graphene layers are electron and hole doped in an alternating fashion, a uniform magnetic field leads to the formation of quantum Hall phases with opposite Chern numbers. It thus has the remarkable property that unlike previously proposed and observed TCIs, the non-trivial topology is generated by an external time-reversal breaking perturbation.


title: Automatic Detection of Node-Replication Attack in Vehicular Ad-hoc Networks

abstract: Recent advances in smart cities applications enforce security threads such as node replication attacks. Such attack is take place when the attacker plants a replicated network node within the network. Vehicular Ad hoc networks are connecting sensors that have limited resources and required the response time to be as low as possible. In this type networks, traditional detection algorithms of node replication attacks are not efficient. In this paper, we propose an initial idea to apply a newly adapted statistical methodology that can detect node replication attacks with high performance as compared to state-of-the-art techniques. We provide a sufficient description of this methodology and a road-map for testing and experiment its performance.


title: A discrete Hopf-Rinow-theorem

abstract: We prove a version of the Hopf-Rinow-theorem with respect to path metrics on discrete spaces. The novel aspect is that we do not a priori assume local finiteness but isolate a local finiteness type condition, called essential local finiteness, that is indeed necessary. As a side product we identify the maximal weight, called the geodesic weight, which generates the path metric in the situation when the space is complete with respect to any of the equivalent notions of completeness proven in the Hopf-Rinow theorem. As an application, we characterize the graphs for which the resistance metric is a path metric induced by the graph structure.


title: Astronomical interferometry with near-IR e-APD at CHARA: characterization, optimization and on-sky operation

abstract: We characterize a near-infrared C-RED ONE camera from First Light Imaging (FLI). This camera uses a SAPHIRA electron avalanche photo-diode array (e-APD) from Leonardo (previously Selex). To do so, we developed a model of the signal distribution. This model allows a measurement of the gain and the Excess Noise Factor (ENF) independently of preexisting calibration such as the system gain. The results of this study show a gain which is 0.53 +/- 0.04 times the gain reported by the manufacturer. The measured ENF is 1.47 +/- 0.03 when we expected 1.25. For an avalanche gain of 60 and a frame rate larger than 100 Hz, the total noise can be lower than 1 e-/frame/pixel. The lowest dark current level is 90e-/s/pixel, in agreement with the expected H-band background passing through the camera window. These performance values provide a significant improvement compared to earlier-generation PICNIC camera and allowed us to improve the performance of the Michigan infrared combiner (MIRC) instrument at the Center for High Angular Resolution Astronomy (CHARA), as part of our MIRC-X instrumentation project.


title: Computing $\mathcal{L}$-invariants via the Greenberg-Stevens formula

abstract: In this article, we describe how to compute slopes of $p$-adic $\mathcal{L}$-invariants of arbitrary weight and level by means of the Greenberg-Stevens formula. Our method is based on work of Lauder and Vonk on computing the reverse characteristic series of the $U_p$ operator on overconvergent modular forms. Using higher derivatives of this characteristic series, we construct a polynomial whose zeros are precisely the $\mathcal{L}$-invariants appearing in the corresponding space of modular forms with fixed sign of the Atkin-Lehner involution at $p$. In addition, we describe how to compute this polynomial efficiently. In the final section, we give computational evidence for relations between slopes of $\mathcal{L}$-invariants for small primes.


title: Optimal Design in Hierarchical Models with application in Multi-center Trials

abstract: Hierarchical random effect models are used for different purposes in clinical research and other areas. In general, the main focus is on population parameters related to the expected treatment effects or group differences among all units of an upper level (e.g. subjects in many settings). Optimal design for estimation of population parameters are well established for many models. However, optimal designs for the prediction for the individual units may be different. Several settings are identiffed in which individual prediction may be of interest. In this paper we determine optimal designs for the individual predictions, e.g. in multi-center trials, and compare them to a conventional balanced design with respect to treatment allocation. Our investigations show, that balanced designs are far from optimal if the treatment effects vary strongly as compared to the residual error and more subjects should be recruited to the active (new) treatment in multi-center trials. Nevertheless, effciency loss may be limited resulting in a moderate sample size increase when individual predictions are foreseen with a balanced allocation.


title: Nonreciprocal Photon Blockade

abstract: We propose how to create and manipulate one-way nonclassical light via photon blockade in rotating nonlinear devices. We refer to this effect as nonreciprocal photon blockade (PB). Specifically, we show that in a spinning Kerr resonator, PB happens when the resonator is driven in one direction but not the other. This occurs because of the {Fizeau drag,} leading to a full split of the resonance frequencies of the countercirculating modes. Different types of purely quantum correlations, such as single- and two-photon blockades, can emerge in different directions in a well-controlled manner, and the transition from PB to photon-induced tunneling is revealed as well. Our work opens up a new route to achieve quantum nonreciprocal devices, which are crucial elements in chiral quantum technologies or topological photonics.


title: Black hole geodesic parallel transport and the Marck recipe for isolating cumulative precession effects

abstract: The Wigner rotations arising from the combination of boosts along two different directions are rederived from a relative boost point of view and applied to gyroscope spin precession along timelike geodesics in a Kerr spacetime, clarifying the geometrical properties of Marck's recipe for describing parallel transport along such world lines expressed in terms of the constants of the motion. His final angular velocity isolates the cumulative spin precession angular velocity independent of the spacetime tilting required to keep the spin 4-vector orthogonal to the gyro 4-velocity. As an explicit example the cumulative precession effects are computed for a test gyroscope moving along both bound and unbound equatorial plane geodesic orbits.


title: Rational approximations to fractional powers of self-adjoint positive operators

abstract: We investigate the rational approximation of fractional powers of unbounded positive operators attainable with a specific integral representation of the operator function. We provide accurate error bounds by exploiting classical results in approximation theory involving Pad\'{e} approximants. The analysis improves some existing results and the numerical experiments proves its accuracy.


title: Observations of meteors in the Earth's atmosphere: Reducing data from dedicated double-station wide-angle cameras

abstract: Meteoroids entering the Earth's atmosphere can be observed as meteors, thereby providing useful information on their formation and hence on their parent bodies. We developed a data reduction software package for double station meteor data from the SPOSH camera, which includes event detection, image geometric and radiometric calibration, radiant and speed estimates, trajectory and orbit determination, and meteor light curve recovery. The software package is designed to fully utilise the high photometric quality of SPOSH images. This will facilitate the detection of meteor streams and studies of their trajectories. We have run simulations to assess the performance of the software by estimating the radiants, speeds, and magnitudes of synthetic meteors and comparing them with the a priori values. The estimated uncertainties in radiant location had a zero mean with a median deviation between 0.03$^{\circ}$ and 0.11$^{\circ}$ for the right ascension and 0.02$^{\circ}$ and 0.07$^{\circ}$ for the declination. The estimated uncertainties for the speeds had a median deviation between 0.40 and 0.45 km s$^{-1}$. The brightness of synthetic meteors was estimated to within +0.01m. We have applied the software package to 177 real meteors acquired by the SPOSH camera. The median propagated uncertainties in geocentric right ascension and declination were found to be of 0.64$^{\circ}$ and 0.29$^{\circ}$, while the median propagated error in geocentric speed was 1.21 km $s^{-1}$.


title: AlphaGAN: Generative adversarial networks for natural image matting

abstract: We present the first generative adversarial network (GAN) for natural image matting. Our novel generator network is trained to predict visually appealing alphas with the addition of the adversarial loss from the discriminator that is trained to classify well-composited images. Further, we improve existing encoder-decoder architectures to better deal with the spatial localization issues inherited in convolutional neural networks (CNN) by using dilated convolutions to capture global context information without downscaling feature maps and losing spatial information. We present state-of-the-art results on the alphamatting online benchmark for the gradient error and give comparable results in others. Our method is particularly well suited for fine structures like hair, which is of great importance in practical matting applications, e.g. in film/TV production.


title: Discovery of Kolmogorov-like Magnetic Energy Spectrum in Tycho's Supernova Remnant by Two-point Correlations of Synchrotron Intensity

abstract: The spectral slope of the magnetic energy in supernova remnants (SNRs) can be obtained by analysis of spatial two-point correlation functions of synchrotron intensities. This method has been originally developed for the analysis of magnetic field structure in diffuse interstellar medium and applied when the geometry of the emission region is simple and known. In this paper, applying this correlation analysis to Tycho's SNR, for which the synchrotron emission region is known to be a spherical shell, we find that the magnetic energy spectrum shows the Kolmogorov-like scaling. Our results can be explained by turbulence developed downstream of the shock front via Richtmyer-Meshkov instability or an amplification of upstream magnetic field induced by cosmic rays. They could be discriminated by future observations with a sub arcsecond resolution such as Square Kilometer Array.


title: A variational model for data fitting on manifolds by minimizing the acceleration of a B\'ezier curve

abstract: We derive a variational model to fit a composite B\'ezier curve to a set of data points on a Riemannian manifold. The resulting curve is obtained in such a way that its mean squared acceleration is minimal in addition to remaining close the data points. We approximate the acceleration by discretizing the squared second order derivative along the curve. We derive a closed-form, numerically stable and efficient algorithm to compute the gradient of a B\'ezier curve on manifolds with respect to its control points, expressed as a concatenation of so-called adjoint Jacobi fields. Several examples illustrate the capabilites and validity of this approach both for interpolation and approximation. The examples also illustrate that the approach outperforms previous works tackling this problem.


title: Initiation of fluid-induced fracture in a thick-walled hollow permeable sphere

abstract: The initiation of fluid-induced fracture in formations of permeable geomaterials subjected to quasi-stationary flow processes (drained response) can be strongly affected by Biot's coefficient and the size of the formation. The aim of this paper is to analyse the influence of these parameters on the initial fracture process of a thick-walled hollow permeable sphere subjected to fluid injection in the hole. Assuming that fracture patterns are distributed uniformly during the hardening stage of the fracture initiation process, the coupled fluid-solid problem is described by a nonlinear ordinary differential equation, which is solved numerically by means of finite differences combined with shooting and Newton methods. The finite difference code has also been validated in the elastic range, i.e., before initiation of fracture, against an original closed-form analytical solution of the above differential equation. The results show that the nominal strength of the sphere increases with increasing Biot's coefficient and decreases with increasing size.


title: High-performance source of spectrally pure, polarization entangled photon pairs based on hybrid integrated-bulk optics

abstract: Entangled photon pair sources based on bulk optics are approaching optimal design and implementation, with high state fidelities, spectral purities and heralding efficiencies, but generally low brightness. Integrated entanglement sources, while providing higher brightness and low-power operation, often sacrifice performance in output state quality and coupling efficiency. Here we present a polarization-entangled pair source based on a hybrid approach of waveguiding and bulk optics, addressing every metric simultaneously. We show 96% fidelity to the singlet state, 82% Hong-Ou-Mandel interference visibility, 43% average Klyshko efficiency, and a high brightness of $2.9\times10^6$ pairs/(mode$\cdot$s$\cdot$mW), while requiring only microwatts of pump power.


title: Computing optimal shortcuts for networks

abstract: We study augmenting a plane Euclidean network with a segment, called a shortcut, to minimize the largest distance between any two points along the edges of the resulting network. Problems of this type have received considerable attention recently, mostly for discrete variants of the problem. We consider a fully continuous setting, where the problem of computing distances and placing a shortcut is much harder as all points on the network, instead of only the vertices, must be taken into account. We present the first results on the computation of optimal shortcuts for general networks in this model: a polynomial time algorithm and a discretization of the problem that leads to an approximation algorithm. We also improve the general method for networks that are paths, restricted to two types of shortcuts: those with a fixed orientation and simple shortcuts.


title: Semi-regular Dubuc-Deslauriers wavelet tight frames

abstract: In this paper, we construct wavelet tight frames with n vanishing moments for Dubuc-Deslauriers 2npoint semi-regular interpolatory subdivision schemes. Our motivation for this construction is its practical use for further regularity analysis of wide classes of semi-regular subdivision. Our constructive tools are local eigenvalue convergence analysis for semi-regular Dubuc-Deslauriers subdivision, the Unitary Extension Principle and the generalization of the Oblique Extension Principle to the irregular setting by Chui, He and St\"ockler. This group of authors derives suitable approximation of the inverse Gramian for irregular Bspline subdivision. Our main contribution is the derivation of the appropriate approximation of the inverse Gramian for the semi-regular Dubuc-Deslauriers scaling functions ensuring n vanishing moments of the corresponding framelets.


title: Some results related to finiteness properties of groups for families of subgroups

abstract: For a group $G$ we consider the classifying space $E_{\mathcal{VC}yc}(G)$ for the family of virtually cyclic subgroups. We show that an Artin group admits a finite model for $E_{\mathcal{VC}yc}(G)$ if and only if it is virtually cyclic. This solves a conjecture of Juan-Pineda and Leary and a question of L\"uck-Reich-Rognes-Varisco for Artin groups. We then study the conjugacy growth of CAT(0) groups and show that if a CAT(0) group contains a free abelian group of rank two, its conjugacy growth is strictly faster than linear. This also yields an alternative proof for the fact that a CAT(0) cube group admits a finite model for $E_{\mathcal{VC}yc}(G)$ if and only if it is virtually cyclic. Our last result deals with the homotopy type of the quotient space $B_{\mathcal{VC}yc}(G) = E_{\mathcal{VC}yc}(G)/G$. We show for a poly-$\mathbb Z$-group $G$, that $B_{\mathcal{VC}yc}(G)$ is homotopy equivalent to a finite CW-complex if and only if $G$ is cyclic.


title: Towards a Deep Unified Framework for Nuclear Reactor Perturbation Analysis

abstract: In this paper, we take the first steps towards a novel unified framework for the analysis of perturbations in both the Time and Frequency domains. The identification of type and source of such perturbations is fundamental for monitoring reactor cores and guarantee safety while running at nominal conditions. A 3D Convolutional Neural Network (3D-CNN) was employed to analyse perturbations happening in the frequency domain, such as an absorber of variable strength or propagating perturbation. Recurrent neural networks (RNN), specifically Long Short-Term Memory (LSTM) networks were used to study signal sequences related to perturbations induced in the time domain, including the vibrations of fuel assemblies and the fluctuations of thermal-hydraulic parameters at the inlet of the reactor coolant loops. 512 dimensional representations were extracted from the 3D-CNN and LSTM architectures, and used as input to a fused multi-sigmoid classification layer to recognise the perturbation type. If the perturbation is in the frequency domain, a separate fully-connected layer utilises said representations to regress the coordinates of its source. The results showed that the perturbation type can be recognised with high accuracy in all cases, and frequency domain scenario sources can be localised with high precision.


title: Learning to predict crisp boundaries

abstract: Recent methods for boundary or edge detection built on Deep Convolutional Neural Networks (CNNs) typically suffer from the issue of predicted edges being thick and need post-processing to obtain crisp boundaries. Highly imbalanced categories of boundary versus background in training data is one of main reasons for the above problem. In this work, the aim is to make CNNs produce sharp boundaries without post-processing. We introduce a novel loss for boundary detection, which is very effective for classifying imbalanced data and allows CNNs to produce crisp boundaries. Moreover, we propose an end-to-end network which adopts the bottom-up/top-down architecture to tackle the task. The proposed network effectively leverages hierarchical features and produces pixel-accurate boundary mask, which is critical to reconstruct the edge map. Our experiments illustrate that directly making crisp prediction not only promotes the visual results of CNNs, but also achieves better results against the state-of-the-art on the BSDS500 dataset (ODS F-score of .815) and the NYU Depth dataset (ODS F-score of .762).


title: Glueing a peak to a non-zero limiting profile for a critical Moser-Trudinger equation

abstract: Druet [6] proved that if $(f_\gamma)_\gamma$ is a sequence of Moser-Trudinger type nonlinearities with critical growth, and if $(u_\gamma)_\gamma$ solves $$ \begin{cases} &\Delta u =f_\gamma(x,u)\,,~~ u>0\text{ in }\Omega\,,\\ &u =0\text{ on }\partial\Omega\,, \end{cases} $$ and converges weakly in $H^1_0$ to some $u_\infty$, then the Dirichlet energy is quantified, namely there exists an integer $N\ge 0$ such that the energy of $u_\gamma$ converges to $4\pi N$ plus the Dirichlet energy of $u_\infty$. As a crucial step to get the general existence results of [7], it was more recently proved in [8] that, for a specific class of nonlinearities, the loss of compactness (i.e. $N>0$) implies that $u_\infty\equiv 0$. In contrast, we prove here that there exist sequences $(f_\gamma)_\gamma$ of Moser-Trudinger type nonlinearities which admit a noncompact sequence of solutions $(u_\gamma)_\gamma$ having a nontrivial weak limit.


title: Scattering due to geometry: The case of a spinless particle moving on an asymptotically flat embedded surface

abstract: A nonrelativistic quantum mechanical particle moving freely on a curved surface feels the effect of the nontrivial geometry of the surface through the kinetic part of the Hamiltonian, which is proportional to the Laplace-Beltrami operator, and a geometric potential, which is a linear combination of the mean and Gaussian curvatures of the surface. The coefficients of these terms cannot be uniquely determined by general principles of quantum mechanics but enter the calculation of various physical quantities. We examine their contribution to the geometric scattering of a scalar particle moving on an asymptotically flat embedded surface. In particular, having in mind the possibility of an experimental realization of the geometric scattering in a low density electron gas formed on a bumped surface, we determine the scattering amplitude for arbitrary choices of the curvature coefficients for a surface with global or local cylindrical symmetry. We also examine the effect of perturbations that violate this symmetry and consider surfaces involving bumps that form a lattice.


title: Two-Step Estimation and Inference with Possibly Many Included Covariates

abstract: We study the implications of including many covariates in a first-step estimate entering a two-step estimation procedure. We find that a first order bias emerges when the number of \textit{included} covariates is "large" relative to the square-root of sample size, rendering standard inference procedures invalid. We show that the jackknife is able to estimate this "many covariates" bias consistently, thereby delivering a new automatic bias-corrected two-step point estimator. The jackknife also consistently estimates the standard error of the original two-step point estimator. For inference, we develop a valid post-bias-correction bootstrap approximation that accounts for the additional variability introduced by the jackknife bias-correction. We find that the jackknife bias-corrected point estimator and the bootstrap post-bias-correction inference perform excellent in simulations, offering important improvements over conventional two-step point estimators and inference procedures, which are not robust to including many covariates. We apply our results to an array of distinct treatment effect, policy evaluation, and other applied microeconomics settings. In particular, we discuss production function and marginal treatment effect estimation in detail.


title: Series representations for fractional-calculus operators involving generalised Mittag-Leffler functions

abstract: We consider an integral transform introduced by Prabhakar, involving generalised multi-parameter Mittag-Leffler functions, which can be used to introduce and investigate several different models of fractional calculus. We derive a new series expression for this transform, in terms of classical Riemann-Liouville fractional integrals, and use it to obtain or verify series formulae in various specific cases corresponding to different fractional-calculus models. We demonstrate the power of our result by applying the series formula to derive analogues of the product and chain rules in more general fractional contexts. We also discuss how the Prabhakar model can be used to explore the idea of fractional iteration in connection with semigroup properties.


title: Tunneling blockade and single-photon emission in GaAs double quantum wells

abstract: We report on the selective excitation of single impurity-bound exciton states in a GaAs double quantum well (DQW). The structure consists of two quantum wells (QWs) coupled by a thin tunnel barrier. The DQW is subject to a transverse electric field to create spatially indirect inter-QW excitons with electrons and holes located in different QWs. We show that the presence of intra-QW charged excitons (trions) blocks carrier tunneling across the barrier to form indirect excitons, thus opening a gap in their emission spectrum. This behavior is attributed to the low binding energy of the trions. Within the tunneling blockade regime, emission becomes dominated by processes involving excitons bound to single shallow impurities, which behave as two-level centers activated by resonant tunneling. The quantum nature of the emission is confirmed by the anti-bunched photon emission statistics. The narrow distribution of emission energies ($\sim 10$~meV) and the electrical connection to the QWs make these single-exciton centers interesting candidates for applications in single-photon sources.


title: Basins of Attraction, Commitment Sets and Phenotypes of Boolean Networks

abstract: The attractors of Boolean networks and their basins have been shown to be highly relevant for model validation and predictive modelling, e.g., in systems biology. Yet there are currently very few tools available that are able to compute and visualise not only attractors but also their basins. In the realm of asynchronous, non-deterministic modeling not only is the repertoire of software even more limited, but also the formal notions for basins of attraction are often lacking. In this setting, the difficulty both for theory and computation arises from the fact that states may be ele- ments of several distinct basins. In this paper we address this topic by partitioning the state space into sets that are committed to the same attractors. These commitment sets can easily be generalised to sets that are equivalent w.r.t. the long-term behaviours of pre-selected nodes which leads us to the notions of markers and phenotypes which we illustrate in a case study on bladder tumorigenesis. For every concept we propose equivalent CTL model checking queries and an extension of the state of the art model checking software NuSMV is made available that is capa- ble of computing the respective sets. All notions are fully integrated as three new modules in our Python package PyBoolNet, including functions for visualising the basins, commitment sets and phenotypes as quotient graphs and pie charts.


title: Term Set Expansion based on Multi-Context Term Embeddings: an End-to-end Workflow

abstract: We present SetExpander, a corpus-based system for expanding a seed set of terms into a more complete set of terms that belong to the same semantic class. SetExpander implements an iterative end-to end workflow for term set expansion. It enables users to easily select a seed set of terms, expand it, view the expanded set, validate it, re-expand the validated set and store it, thus simplifying the extraction of domain-specific fine-grained semantic classes. SetExpander has been used for solving real-life use cases including integration in an automated recruitment system and an issues and defects resolution system. A video demo of SetExpander is available at https://drive.google.com/open?id=1e545bB87Autsch36DjnJHmq3HWfSd1Rv (some images were blurred for privacy reasons).


title: Global Existence and Ulam--Hyers Stability of $\Psi$--Hilfer Fractional Differential Equations

abstract: In this paper, we consider the Cauchy-type problem for a nonlinear differential equation involving $\Psi$-Hilfer fractional derivative and prove the existence and uniqueness of solutions in the weighted space of functions. The Ulam--Hyers and Ulam--Hyers--Rassias stability of Cauchy--type problem is investigated via successive approximation method. Further, we investigate the dependence of solutions on the initial conditions and uniqueness via $\epsilon$-approximated solution. An example is provided to illustrate the results we obtained.


title: Trend Analysis on the Metadata of Program Comprehension Papers

abstract: As program comprehension is a vast research area, it is necessary to get an overview of its rising and falling trends. We performed an n-gram frequency analysis on titles, abstracts and keywords of 1885 articles about program comprehension from the years 2000-2014. According to this analysis, the most rising trends are feature location and open source systems, the most falling ones are program slicing and legacy systems.


title: Transiting Exoplanet Monitoring Project (TEMP). I. Refined System Parameters and Transit Timing Variations of HAT-P-29b

abstract: We report the photometry of six transits of the hot Jupiter HAT-P-29b obtained from 2013 October to 2015 January. We analyze the new light curves, in combination with the published photometric, and Doppler velocimetric, and spectroscopic measurements, finding an updated orbital ephemeris for the HAT-P-29 system, $T_{\rm C}[0]= 2456170.5494(15)\,[\rm{BJD_{TDB}}]$ and $P=5.723390(13) \,{\rm days}$. It is $17.63\,{\rm s}$ ($4.0\,\sigma$) longer than the previously published value, amounting to errors exceeding $2.5\,\rm{hrs}$ at the time of writing (on UTC 2018 June 1). The measured transit mid-times for HAT-P-29b show no compelling evidence of timing anomalies from a linear model, which rules out the presence of a perturbers with masses greater than $0.6$, $0.7$, $0.5$, and $0.4\,{\rm M_\oplus}$ near the $1:2$, $2:3$, $3:2$, and $2:1$ resonances with HAT-P-29b, respectively.


title: Effects of Degradations on Deep Neural Network Architectures

abstract: Deep convolutional neural networks (CNN) have massively influenced recent advances in large-scale image classification. More recently, a dynamic routing algorithm with capsules (groups of neurons) has shown state-of-the-art recognition performance. However, the behavior of such networks in the presence of a degrading signal (noise) is mostly unexplored. An analytical study on different network architectures toward noise robustness is essential for selecting the appropriate model in a specific application scenario. This paper presents an extensive performance analysis of six deep architectures for image classification on six most common image degradation models. In this study, we have compared VGG-16, VGG-19, ResNet-50, Inception-v3, MobileNet and CapsuleNet architectures on Gaussian white, Gaussian color, salt-and-pepper, Gaussian blur, motion blur and JPEG compression noise models.


title: The Statistical Atom

abstract: This review is the updated and enlarged version of a talk delivered by J. S. on the occasion of the 1982 meeting of Nobel laureates at Lindau, and of talks given by B.-G. E. at several West German universities and Max Planck institutes in 1984.


title: ToriLLE: Learning Environment for Hand-to-Hand Combat

abstract: We present Toribash Learning Environment (ToriLLE), a learning environment for machine learning agents based on the video game Toribash. Toribash is a MuJoCo-like environment of two humanoid character fighting each other hand-to-hand, controlled by changing actuation modes of the joints. Competitive nature of Toribash as well its focused domain provide a platform for evaluating self-play methods, and evaluating machine learning agents against human players. In this paper we describe the environment with ToriLLE's capabilities and limitations, and experimentally show its applicability as a learning environment. The source code of the environment and conducted experiments can be found at https://github.com/Miffyli/ToriLLE.


title: MRI to FDG-PET: Cross-Modal Synthesis Using 3D U-Net For Multi-Modal Alzheimer's Classification

abstract: Recent studies suggest that combined analysis of Magnetic resonance imaging~(MRI) that measures brain atrophy and positron emission tomography~(PET) that quantifies hypo-metabolism provides improved accuracy in diagnosing Alzheimer's disease. However, such techniques are limited by the availability of corresponding scans of each modality. Current work focuses on a cross-modal approach to estimate FDG-PET scans for the given MR scans using a 3D U-Net architecture. The use of the complete MR image instead of a local patch based approach helps in capturing non-local and non-linear correlations between MRI and PET modalities. The quality of the estimated PET scans is measured using quantitative metrics such as MAE, PSNR and SSIM. The efficacy of the proposed method is evaluated in the context of Alzheimer's disease classification. The accuracy using only MRI is 70.18% while joint classification using synthesized PET and MRI is 74.43% with a p-value of $0.06$. The significant improvement in diagnosis demonstrates the utility of the synthesized PET scans for multi-modal analysis.


title: Spectra of Adjacency and Laplacian Matrices of Inhomogeneous Erd\H{o}s-R\'enyi Random Graphs

abstract: Inhomogeneous Erd\H{o}s-R\'enyi random graphs $\mathbb G_N$ on $N$ vertices in the non-dense regime are considered in this paper. The edge between the pair of vertices $\{i,j\}$ is retained with probability $\varepsilon_N\,f(\frac{i}{N},\frac{j}{N})$, $1 \leq i \neq j \leq N$, independently of other edges, where $f\colon\,[0,1] \times [0,1] \to [0,\infty)$ is a continuous function such that $f(x,y)=f(y,x)$ for all $x,y \in [0,1]$. We study the empirical distribution of both the adjacency matrix $A_N$ and the Laplacian matrix $\Delta_N$ associated with $\mathbb G_N$ in the limit as $N \to \infty$ when $\lim_{N\to\infty} \varepsilon_N = 0$ and $\lim_{N\to\infty} N\varepsilon_N = \infty$. In particular, it is shown that the empirical spectral distributions of $A_N$ and $\Delta_N$, after appropriate scaling and centering, converge to deterministic limits weakly in probability. For the special case where $f(x,y) = r(x)r(y)$ with $r\colon\,[0,1] \to [0,\infty)$ a continuous function, we give an explicit characterization of the limiting distributions. Furthermore, applications of the results to constrained random graphs, Chung-Lu random graphs and social networks are shown.


title: On compactifications and product-free sets

abstract: A subset of a group is said to be product-free if it does not contain three elements satisfying the equation $xy=z$. We give a negative answer to a question of Babai and S\'os on the existence of large product-free sets by model theoretic means. This question was originally answered by Gowers. Furthermore, we give a natural and sufficient model theoretic condition for a group to have a large product-free subset, as well as a model theoretic account of a result of Nikolov and Pyber on triple products.


title: The Evolution of Security Prices Is Not Stochastic but Governed by a Physicomathematical Law

abstract: Since Bachelier's thesis in 1900 (laying the foundation of the stochastic process, or Brownian motion, as a model of stock price changes), attempts at understanding the nature of prices and at predicting them have failed. Statistical methods have only found minor regularities/anomalies, and other mathematical and physical approaches do not work. This leads researchers to consider that the evolution of security prices is basically random, and, thus, inherently not predictable. We show that the evolution of security prices is not a stochastic process but largely deterministic and governed by a physical law. The law takes the form of a physicomathematical theory centered around a purely mathematical function, unrelated to models and statistical methods. It can be described as an "isodense" network of moving regression curves of an order greater than or equal to 1. The salient aspect of the function is that, when inputting a time series of any security into the function, new mathematical objects emerge spontaneously, and these objects exhibit the unique property of attracting and repelling the quantity. The graphical representation of the function is called a "topological network" due to the preeminence of shapes over metrics, and the emergent objects are called "characteristic figures" (mainly "cords"). The attraction and repulsion of the price by the cords results in the price bouncing from cord to cord. Thus, the price has to be considered as driven by the cords in a semi-deterministic manner (leaning towards deterministic). With a function that describes the evolution of the price, we now understand the reason behind each price movement and can also predict prices both qualitatively and quantitatively. The function is universal, does not rely on any fitting, and, due to its extreme sensitivity, reveals the hidden order in financial time series data that existing research never uncovered.


title: Key Borrowers Detection by Long-Range Interactions

abstract: We propose a new method for assessing agents' influence in financial network structures, which takes into consideration the intensity of interactions. A distinctive feature of this approach is that it considers not only direct interactions of agents of the first level and indirect interactions of the second level, but also long-range indirect interactions. At the same time we take into account the attributes of agents as well as the possibility of impact to a single agent from a group of other agents. This approach helps us to identify systemically important elements which cannot be detected by classical centrality measures or other indices. The proposed method was used to analyze the banking foreign claims for the end of 1Q 2015. Under the approach, two types of key borrowers were detected: a) major players with high ratings and positive credit history; b) intermediary players, which have a great scale of financial activities through the organization of favorable investment conditions and positive business climate.


title: Lattice sums for polyanalytic functions

abstract: In 1892, Lord Rayleigh estimated the effective conductivity of rectangular arrays of disks and proved, by means of the Eisenstein summation, that the lattice sum $S_2$ is equal to $\pi$ for the square array. Further, it became clear that such an equality can be treated as a necessary condition of the macroscopic isotropy of composites governed by the Laplace equation. This yielded the description of two-dimensional conducting composites by the classic elliptic functions including the conditionally convergent Eisenstein series. This paper is devoted to extension of the lattice sums to double periodic polyanalytic functions. The exact relations and computationally effective formulae between the polyanalytic and classic lattice sums are established. Polynomial representations of the lattice sums are obtained. They are a source of new exact formulae for the lattice sums where the number $\pi$ persists for macroscopically isotropic composites by our suggestion.


title: Effectiveness of Scaled Exponentially-Regularized Linear Units (SERLUs)

abstract: Recently, self-normalizing neural networks (SNNs) have been proposed with the intention to avoid batch or weight normalization. The key step in SNNs is to properly scale the exponential linear unit (referred to as SELU) to inherently incorporate normalization based on central limit theory. SELU is a monotonically increasing function, where it has an approximately constant negative output for large negative input. In this work, we propose a new activation function to break the monotonicity property of SELU while still preserving the self-normalizing property. Differently from SELU, the new function introduces a bump-shaped function in the region of negative input by regularizing a linear function with a scaled exponential function, which is referred to as a scaled exponentially-regularized linear unit (SERLU). The bump-shaped function has approximately zero response to large negative input while being able to push the output of SERLU towards zero mean statistically. To effectively combat over-fitting, we develop a so-called shift-dropout for SERLU, which includes standard dropout as a special case. Experimental results on MNIST, CIFAR10 and CIFAR100 show that SERLU-based neural networks provide consistently promising results in comparison to other 5 activation functions including ELU, SELU, Swish, Leakly ReLU and ReLU.


title: Interplay between low-lying isoscalar and isovector dipole modes: a comparative analysis between semi-classical and quantum approaches

abstract: We perform Time Dependent Hartree-Fock (TDHF) calculations to investigate the small amplitude dipole response of selected neutron-rich nuclei and Sn isotopes. A detailed comparison with the dipole strength predicted by Random-Phase Approximation (RPA) calculations is presented for the first time. TDHF results are also confronted to Vlasov calculations, to explore up to which extent a semi-classical picture can explain the properties of the nuclear response. The focus is on the low-energy response, below the Giant Dipole Resonance region, where different modes of non negligible strength are identified. We show that the relative weight of these excitations evolves with nuclear global features, such as density profile and neutron skin, which in turn reflect impor tant properties of the nuclear effective interaction. A thorough analysis of the associated transition densities turns out to be quite useful to better characterize the mixed isoscalar(IS)/isovector(IV) nature of the different modes and their surface/volume components. In particular, we show that the dipole response in the so-called Pygmy Dipole Resonance region corresponds to isoscalar-like surface oscillations, of larger strength in nuclei with a more diffuse surface. The ratio between the IV and IS Energy Weighted Sum Rule fractions exhausted in this region is shown to almost linearly increase with the neutron skin thickness in Sn isotopes.


title: A Unified Approximation Framework for Compressing and Accelerating Deep Neural Networks

abstract: Deep neural networks (DNNs) have achieved significant success in a variety of real world applications, i.e., image classification. However, tons of parameters in the networks restrict the efficiency of neural networks due to the large model size and the intensive computation. To address this issue, various approximation techniques have been investigated, which seek for a light weighted network with little performance degradation in exchange of smaller model size or faster inference. Both low-rankness and sparsity are appealing properties for the network approximation. In this paper we propose a unified framework to compress the convolutional neural networks (CNNs) by combining these two properties, while taking the nonlinear activation into consideration. Each layer in the network is approximated by the sum of a structured sparse component and a low-rank component, which is formulated as an optimization problem. Then, an extended version of alternating direction method of multipliers (ADMM) with guaranteed convergence is presented to solve the relaxed optimization problem. Experiments are carried out on VGG-16, AlexNet and GoogLeNet with large image classification datasets. The results outperform previous work in terms of accuracy degradation, compression rate and speedup ratio. The proposed method is able to remarkably compress the model (with up to 4.9x reduction of parameters) at a cost of little loss or without loss on accuracy.


title: SVD-based Visualisation and Approximation for Time Series Data in Smart Energy Systems

abstract: Many time series in smart energy systems exhibit two different timescales. On the one hand there are patterns linked to daily human activities. On the other hand, there are relatively slow trends linked to seasonal variations. In this paper we interpret these time series as matrices, to be visualized as images. This approach has two advantages: First of all, interpreting such time series as images enables one to visually integrate across the image and makes it therefore easier to spot subtle or faint features. Second, the matrix interpretation also grants elucidation of the underlying structure using well-established matrix decomposition methods. We will illustrate both these aspects for data obtained from the German day-ahead market.


title: A new criterion on k-normal elements over finite fields

abstract: The notion of normal elements for finite fields extension has been generalized as k-normal elements by Huczynska et al. [3]. The number of k-normal elements for a fixed finite field extension has been calculated and estimated [3], and several methods to construct k-normal elements have been presented [1,3]. Several criteria on k-normal element have been given [1,2]. In this paper we present a new criterion on k-normal elements by using idempotents and show some examples. Such criterion has been given for usual normal element before [6].


title: On the number of monic admissible polynomials in the ring $\mathbb{Z}[x]$

abstract: In this paper we study admissible polynomials. We establish an estimate for the number of admissible polynomials of degree $n$ with coeffients $a_i$ satisfying $0\leq a_i\leq H$ for a fixed $H$, for $i=0,1,2, \ldots, n-1$. In particular, letting $\mathcal{N}(H)$ denotes the number of monic admissible polynomials of degree $n\geq 3$ with coefficients satisfying the inequality $0\leq a_i\leq H$, we show that \begin{align}\frac{H^{n-1}}{(n-1)!}+O(H^{n-2})\leq \mathcal{N}(H) \leq \frac{n^{n-1}H^{n-1}}{(n-1)!}+O(H^{n-2}).\nonumber \end{align} Also letting $\mathcal{A}(H)$ denotes the number of monic irreducible admissible polynomials, with coefficients satisfying the same condition , we show that \begin{align}\mathcal{A}(H)\geq \frac{H^{n-1}}{(n-1)!}+O\bigg( H^{n-4/3}(\log H)^{2/3}\bigg).\nonumber \end{align}


title: Well-Posedness for the Two-dimensional Zakharov-Kuznetsov Equation

abstract: We show the global well-posedness for the two-dimensional Zakharov-Kuznetsov equation in $H^{s}({\mathbb{R}^2})$ when $\frac{11}{13}<s<1$ via the I-method. Additionally, local well-posedness for the symmetrized ZK equation in $ B^\frac{1}{2}_{2,1}(\mathbb{R}^2)$ is established by using atomic spaces.


title: Interacting particles with L\'{e}vy strategies: limits of transport equations for swarm robotic systems

abstract: L\'{e}vy robotic systems combine superdiffusive random movement with emergent collective behaviour from local communication and alignment in order to find rare targets or track objects. In this article we derive macroscopic fractional PDE descriptions from the movement strategies of the individual robots. Starting from a kinetic equation which describes the movement of robots based on alignment, collisions and occasional long distance runs according to a L\'{e}vy distribution, we obtain a system of evolution equations for the fractional diffusion for long times. We show that the system allows efficient parameter studies for a search problem, addressing basic questions like the optimal number of robots needed to cover an area in a certain time. For shorter times, in the hyperbolic limit of the kinetic equation, the PDE model is dominated by alignment, irrespective of the long range movement. This is in agreement with previous results in swarming of self-propelled particles. The article indicates the novel and quantitative modeling opportunities which swarm robotic systems provide for the study of both emergent collective behaviour and anomalous diffusion, on the respective time scales.


title: Chudnovsky's formula for $1/\pi$ revisited

abstract: The document contains an outline of a modular proof for Ramanujan-Chudnovsky identity.


title: Time dependent simulation of the flow reduction of D$_2$ and T$_2$ in the KATRIN experiment

abstract: The KArlsruhe TRItium Neutrino experiment (KATRIN) aims to measure the effective electron anti-neutrino mass with an unprecedented sensitivity of 0.2 eV/c$^2$, using beta-electrons from tritium decay. Super-conducting magnets will guide the electrons through a vacuum beam-line from the windowless gaseous tritium source through differential and cryogenic pumping sections to a high resolution spectrometer. At the same time tritium gas has to be prevented from entering the spectrometer. Therefore, the pumping sections have to reduce the tritium flow by at least 14 orders of magnitude. This paper describes various simulation methods in the molecular flow regime used to determine the expected gas flow reduction in the pumping sections for deuterium (commissioning runs) and for radioactive tritium. Simulations with MolFlow+ and with an analytical model are compared with each other, and with the stringent requirements of the KATRIN experiment.


title: Do quantum propositions obey the principle of excluded middle?

abstract: The present paper demonstrates the failure of the principle of excluded middle (PEM) in the lattice of all closed linear subspaces of a Hilbert space (usually defined as quantum logic). Namely, it is shown that for a qubit, a proposition and its negation can be both false. Since PEM is the assumed theorem of quantum logic, this raises the question: If PEM holds in the orthocomplemented lattice of all propositions of the quantum system, then how the failure of PEM in quantum logic can be explained? Alternatively, if the propositions relating to the quantum system do not obey PEM, then what is the semantics of those propositions? Possible answers to these questions are analyzed in the present paper.


title: Joint Queue-Aware and Channel-Aware Delay Optimal Scheduling of Arbitrarily Bursty Traffic over Multi-State Time-Varying Channels

abstract: This paper is motivated by the observation that the average queueing delay can be decreased by sacrificing power efficiency in wireless communications. In this sense, we naturally wonder what is the minimum queueing delay when the available power is limited and how to achieve the minimum queueing delay. To answer these two questions in the scenario where randomly arriving packets are transmitted over multi-state wireless fading channel, a probabilistic cross-layer scheduling policy is proposed in this paper, and characterized by a constrained Markov Decision Process (MDP). Using the steady-state probability of the underlying Markov chain, we are able to derive the mathematical expressions of the concerned metrics, namely, the average queueing delay and the average power consumption. To describe the delay-power tradeoff, we formulate a non-linear programming problem, which, however, is very challenging to solve. By analyzing its structure, this optimization problem can be converted into an equivalent Linear Programming (LP) problem via variable substitution, which allows us to derive the optimal delay-power tradeoff as well as the optimal scheduling policy. The optimal scheduling policy turns out to be dual-threshold-based, which means transmission decisions should be made based on the optimal thresholds imposed on the queue length and the channel state.


title: A Benchmark of Selected Algorithmic Differentiation Tools on Some Problems in Computer Vision and Machine Learning

abstract: Algorithmic differentiation (AD) allows exact computation of derivatives given only an implementation of an objective function. Although many AD tools are available, a proper and efficient implementation of AD methods is not straightforward. The existing tools are often too different to allow for a general test suite. In this paper, we compare fifteen ways of computing derivatives including eleven automatic differentiation tools implementing various methods and written in various languages (C++, F#, MATLAB, Julia and Python), two symbolic differentiation tools, finite differences, and hand-derived computation. We look at three objective functions from computer vision and machine learning. These objectives are for the most part simple, in the sense that no iterative loops are involved, and conditional statements are encapsulated in functions such as {\tt abs} or {\tt logsumexp}. However, it is important for the success of algorithmic differentiation that such `simple' objective functions are handled efficiently, as so many problems in computer vision and machine learning are of this form. Of course, our results depend on programmer skill, and familiarity with the tools. However, we contend that this paper presents an important datapoint: a skilled programmer devoting roughly a week to each tool produced the timings we present. We have made our implementations available as open source to allow the community to replicate and update these benchmarks.


title: Bestow and Atomic: Concurrent Programming using Isolation, Delegation and Grouping

abstract: Any non-trivial concurrent system warrants synchronisation, regardless of the concurrency model. Actor-based concurrency serialises all computations in an actor through asynchronous message passing. In contrast, lock-based concurrency serialises some computations by following a lock--unlock protocol for accessing certain data. Both systems require sound reasoning about pointers and aliasing to exclude data-races. If actor isolation is broken, so is the single-thread-of-control abstraction. Similarly for locks, if a datum is accessible outside of the scope of the lock, the datum is not governed by the lock. In this paper we discuss how to balance aliasing and synchronisation. In previous work, we defined a type system that guarantees data-race freedom of actor-based concurrency and lock-based concurrency. This paper extends this work by the introduction of two programming constructs; one for decoupling isolation and synchronisation and one for constructing higher-level atomicity guarantees from lower-level synchronisation. We focus predominantly on actors, and in particular the Encore programming language, but our ultimate goal is to define our constructs in such a way that they can be used both with locks and actors, given that combinations of both models occur frequently in actual systems. We discuss the design space, provide several formalisations of different semantics and discuss their properties, and connect them to case studies showing how our proposed constructs can be useful. We also report on an on-going implementation of our proposed constructs in Encore.


title: Bridge trisections in $\mathbb{CP}^2$ and the Thom conjecture

abstract: In this paper, we develop new techniques for understanding surfaces in $\mathbb{CP}^2$ via bridge trisections. Trisections are a novel approach to smooth 4-manifold topology, introduced by Gay and Kirby, that provide an avenue to apply 3-dimensional tools to 4-dimensional problems. Meier and Zupan subsequently developed the theory of bridge trisections for smoothly embedded surfaces in 4-manifolds. The main application of these techniques is a new proof of the Thom conjecture, which posits that algebraic curves in $\mathbb{CP}^2$ have minimal genus among all smoothly embedded, oriented surfaces in their homology class. This new proof is notable as it completely avoids any gauge theory or pseudoholomorphic curve techniques.


title: Proper effective temperature of nonequilibrium steady state

abstract: We define a proper effective temperature for relativistic nonequilibrium steady states (NESSs). A conventional effective temperature of NESSs is defined from the ratio of the fluctuation to the dissipation. However, NESSs have relative velocities to the heat bath in general, and hence the conventional effective temperature can be frame dependent in relativistic systems. The proper effective temperature is introduced as a frame-independent (Lorentz invariant) quantity that characterizes NESSs. We find that the proper effective temperature of NESSs is higher than the proper temperature of the heat bath in a wide range of holographic models even when the conventional effective temperature is lower than the temperature of the heat bath.

categories: hep-th cond-mat.stat-mech gr-qc hep-ph nucl-th


title: Positive Integer Solutions of Certain Diophantine Equations Involving Lucas-Balancing numbers

abstract: In this article the Diophantine equations of the form $x^{2}-8C_{n}xy+ 16y^{2}=\pm2^{r}$, $x^{4}-8C_{n}x^{2}y+ 16y^{2}=\pm2^{r}$ and $x^{2}-8C_{n}xy^{2}+ 16y^{4}=\pm2^{r}$ are taken into consideration for the investigations of existence of any positive integer solutions. Furthermore, the solutions are manifested in terms of the balancing, Pell and Lucas-Pell numbers.


title: Analytic Geometry of Homogeneous Spaces

abstract: The theory uses methods and language of linear algebra to study nonlinear spaces. These techniques can be used particularly to describe analytic geometry of non-linear elliptic, hyperbolic, De Sitter and Anti de Sitter spaces. The main innovation of elaborated theory is space parameterization by introduction of space signature. This parameterization allows studying of different homogeneous spaces in one global framework. When the parameters are used as variables in definitions, axioms, equations, theorems, proofs, all these have exactly the same form that describes the reality of all homogeneous spaces simultaneously. When it is necessary to describe some space particularities or to see the difference between two concrete spaces, the concrete values can be put in parameters of each definition, axiom, equation, theorem and proof.


title: On the nodal set of solutions to a class of nonlocal parabolic equations

abstract: We investigate the local properties, including the nodal set and the nodal properties of solutions to the following parabolic problem of Muckenhoupt-Neumann type: \begin{equation*} \begin{cases} \partial_t \overline{u} - y^{-a} \nabla \cdot(y^a \nabla \overline{u}) = 0 \quad &\text{ in } \mathbb{B}_1^+ \times (-1,0) \\ -\partial_y^a \overline{u} = q(x,t)u \quad &\text{ on } B_1 \times \{0\} \times (-1,0), \end{cases} \end{equation*} where $a\in(-1,1)$, is a fixed parameter $\mathbb{B}_1^+\subset \mathbb{R}^{N+1}$ is the upper unit half ball and $B_1$ is the unit ball in $\mathbb{R}^N$. Our main motivation comes from its relation with a class of nonlocal parabolic equations involving the fractional power of the heat operator \begin{equation*} H^su(x,t) = \frac{1}{|\Gamma(-s)|} \int_{-\infty}^t \int_{\mathbb{R}^N} \left[u(x,t) - u(z,\tau)\right] \frac{G_N(x-z,t-\tau)}{(t-\tau)^{1+s}} dzd\tau. \end{equation*} We characterise the possible blow-ups and we examine the structure of the nodal set of solutions vanishing with a finite order. More precisely, we prove that the nodal set has at least parabolic Hausdorff codimension one in $\mathbb{R}^N\times\mathbb{R}$, and can be written as the union of a locally smooth part and a singular part, which turns out to possess remarkable stratification properties. Moreover, the asymptotic behaviour of general solutions near their nodal points is classified in terms of a class of explicit polynomials of Hermite and Laguerre type, obtained as eigenfunctions to an Ornstein-Uhlenbeck type operator. Our main results are obtained through a fine blow-up analysis which relies on the monotonicity of an Almgren-Poon type quotient and some new Liouville type results for parabolic equations, combined with more classical results including Federer's reduction principle and the parabolic Whitney's extension.


title: Semiclassical calculation of spectral correlation functions of chaotic systems

abstract: We present a semiclassical approach to n-point spectral correlation functions of quantum systems whose classical dynamics is chaotic, for arbitrary n. The basic ingredients are sets of periodic orbits that have nearly the same action and therefore provide constructive interference. We calculate explicitly the first correlation functions, to leading orders in their energy arguments, for both unitary and orthogonal symmetry classes. The results agree with corresponding predictions from random matrix theory, thereby giving solid support to the conjecture of universality.


title: A simple device for microinjections, manipulations and measurements using an electromorphological chip under microinterferometric control of the interface and membrane processes at the thickness range of 5-1000 nm at different angles

abstract: Micromanipulations, perfusions and measurements performed using glass microelectrodes filled with an electrolyte is a conventional technique for experimental morphological and membrane electrophysiological studies at a single cell and membrane surface level. The typical (effective) diameter of the end of the glass microelectrode is from 500 up to less than 100 nm, which prevents one from observing it using a standard optical microscope in accordance with the optical resolution criteria, since the diameter less than 500 nm is indistinguishable within the interference zone. Microprocessor programming of the puller (microforge) that provides pulling and tearing allows to obtain in certain regimes the adjusted diameter and shape of the micropipette tip, although this result is not fully controlled due to the above limitations. In this connection it is necessary to design the control devices for the micropipette tips both at the preparation and operation stages (intracellular or extracellular insertion). This method also should provide visualization of the processes occurring upon interaction of the microelectrode tip with the cell in real time, depending on the electrode type and state, which allows to level the artifacts arising with the systematic error frequency from the uncontrolled operation of the micropipette tip after different ways of the microelectrode filling with the electrolyte. We propose an installation scheme that solves the above problems by means of introducing an interferometric device for microscopic control of the microelectrode and micromanipulator or microperfusor, for the first time for a given type of optical instruments combined with the interferometric optical scheme.

categories: q-bio.SC cond-mat.soft q-bio.CB q-bio.NC


title: Block models for multipartite networks.Applications in ecology and ethnobiology

abstract: Modeling relations between individuals is a classical question in social sciences, ecology, etc. In order to uncover a latent structure in the data, a popular approach consists in clustering individuals according to the observed patterns of interactions. To do so, Stochastic block models (SBM) and Latent Block models (LBM) are standard tools for clustering the individuals with respect to their comportment in a unique network. However, when adopting an integrative point of view, individuals are not involved in a unique network but are part of several networks, resulting into a potentially complex multipartite network. In this contribution, we propose a stochastic block model able to handle multipartite networks, thus supplying a clustering of the individuals based on their connection behavior in more than one network. Our model is an extension of the latent block models (LBM) and stochastic block model (SBM). The parameters -- such as the marginal probabilities of assignment to blocks and the matrix of probabilities of connections between blocks -- are estimated through a variational Expectation-Maximization procedure. The numbers of blocks are chosen with the Integrated Completed Likelihood criterion, a penalized likelihood criterion. The pertinence of our methodology is illustrated on two datasets issued from ecology and ethnobiology.


title: Can graphene bilayers be the membrane mimetic materials? "Ion channels" in graphene-based nanostructures

abstract: The prospects of application of graphene and related structures as the membrane mimetic materials, capable of reproducing several biomembrane functions up to the certain limit, are analyzed in the series of our papers. This paper considers the possibility of the ion channel function modeling using graphene and its derivatives. The physical mechanisms providing selective permeability for different membrane mimetic materials, as well as the limits of the adequate simulation of the transport, catalytic, sensing and electrogenic properties of the cell membrane ion channels using bilayered graphene-based structures are discussed.


title: Electromagnetic effect on anisotropic scalar field collapse in higher curvature gravity

abstract: We consider a `Scalar-Maxwell-Einstein-Gauss-Bonnet' theory in four dimension, where the scalar field couples non-minimally with the Gauss-Bonnet (GB) term. This coupling with the scalar field ensures the non topological character of the GB term. In such higher curvature scenario, we explore the effect of electromagnetic field on scalar field collapse. Our results reveal that the presence of a time dependent electromagnetic field requires an anisotropy in the background spacetime geometry and such anisotropic spacetime allows a collapsing solution for the scalar field. The singularity formed as a result of the collapse is found to be a curvature singularity which may be point like or line like depending on the strength of the anisotropy. We also show that the singularity is always hidden from exterior by an apparent horizon.