title: Exceeding the asymptotic limit of polymer drag reduction

abstract: The drag of turbulent flows can be drastically decreased by addition of small amounts of high molecular weight polymers. While drag reduction initially increases with polymer concentration, it eventually saturates to what is known as the maximum drag reduction (MDR) asymptote; this asymptote is generally attributed to the dynamics being reduced to a marginal yet persistent state of subdued turbulent motion. Contrary to this accepted view we will show in the following that for an appropriate choice of parameters polymers can reduce the drag beyond the suggested asymptotic limit, eliminating turbulence and giving way to laminar flow. However at higher polymer concentrations the laminar state becomes unstable, resulting in a fluctuating flow with the characteristic drag of the MDR asymptote. The asymptotic state is hence dynamically disconnected from ordinary turbulence.


title: An Automated Auto-encoder Correlation-based Health-Monitoring and Prognostic Method for Machine Bearings

abstract: This paper studies an intelligent ultimate technique for health-monitoring and prognostic of common rotary machine components, particularly bearings. During a run-to-failure experiment, rich unsupervised features from vibration sensory data are extracted by a trained sparse auto-encoder. Then, the correlation of the extracted attributes of the initial samples (presumably healthy at the beginning of the test) with the succeeding samples is calculated and passed through a moving-average filter. The normalized output is named auto-encoder correlation-based (AEC) rate which stands for an informative attribute of the system depicting its health status and precisely identifying the degradation starting point. We show that AEC technique well-generalizes in several run-to-failure tests. AEC collects rich unsupervised features form the vibration data fully autonomous. We demonstrate the superiority of the AEC over many other state-of-the-art approaches for the health monitoring and prognostic of machine bearings.


title: Distributed Stochastic Model Predictive Control for Large-Scale Linear Systems with Private and Common Uncertainty Sources

abstract: This paper presents a distributed stochastic model predictive control (SMPC) approach for large-scale linear systems with private and common uncertainties in a plug-and-play framework. Using the so-called scenario approach, the centralized SMPC involves formulating a large-scale finite-horizon scenario optimization problem at each sampling time, which is in general computationally demanding, due to the large number of required scenarios. We present two novel ideas in this paper to address this issue. We first develop a technique to decompose the large-scale scenario program into distributed scenario programs that exchange a certain number of scenarios with each other in order to compute local decisions using the alternating direction method of multipliers (ADMM). We show the exactness of the decomposition with a-priori probabilistic guarantees for the desired level of constraint fulfillment for both uncertainty sources. As our second contribution, we develop an inter-agent soft communication scheme based on a set parametrization technique together with the notion of probabilistically reliable sets to reduce the required communication between the subproblems. We show how to incorporate the probabilistic reliability notion into existing results and provide new guarantees for the desired level of constraint violations. Two different simulation studies of two types of systems interactions, dynamically coupled and coupling constraints, are presented to illustrate the advantages of the proposed framework.


title: Hybrid correlation holography with a single pixel detector

abstract: The correlation holography reconstructs 3D objects as a distribution of two-point correlation of the random field detected by two dimensional detector arrays. Here, we describe a hybrid method, a combination of optical and computational channels, to reconstruct the objects from only a single pixel detector. An experimental arrangement is proposed and as a first step to realize the technique, we have simulated the experimental model for both scalar and vectorial objects. Proposed technique affords background free imaging with 3D capability.


title: Evolving Game Skill-Depth using General Video Game AI Agents

abstract: Most games have, or can be generalised to have, a number of parameters that may be varied in order to provide instances of games that lead to very different player experiences. The space of possible parameter settings can be seen as a search space, and we can therefore use a Random Mutation Hill Climbing algorithm or other search methods to find the parameter settings that induce the best games. One of the hardest parts of this approach is defining a suitable fitness function. In this paper we explore the possibility of using one of a growing set of General Video Game AI agents to perform automatic play-testing. This enables a very general approach to game evaluation based on estimating the skill-depth of a game. Agent-based play-testing is computationally expensive, so we compare two simple but efficient optimisation algorithms: the Random Mutation Hill-Climber and the Multi-Armed Bandit Random Mutation Hill-Climber. For the test game we use a space-battle game in order to provide a suitable balance between simulation speed and potential skill-depth. Results show that both algorithms are able to rapidly evolve game versions with significant skill-depth, but that choosing a suitable resampling number is essential in order to combat the effects of noise.


title: Droplet Ripening in Concentration Gradients

abstract: Living cells use phase separation and concentration gradients to organize chemical compartments in space. Here, we present a theoretical study of droplet dynamics in gradient systems. We derive the corresponding growth law of droplets and find that droplets exhibit a drift velocity and position dependent growth. As a consequence, the dissolution boundary moves through the system, thereby segregating droplets to one end. We show that for steep enough gradients, the ripening leads to a transient arrest of droplet growth that is induced by an narrowing of the droplet size distribution.


title: Unsupervised Learning of Mixture Regression Models for Longitudinal Data

abstract: This paper is concerned with learning of mixture regression models for individuals that are measured repeatedly. The adjective "unsupervised" implies that the number of mixing components is unknown and has to be determined, ideally by data driven tools. For this purpose, a novel penalized method is proposed to simultaneously select the number of mixing components and to estimate the mixing proportions and unknown parameters in the models. The proposed method is capable of handling both continuous and discrete responses by only requiring the first two moment conditions of the model distribution. It is shown to be consistent in both selecting the number of components and estimating the mixing proportions and unknown regression parameters. Further, a modified EM algorithm is developed to seamlessly integrate model selection and estimation. Simulation studies are conducted to evaluate the finite sample performance of the proposed procedure. And it is further illustrated via an analysis of a primary biliary cirrhosis data set.


title: A Gaia-PS1-SDSS (GPS1) Proper Motion Catalog Covering 3/4 of the Sky

abstract: We combine Gaia DR1, PS1, SDSS and 2MASS astrometry to measure proper motions for 350 million sources across three-fourths of the sky down to a magnitude of $m_r\sim20$\,. Using positions of galaxies from PS1, we build a common reference frame for the multi-epoch PS1, single-epoch SDSS and 2MASS data, and calibrate the data in small angular patches to this frame. As the Gaia DR1 excludes resolved galaxy images, we choose a different approach to calibrate its positions to this reference frame: we exploit the fact that the proper motions of stars in these patches are {\it linear}. By simultaneously fitting the positions of stars at different epochs of -- Gaia DR1, PS1, SDSS, and 2MASS -- we construct an extensive catalog of proper motions dubbed GPS1. GPS1 has a characteristic systematic error of less than 0.3 \masyr\, and a typical precision of $ 1.5-2.0$\masyr. The proper motions have been validated using galaxies, open clusters, distant giant stars and QSOs. In comparison with other published faint proper motion catalogs, GPS1's systematic error ($<0.3$ \masyr) should be nearly an order of magnitude better than that of PPMXL and UCAC4 ($>2.0$ \masyr). Similarly, its precision ($\sim 1.5$ \masyr) is a four-fold improvement relative to PPMXL and UCAC4 ($\sim 6.0$ \masyr). For QSOs, the precision of GPS1 is found to be worse ($\sim 2.0-3.0$\masyr), possibly due to their particular differential chromatic refraction (DCR). The GPS1 catalog will be released on-line and available via the VizieR Service and VO Service. (===GPS1 is available with VO TAP Query now, see http://www2.mpia-hd.mpg.de/~tian/GPS1/ for details=== )


title: First- and Second-Order Hypothesis Testing for Mixed Memoryless Sources with General Mixture

abstract: The first- and second-order optimum achievable exponents in the simple hypothesis testing problem are investigated. The optimum achievable exponent for type II error probability, under the constraint that the type I error probability is allowed asymptotically up to epsilon, is called the epsilon-optimum exponent. In this paper, we first give the second-order epsilon-exponent in the case where the null hypothesis and the alternative hypothesis are a mixed memoryless source and a stationary memoryless source, respectively. We next generalize this setting to the case where the alternative hypothesis is also a mixed memoryless source. We address the first-order epsilon-optimum exponent in this setting. In addition, an extension of our results to more general setting such as the hypothesis testing with mixed general source and the relationship with the general compound hypothesis testing problem are also discussed.


title: Coexistence of non-Fermi liquid and Fermi liquid self-energies at all dopings in cuprates

abstract: Various angle-dependent measurements in hole-doped cuprates suggested that Non-Fermi liquid (NFL) and Fermi-liquid (FL) self-energies coexist in the Brillouin zone. Moreover, it is also found that NFL self-energies survive up to the overdoped region where the resistivity features a global FL-behavior. To address this problem, here we compute the momentum dependent self-energy from a single band Hubbard model. The self-energy is calculated self-consistently by using a momentum-dependent density-fluctuation (MRDF) method. One of our main result is that the computed self-energy exhibits a NFL-like frequency dependence only in the antinodal region, and FL-like behavior elsewhere, and retains its analytic form at all momenta and dopings. The dominant source of NFL self-energy in the antinodal region stems from the self-energy-dressed fluctuations between the itinerant and localized densities as self-consistency is invoked. We also calculate the DC conductivity by including the full momentum dependent self-energy. We find that the resistivity-temperature exponent n becomes 1 near the optimal doping, while the NFL self-energy occupies largest momentum-space volume. Surprisingly, even in the NFL state near the optimal doping, the nodal region contains FL-like self-energies; while in the under- and over-dopings (n ~ 2), the antinodal region remains NFL-like. These results highlight the non-local correlation physics in cuprates and in other similar intermediately correlated materials, where a direct link between the microscopic single-particle spectral properties and the macroscopic transport behavior can not be well established.


title: $\beta$-decay properties of waiting point nuclei for astrophysical applications

abstract: We report microscopic calculation of key $\beta$-decay properties for some of the crucial waiting point species having neutron closed magic shells 50 and 82. Our calculation bear astrophysical significance vis-\'{a}-vis speeding of the $r$-process. The $\beta$-decay properties include half-lives, energy rates of $\beta$-delayed neutrons and their emission probabilities, both under terrestrial and stellar conditions. We perform a pn-QRPA calculation with a separable multi-shell interaction and include both allowed and unique first-forbidden (U1F) transitions in our calculation. We compare our results with previous calculations and measured data. Our calculation is in good agreement with the experimental data. For certain cases, we noted a significant decrease in the half-life calculation with the contribution of U1F transitions. This is in contradiction to the shell model study where only for $N$ = 126 waiting-point nuclei, the forbidden transitions were reported to significantly reduce the calculated half-lives. Our model fulfills the Ikeda sum rule for even-even cases. For odd-A cases the rule is violated by up to 15$\%$ for $^{125}$Tc.


title: Comparing the dark matter models, modified Newtonian dynamics and modified gravity in accounting for the galaxy rotation curves

abstract: We compare six models (including the baryonic model, two dark matter models, two modified Newtonian dynamics models and one modified gravity model) in accounting for the galaxy rotation curves. For the dark matter models, we assume NFW profile and core-modified profile for the dark halo, respectively. For the modified Newtonian dynamics models, we discuss Milgrom's MOND theory with two different interpolation functions, i.e. the standard and the simple interpolation functions. As for the modified gravity, we focus on Moffat's MSTG theory. We fit these models to the observed rotation curves of 9 high-surface brightness and 9 low-surface brightness galaxies. We apply the Bayesian Information Criterion and the Akaike Information Criterion to test the goodness-of-fit of each model. It is found that non of the six models can well fit all the galaxy rotation curves. Two galaxies can be best fitted by the baryonic model without involving the nonluminous dark matter. MOND can fit the largest number of galaxies, and only one galaxy can be best fitted by MSTG model. Core-modified model can well fit about one half LSB galaxies but no HSB galaxy, while NFW model can fit only a small fraction of HSB galaxies but no LSB galaxy. This may imply that the oversimplified NFW and Core-modified profiles couldn't well mimic the postulated dark matter halo.


title: Expecting the Unexpected: Training Detectors for Unusual Pedestrians with Adversarial Imposters

abstract: As autonomous vehicles become an every-day reality, high-accuracy pedestrian detection is of paramount practical importance. Pedestrian detection is a highly researched topic with mature methods, but most datasets focus on common scenes of people engaged in typical walking poses on sidewalks. But performance is most crucial for dangerous scenarios, such as children playing in the street or people using bicycles/skateboards in unexpected ways. Such "in-the-tail" data is notoriously hard to observe, making both training and testing difficult. To analyze this problem, we have collected a novel annotated dataset of dangerous scenarios called the Precarious Pedestrian dataset. Even given a dedicated collection effort, it is relatively small by contemporary standards (around 1000 images). To allow for large-scale data-driven learning, we explore the use of synthetic data generated by a game engine. A significant challenge is selected the right "priors" or parameters for synthesis: we would like realistic data with poses and object configurations that mimic true Precarious Pedestrians. Inspired by Generative Adversarial Networks (GANs), we generate a massive amount of synthetic data and train a discriminative classifier to select a realistic subset, which we deem the Adversarial Imposters. We demonstrate that this simple pipeline allows one to synthesize realistic training data by making use of rendering/animation engines within a GAN framework. Interestingly, we also demonstrate that such data can be used to rank algorithms, suggesting that Adversarial Imposters can also be used for "in-the-tail" validation at test-time, a notoriously difficult challenge for real-world deployment.


title: Multi-talker Speech Separation with Utterance-level Permutation Invariant Training of Deep Recurrent Neural Networks

abstract: In this paper we propose the utterance-level Permutation Invariant Training (uPIT) technique. uPIT is a practically applicable, end-to-end, deep learning based solution for speaker independent multi-talker speech separation. Specifically, uPIT extends the recently proposed Permutation Invariant Training (PIT) technique with an utterance-level cost function, hence eliminating the need for solving an additional permutation problem during inference, which is otherwise required by frame-level PIT. We achieve this using Recurrent Neural Networks (RNNs) that, during training, minimize the utterance-level separation error, hence forcing separated frames belonging to the same speaker to be aligned to the same output stream. In practice, this allows RNNs, trained with uPIT, to separate multi-talker mixed speech without any prior knowledge of signal duration, number of speakers, speaker identity or gender. We evaluated uPIT on the WSJ0 and Danish two- and three-talker mixed-speech separation tasks and found that uPIT outperforms techniques based on Non-negative Matrix Factorization (NMF) and Computational Auditory Scene Analysis (CASA), and compares favorably with Deep Clustering (DPCL) and the Deep Attractor Network (DANet). Furthermore, we found that models trained with uPIT generalize well to unseen speakers and languages. Finally, we found that a single model, trained with uPIT, can handle both two-speaker, and three-speaker speech mixtures.


title: Symmetric powers of permutation representations of finite groups and primitive colorings on polyhedrons

abstract: In this paper, we define a set which has a finite group action and is generated by a finite color set, a set which has a finite group action, and a subset of the set of non negative integers. we state its properties to apply one of solution of the following two problems, respectively. First, we calculate the generating function of the character of symmetric powers of permutation representation associated with a set which has a finite group action. Second, we calculate the number of primitive colorings on some objects of polyhedrons. It is a generalization of the calculation of the number of primitive necklaces by N.Metropolis and G-C.Rota.


title: A Characterization of the Graphs to Induce Periodic Grover Walk

abstract: This paper explains the periodicity of the Grover walk on finite graphs. We characterize the graphs to induce 2, 3, 4, 5-periodic Grover walk and obtain a necessary condition of the graphs to induce an odd-periodic Grover walk.


title: Differential operators and contravariant derivatives in Poisson geometry

abstract: We inquire into the relation between the curl operators, the Poisson coboundary operators and contravariant derivatives on Poisson manifolds to study the theory of differential operators in Poisson geometry. Given an oriented Poisson manifold, we describe locally those two differential operators in terms of Poisson connection whose torsion is vanishing. Moreover, we introduce the notion of the modular operator for an oriented Poisson manifold. For a symplectic manifold, we describe explicitly the modular operator in terms of the curvature 2-section of Poisson connection, analogously to the Weitzenb$\ddot{\rm o}$ck formula in Riemannian geometry.


title: Gender Matters! Analyzing Global Cultural Gender Preferences for Venues Using Social Sensing

abstract: Gender differences is a phenomenon around the world actively researched by social scientists. Traditionally, the data used to support such studies is manually obtained, often through surveys with volunteers. However, due to their inherent high costs because of manual steps, such traditional methods do not quickly scale to large-size studies. We here investigate a particular aspect of gender differences: preferences for venues. To that end we explore the use of check-in data collected from Foursquare to estimate cultural gender preferences for venues in the physical world. For that, we first demonstrate that by analyzing the check-in data in various regions of the world we can find significant differences in preferences for specific venues between gender groups. Some of these significant differences reflect well-known cultural patterns. Moreover, we also gathered evidence that our methodology offers useful information about gender preference for venues in a given region in the real world. This suggests that gender and venue preferences observed may not be independent. Our results suggests that our proposed methodology could be a promising tool to support studies on gender preferences for venues at different spatial granularities around the world, being faster and cheaper than traditional methods, besides quickly capturing changes in the real world.


title: Constraints and Implications on Higgs FCNC Couplings from Precision Measurement of $B_s \to \mu^+\mu^-$ Decay

abstract: We study constraints and implications of the recent LHCb measurement of ${\cal B}(B_s \to \mu^+\mu^-)$ for tree-level Higgs-mediated flavor-changing neutral current (FCNC) interactions. Combined with experimental data on $B_s$ mass difference $\Delta m_s$, the $h \to \mu \tau$, and the $h \to \tau^+\tau^-$ decay branching ratios from the LHC, we find that the Higgs FCNC couplings are severely constrained. The allowed regions for $B_s \to \mu \tau$, $\tau\tau$ and $h \to sb$ decays are obtained. Current data allow large CP violation in the $h \to \tau^+ \tau^-$ decay. Consequences of the Cheng-Sher ansatz for the Higgs Yukawa couplings are discussed in some detail.


title: A wake-sleep algorithm for recurrent, spiking neural networks

abstract: We investigate a recently proposed model for cortical computation which performs relational inference. It consists of several interconnected, structurally equivalent populations of leaky integrate-and-fire (LIF) neurons, which are trained in a self-organized fashion with spike-timing dependent plasticity (STDP). Despite its robust learning dynamics, the model is susceptible to a problem typical for recurrent networks which use a correlation based (Hebbian) learning rule: if trained with high learning rates, the recurrent connections can cause strong feedback loops in the network dynamics, which lead to the emergence of attractor states. This causes a strong reduction in the number of representable patterns and a decay in the inference ability of the network. As a solution, we introduce a conceptually very simple "wake-sleep" algorithm: during the wake phase, training is executed normally, while during the sleep phase, the network "dreams" samples from its generative model, which are induced by random input. This process allows us to activate the attractor states in the network, which can then be unlearned effectively by an anti-Hebbian mechanism. The algorithm allows us to increase learning rates up to a factor of ten while avoiding clustering, which allows the network to learn several times faster. Also for low learning rates, where clustering is not an issue, it improves convergence speed and reduces the final inference error.


title: The effect of coherent stirring on the advection-condensation of water vapour

abstract: Atmospheric water vapour is an essential ingredient of weather and climate. Key features of its distribution can be represented by kinematic models which treat it as a passive scalar advected by a prescribed flow and reacting through condensation. Condensation acts as a sink that maintains specific humidity below a prescribed, space-dependent saturation value. In order to investigate how the interplay between large-scale advection, small-scale turbulence and condensation controls the moisture distribution, we develop simple kinematic models which combine a single circulating flow with a Brownian-motion representation of turbulence. We first study the drying mechanism of a water-vapour anomaly released inside a vortex at an initial time. Next, we consider a cellular flow with a moisture source at a boundary. The statistically steady state attained shows features reminiscent of the Hadley cell such as boundary layers, a region of intense precipitation and a relative humidity minimum. Explicit results provide a detailed characterisation of these features in the limit of strong flow.


title: Hydrodynamic limit for the Ginzburg-Landau $\nabla\phi$ interface model with non-convex potential

abstract: Hydrodynamic limit for the Ginzburg-Landau $\nabla\phi$ interface model was established in [Nishikawa, 2003] under the Dirichlet boundary conditions. This paper studies the similar problem, but with non-convex potentials. Because of the lack of strict convexity, a lot of difficulties arise, especially, on the identification of equilibrium states. We give a proof of the equivalence between the stationarity and the Gibbs property under quite general settings, and as its conclusion, we complete the identification of equilibrium states under the high temparature regime in [Deuschel and Cotar, 2008]. We also establish some uniform estimates for variances of extremal Gibbs measures under quite general settings.


title: Software Metrics in Boa Large-Scale Software Mining Infrastructure: Challenges and Solutions

abstract: In this paper, we describe our experience implementing some of classic software engineering metrics using Boa - a large-scale software repository mining platform - and its dedicated language. We also aim to take an advantage of the Boa infrastructure to propose new software metrics and to characterize open source projects by software metrics to provide reference values of software metrics based on large number of open source projects. Presented software metrics, well known and proposed in this paper, can be used to build large-scale software defect prediction models. Additionally, we present the obstacles we met while developing metrics, and our analysis can be used to improve Boa in its future releases. The implemented metrics can also be used as a foundation for more complex explorations of open source projects and serve as a guide how to implement software metrics using Boa as the source code of the metrics is freely available to support reproducible research.


title: Collective atomic-population-inversion and stimulated radiation for two-componentBose-Einstein condensate in an optical cavity

abstract: In this paper we investigate the ground-state properties and related quantum phase transitions for the two-component Bose-Einstein condensate in a single-mode optical cavity. Apart from the usual normal and superradiant phases multi-stable macroscopic quantum states are realized by means of the spin-coherent-state variational method. We demonstrate analytically the stimulated radiation from collective state of atomic population inversion, which does not exist in the normal Dicke model with single-component atoms. It is also revealed that the stimulated radiation can be generated only from one component of atoms and the other remains in the ordinary superradiant state. However the order of superradiant and stimulatedradiation states is interchangeable between two components of atoms by tuning the relative atom-field couplings and the frequency detuning as well.


title: An almost complex Chern-Ricci flow

abstract: We consider the evolution of an almost Hermitian metric by the $(1,1)$ part of its Chern-Ricci form on almost complex manifolds. This is an evolution equation first studied by Chu and coincides with the Chern-Ricci flow if the complex structure is integrable and with the K\"ahler-Ricci flow if moreover the initial metric is K\"ahler. We find the maximal existence time for the flow in term of the initial data and also give some convergence results. As an example, we study this flow on the (locally) homogeneous manifolds in more detail.


title: Geometric RTT realization of $U_v(gl_n)^+$

abstract: In this paper, we give a BLM realization of the positive part of the quantum group of $U_v(gl_n)$ with respect to RTT relations.


title: Testing Holographic Conjectures of Complexity with Born-Infeld Black Holes

abstract: In this paper, we use Born-Infeld black holes to test two recent holographic conjectures of complexity, the "Complexity = Action" (CA) duality and "Complexity = Volume 2.0" (CV) duality. The complexity of a boundary state is identified with the action of the Wheeler-deWitt patch in CA duality, while this complexity is identified with the spacetime volume of the WdW patch in CV duality. In particular, we check whether the Born-Infeld black holes violate the Lloyd bound: $\mathcal{\dot{C}\leq}\frac{2}{\pi\hbar}\left[ \left( M-Q\Phi\right) -\left( M-Q\Phi\right) _{\text{gs}}\right] $, where gs stands for the ground state for a given electrostatic potential. We find that the ground states are either some extremal black hole or regular spacetime with nonvanishing charges. Near extremality, the Lloyd bound is violated in both dualities. Near the charged regular spacetime, this bound is satisfied in CV duality but violated in CA duality. When moving away from the ground state on a constant potential curve, the Lloyd bound tend to be saturated from below in CA duality, while $\mathcal{\dot{C}}$ is $\pi/2$ times as large as the Lloyd bound in CV duality.


title: Search for High-energy Neutrinos from Gravitational Wave Event GW151226 and Candidate LVT151012 with ANTARES and IceCube

abstract: The Advanced LIGO observatories detected gravitational waves from two binary black hole mergers during their first observation run (O1). We present a high-energy neutrino follow-up search for the second gravitational wave event, GW151226, as well as for gravitational wave candidate LVT151012. We find 2 and 4 neutrino candidates detected by IceCube, and 1 and 0 detected by ANTARES, within $\pm500$ s around the respective gravitational wave signals, consistent with the expected background rate. None of these neutrino candidates are found to be directionally coincident with GW151226 or LVT151012. We use non-detection to constrain isotropic-equivalent high-energy neutrino emission from GW151226 adopting the GW event's 3D localization, to less than $2\times 10^{51}-2\times10^{54}$ erg.


title: Extension of local smooth maps of Banach spaces

abstract: It is known that smooth bump functions are absent in the majority of infinite-dimensional Banach spaces. This is an obstacle in the development of local analysis, in particular in the questions of extending local maps onto the whole space. We suggest an approach that substitutes bump functions with special maps, which we call K-maps. It allows us to extend smooth local maps from non-smooth spaces, such as $C^q[0,1], q=0,1,...$. We also prove the Borel lemma for spaces possessing K-maps.


title: Defect prediction with bad smells in code

abstract: Background: Defect prediction in software can be highly beneficial for development projects, when prediction is highly effective and defect-prone areas are predicted correctly. One of the key elements to gain effective software defect prediction is proper selection of metrics used for dataset preparation. Objective: The purpose of this research is to verify, whether code smells metrics, collected using Microsoft CodeAnalysis tool, added to basic metric set, can improve defect prediction in industrial software development project. Results: We verified, if dataset extension by the code smells sourced metrics, change the effectiveness of the defect prediction by comparing prediction results for datasets with and without code smells-oriented metrics. In a result, we observed only small improvement of effectiveness of defect prediction when dataset extended with bad smells metrics was used: average accuracy value increased by 0.0091 and stayed within the margin of error. However, when only use of code smells based metrics were used for prediction (without basic set of metrics), such process resulted with surprisingly high accuracy (0.8249) and F-measure (0.8286) results. We also elaborated data anomalies and problems we observed when two different metric sources were used to prepare one, consistent set of data. Conclusion: Extending the dataset by the code smells sourced metric does not significantly improve the prediction effectiveness. Achieved result did not compensate effort needed to collect additional metrics. However, we observed that defect prediction based on the code smells only is still highly effective and can be used especially where other metrics hardly be used.


title: The rate of period change in DAV stars

abstract: Grids of DAV star models are evolved by \texttt{WDEC}, taking the element diffusion effect into account. The grid parameters are hydrogen mass log($M_{H}/M_{*}$), helium mass log($M_{He}/M_{*}$), stellar mass $M_{\rm *}$, and effective temperature $T_{\rm eff}$ for DAV stars. The core compositions are from white dwarf models evolved by \texttt{MESA}. Therefore, those DAV star models evolved by \texttt{WDEC} have historically viable core compositions. Based on those DAV star models, we studied the rate of period change ($\dot{P}(k)$) for different values of H, He, $M_{\rm *}$, and $T_{\rm eff}$. The results are consistent with previous work. Two DAV stars G117-B15A and R548 have been observed around forty years. The rates of period change of two large-amplitude modes were obtained through O-C method. We did asteroseismological study on the two DAV stars and then obtained a best-fitting model for each star. Based on the two best-fitting models, the mode identifications ($l$, $k$) of the observed modes for G117-B15A and R548 are consistent with previous work. Both the observed modes and the observed $\dot{P}$s can be fitted by calculated ones. The results indicate that our method of evolving DAV star models is feasible.


title: Analytical expressions for thermophysical properties of solid and liquid tungsten relevant for fusion applications

abstract: The status of the literature is reviewed for several thermophysical properties of pure solid and liquid tungsten which constitute input for the modelling of intense plasma-surface interaction phenomena that are important for fusion applications. Reliable experimental data are analyzed for the latent heat of fusion, the electrical resistivity, the specific isobaric heat capacity, the thermal conductivity and the mass density from the room temperature up to the boiling point of tungsten as well as for the surface tension and the dynamic viscosity across the liquid state. Analytical expressions of high accuracy are recommended for these thermophysical properties that involved a minimum degree of extrapolations. In particular, extrapolations were only required for the surface tension and viscosity.


title: A Fast Algorithm for a Weighted Low Rank Approximation

abstract: Matrix low rank approximation including the classical PCA and the robust PCA (RPCA) method have been applied to solve the background modeling problem in video analysis. Recently, it has been demonstrated that a special weighted low rank approximation of matrices can be made robust to the outliers similar to the $\ell_1$-norm in RPCA method. In this work, we propose a new algorithm that can speed up the existing algorithm for solving the special weighted low rank approximation and demonstrate its use in background estimation problem.


title: Extended Calculations of Spectroscopic Data: Energy Levels, Lifetimes and Transition rates for O-like ions from Cr XVII to Zn XXIII

abstract: Employing two state-of-the-art methods, multiconfiguration Dirac--Hartree--Fock and second-order many-body perturbation theory, the excitation energies and lifetimes for the lowest 200 states of the $2s^2 2p^4$, $2s 2p^5$, $2p^6$, $2s^2 2p^3 3s$, $2s^2 2p^3 3p$, $2s^2 2p^3 3d$, $2s 2p^4 3s$, $2s 2p^4 3p$, and $2s 2p^4 3d$ configurations, and multipole (electric dipole (E1), magnetic dipole (M1), and electric quadrupole (E2)) transition rates, line strengths, and oscillator strengths among these states are calculated for each O-like ion from Cr XVII to Zn XXIII. Our two data sets are compared with the NIST and CHIANTI compiled values, and previous calculations. The data are accurate enough for identification and deblending of new emission lines from the sun and other astrophysical sources. The amount of data of high accuracy is significantly increased for the $n = 3$ states of several O-like ions of astrophysics interest, where experimental data are very scarce.


title: Hardness of almost embedding simplicial complexes in $\mathbb R^d$

abstract: A map $f\colon K\to \mathbb R^d$ of a simplicial complex is an almost embedding if $f(\sigma)\cap f(\tau)=\emptyset$ whenever $\sigma,\tau$ are disjoint simplices of $K$. Theorem. Fix integers $d,k\ge2$ such that $d=\frac{3k}2+1$. (a) Assume that $P\ne NP$. Then there exists a finite $k$-dimensional complex $K$ that does not admit an almost embedding in $\mathbb R^d$ but for which there exists an equivariant map $\tilde K\to S^{d-1}$. (b) The algorithmic problem of recognition almost embeddability of finite $k$-dimensional complexes in $\mathbb R^d$ is NP hard. The proof is based on the technique from the Matou\v{s}ek-Tancer-Wagner paper (proving an analogous result for embeddings), and on singular versions of the higher-dimensional Borromean rings lemma and a generalized van Kampen--Flores theorem.


title: The ladder physics in the Spin Fermion model

abstract: A link is established between the spin-fermion (SF) model of the cuprates and the approach based on the analogy between the physics of doped Mott insulators in two dimensions and the physics of fermionic ladders. This enables one to use nonperturbative results derived for fermionic ladders to move beyond the large-N approximation in the SF model. It is shown that the paramagnon exchange postulated in the SF model has exactly the right form to facilitate the emergence of the fully gapped d-Mott state in the region of the Brillouin zone at the hot spots of the Fermi surface. Hence the SF model provides an adequate description of the pseudogap.


title: Definition of geometric space around analytic fractal trees using derivative coordinate funtions

abstract: The concept of derivative coordinate functions proved useful in the formulation of analytic fractal functions to represent smooth symmetric binary fractal trees [1]. In this paper we introduce a new geometry that defines the fractal space around these fractal trees. We present the canonical and degenerate form of this fractal space and extend the fractal geometrical space to R3 specifically and Rn by a recurrence relation. We also discuss the usage of such fractal geometry.


title: On the construction of Wannier functions in topological insulators: the 3D case

abstract: We investigate the possibility of constructing exponentially localized composite Wannier bases, or equivalently smooth periodic Bloch frames, for 3-dimensional time-reversal symmetric topological insulators, both of bosonic and of fermionic type, so that the bases in question are also compatible with time-reversal symmetry. This problem is translated in the study, of independent interest, of homotopy classes of continuous, periodic, and time-reversal symmetric families of unitary matrices. We identify three $\mathbb{Z}_2$-valued complete invariants for these homotopy classes. When these invariants vanish, we provide an algorithm which constructs a "multi-step" logarithm that is employed to continuously deform the given family into a constant one, identically equal to the identity matrix. This algorithm leads to a constructive procedure to produce the composite Wannier bases mentioned above.

categories: math-ph cond-mat.mes-hall cond-mat.mtrl-sci math.MP


title: Interference Between Source-Free Radiation and Radiation from Sources: Particle-Like Behavior for Classical Radiation

abstract: A simple junior-level electrodynamics problem is used to illustrate the interference between a source-free standing plane wave and a wave generated by a pulse in a current sheet. Depending upon the relative phases between the standing wave and the current pulse and also upon the relative magnitudes, we can find quite different patterns of emitted energy and momentum. If the source gives a large radiation pulse so that the source-free plane wave can be neglected, then the radiation spreads out symmetrically on either side of the current sheet. However, if the radiation sheet gives a pulse with fields comparable to those of the standing wave, then we can find a single radiation pulse moving to the right while the current sheet recoils to the left, or the situation with the directions reversed. The example is a crude illustration of particle-like behavior arising from conventional classical electromagnetic behavior in the presence of source-free radiation. The discussion makes contact with the ideas of photons in modern physics.


title: Nitrogen-vacancy centers created by N$^+$ ion implantation through screening SiO$_2$ layers on diamond

abstract: We report on an ion implantation technique utilizing a screening mask made of SiO$_2$ to control both the depth profile and the dose. By appropriately selecting the thickness of the screening layer, this method fully suppresses the ion channeling, brings the location of the highest NV density to the surface, and effectively reduces the dose by more than three orders of magnitude. With a standard ion implantation system operating at the energy of 10 keV and the dose of 10$^{11}$ cm$^2$ and without an additional etching process, we create single NV centers close to the surface with coherence times of a few tens of $\mu$s.


title: Imitating the winner leads to discrimination in spatial prisoner's dilemma model

abstract: The occurrence of discrimination is an important problem in the social and economical sciences. Much of the discrimination observed in empirical studies can be explained by the theory of in-group favoritism, which states that people tend to act more positively towards peers whose appearances are more similar to their own. Some studies, however, find hierarchical structures in inter-group relations, where members of low-status groups also favor the high-status group members. These observations cannot be understood in the light of in-group favoritism. Here we present an agent based model in which evolutionary dynamics can result in a hierarchical discrimination between two groups characterized by a meaningless, but observable binary label. We find that discriminating strategies end up dominating the system when the selection pressure is high, i.e. when agents have a much higher probability of imitating their neighbor with the highest payoff. These findings suggest that the puzzling persistence of hierarchical discrimination may result from the evolutionary dynamics of the social system itself, namely the social imitation dynamics. It also predicts that discrimination will occur more often in highly competitive societies.


title: Construction of constant scalar curvature K\"ahler cone metrics

abstract: Over a compact K\"ahler manifold, we provide a Fredholm alternative result for the Lichnerowicz operator associated to a K\"ahler metric with conic singularities along a divisor. We deduce several existence results of constant scalar curvature K\"ahler metrics with conic singularities: existence result under small deformations of K\"ahler classes, existence result over a Fano manifold, existence result over certain ruled manifolds. In this last case, we consider the projectivisation of a parabolic stable holomorphic bundle. This leads us to prove that the existing Hermitian-Einstein metric on this bundle enjoys a regularity property along the divisor on the base.


title: Constrained Spacecraft Relative Motion Planning Exploiting Periodic Natural Motion Trajectories and Invariance

abstract: Spacecraft relative motion planning is concerned with the design and execution of maneuvers relative to a nominal target. These types of maneuvers are frequently utilized in missions such as rendezvous and docking, satellite inspection and formation flight where exclusion zones representing spacecraft or other obstacles must be avoided. The presence of these exclusion zones leads to non-linear and non-convex constraints which must be satisfied. In this paper, a novel approach to spacecraft relative motion planning with obstacle avoidance and thrust constraints is developed. This approach is based on a graph search applied to a virtual net of closed (periodic) natural motion trajectories, where the natural motion trajectories represent virtual net nodes (vertices), and adjacency and connection information is determined by conditions defined in terms of safe, positively-invariant tubes built around each trajectory. These conditions guarantee that transitions from one natural motion trajectory to another natural motion trajectory can be completed without constraint violations. The proposed approach improves the flexibility of a previous approach based on the use of forced equilibria, and has other advantages in terms of reduced fuel consumption and passive safety. The resulting maneuvers, if planned on-board, can be executed directly or, if planned off board, can be used to warm start trajectory optimizers to generate further improvements.


title: Representability of Lyndon-Maddux relation algebras

abstract: In Alm-Hirsch-Maddux (2016), relation algebras $\mathfrak{L}(q,n)$ were defined that generalize Roger Lyndon's relation algebras from projective lines, so that $\mathfrak{L}(q,0)$ is a Lyndon algebra. In that paper, it was shown that if $q>2304n^2+1$, $\mathfrak{L}(q,n)$ is representable, and if $q<2n$, $\mathfrak{L}(q,n)$ is not representable. In the present paper, we reduced this gap by proving that if $q\geq n(\log n)^{1+\varepsilon}$, $\mathfrak{L}(q,n)$ is representable.


title: Resonance a_1(1420) and the Three-Pion Decays of the Tauon

abstract: The role of the a_1(1420) resonance in the three-pion decays of the tau lepton is investigated using a phenomenological model. For all data before 2008, roughly equal fit quality is achieved when the basic a_1(1260) resonance is supplemented with either a_1(1640) or a_1(1420). However, two recent and more precise data sets require resonances with masses that are not very far from that of a_1(1420). This suggests that the axial-vector resonance that accompanies the a_1(1260) in the three-pion decays of the tauon is a_1(1420), not a_1(1640), as believed up to now. More data are needed to demonstrate this definitely.


title: On the linear polarization constants of finite dimensional spaces

abstract: We study the linear polarization constants of finite dimensional Banach spaces. We obtain the correct asymptotic behaviour of these constants for the spaces $\ell_p^d$: they behave as $\sqrt[p]{d}$ if $1\le p\le 2$ and as $\sqrt{d}$ if $2\le p<\infty$. For $p=\infty$ we get the asymptotic behavior up to a logarithmic factor


title: Designing for older adults: review of touchscreen design guidelines

abstract: The distinct abilities of older adults to interact with computers has motivated a wide range of contributions in the the form of design guidelines for making technologies usable and accessible for the elderly population. However, despite the growing effort by the research community, the adoption of guidelines by developers and designers has been scant or not properly translated into more accessible interaction systems. In this paper we explore this issue by reporting on a qualitative outcomes of a systematic review of 204 research-derived design guidelines for touchscreen applications. We report first on the different definitions of "elderly" and assess the reliability, organization and accessibility of the guidelines. Then we present our early attempt at facilitating the reporting and access of such guidelines to researchers and practitioners.


title: Can one determine the neutrino mass by electron capture?

abstract: There are three different methods used to search the neutrino mass: - The electron antineutrino mass can probably best be determined by the Triton decay. - The neutrinoless Double Beta Decay yields information, if the neutrino is a Dirac or a Majorana particle. It can also determine the Majorana neutrino mass. - Electron capture of an atomic bound electron by a proton in a nucleus bound electron plus proton to neutron plus electron-neutrino can give the mass of the electron neutrino. This contribution summarizes our theoretical work on the possibility to determine the electron neutrino mass by electron capture. One expects the largest influence of the neutrino mass on this decay for a small Q = 2.8 keV for electron capture in Holmium. The energy of the Q value is distributed to the emitted neutrino and the excitation of the Dy atom. Thus the energy difference between the Q value and the upper end of the deexcitation spectrum is the electron neutrino mass. The excitation spectrum of Dy is calculate by one-, two- and three-electron hole excitations, and by the shake-off process. The electron wave functions are calculated selfconsistently by the Dirac-Hartree-Fock approach for the bound and the continuum states. To extract the neutrino mass from the spectrum one must adjust simultaneously the neutrino mass, the Q value, the position, the relative strength and the width of the highest resonance. This fit is only possible, if the background is reduced relative to the present situation. In case of a drastically reduced background a fit of the Q-value and the neutrino mass only seems also to be possible. The analysis presented here shows, that the determination of the electron neutrino mass by electron capture is difficult, but seems not to be impossible.


title: Precise LIGO Lensing Rate Predictions for Binary Black Holes

abstract: We show how LIGO is expected to detect coalescing binary black holes at $z>1$, that are lensed by the intervening galaxy population. Gravitational magnification, $\mu$, strengthens gravitational wave signals by $\sqrt{\mu}$, without altering their frequencies, which if unrecognised leads to an underestimate of the event redshift and hence an overestimate of the binary mass. High magnifications can be reached for coalescing binaries because the region of intense gravitational wave emission during coalescence is so small ($\sim$100km), permitting very close projections between lensing caustics and gravitational-wave events. Our simulations incorporate accurate waveforms convolved with the LIGO power spectral density. Importantly, we include the detection dependence on sky position and orbital orientation, which for the LIGO configuration translates into a wide spread in observed redshifts and chirp masses. Currently we estimate a detectable rate of lensed events \rateEarly{}, that rises to \rateDesign{}, at LIGO's design sensitivity limit, depending on the high redshift rate of black hole coalescence.


title: Hardware-Efficient Schemes of Quaternion Multiplying Units for 2D Discrete Quaternion Fourier Transform Processors

abstract: In this paper, we offer and discuss three efficient structural solutions for the hardware-oriented implementation of discrete quaternion Fourier transform basic operations with reduced implementation complexities. The first solution: a scheme for calculating sq product, the second solution: a scheme for calculating qt product, and the third solution: a scheme for calculating sqt product, where s is a so-called i-quaternion, t is an j-quaternion, and q is an usual quaternion. The direct multiplication of two usual quaternions requires 16 real multiplications (or two-operand multipliers in the case of fully parallel hardware implementation) and 12 real additions (or binary adders). At the same time, our solutions allow to design the computation units, which consume only 6 multipliers plus 6 two input adders for implementation of sq or qt basic operations and 9 binary multipliers plus 6 two-input adders and 4 four-input adders for implementation of sqt basic operation.


title: Solving the Goddard problem by an influence diagram

abstract: Influence diagrams are a decision-theoretic extension of probabilistic graphical models. In this paper we show how they can be used to solve the Goddard problem. We present results of numerical experiments with this problem and compare the solutions provided by influence diagrams with the optimal solution.


title: An empirical analysis of smart contracts: platforms, applications, and design patterns

abstract: Smart contracts are computer programs that can be consistently executed by a network of mutually distrusting nodes, without the arbitration of a trusted authority. Because of their resilience to tampering, smart contracts are appealing in many scenarios, especially in those which require transfers of money to respect certain agreed rules (like in financial services and in games). Over the last few years many platforms for smart contracts have been proposed, and some of them have been actually implemented and used. We study how the notion of smart contract is interpreted in some of these platforms. Focussing on the two most widespread ones, Bitcoin and Ethereum, we quantify the usage of smart contracts in relation to their application domain. We also analyse the most common programming patterns in Ethereum, where the source code of smart contracts is available.


title: Robust and scalable domain decomposition solvers for unfitted finite element methods

abstract: Unfitted finite element methods, e.g., extended finite element techniques or the so-called finite cell method, have a great potential for large scale simulations, since they avoid the generation of body-fitted meshes and the use of graph partitioning techniques, two main bottlenecks for problems with non-trivial geometries. However, the linear systems that arise from these discretizations can be much more ill-conditioned, due to the so-called small cut cell problem. The state-of-the-art approach is to rely on sparse direct methods, which have quadratic complexity and are thus not well suited for large scale simulations. In order to solve this situation, in this work we investigate the use of domain decomposition preconditioners (balancing domain decomposition by constraints) for unfitted methods. We observe that a straightforward application of these preconditioners to the unfitted case has a very poor behavior. As a result, we propose a {customization} of the classical BDDC methods based on the stiffness weighting operator and an improved definition of the coarse degrees of freedom in the definition of the preconditioner. These changes lead to a robust and algorithmically scalable solver able to deal with unfitted grids. A complete set of complex 3D numerical experiments show the good performance of the proposed preconditioners.


title: Deep Tensor Encoding

abstract: Learning an encoding of feature vectors in terms of an over-complete dictionary or a information geometric (Fisher vectors) construct is wide-spread in statistical signal processing and computer vision. In content based information retrieval using deep-learning classifiers, such encodings are learnt on the flattened last layer, without adherence to the multi-linear structure of the underlying feature tensor. We illustrate a variety of feature encodings incl. sparse dictionary coding and Fisher vectors along with proposing that a structured tensor factorization scheme enables us to perform retrieval that can be at par, in terms of average precision, with Fisher vector encoded image signatures. In short, we illustrate how structural constraints increase retrieval fidelity.


title: Adaptive Local (AL) Basis for Elliptic Problems with $L^\infty$-Coefficients

abstract: We define a generalized finite element method for the discretization of elliptic partial differential equations in heterogeneous media. An adaptive local finite element basis (AL basis) on a coarse mesh which does not resolve the matrix of the media is constructed by solving finite-dimensional localized problems. The method requires $O(log(1/H)^{d+1})$ basis functions per mesh point. We prove that the optimal finite element convergence rates are preserved.


title: Static trace free Einstein equations and stellar distributions

abstract: We construct models of static spherical distributions of perfect fluid in trace--free Einstein gravity theory. The equations governing the gravitational field are equivalent to the standard Einstein's equations however, their presentation is manifestly different which motivates the question whether new information would emerge due to the nonlinearity of the field equations. The incompressible fluid assumption does not lead to the well known Schwarzschild interior metric of Einstein gravity and a term denoting the presence of a cosmological constant is present on account of the integration process. The Schwarzschild interior is regained as a special case of a richer geometry. On the other hand, when the Schwarzschild geometry is prescribed, a constant density fluid emerges consistent with the standard equations. A complete model of an isothermal fluid sphere with pressure and density obeying the inverse square law is obtained. Corrections to the model previously presented in the literature by Saslaw {\it {et al}} are exhibited. The isothermal ansatz does not yield a constant gravitational potential in general but both potentials are position dependent. Conversely, it is shown that assuming a constant $g_{rr}$ gravitational potential does not yield an isothermal fluid in general as is the case in standard general relativity. The results of the standard Einstein equations are special cases of the models reported here. Noteworthy is the fact that whereas the previously reported isothermal solution was only of cosmological interest, the solution reported herein admit compact objects by virtue of the fact that a pressure-free hypersurface exists. Finally we analyze the consequences of selecting the Finch--Skea metric as the seed solution.


title: Spectrum Estimation from a Few Entries

abstract: Singular values of a data in a matrix form provide insights on the structure of the data, the effective dimensionality, and the choice of hyper-parameters on higher-level data analysis tools. However, in many practical applications such as collaborative filtering and network analysis, we only get a partial observation. Under such scenarios, we consider the fundamental problem of recovering spectral properties of the underlying matrix from a sampling of its entries. We are particularly interested in directly recovering the spectrum, which is the set of singular values, and also in sample-efficient approaches for recovering a spectral sum function, which is an aggregate sum of the same function applied to each of the singular values. We propose first estimating the Schatten $k$-norms of a matrix, and then applying Chebyshev approximation to the spectral sum function or applying moment matching in Wasserstein distance to recover the singular values. The main technical challenge is in accurately estimating the Schatten norms from a sampling of a matrix. We introduce a novel unbiased estimator based on counting small structures in a graph and provide guarantees that match its empirical performance. Our theoretical analysis shows that Schatten norms can be recovered accurately from strictly smaller number of samples compared to what is needed to recover the underlying low-rank matrix. Numerical experiments suggest that we significantly improve upon a competing approach of using matrix completion methods.


title: Functional Central Limit Theorem For Susceptible-Infected Process On Configuration Model Graphs

abstract: We study a stochastic compartmental susceptible-infected (SI) epidemic process on a configuration model random graph with a given degree distribution over a finite time interval $[0,T],$ for some $ T>0$. In this setting, we split the population of graph nodes into two compartments, namely, $S$ and $I$, denoting the susceptible and infected nodes, respectively. In addition to the sizes of these two compartments, we study counts of $SI$-edges (those connecting a susceptible and an infected node) and $SS$-edges (those connecting two susceptible nodes). We describe the dynamical process in terms of these counts and present a functional central limit theorem (FCLT) for them, a scaling limit of the dynamical process as $n$, the number of nodes in the random graph, grows to infinity. To be precise, we show that these counts, when appropriately scaled, converge weakly to a continuous Gaussian vector martingale process the usual Skorohod space of real 3-dimensional vector-valued \cadlag\, functions on $[0,T]$ endowed with the Skorohod topology. We assume certain technical requirements for this purpose. We discuss applications of our FCLT in percolation theory (from a non-equilibrium statistical mechanics point of view), and in computer science in the context of spread of computer viruses. We also provide simulation results for some common degree distributions.


title: G2 instantons and the Seiberg-Witten monopoles

abstract: I describe a relation (mostly conjectural) between the Seiberg-Witten monopoles, Fueter sections, and G2 instantons. In the last part of this article I gathered some open questions connected with this relation.


title: Measurement of the jet mass in highly boosted t t-bar events from pp collisions at sqrt(s) = 8 TeV

abstract: The first measurement of the jet mass m[jet] of top quark jets produced in t t-bar events from pp collisions at sqrt(s) = 8 TeV is reported for the jet with the largest transverse momentum pt in highly boosted hadronic top quark decays. The data sample, collected with the CMS detector, corresponds to an integrated luminosity of 19.7 inverse femtobarns. The measurement is performed in the lepton+jets channel in which the products of the semileptonic decay t to bW with W to l nu where l is an electron or muon, are used to select t t-bar events with large Lorentz boosts. The products of the fully hadronic decay t to bW with W to q q'-bar are reconstructed using a single Cambridge-Aachen jet with distance parameter R = 1.2, and pt > 400 GeV. The t t-bar cross section as a function of m[jet] is unfolded at the particle level and is used to test the modelling of highly boosted top quark production. The peak position of the m[jet] distribution is sensitive to the top quark mass m[t], and the data are used to extract a value of m[t] to assess this sensitivity.


title: On 3+1 Lorentzian Einstein Manifolds with One Rotational Isometry

abstract: We consider 3+1 rotationally symmetric Lorentzian Einstein spacetime manifolds with $\Lambda >0$ and reduce the equations to 2+1 Einstein equations coupled to `shifted' wave maps. Subsequently, we prove various (explicit) positive mass-energy theorems. No smallness is assumed.


title: On The Uniqueness of SRB Measures for Endomorphisms

abstract: In this paper we improve the results of \cite{MT} and show that a weak hyperbolic transitivity implies the uniqueness of hyperbolic SRB measures. As an important corollary, it arises the ergodicity of the system in a conservative setting. It also arises the condition which implies the stable ergodicity as well as the statistical stability for a general $C^2$-regular map.


title: Sharp estimates for the gradient of the generalized Poisson integral for a half-space

abstract: A representation of the sharp coefficient in a pointwise estimate for the gradient of the generalized Poisson integral of a function $f$ on ${\mathbb R}^n$ is obtained under the assumption that $f$ belongs to $L^p$. The explicit value of the coefficient is found for the cases $p=1$ and $p=2$.


title: Successive Coordinate Search and Component-by-Component Construction of Rank-1 Lattice Rules

abstract: The (fast) component-by-component (CBC) algorithm is an efficient tool for the construction of generating vectors for quasi-Monte Carlo rank-1 lattice rules in weighted reproducing kernel Hilbert spaces. We consider product weights, which assigns a weight to each dimension. These weights encode the effect a certain variable (or a group of variables by the product of the individual weights) has. Smaller weights indicate less importance. Kuo (2003) proved that the CBC algorithm achieves the optimal rate of convergence in the respective function spaces, but this does not imply the algorithm will find the generating vector with the smallest worst-case error. In fact it does not. We investigate a generalization of the component-by-component construction that allows for a general successive coordinate search (SCS), based on an initial generating vector, and with the aim of getting closer to the smallest worst-case error. The proposed method admits the same type of worst-case error bounds as the CBC algorithm, independent of the choice of the initial vector. Under the same summability conditions on the weights as in [Kuo,2003] the error bound of the algorithm can be made independent of the dimension $d$ and we achieve the same optimal order of convergence for the function spaces from [Kuo,2003]. Moreover, a fast version of our method, based on the fast CBC algorithm by Nuyens and Cools, is available, reducing the computational cost of the algorithm to $O(d \, n \log(n))$ operations, where $n$ denotes the number of function evaluations. Numerical experiments seeded by a Korobov-type generating vector show that the new SCS algorithm will find better choices than the CBC algorithm and the effect is better when the weights decay slower.


title: Invariant measures, matching and the frequency of 0 for signed binary expansions

abstract: We introduce a parametrised family of maps $\{S_{\eta}\}_{\eta \in [1,2]}$, called symmetric doubling maps, defined on $[-1,1]$ by $S_\eta (x)=2x-d\eta$, where $d\in \{-1,0,1 \}$. Each map $S_\eta$ generates binary expansions with digits $-1$, 0 and 1. We study the frequency of the digit 0 in typical expansions as a function of the parameter $\eta$. The transformations $S_\eta$ have a natural ergodic invariant measure $\mu_\eta$ that is absolutely continuous with respect to Lebesgue measure. The frequency of the digit 0 is related to the measure $\mu_{\eta}([-\frac12,\frac12])$ by the Ergodic Theorem. We show that the density of $\mu_\eta$ is piecewise smooth except for a set of parameters of zero Lebesgue measure and full Hausdorff dimension and give a full description of the structure of the maximal parameter intervals on which the density is piecewise smooth. We give an explicit formula for the frequency of the digit 0 in typical signed binary expansions on each of these parameter intervals and show that this frequency depends continuously on the parameter $\eta$. Moreover, it takes the value $\frac23$ only on the interval $\big[ \frac65, \frac32\big]$ and it is strictly less than $\frac23$ on the remainder of the parameter space.


title: Analysis of error control in large scale two-stage multiple hypothesis testing

abstract: When dealing with the problem of simultaneously testing a large number of null hypotheses, a natural testing strategy is to first reduce the number of tested hypotheses by some selection (screening or filtering) process, and then to simultaneously test the selected hypotheses. The main advantage of this strategy is to greatly reduce the severe effect of high dimensions. However, the first screening or selection stage must be properly accounted for in order to maintain some type of error control. In this paper, we will introduce a selection rule based on a selection statistic that is independent of the test statistic when the tested hypothesis is true. Combining this selection rule and the conventional Bonferroni procedure, we can develop a powerful and valid two-stage procedure. The introduced procedure has several nice properties: (i) it completely removes the selection effect; (ii) it reduces the multiplicity effect; (iii) it does not "waste" data while carrying out both selection and testing. Asymptotic power analysis and simulation studies illustrate that this proposed method can provide higher power compared to usual multiple testing methods while controlling the Type 1 error rate. Optimal selection thresholds are also derived based on our asymptotic analysis.


title: Rediscovery Datasets: Connecting Duplicate Reports

abstract: The same defect can be rediscovered by multiple clients, causing unplanned outages and leading to reduced customer satisfaction. In the case of popular open source software, high volume of defects is reported on a regular basis. A large number of these reports are actually duplicates / rediscoveries of each other. Researchers have analyzed the factors related to the content of duplicate defect reports in the past. However, some of the other potentially important factors, such as the inter-relationships among duplicate defect reports, are not readily available in defect tracking systems such as Bugzilla. This information may speed up bug fixing, enable efficient triaging, improve customer profiles, etc. In this paper, we present three defect rediscovery datasets mined from Bugzilla. The datasets capture data for three groups of open source software projects: Apache, Eclipse, and KDE. The datasets contain information about approximately 914 thousands of defect reports over a period of 18 years (1999-2017) to capture the inter-relationships among duplicate defects. We believe that sharing these data with the community will help researchers and practitioners to better understand the nature of defect rediscovery and enhance the analysis of defect reports.


title: On Piercing Numbers of Families Satisfying the $(p,q)_r$ Property

abstract: The Hadwiger-Debrunner number $HD_d(p,q)$ is the minimal size of a piercing set that can always be guaranteed for a family of compact convex sets in $\mathbb{R}^d$ that satisfies the $(p,q)$ property. Hadwiger and Debrunner showed that $HD_d(p,q) \geq p-q+1$ for all $q$, and equality is attained for $q > \frac{d-1}{d}p +1$. Almost tight upper bounds for $HD_d(p,q)$ for a `sufficiently large' $q$ were obtained recently using an enhancement of the celebrated Alon-Kleitman theorem, but no sharp upper bounds for a general $q$ are known. In [L. Montejano and P. Sober\'{o}n, Piercing numbers for balanced and unbalanced families, Disc. Comput. Geom., 45(2) (2011), pp. 358-364], Montejano and Sober\'{o}n defined a refinement of the $(p,q)$ property: $F$ satisfies the $(p,q)_r$ property if among any $p$ elements of $F$, at least $r$ of the $q$-tuples intersect. They showed that $HD_d(p,q)_r \leq p-q+1$ holds for all $r>{{p}\choose{q}}-{{p+1-d}\choose{q+1-d}}$; however, this is far from being tight. In this paper we present improved asymptotic upper bounds on $HD_d(p,q)_r$ which hold when only a tiny portion of the $q$-tuples intersect. In particular, we show that for $p,q$ sufficiently large, $HD_d(p,q)_r \leq p-q+1$ holds with $r = \frac{1}{p^{\frac{q}{2d}}}{{p}\choose{q}}$. Our bound misses the known lower bound for the same piercing number by a factor of less than $pq^d$. Our results use Kalai's Upper Bound Theorem for convex sets, along with the Hadwiger-Debrunner theorem and the recent improved upper bound on $HD_d(p,q)$ mentioned above.


title: PatternNet: Visual Pattern Mining with Deep Neural Network

abstract: Visual patterns represent the discernible regularity in the visual world. They capture the essential nature of visual objects or scenes. Understanding and modeling visual patterns is a fundamental problem in visual recognition that has wide ranging applications. In this paper, we study the problem of visual pattern mining and propose a novel deep neural network architecture called PatternNet for discovering these patterns that are both discriminative and representative. The proposed PatternNet leverages the filters in the last convolution layer of a convolutional neural network to find locally consistent visual patches, and by combining these filters we can effectively discover unique visual patterns. In addition, PatternNet can discover visual patterns efficiently without performing expensive image patch sampling, and this advantage provides an order of magnitude speedup compared to most other approaches. We evaluate the proposed PatternNet subjectively by showing randomly selected visual patterns which are discovered by our method and quantitatively by performing image classification with the identified visual patterns and comparing our performance with the current state-of-the-art. We also directly evaluate the quality of the discovered visual patterns by leveraging the identified patterns as proposed objects in an image and compare with other relevant methods. Our proposed network and procedure, PatterNet, is able to outperform competing methods for the tasks described.


title: Weighted spherical means generated by generalized translation and general Euler-Poisson-Darboux equation

abstract: We consider the spherical mean generated by a multidimensional generalized translation and general Euler-Poisson-Darboux equation corresponding to this mean. The Asgeirsson property of solutions of the ultrahyperbolic equation that includes singular differential Bessel operators acting by each variable is provided.


title: Cosmological Evolution and Solar System Consistency of Massive Scalar-Tensor Gravity

abstract: The scalar-tensor theory of Damour and Esposito-Far\`ese recently gained some renewed interest because of its ability to suppress modifications to General Relativity in the weak field, while introducing large corrections in the strong field of compact objects through a process called scalarization. A large sector of this theory that allows for scalarization, however, has been shown to be in conflict with Solar System observations when accounting for the cosmological evolution of the scalar field. We here study an extension of this theory by endowing the scalar field with a mass to determine whether this allows the theory to pass Solar System constraints upon cosmological evolution for a larger sector of coupling parameter space. We show that the cosmological scalar field goes first through a quiescent phase, similar to the behavior of a massless field, but then it enters an oscillatory phase, with an amplitude (and frequency) that decays (and grows) exponentially. We further show that after the field enters the oscillatory phase, its effective energy density and pressure are approximately those of dust, as expected from previous cosmological studies. Due to these oscillations, we show that the scalar field cannot be treated as static today on astrophysical scales, and so we use time-dependent perturbation theory to compute the scalar-field-induced modifications to Solar System observables. We find that these modifications are suppressed when the mass of the scalar field and the coupling parameter of the theory are in a wide range, allowing the theory to pass Solar System constraints, while in principle possibly still allowing for scalarization.


title: Restricted sums of four integral squares

abstract: We give a simple quaternionic proof of a recent result of Goldmakher and Pollack on restricted sums of four integral squares.


title: New conformal mapping for adaptive resolving of the complex singularities of Stokes wave

abstract: A new highly efficient method is developed for computation of traveling periodic waves (Stokes waves) on the free surface of deep water. A convergence of numerical approximation is determined by the complex singularites above the free surface for the analytical continuation of the travelling wave into the complex plane. An auxiliary conformal mapping is introduced which moves singularities away from the free surface thus dramatically speeding up numerical convergence by adapting the numerical grid for resolving singularities while being consistent with the fluid dynamics. The efficiency of that conformal mapping is demonstrated for Stokes wave approaching the limiting Stokes wave (the wave of the greatest height) which significantly expands the family of numerically accessible solutions. It allows to provide a detailed study of the oscillatory approach of these solutions to the limiting wave. Generalizations of the conformal mapping to resolve multiple singularities are also introduced.

categories: physics.flu-dyn math.NA nlin.PS physics.ao-ph physics.comp-ph


title: Essential spectrum of elliptic systems of pseudo-differential operators on $L^2(\mathbb{R}^N)\oplus L^2(\mathbb{R}^N)$

abstract: Inspired by a result of Wong (Commun. Partial Differ. Equ. 13(10):1209-1221, 1988), we establish an analytic description of the essential spectrum of non-self-adjoint mixed-order systems of pseudo-differential operators on $L^2(\mathbb{R}^N) \oplus L^2(\mathbb{R}^N)$ that are uniformly Douglis-Nirenberg elliptic with positive-order diagonal entries. We apply our result to a problem arising in the dynamics of falling liquid films.


title: Transfer Learning for Sequence Tagging with Hierarchical Recurrent Networks

abstract: Recent papers have shown that neural networks obtain state-of-the-art performance on several different sequence tagging tasks. One appealing property of such systems is their generality, as excellent performance can be achieved with a unified architecture and without task-specific feature engineering. However, it is unclear if such systems can be used for tasks without large amounts of training data. In this paper we explore the problem of transfer learning for neural sequence taggers, where a source task with plentiful annotations (e.g., POS tagging on Penn Treebank) is used to improve performance on a target task with fewer available annotations (e.g., POS tagging for microblogs). We examine the effects of transfer learning for deep hierarchical recurrent networks across domains, applications, and languages, and show that significant improvement can often be obtained. These improvements lead to improvements over the current state-of-the-art on several well-studied tasks.


title: Roughness as a measurement of random disorder

abstract: We have studied the quenched random disorder (QRD) effects created by aerosil dispersion in octylcyanobiphenyl (8CB) liquid crystal (LC) using Atomic Force Microscopy (AFM) technique. Gelation process in the 8CB+aerosil gels yields a QRD network which also changes the surface topography. By increasing the aerosil concentration, the original smooth pattern of LC sample surfaces is suppressed by creating a fractal aerosil surface effect, these surfaces become more porous, rougher with more and bigger crevices. The dispersed aerosil also serves as pinning centers for the liquid crystal molecules. It has been observed that via the diffusion-limited-aggregation process, aerosil nano-particles yield a fractal like surface pattern for the less disordered samples. As the aerosil dispersion increases, the surface can be described by more aggregated regions which also introduce more roughness. Using this fact, we have shown that there is a net correlation between the short-ranged x-ray peak widths (-the results of previous x-ray diffraction experiments-) and the calculated surface roughness. In other words, we have shown that these QRD gels can also be characterized by their surface roughness values.


title: Triangle-free induced subgraphs of polarity graphs

abstract: Given a finite projective plane $\Pi$ and a polarity $\theta$ of $\Pi$, the corresponding polarity graph is the graph whose vertices are the points of $\Pi$. Two distinct vertices $p$ and $p'$ are adjacent if $p$ is incident to $\theta (p')$. Polarity graphs have been used in a variety of extremal problems, perhaps the most well-known being the Tur\'{a}n number of the cycle of length four. We investigate the problem of finding the maximum number of vertices in an induced triangle-free subgraph of a polarity graph. Mubayi and Williford showed that when $\Pi$ is the projective geometry $PG(2,q)$ and $\theta$ is the orthogonal polarity, an induced triangle-free subgraph has at most $\frac{1}{2}q^2 + O(q^{3/2})$ vertices. We generalize this result to all polarity graphs, and provide some interesting computational results that are relevant to an unresolved conjecture of Mubayi and Williford.


title: Application of Information Centric Networking to NoSQL Databases: the Spatio-Temporal use case

abstract: This paper explores methodologies, advantages and challenges related to the use of the Information Centric Network technology for developing NoSQL distributed databases, which are expected to play a central role in the forthcoming IoT and BigData era. ICN services make possible to simplify the development of the database software, improve performance, and provide data-level access control. We use our findings for devising a NoSQL spatio-temporal database, named OpenGeoBase, and evaluate its performance with a real data set related to Intelligent Transport System applications.


title: Stability for hyperplane complements of type B/C and statistics on squarefree polynomials over finite fields

abstract: In this paper we explore a relationship between the topology of the complex hyperplane complements $\mathcal{M}_{BC_n} (\mathbb{C})$ in type B/C and the combinatorics of certain spaces of degree-$n$ polynomials over a finite field $\mathbb{F}_q$. This relationship is a consequence of the Grothendieck trace formula and work of Lehrer and Kim. We use it to prove a correspondence between a representation-theoretic convergence result on the cohomology algebras $H^*(\mathcal{M}_{BC_n} (\mathbb{C});\mathbb{C})$, and an asymptotic stability result for certain polynomial statistics on monic squarefree polynomials over $\mathbb{F}_q$ with nonzero constant term. This result is the type B/C analogue of a theorem due to Church, Ellenberg, and Farb in type A, and we include a new proof of their theorem. To establish these convergence results, we realize the sequences of cohomology algebras of the hyperplane complements as FI$_\mathcal{W}$-algebras finitely generated in FI$_\mathcal{W}$- degree $2$, and we investigate the asymptotic behaviour of general families of algebras with this structure. We prove a negative result implying that this structure alone is not sufficient to prove the necessary convergence conditions. Our proof of convergence for the cohomology algebras involves the combinatorics of their relators.


title: Engineering Trustworthy Self-Adaptive Software with Dynamic Assurance Cases

abstract: Building on concepts drawn from control theory, self-adaptive software handles environmental and internal uncertainties by dynamically adjusting its architecture and parameters in response to events such as workload changes and component failures. Self-adaptive software is increasingly expected to meet strict functional and non-functional requirements in applications from areas as diverse as manufacturing, healthcare and finance. To address this need, we introduce a methodology for the systematic ENgineering of TRUstworthy Self-adaptive sofTware (ENTRUST). ENTRUST uses a combination of (1) design-time and runtime modelling and verification, and (2) industry-adopted assurance processes to develop trustworthy self-adaptive software and assurance cases arguing the suitability of the software for its intended application. To evaluate the effectiveness of our methodology, we present a tool-supported instance of ENTRUST and its use to develop proof-of-concept self-adaptive software for embedded and service-based systems from the oceanic monitoring and e-finance domains, respectively. The experimental results show that ENTRUST can be used to engineer self-adaptive software systems in different application domains and to generate dynamic assurance cases for these systems.


title: Election Predictions as Martingales: An Arbitrage Approach

abstract: We consider the estimation of binary election outcomes as martingales and propose an arbitrage pricing when one continuously updates estimates. We argue that the estimator needs to be priced as a binary option as the arbitrage valuation minimizes the conventionally used Brier score for tracking the accuracy of probability assessors. We create a dual martingale process $Y$, in $[L,H]$ from the standard arithmetic Brownian motion, $X$ in $(-\infty, \infty)$ and price elections accordingly. The dual process $Y$ can represent the numerical votes needed for success. We show the relationship between the volatility of the estimator in relation to that of the underlying variable. When there is a high uncertainty about the final outcome, 1) the arbitrage value of the binary gets closer to 50\%, 2) the estimate should not undergo large changes even if polls or other bases show significant variations. There are arbitrage relationships between 1) the binary value, 2) the estimation of $Y$, 3) the volatility of the estimation of $Y$ over the remaining time to expiration. We note that these arbitrage relationships were often violated by the various forecasting groups in the U.S. presidential elections of 2016, as well as the notion that all intermediate assessments of the success of a candidate need to be considered, not just the final one.


title: Search for a heavy resonance decaying to a top quark and a vector-like top quark at sqrt(s) = 13 TeV

abstract: A search is presented for massive spin-1 Z' resonances decaying to a top quark and a heavy vector-like top quark partner T. The search is based on a 2.6 inverse femtobarns sample of proton-proton collisions at 13 TeV collected with the CMS detector at the LHC. The analysis is optimized for final states in which the T quark decays to a W boson and a bottom quark. The focus is on all-jet final states in which both the W boson and the top quark decay into quarks that evolve into jets. The decay products of the top quark and of the W boson are assumed to be highly Lorentz-boosted and cannot be reconstructed as separate jets, but are instead reconstructed as merged, wide jets. Techniques for the identification of jet substructure and jet flavour are used to distinguish signal from background events. Several models for Z' bosons decaying to T quarks are excluded at 95% confidence level, with upper limits on the cross section ranging from 0.13 to 10 pb, depending on the chosen hypotheses. This is the first search for a neutral spin-1 heavy resonance decaying to a top quark and a vector-like T quark in the all-hadronic final state.


title: Systematic Mapping Study of Template-based Code Generation

abstract: Template-based code generation (TBCG) is a synthesis technique that produces code from high-level specifications, called templates. TBCG is a popular technique in model-driven engineering (MDE) given that they both emphasize abstraction and automation. Given the diversity of tools and approaches, it is necessary to classify existing TBCG techniques to better guide developers in their choices. The goal of this article is to better understand the characteristics of TBCG techniques and associated tools, identify research trends, and assess the importance of the role of MDE in this code synthesis approach. We conducted a systematic mapping study of the literature to paint an interesting picture about the trends and uses of TBCG. Our study shows that the community has been diversely using TBCG over the past 15 years. TBCG has greatly benefited from MDE. It has favored a template style that is output-based and high level modeling languages as input. TBCG is mainly used to generate source code and has been applied in a variety of domains. Furthermore, both MDE and non-MDE tools are becoming effective development resources in industry.


title: Goal Conflict in Designing an Autonomous Artificial System

abstract: Research on human self-regulation has shown that people hold many goals simultaneously and have complex self-regulation mechanisms to deal with this goal conflict. Artificial autonomous systems may also need to find ways to cope with conflicting goals. Indeed, the intricate interplay among different goals may be critical to the design as well as long-term safety and stability of artificial autonomous systems. I discuss some of the critical features of the human self-regulation system and how it might be applied to an artificial system. Furthermore, the implications of goal conflict for the reliability and stability of artificial autonomous systems and ensuring their alignment with human goals and ethics is examined.


title: Weyl relativity: A novel approach to Weyl's ideas

abstract: In this paper we revisit the motivation and construction of a unified theory of gravity and electromagnetism, following Weyl's insights regarding the appealing potential connection between the gauge invariance of electromagnetism and the conformal invariance of the gravitational field. We highlight that changing the local symmetry group of spacetime permits to construct a theory in which these two symmetries are combined into a putative gauge symmetry but with second-order field equations and non-trivial mass scales, unlike the original higher-order construction by Weyl. We prove that the gravitational field equations are equivalent to the (trace-free) Einstein field equations, ensuring their compatibility with known tests of general relativity. As a corollary, the effective cosmological constant is rendered radiatively stable due to Weyl invariance. A novel phenomenological consequence characteristic of this construction, potentially relevant for cosmological observations, is the existence of an energy scale below which effects associated with the non-integrability of spacetime distances, and an effective mass for the electromagnetic field, appear simultaneously (as dual manifestations of the use of Weyl connections). We explain how former criticisms against Weyl's ideas lose most of their power in its present reincarnation, which we refer to as Weyl relativity, as it represents a Weyl-invariant, unified description of both the Einstein and Maxwell field equations.


title: Subset Synchronization in Monotonic Automata

abstract: We study extremal and algorithmic questions of subset and careful synchronization in monotonic automata. We show that several synchronization problems that are hard in general automata can be solved in polynomial time in monotonic automata, even without knowing a linear order of the states preserved by the transitions. We provide asymptotically tight bounds on the maximum length of a shortest word synchronizing a subset of states in a monotonic automaton and a shortest word carefully synchronizing a partial monotonic automaton. We provide a complexity framework for dealing with problems for monotonic weakly acyclic automata over a three-letter alphabet, and use it to prove NP-completeness and inapproximability of problems such as {\sc Finite Automata Intersection} and the problem of computing the rank of a subset of states in this class. We also show that checking whether a monotonic partial automaton over a four-letter alphabet is carefully synchronizing is NP-hard. Finally, we give a simple necessary and sufficient condition when a strongly connected digraph with a selected subset of vertices can be transformed into a deterministic automaton where the corresponding subset of states is synchronizing.


title: Thin obstacle problem: estimates of the distance to the exact solution

abstract: We consider elliptic variational inequalities generated by obstacle type problems with thin obstacles. For this class of problems, we deduce estimates of the distance (measured in terms of the natural energy norm) between the exact solution and any function that satisfies the boundary condition and is admissible with respect to the obstacle condition (i.e., it is valid for any approximation regardless of the method by which it was found). Computation of the estimates does not require knowledge of the exact solution and uses only the problem data and an approximation. The estimates provide guaranteed upper bounds of the error (error majorants) and vanish if and only if the approximation coincides with the exact solution. In the last section, the efficiency of error majorants is confirmed by an example, where the exact solution is known.


title: The Extreme Energy Events HECR array: status and perspectives

abstract: The Extreme Energy Events Project is a synchronous sparse array of 52 tracking detectors for studying High Energy Cosmic Rays (HECR) and Cosmic Rays-related phenomena. The observatory is also meant to address Long Distance Correlation (LDC) phenomena: the network is deployed over a broad area covering 10 degrees in latitude and 11 in longitude. An overview of a set of preliminary results is given, extending from the study of local muon flux dependance on solar activity to the investigation of the upward-going component of muon flux traversing the EEE stations; from the search for anisotropies at the sub-TeV scale to the hints for observations of km-scale Extensive Air Shower (EAS).


title: Fully symmetric kernel quadrature

abstract: Kernel quadratures and other kernel-based approximation methods typically suffer from prohibitive cubic time and quadratic space complexity in the number of function evaluations. The problem arises because a system of linear equations needs to be solved. In this article we show that the weights of a kernel quadrature rule can be computed efficiently and exactly for up to tens of millions of nodes if the kernel, integration domain, and measure are fully symmetric and the node set is a union of fully symmetric sets. This is based on the observations that in such a setting there are only as many distinct weights as there are fully symmetric sets and that these weights can be solved from a linear system of equations constructed out of row sums of certain submatrices of the full kernel matrix. We present several numerical examples that show feasibility, both for a large number of nodes and in high dimensions, of the developed fully symmetric kernel quadrature rules. Most prominent of the fully symmetric kernel quadrature rules we propose are those that use sparse grids.


title: Casting the Coronal Magnetic Field Reconstruction Tools in 3D Using MHD Bifrost Model

abstract: Quantifying coronal magnetic field remains a central problem in solar physics. Nowadays the coronal magnetic field is often modelled using nonlinear force-free field (NLFFF) reconstructions, whose accuracy has not yet been comprehensively assessed. Here we perform a detailed casting of the NLFFF reconstruction tools, such as pi-disambiguation, photospheric field preprocessing, and volume reconstruction methods using a 3D snapshot of the publicly available full-fledged radiative MHD model. Specifically, from the MHD model we know the magnetic field vector in the entire 3D domain, which enables us to perform "voxel-by-voxel" comparison of the restored and the true magnetic field in the 3D model volume. Our tests show that the available pi-disambiguation methods often fail at the quiet sun areas dominated by small-scale magnetic elements, while they work well at the AR photosphere and (even better) chromosphere. The preprocessing of the photospheric magnetic field, although does produce a more force-free boundary condition, also results in some effective `elevation' of the magnetic field components. This `elevation' height is different for the longitudinal and transverse components, which results in a systematic error in absolute heights in the reconstructed magnetic data cube. The extrapolations performed starting from actual AR photospheric magnetogram are free from this systematic error, while have other metrics comparable with those for extrapolations from the preprocessed magnetograms. This finding favors the use of extrapolations from the original photospheric magnetogram without preprocessing. Our tests further suggest that extrapolations from a force-free chromospheric boundary produce measurably better results, than those from the photospheric boundary.


title: Which friends are more popular than you? Contact strength and the friendship paradox in social networks

abstract: The friendship paradox states that in a social network, egos tend to have lower degree than their alters, or, "your friends have more friends than you do". Most research has focused on the friendship paradox and its implications for information transmission, but treating the network as static and unweighted. Yet, people can dedicate only a finite fraction of their attention budget to each social interaction: a high-degree individual may have less time to dedicate to individual social links, forcing them to modulate the quantities of contact made to their different social ties. Here we study the friendship paradox in the context of differing contact volumes between egos and alters, finding a connection between contact volume and the strength of the friendship paradox. The most frequently contacted alters exhibit a less pronounced friendship paradox compared with the ego, whereas less-frequently contacted alters are more likely to be high degree and give rise to the paradox. We argue therefore for a more nuanced version of the friendship paradox: "your closest friends have slightly more friends than you do", and in certain networks even: "your best friend has no more friends than you do". We demonstrate that this relationship is robust, holding in both a social media and a mobile phone dataset. These results have implications for information transfer and influence in social networks, which we explore using a simple dynamical model.


title: Resonance tongues for the Hill equation with Duffing coefficients and instabilities in a nonlinear beam equation

abstract: We consider a class of Hill equations where the periodic coefficient is the squared solution of some Duffing equation plus a constant. We study the stability of the trivial solution of this Hill equation and we show that a criterion due to Burdina (V.I. Burdina, Boundedness of solutions of a system of differential equations) is very helpful for this analysis. In some cases, we are also able to determine exact solutions in terms of Jacobi elliptic functions. Overall, we obtain a fairly complete picture of the stability and instability regions. These results are then used to study the stability of nonlinear modes in some beam equations.