title: Optimal control of membrane filtration systems

abstract: This paper studies an optimal control problem related to membrane filtration processes. A simple mathematical model of membrane fouling is used to capture the dynamic behavior of the filtration process which consists in the attachment of matter onto the membrane during the filtration period and the detachment of matter during the cleaning period. We consider the maximization of the net production of a membrane filtration system (i.e. the filtrate) over a finite time horizon, where control variable is the sequence of filtration/backwashing cycles over the operation time of process. Based on the Pontryagin Maximum Principle, we characterize the optimal control strategy and show that it presents a singular arc. Moreover we prove the existence of an additional switching curve before reaching the terminal state, and also the possibility of having a dispersal curve as a locus where two different strategies are both optimal.


title: Optically bistable driven-dissipative Bose-Hubbard dimer: Gutzwiller approaches and entanglement

abstract: We theoretically examine the driven-dissipative Bose-Hubbard dimer in the optical bistable regime. Various approximation schemes based on a Gutzwiller mean field decoupling are applied and compared. Depending on the system parameters we show that a decoupling with respect to the real space or to the reciprocal space can be more accurate. The Gutzwiller decoupling is applied both at the level of the density matrix and for the wavefunction during a quantum trajectory simulation. The latter is shown to be a more accurate approximation. A Gaussian approximation for the non-homogeneous anti-bonding mode is also explored. We also show that entanglement in this system is witnessed by squeezing in reciprocal space.


title: Scene-adapted plug-and-play algorithm with convergence guarantees

abstract: Recent frameworks, such as the so-called plug-and-play, allow us to leverage the developments in image denoising to tackle other, and more involved, problems in image processing. As the name suggests, state-of-the-art denoisers are plugged into an iterative algorithm that alternates between a denoising step and the inversion of the observation operator. While these tools offer flexibility, the convergence of the resulting algorithm may be difficult to analyse. In this paper, we plug a state-of-the-art denoiser, based on a Gaussian mixture model, in the iterations of an alternating direction method of multipliers and prove the algorithm is guaranteed to converge. Moreover, we build upon the concept of scene-adapted priors where we learn a model targeted to a specific scene being imaged, and apply the proposed method to address the hyperspectral sharpening problem.


title: Tiered trees, weights, and q-Eulerian numbers

abstract: Maxmin trees are labeled trees with the property that each vertex is either a local maximum or a local minimum. Such trees were originally introduced by Postnikov, who gave a formula to count them and different combinatorial interpretations for their number. In this paper we generalize this construction and define tiered trees by allowing more than two classes of vertices. Tiered trees arise naturally when counting the absolutely indecomposable representations of certain quivers, and also when one enumerates torus orbits on certain homogeneous varieties. We define a notion of weight for tiered trees and prove bijections between various weight 0 tiered trees and other combinatorial objects; in particular order n weight 0 maxmin trees are naturally in bijection with permutations on n-1 letters. We conclude by using our weight function to define a new q-analogue of the Eulerian numbers.


title: Region Ensemble Network: Improving Convolutional Network for Hand Pose Estimation

abstract: Hand pose estimation from monocular depth images is an important and challenging problem for human-computer interaction. Recently deep convolutional networks (ConvNet) with sophisticated design have been employed to address it, but the improvement over traditional methods is not so apparent. To promote the performance of directly 3D coordinate regression, we propose a tree-structured Region Ensemble Network (REN), which partitions the convolution outputs into regions and integrates the results from multiple regressors on each regions. Compared with multi-model ensemble, our model is completely end-to-end training. The experimental results demonstrate that our approach achieves the best performance among state-of-the-arts on two public datasets.


title: Nuclear matrix element of neutrinoless double-$\beta$ decay: Relativity and short-range correlations

abstract: Background: The discovery of neutrinoless double-beta ($0\nu\beta\beta$) decay would demonstrate the nature of neutrinos, have profound implications for our understanding of matter-antimatter mystery, and solve the mass hierarchy problem of neutrinos. The calculations for the nuclear matrix elements $M^{0\nu}$ of $0\nu\beta\beta$ decay are crucial for the interpretation of this process. Purpose: We study the effects of relativity and nucleon-nucleon short-range correlations on the nuclear matrix elements $M^{0\nu}$ by assuming the mechanism of exchanging light or heavy neutrinos for the $0\nu\beta\beta$ decay. Methods: The nuclear matrix elements $M^{0\nu}$ are calculated within the framework of covariant density functional theory, where the beyond-mean-field correlations are included in the nuclear wave functions by configuration mixing of both angular-momentum and particle-number projected quadrupole deformed mean-field states. Results: The nuclear matrix elements $M^{0\nu}$ are obtained for ten $0\nu\beta\beta$-decay candidate nuclei. The impact of relativity is illustrated by adopting relativistic or nonrelativistic decay operators. The effects of short-range correlations are evaluated. Conclusions: The effects of relativity and short-range correlations play an important role in the mechanism of exchanging heavy neutrinos though the influences are marginal for light neutrinos. Combining the nuclear matrix elements $M^{0\nu}$ with the observed lower limits on the $0\nu\beta\beta$-decay half-lives, the predicted strongest limits on the effective masses are $|\langle m_\nu\rangle|<0.06~\mathrm{eV}$ for light neutrinos and $|\langle m_{\nu_h}^{-1}\rangle|^{-1}>3.065\times 10^8~\mathrm{GeV}$ for heavy neutrinos.


title: Nonparametric mean curvature type flows of graphs with contact angle conditions

abstract: In this paper we study nonparametric mean curvature type flows in $M\times\mathbb{R}$ which are represented as graphs $(x,u(x,t))$ over a domain in a Riemannian manifold $M$ with prescribed contact angle. The speed of $u$ is the mean curvature speed minus an admissible function $\psi(x,u,Du)$. Long time existence and uniformly convergence are established if $\psi(x,u, Du)\equiv 0$ with vertical contact angle and $\psi(x,u,Du)=h(x,u)\omega$ with $h_u(x,u)\geq h_0>0$ and $\omega=\sqrt{1+|Du|^2}$. Their applications include mean curvature type equations with prescribed contact angle boundary condition and the asymptotic behavior of nonparametric mean curvature flows of graphs over a convex domain in $M^2$ which is a surface with nonnegative Ricci curvature.


title: Ironwood Meta Key Agreement and Authentication Protocol

abstract: Number theoretic public-key solutions, currently used in many applications worldwide, will be subject to various quantum attacks, making them less attractive for longer-term use. Certain group theoretic constructs are now showing promise in providing quantum-resistant cryptographic primitives, and may provide suitable alternatives for those looking to address known quantum attacks. In this paper, we introduce a new protocol called a Meta Key Agreement and Authentication Protocol (MKAAP) that has some characteristics of a public-key solution and some of a shared-key solution. Specifically it has the deployment benefits of a public-key system, allowing two entities that have never met before to authenticate without requiring real-time access to a third-party, but does require secure provisioning of key material from a trusted key distribution system (similar to a symmetric system) prior to deployment. We then describe a specific MKAAP instance, the Ironwood MKAAP, discuss its security, and show how it resists certain quantum attacks such as Shor's algorithm or Grover's quantum search algorithm. We also show Ironwood implemented on several ``internet of things'' (IoT devices), measure its performance, and show how it performs significantly better than ECC using fewer device resources.


title: Magnetic vortices as localized mesoscopic domain wall pinning sites

abstract: We report on the controllable pinning of domain walls in stripes with perpendicular magnetic anisotropy by magnetostatic coupling to magnetic vortices in disks located above the stripe. Pinning mechanisms and depinning fields are reported. This novel pinning strategy, which can be realized by current nanofabrication techniques, opens up new possibilities for the non-destructive control of domain wall mobility in domain wall based spintronic devices.


title: kappa-Poincare as a symmetry of flat quantum spacetime

abstract: In this short note, based on the talk given at the 3rd Conference of the Polish Society on Relativity, I present the basic points of our recent paper "Symmetries of quantum spacetime in three dimensions", stressing their physical meaning, and avoiding technical details.


title: Preparing for the Unknown: Learning a Universal Policy with Online System Identification

abstract: We present a new method of learning control policies that successfully operate under unknown dynamic models. We create such policies by leveraging a large number of training examples that are generated using a physical simulator. Our system is made of two components: a Universal Policy (UP) and a function for Online System Identification (OSI). We describe our control policy as universal because it is trained over a wide array of dynamic models. These variations in the dynamic model may include differences in mass and inertia of the robots' components, variable friction coefficients, or unknown mass of an object to be manipulated. By training the Universal Policy with this variation, the control policy is prepared for a wider array of possible conditions when executed in an unknown environment. The second part of our system uses the recent state and action history of the system to predict the dynamics model parameters mu. The value of mu from the Online System Identification is then provided as input to the control policy (along with the system state). Together, UP-OSI is a robust control policy that can be used across a wide range of dynamic models, and that is also responsive to sudden changes in the environment. We have evaluated the performance of this system on a variety of tasks, including the problem of cart-pole swing-up, the double inverted pendulum, locomotion of a hopper, and block-throwing of a manipulator. UP-OSI is effective at these tasks across a wide range of dynamic models. Moreover, when tested with dynamic models outside of the training range, UP-OSI outperforms the Universal Policy alone, even when UP is given the actual value of the model dynamics. In addition to the benefits of creating more robust controllers, UP-OSI also holds out promise of narrowing the Reality Gap between simulated and real physical systems.


title: Coloring Cantor sets and resolvability of pseudocompact spaces

abstract: Let us denote by $\Phi(\lambda,\mu)$ the statement that $\mathbb{B}(\lambda) = D(\lambda)^\omega$, i.e. the Baire space of weight $\lambda$, has a coloring with $\mu$ colors such that every homeomorphic copy of the Cantor set $\mathbb{C}$ in $\mathbb{B}(\lambda)$ picks up all the $\mu$ colors. We call a space $X\,$ {\em $\pi$-regular} if it is Hausdorff and for every non-empty open set $U$ in $X$ there is a non-empty open set $V$ such that $\overline{V} \subset U$. We recall that a space $X$ is called {\em feebly compact} if every locally finite collection of open sets in $X$ is finite. A Tychonov space is pseudocompact iff it is feebly compact. The main result of this paper is the following. Theorem. Let $X$ be a crowded feebly compact $\pi$-regular space and $\mu$ be a fixed (finite or infinite) cardinal. If $\Phi(\lambda,\mu)$ holds for all $\lambda < \widehat{c}(X)$ then $X$ is $\mu$-resolvable, i.e. contains $\mu$ pairwise disjoint dense subsets. (Here $\widehat{c}(X)$ is the smallest cardinal $\kappa$ such that $X$ does not contain $\kappa$ many pairwise disjoint open sets.) This significantly improves earlier results of van Mill , resp. Ortiz-Castillo and Tomita.


title: Deterministic Protocols in the SINR Model without Knowledge of Coordinates

abstract: Much work has been developed for studying the classical broadcasting problem in the SINR (Signal-to-Interference-plus-Noise-Ratio) model for wireless device transmission. The setting typically studied is when all radio nodes transmit a signal of the same strength. This work studies the challenging problem of devising a distributed algorithm for multi-broadcasting, assuming a subset of nodes are initially awake, for the SINR model when each device only has access to knowledge about the total number of nodes in the network $n$, the range from which each node's label is taken $\lbrace 1,\dots,N \rbrace$, and the label of the device itself. Specifically, we assume no knowledge of the physical coordinates of devices and also no knowledge of the neighborhood of each node. We present a deterministic protocol for this problem in $O(n \lg N \lg n)$ rounds. There is no known polynomial time deterministic algorithm in literature for this setting, and it remains the principle open problem in this domain. A lower bound of $\Omega(n \lg N)$ rounds is known for deterministic broadcasting without local knowledge. In addition to the above result, we present algorithms to achieve multi-broadcast in $O(n \lg N)$ rounds and create a backbone in $O(n \lg N)$ rounds, assuming that all nodes are initially awake. For a given backbone, messages can be exchanged between every pair of connected nodes in the backbone in $O(\lg N)$ rounds and between any node and its designated contact node in the backbone in $O(\Delta \lg N)$ rounds.


title: Understanding the Spatial and Temporal Activity Patterns of Subway Mobility Flows

abstract: In urban transportation systems, mobility flows in the subway system reflect the spatial and temporal dynamics of working days. To investigate the variability of mobility flows, we analyse the spatial community through a series of snapshots of subway stations over sequential periods. Using Shanghai as a case study, we find that the spatial community snapshots reveal dynamic passenger activities. Adopting a dual-perspective, we apply spatial and temporal models separately to explore where and when individuals travel for entertainment. In the two models, microblog topics and spatial facilities such as food venues and entertainment businesses are used to characterise the spatial popularity of each station and people's travelling perceptions. In the studied case, the city centre is characterised by greater social influence, and it is better described by the spatial model. In the temporal model, shorter travel distances motivate individuals to start their trips earlier. Interestingly, as the number of food-related facilities near the starting station increases, until it exceeds 1563, the speed of people's journeys slows down. This study provides a method for modelling the effects of social features on mobility flows and for predicting the spatial-temporal mobility flows of newly built subway stations.


title: PHz-wide spectral interference through coherent plasma-induced fission of higher-order solitons

abstract: We identify a novel regime of soliton-plasma interactions in which high-intensity ultrashort pulses of intermediate soliton order undergo coherent plasma-induced fission. Experimental results obtained in gas-filled hollow-core photonic crystal fibers are supported by rigorous numerical simulations. The cumulative blueshift of higher-order input solitons with ionizing intensities results in pulse splitting before the ultimate self-compression point, leading to the generation of robust pulse pairs with PHz bandwidths. The novel dynamics closes the gap between plasma-induced adiabatic soliton compression and modulational instability.


title: A magnetic phase-transition graphene transistor with tunable spin polarization

abstract: Graphene nanoribbons (GNRs) have been proposed as potential building blocks for field effect transistor (FET) devices due to their quantum confinement bandgap. Here, we propose a novel GNR device concept, enabling the control of both charge and spin signals, integrated within the simplest three-terminal device configuration. In a conventional FET device, a gate electrode is employed to tune the Fermi level of the system in and out of a static bandgap. By contrast, in the switching mechanism proposed here, the applied gate voltage can dynamically open and close an interaction gap, with only a minor shift of the Fermi level. Furthermore, the strong interplay of the band structure and edge spin configuration in zigzag ribbons enables such transistors to carry spin polarized current without employing an external magnetic field or ferromagnetic contacts. Using an experimentally validated theoretical model, we show that such transistors can switch at low voltages and high speed, and the spin polarization of the current can be tuned from 0% to 50% by using the same back gate electrode. Furthermore, such devices are expected to be robust against edge irregularities and can operate at room temperature. Controlling both charge and spin signal within the simplest FET device configuration could open up new routes in data processing with graphene based devices.


title: On the Possibility of Singularities on the Ambient Boundary

abstract: The order horismos induces the Zeeman $Z$ topology, which is coarser than the Fine Zeeman Topology $F$. The causal curves in a spacetime under $Z$ are piecewise null. $F$ is considered to be the most physical topology in a spacetime manifold $M$, as the group of homeomorphisms of $M$ is isomorphic to the group of homothetic transformations of $M$. $Z$ was used in the Ambient Boundary-Ambient Space cosmological model, in order to show that there is no possibility of formation of spacetime singularities. In this article we question this result, by reviewing the corresponding articles, and we propose new questions towards the improvement of this model.


title: Deterministic Backbone Creation in an SINR Network without Knowledge of Location

abstract: For a given network, a backbone is an overlay network consisting of a connected dominating set with additional accessibility properties. Once a backbone is created for a network, it can be utilized for fast communication amongst the nodes of the network. The Signal-to-Interference-plus-Noise-Ratio (SINR) model has become the standard for modeling communication among devices in wireless networks. For this model, the community has pondered what the most realistic solutions for communication problems in wireless networks would look like. Such solutions would have the characteristic that they would make the least number of assumptions about the availability of information about the participating nodes. Solving problems when nothing at all is known about the network and having nodes just start participating would be ideal. However, this is quite challenging and most likely not feasible. The pragmatic approach is then to make meaningful assumptions about the available information and present efficient solutions based on this information. We present a solution for creation of backbone in the SINR model, when nodes do not have access to their physical coordinates or the coordinates of other nodes in the network. This restriction models the deployment of nodes in various situations for sensing hurricanes, cyclones, and so on, where only information about nodes prior to their deployment may be known but not their actual locations post deployment. We assume that nodes have access to knowledge of their label, the labels of nodes within their neighborhood, the range from which labels are taken $[N]$ and the total number of participating nodes $n$. We also assume that nodes wake up spontaneously. We present an efficient deterministic protocol to create a backbone with a round complexity of $O(\Delta \lg^2 N)$.


title: Who Will Win Practical Artificial Intelligence? AI Engineerings in China

abstract: Currently, Artificial Intelligence (AI) has won unprecedented attention and is becoming the increasingly popular focus in China. This change can be judged by the impressive record of academic publications, the amount of state-level investment and the presence of nation-wide participation and devotion. In this paper, we place emphasis on discussing the progress of artificial intelligence engineerings in China. We first introduce the focus on AI in Chinese academia, including the supercomputing brain system, Cambrian Period supercomputer of neural networks, and biometric systems. Then, the development of AI in industrial circles and the latest layout of AI products in companies, such as Baidu, Tencent, and Alibaba, are introduced. Last, we bring in the opinions and arguments of the main intelligentsia of China about the future development of AI, including how to examine the relationship between humanity on one side and science and technology on the other.


title: Integrated Information as a Metric for Group Interaction: Analyzing Human and Computer Groups Using a Technique Developed to Measure Consciousness

abstract: Researchers in many disciplines have previously used a variety of mathematical techniques for analyzing group interactions. Here we use a new metric for this purpose, called 'integrated information' or 'phi.' Phi was originally developed by neuroscientists as a measure of consciousness in brains, but it captures, in a single mathematical quantity, two properties that are important in many other kinds of groups as well: differentiated information and integration. Here we apply this metric to the activity of three types of groups that involve people and computers. First, we find that 4-person work groups with higher measured phi perform a wide range of tasks more effectively, as measured by their collective intelligence. Next, we find that groups of Wikipedia editors with higher measured phi create higher quality articles. Last, we find that the measured phi of the collection of people and computers communicating on the Internet increased over a recent six-year period. Together, these results suggest that integrated information can be a useful way of characterizing a certain kind of interactional complexity that, at least sometimes, predicts group performance. In this sense, phi can be viewed as a potential metric of effective group collaboration. Since the metric was originally developed as a measure of consciousness, the results also raise intriguing questions about the conditions under which it might be useful to regard groups as having a kind of consciousness.


title: Video Frame Synthesis using Deep Voxel Flow

abstract: We address the problem of synthesizing new video frames in an existing video, either in-between existing frames (interpolation), or subsequent to them (extrapolation). This problem is challenging because video appearance and motion can be highly complex. Traditional optical-flow-based solutions often fail where flow estimation is challenging, while newer neural-network-based methods that hallucinate pixel values directly often produce blurry results. We combine the advantages of these two methods by training a deep network that learns to synthesize video frames by flowing pixel values from existing ones, which we call deep voxel flow. Our method requires no human supervision, and any video can be used as training data by dropping, and then learning to predict, existing frames. The technique is efficient, and can be applied at any video resolution. We demonstrate that our method produces results that both quantitatively and qualitatively improve upon the state-of-the-art.


title: SDSS IV MaNGA: Discovery of an $H_{\alpha}$ blob associated with a dry galaxy pair -- ejected gas or a `dark' galaxy candidate?

abstract: We report the discovery of a mysterious giant $H_{\alpha}$ blob that is $\sim 8$ kpc away from the main MaNGA target 1-24145, one component of a dry galaxy merger, identified in the first-year SDSS-IV MaNGA data. The size of the $H_{\alpha}$ blob is $\sim$ 3-4 kpc in radius, and the $H_{\alpha}$ distribution is centrally concentrated. However, there is no optical continuum counterpart in deep broadband images reaching $\sim$26.9 mag arcsec$^{-2}$ in surface brightness. We estimate that the masses of ionized and cold gases are $3.3 \times 10^{5}$ $\rm M_{\odot}$ and $< 1.3 \times 10^{9}$ $\rm M_{\odot}$, respectively. The emission-line ratios indicate that the $H_{\alpha}$ blob is photoionized by a combination of massive young stars and AGN. Furthermore, the ionization line ratio decreases from MaNGA 1-24145 to the $H_{\alpha}$ blob, suggesting that the primary ionizing source may come from MaNGA 1-24145, likely a low-activity AGN. Possible explanations of this $H_{\alpha}$ blob include AGN outflow, the gas remnant being tidally or ram-pressure stripped from MaNGA 1-24145, or an extremely low surface brightness (LSB) galaxy. However, the stripping scenario is less favoured according to galaxy merger simulations and the morphology of the $H_{\alpha}$ blob. With the current data, we can not distinguish whether this $H_{\alpha}$ blob is ejected gas due to a past AGN outburst, or a special category of `ultra-diffuse galaxy' (UDG) interacting with MaNGA 1-24145 that further induces the gas inflow to fuel the AGN in MaNGA 1-24145.


title: Radiogenic Neutron Yield Calculations for Low-Background Experiments

abstract: Nuclear recoil backgrounds are one of the most dangerous backgrounds for many dark matter experiments. A primary source of nuclear recoils is radiogenic neutrons produced in the detector material itself. These neutrons result from fission and $(\alpha,n)$ reactions originating from uranium and thorium contamination. In this paper, we discuss neutron yields from these sources. We compile a list of $(\alpha,n)$ yields for many materials common in low-background detectors, calculated using NeuCBOT, a new tool introduced in this paper, available at https://github.com/shawest/neucbot. These calculations are compared to computations made using data compilations and SOURCES-4A


title: The Majorana Demonstrator calibration system

abstract: The MAJORANA Collaboration is searching for the neutrinoless double-beta decay of the nucleus $^{76}$Ge. The MAJORANA DEMONSTRATOR is an array of germanium detectors deployed with the aim of implementing background reduction techniques suitable for a 1-tonne $^{76}$Ge-based search. The ultra low-background conditions require regular calibrations to verify proper function of the detectors. Radioactive line sources can be deployed around the cryostats containing the detectors for regular energy calibrations. When measuring in low-background mode, these line sources have to be stored outside the shielding so they do not contribute to the background. The deployment and the retraction of the source are designed to be controlled by the data acquisition system and do not require any direct human interaction. In this paper, we detail the design requirements and implementation of the calibration apparatus, which provides the event rates needed to define the pulse-shape cuts and energy calibration used in the final analysis as well as data that can be compared to simulations.


title: Charged gravastar model in $f(T)$ gravity admitting conformal motion

abstract: In this paper, model of charged gravastar under $f(T)$ modified gravity is obtained. The model has been explored by taking the diagonal tetrad field of static spacetime together with electric charge. To solve the Einstein-Maxwell field equations, along with $f(T)$ gravity, we assume the existence of a conformal Killing vector which relates between geometry and matter through the Einstein-Maxwell field equations by an inheritance symmetry. We study several cases of interest to explore physically valid features of the solutions. Some physical properties of the model are discussed and we match our interior spacetime to the exterior Reissner-Nordstr\"{o}m spacetime in presence of thin shell.


title: The phylogeny of quasars and the ontogeny of their central black holes

abstract: The connection between multifrequency quasar observational and physical parameters related to accretion processes is still open to debate. In the last 20 year, Eigenvector 1-based approaches developed since the early papers by Boroson and Green (1992) and Sulentic et al. (2000b) have been proven to be a remarkably powerful tool to investigate this issue, and have led to the definition of a quasar "main sequence". In this paper we perform a cladistic analysis on two samples of 215 and 85 low-z quasars (z 0.7) which were studied in several previous works and which offer a satisfactory coverage of the Eigenvector 1-derived main sequence. The data encompass accurate measurements of observational parameters which represent key aspects associated with the structural diversity of quasars. Cladistics is able to group sources radiating at higher Eddington ratios, as well as to separate radio-quiet (RQ) and radio-loud (RL) quasars. The analysis suggests a black hole mass threshold for powerful radio emission and also properly distinguishes core-dominated and lobe-dominated quasars, in accordance with the basic tenet of RL unification schemes. Considering that black hole mass provides a sort of "arrow of time" of nuclear activity, a phylogenetic interpretation becomes possible if cladistic trees are rooted on black hole mass: the ontogeny of black holes is represented by their monotonic increase in mass. More massive radio-quiet Population B sources at low-z become a more evolved counterpart of Population A i.e., wind dominated sources to which the "local" Narrow-Line Seyfert 1s belong.


title: Anisotropic cosmological models in $f(R,T)$ gravity with variable deceleration parameter

abstract: The objective of this work enclosed with the study of spatially homogeneous anisotropic Bianchi type-I universe in $f(R,T)$ gravity (where $R$ is the Ricci scalar and $T$ is the trace of stress energy momentum tensor) in two different cases viz. $f(R,T)=R+2f(T)$ and $f(R,T)=f_1(R)+f_2(T)$ with bulk viscosity matter content. In this study, we consider a time varying deceleration parameter, which generates an accelerating universe to obtain the exact solution of the field equations. The physical and kinematical properties of both the models are discussed in detail for the future evolution of the universe. We have explored the nature of WEC, DEC, SEC and energy density for both the cases. We have found that both the models, with bulk viscosity matter component, show an acceleration of the universe. We have also shown that the cosmic jerk parameter is compatible with the three kinematical data sets.


title: Propagation via Kernelization: The Vertex Cover Constraint

abstract: The technique of kernelization consists in extracting, from an instance of a problem, an essentially equivalent instance whose size is bounded in a parameter k. Besides being the basis for efficient param-eterized algorithms, this method also provides a wealth of information to reason about in the context of constraint programming. We study the use of kernelization for designing propagators through the example of the Vertex Cover constraint. Since the classic kernelization rules often correspond to dominance rather than consistency, we introduce the notion of "loss-less" kernel. While our preliminary experimental results show the potential of the approach, they also show some of its limits. In particular, this method is more effective for vertex covers of large and sparse graphs, as they tend to have, relatively, smaller kernels.


title: Identifiability and parameter estimation of the single particle lithium-ion battery model

abstract: This paper investigates the identifiability and estimation of the parameters of the single particle model (SPM) for lithium-ion battery simulation. Identifiability is addressed both in principle and in practice. The approach begins by grouping parameters and partially non-dimensionalising the SPM to determine the maximum expected degrees of freedom in the problem. We discover that, excluding open circuit voltage, there are only six independent parameters. We then examine the structural identifiability by considering whether the transfer function of the linearised SPM is unique. It is found that the model is unique provided that the electrode open circuit voltage functions have a known non-zero gradient, the parameters are ordered, and the electrode kinetics are lumped into a single charge transfer resistance parameter. We then demonstrate the practical estimation of model parameters from measured frequency-domain experimental electrochemical impedance spectroscopy (EIS) data, and show additionally that the parametrised model provides good predictive capabilities in the time domain, exhibiting a maximum voltage error of 20 mV between model and experiment over a 10 minute dynamic discharge.


title: Semiclassical Newtonian Field Theories Based On Stochastic Mechanics II

abstract: Continuing the development of the ZSM-Newton/Coulomb approach to semiclassical Newtonian gravity/electrodynamics [1], we formulate a ZSM-Newton/Coulomb version of the large N approximation scheme proposed by Oriols et al. [2]. We show that this new large N scheme makes it possible to self-consistently describe the center-of-mass evolution of a large number of gravitationally/electrostatically interacting, identical, \emph{zbw} particles, without assuming that the particles are weakly coupled, and without entailing the problematic macroscopic semiclassical gravitational/electrostatic cat states characteristic of the mean-field Schr{\"o}dinger-Newton/Coulomb equations. We also show how to recover N-particle classical Newtonian gravity/electrodynamics for many gravitationally/electrostatically interacting macroscopic particles (composed of many interacting \emph{zbw} particles), as well as classical Vlasov-Poisson mean-field theory for macroscopic particles weakly interacting gravitationally/electrostatically. Finally, we outline an explicit model of environmental decoherence that can be incorporated into Oriols et al.'s scheme as applied to ZSM-Newton/Coulomb.


title: CutFEM topology optimization of 3D laminar incompressible flow problems

abstract: This paper studies the characteristics and applicability of the CutFEM approach as the core of a robust topology optimization framework for 3D laminar incompressible flow and species transport problems at low Reynolds number (Re < 200). CutFEM is a methodology for discretizing partial differential equations on complex geometries by immersed boundary techniques. In this study, the geometry of the fluid domain is described by an explicit level set method, where the parameters of a level set function are defined as functions of the optimization variables. The fluid behavior is modeled by the incompressible Navier-Stokes equations. Species transport is modeled by an advection-diffusion equation. The governing equations are discretized in space by a generalized extended finite element method. Face-oriented ghost-penalty terms are added for stability reasons and to improve the conditioning of the system. The boundary conditions are enforced weakly via Nit\-sc\-he's method. The emergence of isolated volumes of fluid surrounded by solid during the optimization process leads to a singular analysis problem. An auxiliary indicator field is modeled to identify these volumes and to impose a constraint on the average pressure. Numerical results for 3D, steady-state and transient problems demonstrate that the CutFEM analyses are sufficiently accurate, and the optimized designs agree well with results from prior studies solved in 2D or by density approaches.


title: A simple field function for solving complex and dynamic fluid-solid system on Cartesian grid

abstract: In this paper, a simple field function is presented for facilitating the solution of complex and dynamic fluid-solid systems on Cartesian grids with interface-resolved fluid-fluid, fluid-solid, and solid-solid interactions. For a Cartesian-grid-discretized computational domain segmented by a set of solid bodies, this field function explicitly tracks each subdomain with multiple resolved interfacial node layers. As a result, the presented field function enables low-memory-cost multidomain node mapping, efficient node remapping, fast collision detection, and expedient surface force integration. Implementation algorithms for the field function and its described functionalities are also presented. Equipped with a deterministic multibody collision model, numerical experiments involving fluid-solid systems with flow conditions ranging from subsonic to supersonic states are conducted to validate and illustrate the applicability of the proposed field function.


title: A numerical stabilization framework for viscoelastic fluid flow using the finite volume method on general unstructured meshes

abstract: A robust finite volume method for viscoelastic flow analysis on general unstructured meshes is developed. It is built upon a general-purpose stabilization framework for high Weissenberg number flows. The numerical framework provides full combinatorial flexibility between different kinds of rheological models on the one hand, and effective stabilization methods on the other hand. A special emphasis is put on the velocity-stress-coupling on co-located computational grids. Using special face interpolation techniques, a semi-implicit stress interpolation correction is proposed to correct the cell-face interpolation of the stress in the divergence operator of the momentum balance. Investigating the entry-flow problem of the 4:1 contraction benchmark, we demonstrate that the numerical methods are robust over a wide range of Weissenberg numbers and significantly alleviate the high Weissenberg number problem. The accuracy of the results is evaluated in a detailed mesh convergence study.


title: Nonlinear effects in photoionization over a broad photon-energy range within the TDCIS scheme

abstract: The present tutorial provides an overview of the time-dependent configuration interaction singles (\acs{TDCIS}) scheme applied to nonlinear ionization over a broad photon-energy range. The efficient propagation of the wave function and the calculation of photoelectron spectra within this approach are described and demonstrated in various applications. Above-threshold ionization of argon and xenon in the extreme ultraviolet energy range is investigated as an example. A particular focus is put on the xenon $4d$ giant dipole resonance and the information that nonlinear ionization can provide about resonance substructure. Furthermore, above-threshold ionization is studied in the x-ray regime and the intensity regime, at which multiphoton ionization starts to play a role at hard x-ray photon energies, is identified.


title: The transition from eruptive to confined flares in the same active region

abstract: Solar flares are sudden and violent releases of magnetic energy in the solar atmosphere that can be divided in eruptive flares, when plasma is ejected from the solar atmosphere, resulting in a coronal mass ejection (CME), and confined flares when no CME is associated with the flare. We present a case-study showing the evolution of key topological structures, such as spines and fans which may determine the eruptive versus non-eruptive behavior of the series of eruptive flares, followed by confined flares, which are all originating from the same site. To study the connectivity of the different flux domains and their evolution, we compute a potential magnetic field model of the active region. Quasi-separatrix layers are retrieved from the magnetic field extrapolation. The change of behavior of the flares from one day to the next -eruptive to confined- can be attributed to the change of orientation of the magnetic field below the fan with respect to the orientation of the overlaying spine, rather than an overall change in the stability of the large scale field. Flares tend to be more-and-more confined when the field that supports the filament and the overlying field gradually become less-and-less anti-parallel, as a direct result of changes in the photospheric flux distribution, being themselves driven by continuous shearing motions of the different magnetic flux concentrations.


title: An Interface-Tracking Technique for Multiphase Flow with Soluble Surfactant

abstract: We adapt and extend a formulation for soluble surfactant transport in multiphase flows recently presented by Muradoglu & Tryggvason (JCP 274 (2014) 737-757) to the context of the Level Contour Reconstruction Method (Shin et al. IJNMF 60 (2009) 753-778) which is a hybrid method that combines the advantages of the Front-tracking and Level Set methods. Particularly close attention is paid to the formulation and numerical implementation of the surface gradients of surfactant concentration and surface tension. Various benchmark tests are performed to demonstrate the accuracy of different elements of the algorithm. To verify surfactant mass conservation, values for surfactant diffusion along the interface are compared with the exact solution for the problem of uniform expansion of a sphere. The numerical implementation of the discontinuous boundary condition for the source term in the bulk concentration is compared with the approximate solution. Surface tension forces are tested for Marangoni drop translation. Our numerical results for drop deformation in simple shear are compared with experiments and results from previous simulations. All benchmarking tests compare well with existing data thus providing confidence that our adapted LCRM formulation for surfactant advection and diffusion is accurate and effective in three-dimensional multiphase flows. We also demonstrate that this approach applies easily to massively parallel simulations.


title: On two chain models for the gravity operad

abstract: In this note we recall the construction of two chain level lifts of the gravity operad, one due to Getzler-Kapranov and one due to Westerland. We prove that these two operads are formal and that they indeed have isomorphic homology.


title: Femtosecond Self-Reconfiguration of Laser-Induced Plasma Patterns in Dielectrics

abstract: Laser-induced modification of transparent solids by intense femtosecond laser pulses allows fast integration of nanophotonic and nanofluidic devices with controlled optical properties. So far, the local and dynamic nature of the interactions between plasma and light needed to correctly explain nanograting fabrication on dielectric surfaces has been missing in the theoretical models. With our numerical approach, we show that a self-consistent dynamic treatment of the plasma formation and its interaction with light triggers an ultrafast reconfiguration of the periodic plasma patterns on a field-cycle time scale. Within this framework, a simple stability analysis of the local interactions explains how the laser-induced plasma patterns change their orientation with respect to the incident light polarization, when a certain energy density threshold is reached. Moreover, the reconfigured sub-wavelength plasma structures grow into the bulk of the sample and agree with the experimental findings of self-organized volume nanogratings. Mode coupling of the incident and transversally scattered light with the periodic plasma structures is sufficient to initiate the growth and the self-organization of the characteristic pattern with a periodicity of a half-wavelength in the medium.


title: Existence of weak solutions for a pseudo-parabolic system coupling chemical reactions, diffusion and momentum equations

abstract: We study the weak solvability of a nonlinearly coupled system of parabolic and pseudo-parabolic equations describing the interplay between mechanics, chemical reactions, diffusion and flow in a mixture theory framework. Our approach relies on suitable discrete-in-time energy-like estimates and discrete Gronwall inequalities. In selected parameter regimes, these estimates ensure the convergence of the Rothe method for the discretized partial differential equations.


title: A study on implementing a multithreaded version of the SIRENE detector simulation software for high energy neutrinos

abstract: The primary objective of SIRENE is to simulate the response to neutrino events of any type of high energy neutrino telescope. Additionally, it implements different geometries for a neutrino detector and different configurations and characteristics of photo-multiplier tubes (PMTs) inside the optical modules of the detector through a library of C+ + classes. This could be considered a massive statistical analysis of photo-electrons. Aim of this work is the development of a multithreaded version of the SIRENE detector simulation software for high energy neutrinos. This approach allows utilization of multiple CPU cores leading to a potentially significant decrease in the required execution time compared to the sequential code. We are making use of the OpenMP framework for the production of multithreaded code running on the CPU. Finally, we analyze the feasibility of a GPU-accelerated implementation.


title: The Dipole Repeller

abstract: In the standard (LCDM) model of cosmology the universe has emerged out of an early homogeneous and isotropic phase. Structure formation is associated with the growth of density irregularities and peculiar velocities. Our Local Group is moving with respect to the cosmic microwave background (CMB) with a velocity 631+/-20 km s-1 and participates in a bulk flow that extends out to distances of at least 20,000 km s-1. Since the discovery of the CMB dipole, the implicit assumption was that excesses in the abundance of galaxies induce the Local Group motion. Yet, underdense regions push as much as overdensities attract but they are deficient of light and consequently difficult to chart. It was suggested a decade ago that an underdensity in the northern hemisphere roughly 15,000 km s-1 away is a significant actor in the local flow. Here we report on kinematic evidence for such an underdensity. We map the large scale 3D velocity field using a Wiener filter reconstruction from the Cosmicflows-2 dataset of peculiar velocities, and identify the attractors and repellers that dominate the local dynamics. We show here that the local flow is dominated by a single attractor -associated with the Shapley Concentration- and a single previously unidentified repeller. Multipole expansion of the local flow provides further support for the existence and role played by the attractor and repeller. The bulk flow (i.e. dipole moment) is closely (anti)aligned with the repeller at a distance of 16,000+/-4,500 km s-1. The expansion eigenvector of the shear tensor (quadrupole moment) is closely aligned with the Shapley Attractor out to 7,000 km s-1. The close alignment of the local bulk flow with the repeller provides further support for its dominant role in shaping the local flow. This Dipole Repeller is predicted to be associated with a void in the distribution of galaxies.


title: Optimization Based Methods for Partially Observed Chaotic Systems

abstract: In this paper we consider filtering and smoothing of partially observed chaotic dynamical systems that are discretely observed, with an additive Gaussian noise in the observation. These models are found in a wide variety of real applications and include the Lorenz 96' model. In the context of a fixed observation interval $T$, observation time step $h$ and Gaussian observation variance $\sigma_Z^2$, we show under assumptions that the filter and smoother are well approximated by a Gaussian with high probability when $h$ and $\sigma^2_Z h$ are sufficiently small. Based on this result we show that the Maximum-a-posteriori (MAP) estimators are asymptotically optimal in mean square error as $\sigma^2_Z h$ tends to $0$. Given these results, we provide a batch algorithm for the smoother and filter, based on Newton's method, to obtain the MAP. In particular, we show that if the initial point is close enough to the MAP, then Newton's method converges to it at a fast rate. We also provide a method for computing such an initial point. These results contribute to the theoretical understanding of widely used 4D-Var data assimilation method. Our approach is illustrated numerically on the Lorenz 96' model with state vector up to 1 million dimensions, with code running in the order of minutes. To our knowledge the results in this paper are the first of their type for this class of models.


title: Biologically-inspired characterization of sparseness in natural images

abstract: Natural images follow statistics inherited by the structure of our physical (visual) environment. In particular, a prominent facet of this structure is that images can be described by a relatively sparse number of features. We designed a sparse coding algorithm biologically-inspired by the architecture of the primary visual cortex. We show here that coefficients of this representation exhibit a heavy-tailed distribution. For each image, the parameters of this distribution characterize sparseness and vary from image to image. To investigate the role of this sparseness, we designed a new class of random textured stimuli with a controlled sparseness value inspired by our measurements on natural images. Then, we provide with a method to synthesize random textures images with a given statistics for sparseness that matches that of some given class of natural images and provide perspectives for their use in neurophysiology.


title: Approximating the Real Structured Stability Radius with Frobenius Norm Bounded Perturbations

abstract: We propose a fast method to approximate the real stability radius of a linear dynamical system with output feedback, where the perturbations are restricted to be real valued and bounded with respect to the Frobenius norm. Our work builds on a number of scalable algorithms that have been proposed in recent years, ranging from methods that approximate the complex or real pseudospectral abscissa and radius of large sparse matrices (and generalizations of these methods for pseudospectra to spectral value sets) to algorithms for approximating the complex stability radius (the reciprocal of the $H_\infty$ norm). Although our algorithm is guaranteed to find only upper bounds to the real stability radius, it seems quite effective in practice. As far as we know, this is the first algorithm that addresses the Frobenius-norm version of this problem. Because the cost mainly consists of computing the eigenvalue with maximal real part for continuous-time systems (or modulus for discrete-time systems) of a sequence of matrices, our algorithm remains very efficient for large-scale systems provided that the system matrices are sparse.


title: Polycyclic Aromatic Hydrocarbons with Aliphatic Sidegroups: Intensity Scaling for the C-H Stretching Modes and Astrophysical Implications

abstract: The so-called unidentified infrared emission (UIE) features at 3.3, 6.2, 7.7, 8.6, and 11.3 $\mu$m ubiquitously seen in a wide variety of astrophysical regions are generally attributed to polycyclic aromatic hydrocarbon (PAH) molecules. Astronomical PAHs may have an aliphatic component as revealed by the detection in many UIE sources of the aliphatic C-H stretching feature at 3.4 $\mu$m. The ratio of the observed intensity of the 3.4 $\mu$m feature to that of the 3.3 $\mu$m aromatic C-H feature allows one to estimate the aliphatic fraction of the UIE carriers. This requires the knowledge of the intrinsic oscillator strengths of the 3.3 $\mu$m aromatic C-H stretch ($A_{3.3}$) and the 3.4 $\mu$m aliphatic C-H stretch ($A_{3.4}$). Lacking experimental data on $A_{3.3}$ and $A_{3.4}$ for the UIE candidate materials, one often has to rely on quantum-chemical computations. Although the second-order Moller-Plesset (MP2) perturbation theory with a large basis set is more accurate than the B3LYP density functional theory, MP2 is computationally very demanding and impractical for large molecules. Based on methylated PAHs, we show here that, by scaling the band strengths computed at an inexpensive level (e.g., B3LYP/6-31G*) we are able to obtain band strengths as accurate as that computed at far more expensive levels (e.g., MP2/6-311+G(3df,3pd)).


title: Local time of Levy random walks: a path integral approach

abstract: Local time of a stochastic process quantifies the amount of time that sample trajectories $x(\tau)$ spend in the vicinity of an arbitrary point $x$. For a generic Hamiltonian, we employ the phase-space path-integral representation of random walk transition probabilities in order to quantify the properties of the local time. For time-independent systems, the resolvent of the Hamiltonian operator proves to be a central tool for this purpose. In particular, we focus on local times of Levy random walks (or Levy flights), which correspond to fractional diffusion equations.


title: Hash functions using chaotic iterations

abstract: In this paper, a novel formulation of discrete chaotic iterations in the field of dynamical systems is given. Their topological properties are studied: it is mathematically proved that, under some conditions, these iterations have a chaotic behavior in the meaning of Devaney. This chaotic behavior allows us to propose a way to generate new hash functions. An illustration example is detailed in order to show how to use our theoretical study in practice.


title: A pointwise bipolar theorem

abstract: We provide a pointwise bipolar theorem for liminf-closed convex sets of positive Borel measurable functions on a sigma-compact metric space without the assumption that the polar is a tight set of measures. As applications we derive a version of the transport duality under non-tight marginals, and a superhedging duality for semistatic hedging in discrete time.


title: Anomalous Doppler Effect Singularities in the Radiative Heat Generation, Interaction Force and Frictional Torque for two Rotating Nanoparticles

abstract: We calculate the quantum heat generation, the interaction force and the frictional torque for two rotating spherical nanoparticles with a radius $R$. In contrast to the static case, when there is an upper limit in the radiative heat transfer between the particles, for two rotating nanoparticles the quantum heat generation rate diverges when the angular velocity becomes equal to the poles in the photon emission rate. These poles arise for the separation $d <d_0= R(3/\varepsilon''(\omega_0))^{1/3}$ (where $\varepsilon''(\omega_0)$ is the imaginary part of the dielectric function for the particle material at the surface phonon or plasmon polariton frequency $\omega_0$ ) due to the anomalous Doppler effect and the mutual polarization of the particles and they exist even for the particles with losses. Similar singularities exist also for the interaction force and the frictional torque. The obtained results can be important for biomedical applications.


title: To do things with words (only): An introduction to the role of noise in coordination dynamics without equations

abstract: Uncertainty, spatial or temporal errors, variability, are classic themes in the study of human and animal behaviors. Several theoretical approaches1 and concepts have been adopted to tackle those issues, often considering the CNS as an observer, using Shannon information and entropy, signal to noise ratio, and recently a Bayesian approach, and free energy minimization. In the coordination dynamics framework, addressing pattern formation processes underlying cognitive functions, and goal directed human movement among others, the tools employed originate from the statistical physics of Brownian motion and stochastic processes. The relations between those theories/ concepts have been drawn for some cases in their original fields, for example between free energy minimization and diffusion model in a force field (Jordan et al., 1998). Here an introductory presentation of the approach of noise in coordination dynamics is proposed, aimed at the experimentalist.


title: Leading twist generalized parton distributions and spin densities in a proton

abstract: We evaluate both chirally even and odd generalized parton distributions(GPDs) in the leading twist in a recently proposed quark-diquark model for the proton where the light front wavefunctions are constructed from the soft-wall AdS/QCD prediction. The GPDs in transverse impact parameter space give the spin densities for different quark and proton polarizations. For longitudinally polarized proton only chiral even GPDs contribute but for transversely polarized proton both chiral even and chiral odd GPDs contribute to the spin densities. We present a detail study of the spin densities in this model.


title: Planets in Mean-Motion Resonances and the System Around HD45364

abstract: In some planetary systems, the orbital periods of two of its members present a commensurability, usually known by mean-motion resonance. These resonances greatly enhance the mutual gravitational influence of the planets. As a consequence, these systems present uncommon behaviors, and their motions need to be studied with specific methods. Some features are unique and allow us a better understanding and characterization of these systems. Moreover, mean-motion resonances are a result of an early migration of the orbits in an accretion disk, so it is possible to derive constraints on their formation. Here we review the dynamics of a pair of resonant planets and explain how their orbits evolve in time. We apply our results to the HD 45365 planetary system.


title: Spin Filter Properties of Armchair Graphene Nanoribbons with Substitutional Fe Atoms

abstract: The spin filter capability of a (0,8) armchair graphene nanoribbon with Fe atoms at substitutional sites is investigated by density functional theory in combination with the non-equilibrium Greens function technique. For specific arrangements, a high degree of spin polarization is achieved. These include a single substitution at an edge position or double substitution in the central sector of the transmission element. The possibility of switching between majority and minority spin polarization by changing the double substitution geometry is predicted. Including the bias dependence of the transmission function proves to be essential for correct representation of the spin-resolved current-voltage profiles.


title: Relative Permittivity in the Electrical Double Layer from Nonlinear Optics

abstract: Second harmonic generation (SHG) spectroscopy has been applied to probe the fused silica/water interface at pH 7 and the uncharged 11bar02 sapphire/water interface at pH 5.2 in contact with aqueous solutions of NaCl, NaBr, NaI, KCl, RbCl, and CsCl as low as several 10 microM. For ionic strengths up to about 0.1 mM, the SHG responses were observed to increase, reversibly for all salts surveyed, when compared to the condition of zero salt added. Further increases in the salt concentration led to monotonic decreases in the SHG response. The SHG increases followed by decreases are found to be consistent with recent reports of phase interference and phase matching in nonlinear optics. By varying the relative permittivity employed in common mean field theories used to describe electrical double layers, and by comparing our results to available literature data, we find that models recapitulating the experimental observations are ones in which 1) the relative permittivity of the diffuse layer is that of bulk water, with other possible values as low as 30, 2) the surface charge density varies with salt concentration, and 3) the charge in the Stern layer or its thickness vary with salt concentration. We also note that the experimental data exhibit sensitivity depending on whether the salt concentration is increased from low to high values or decreased from high to low values, which, however, is not borne out in the fits, at least within the current uncertainties associated with the model point estimates.


title: Two-Dimensional AoD and AoA Acquisition for Wideband mmWave Systems with Cross-Polarized MIMO

abstract: In this paper, a novel two-dimensional super-resolution angle-of-departure (AoD) and angle-of-arrival (AoA) estimation technique is proposed for wideband millimeter-wave multiple-input multiple-output systems with cross-polarized antenna elements. The key ingredient of the proposed method is to form custom designed beam pairs, and devise an invertible function of the AoD/AoA to be estimated from the corresponding beam pairs. Further, a new multi-layer reference signal structure is developed for the proposed method to facilitate angle estimation for wideband channels with cross-polarized antenna elements. To facilitate feedback in closed-loop frequency division duplexing systems, a novel differential feedback strategy is proposed aiming at feedback reduction for the two-dimensional angle estimation. Numerical results demonstrate that by using the proposed method, good azimuth/elevation AoD and AoA estimation performance can be achieved under different levels of signal-to-noise ratio, channel conditions, and antenna array configurations.


title: STEREO: Search for sterile neutrinos at the ILL

abstract: Recent studies have shown that there are discrepancies between observations and the theoretical predictions in some neutrino experiments at short distances. In the so-called "Reactor Antineutrino Anomaly" and in the "Gallium Anomaly", these differences from the expectations are at the $\sim$3 $\sigma$ level in both cases. Oscillations into a light sterile neutrino state ($\Delta m^{2} \sim 1eV^{2}$) could account for the deficits in observed rates. The STEREO experiment has been conceived to confirm or reject the sterile neutrino hypothesis. It will search for an oscillation pattern at short baselines (9-11 m) in the energy spectrum of the antineutrinos emitted by the research nuclear reactor of the Institut Laue Langevin. To this end, the detector is filled with two tons of Gd-loaded liquid scintillator read out by an array of PMTs and is segmented into 6 cells in the direction of the antineutrino's propagation. STEREO should be capable of excluding the best fit parameters region of a candidate sterile neutrino at 5 $\sigma$. Data taking will start at the end of 2016 and the first physics results are expected by the first semester of 2017.


title: Molecular Simulation of Thermo-osmotic slip

abstract: Thermo-osmotic slip -- the flow induced by a thermal gradient along a surface -- is a well-known phenomenon, but curiously there is a lack of robust molecular-simulation techniques to predict its magnitude. Here, we compare three different molecular simulation techniques to compute the thermo-osmotic slip at a simple solid-fluid interface. Although we do not expect the different approaches to be in perfect agreement, we find that the differences are barely significant for a range of different physical conditions, suggesting that practical molecular simulations of thermo-osmotic slip are feasible.


title: Free energy of formation of a crystal nucleus in incongruent solidification: Implication for modeling the crystallization of aqueous nitric acid droplets in type 1 polar stratospheric clouds

abstract: Using the formalism of the classical nucleation theory, we derive an expression for the reversible work $W_*$ of formation of a binary crystal nucleus in a liquid binary solution of non-stoichiometric composition (incongruent crystallization). Applied to the crystallization of aqueous nitric acid (NA) droplets, the new expression more adequately takes account of the effect of nitric acid vapor compared to the conventional expression of MacKenzie, Kulmala, Laaksonen, and Vesala (MKLV) [J.Geophys.Res. 102, 19729 (1997)]. The predictions of both MKLV and modified expressions for the average liquid-solid interfacial tension $\sigma^{ls}$ of nitric acid dihydrate (NAD) crystals are compared by using existing experimental data on the incongruent crystallization of aqueous NA droplets of composition relevant to polar stratospheric clouds (PSCs). The predictions based on the MKLV expression are higher by about 5% compared to predictions based on our modified expression. This results in similar differences between the predictions of both expressions for the solid-vapor interfacial tension $\sigma^{sv}$ of NAD crystal nuclei. The latter can be obtained by analyzing of experimental data on crystal nucleation rates in aqueous NA droplets and exploiting the dominance of the surface-stimulated mode of crystal nucleation in small droplets and its negligibility in large ones. Applying that method, our expression for $W_*$ provides an estimate for $\sigma^{sv}$ of NAD in the range from 92 dyn/cm to 100 dyn/cm, while the MKLV expression predicts it in the range from 95 dyn/cm to 105 dyn/cm. The predictions of both expressions for $W_*$ become identical in the case of congruent crystallization; this was also demonstrated by applying our method to the nucleation of nitric acid trihydrate (NAT) crystals in PSC droplets of stoichiometric composition.


title: Scaling of geometric phase versus band structure in cluster-Ising models

abstract: We study the phase diagram of a class of models in which a generalized cluster interaction can be quenched by Ising exchange interaction and external magnetic field. We characterize the various phases through winding numbers. They may be ordinary phases with local order parameter or exotic ones, known as symmetry protected topologically ordered phases. Quantum phase transitions with dynamical critical exponents z = 1 or z = 2 are found. Quantum phase transitions are analyzed through finite-size scaling of the geometric phase accumulated when the spins of the lattice perform an adiabatic precession. In particular, we quantify the scaling behavior of the geometric phase in relation with the topology and low energy properties of the band structure of the system.


title: A Note on Prediction Markets

abstract: In a prediction market, individuals can sequentially place bets on the outcome of a future event. This leaves a trail of personal probabilities for the event, each being conditional on the current individual's private background knowledge and on the previously announced probabilities of other individuals, which give partial information about their private knowledge. By means of theory and examples, we revisit some results in this area. In particular, we consider the case of two individuals, who start with the same overall probability distribution but different private information, and then take turns in updating their probabilities. We note convergence of the announced probabilities to a limiting value, which may or may not be the same as that based on pooling their private information.


title: Motion of discrete interfaces in low-contrast random environments

abstract: We study the asymptotic behavior of a discrete-in-time minimizing movement scheme for square lattice interfaces when both the lattice spacing and the time step vanish. The motion is assumed to be driven by minimization of a weighted random perimeter functional with an additional deterministic dissipation term. We consider rectangular initial sets and lower order random perturbations of the perimeter functional. In case of stationary, $\alpha$-mixing perturbations we prove a stochastic homogenization result for the interface velocity. We also provide an example which indicates that stationary, ergodic perturbations do not yield a spatially homogenized limit velocity for this minimizing movement scheme.


title: Multi-wavelength variability study of the classical BL Lac object PKS 0735+178 on timescales ranging from decades to minutes

abstract: We present the results of our power spectral analysis for the BL Lac object PKS 0735+178 utilizing the Fermi-LAT survey at high-energy $\gamma$-rays, several ground-based optical telescopes, and single-dish radio telescopes operating at GHz frequencies. The novelty of our approach is that, by combining long-term and densely sampled intra-night light curves in the optical regime, we were able to construct for the first time the optical power spectrum of the blazar for a time domain extending from 23 years down to minutes. Our analysis reveals that: (i) the optical variability is consistent with a pure red noise, for which the power spectral density can well be approximated by a single power-law throughout the entire time domain probed; (ii) the slope of power spectral density at high-energy $\gamma$-rays ($\sim 1$), is significantly flatter than that found at radio and optical frequencies ($\sim 2$) within the corresponding time variability range; (iii) for the derived power spectra we did not detect any low-frequency flattening, nor do we see any evidence for cut-offs at the highest frequencies down to the noise floor levels due to measurement uncertainties. We interpret our findings in terms of a model where the blazar variability is generated by the underlying single stochastic process (at radio and optical frequencies), or a linear superposition of such processes (in the $\gamma$-ray regime). Along with the detailed PSD analysis, we also present the results of our extended (1998-2015) intra-night optical monitoring program and newly acquired optical photo-polarimetric data for the source.


title: Inertial Proximal Alternating Linearized Minimization (iPALM) for Nonconvex and Nonsmooth Problems

abstract: In this paper we study nonconvex and nonsmooth optimization problems with semi-algebraic data, where the variables vector is split into several blocks of variables. The problem consists of one smooth function of the entire variables vector and the sum of nonsmooth functions for each block separately. We analyze an inertial version of the Proximal Alternating Linearized Minimization (PALM) algorithm and prove its global convergence to a critical point of the objective function at hand. We illustrate our theoretical findings by presenting numerical experiments on blind image deconvolution, on sparse non-negative matrix factorization and on dictionary learning, which demonstrate the viability and effectiveness of the proposed method.


title: Nichols algebras that are quantum planes

abstract: We compute all Nichols algebras of rigid vector spaces of dimension 2 that admit a non-trivial quadratic relation.


title: Nuclear cluster structure effect on elliptic and triangular flows in heavy-ion collisions

abstract: The initial geometry effect on collective flows, which are inherited from initial projectile structure, is studied in relativistic heavy-ion collisions of $^{12}$C + Au by using a multi-phase transport model (AMPT). Elliptic flow ($v_2$) and triangular flow ($v_3$) which are significantly resulted from the chain and triangle structure of $^{12}\mathrm{C}$ with three-$\alpha$ clusters, respectively, in central $^{12}\mathrm{C}$+$^{197}\mathrm{Au}$ collisions are compared with the flow from the Woods-Saxon distribution of nucleons in $^{12}\mathrm{C}$. $v_3/v_2$ is proposed as a probe to distinguish the pattern of $\alpha$-clustered $^{12}\mathrm{C}$. This study demonstrates that the initial geometry of the collision zone inherited from nuclear structure can be explored by collective flow at the final stage in heavy-ion collisions.


title: Computational Techniques in Multispectral Image Processing: Application to the Syriac Galen Palimpsest

abstract: Multispectral and hyperspectral image analysis has experienced much development in the last decade. The application of these methods to palimpsests has produced significant results, enabling researchers to recover texts that would be otherwise lost under the visible overtext, by improving the contrast between the undertext and the overtext. In this paper we explore an extended number of multispectral and hyperspectral image analysis methods, consisting of supervised and unsupervised dimensionality reduction techniques, on a part of the Syriac Galen Palimpsest dataset (www.digitalgalen.net). Of this extended set of methods, eight methods gave good results: three were supervised methods Generalized Discriminant Analysis (GDA), Linear Discriminant Analysis (LDA), and Neighborhood Component Analysis (NCA); and the other five methods were unsupervised methods (but still used in a supervised way) Gaussian Process Latent Variable Model (GPLVM), Isomap, Landmark Isomap, Principal Component Analysis (PCA), and Probabilistic Principal Component Analysis (PPCA). The relative success of these methods was determined visually, using color pictures, on the basis of whether the undertext was distinguishable from the overtext, resulting in the following ranking of the methods: LDA, NCA, GDA, Isomap, Landmark Isomap, PPCA, PCA, and GPLVM. These results were compared with those obtained using the Canonical Variates Analysis (CVA) method on the same dataset, which showed remarkably accuracy (LDA is a particular case of CVA where the objects are classified to two classes).


title: One side continuity of meromorphic mappings between real analytic hypersurfaces

abstract: We prove that a meromorphic mapping, which sends a peace of a real analytic strictly pseudoconvex hypersurface in $\cc^2$ to a compact subset of $\cc^N$ which doesn't contain germs of non-constant complex curves is continuous from the concave side of the hypersurface. This implies the analytic continuability along CR-paths of germs of holomorphic mappings from real analytic hypersurfaces with non-vanishing Levi form to the locally spherical ones in all dimensions.


title: AMIGOS: A Dataset for Affect, Personality and Mood Research on Individuals and Groups

abstract: We present AMIGOS-- A dataset for Multimodal research of affect, personality traits and mood on Individuals and GrOupS. Different to other databases, we elicited affect using both short and long videos in two social contexts, one with individual viewers and one with groups of viewers. The database allows the multimodal study of the affective responses, by means of neuro-physiological signals of individuals in relation to their personality and mood, and with respect to the social context and videos' duration. The data is collected in two experimental settings. In the first one, 40 participants watched 16 short emotional videos. In the second one, the participants watched 4 long videos, some of them alone and the rest in groups. The participants' signals, namely, Electroencephalogram (EEG), Electrocardiogram (ECG) and Galvanic Skin Response (GSR), were recorded using wearable sensors. Participants' frontal HD video and both RGB and depth full body videos were also recorded. Participants emotions have been annotated with both self-assessment of affective levels (valence, arousal, control, familiarity, liking and basic emotions) felt during the videos as well as external-assessment of levels of valence and arousal. We present a detailed correlation analysis of the different dimensions as well as baseline methods and results for single-trial classification of valence and arousal, personality traits, mood and social context. The database is made publicly available.


title: Design of Stochastic Robotic Swarms for Target Performance Metrics in Boundary Coverage Tasks

abstract: In this work, we analyze \textit{stochastic coverage schemes} (SCS) for robotic swarms in which the robots randomly attach to a one-dimensional boundary of interest using local communication and sensing, without relying on global position information or a map of the environment. Robotic swarms may be required to perform boundary coverage in a variety of applications, including environmental monitoring, collective transport, disaster response, and nanomedicine. We present a novel analytical approach to computing and designing the statistical properties of the communication and sensing networks that are formed by random robot configurations on a boundary. We are particularly interested in the event that a robot configuration forms a connected communication network or maintains continuous sensor coverage of the boundary. Using tools from order statistics, random geometric graphs, and computational geometry, we derive formulas for properties of the random graphs generated by robots that are independently and identically distributed along a boundary. We also develop order-of-magnitude estimates of these properties based on Poisson approximations and threshold functions. For cases where the SCS generates a uniform distribution of robots along the boundary, we apply our analytical results to develop a procedure for computing the robot population size, diameter, sensing range, or communication range that yields a random communication network or sensor network with desired properties.


title: Semi-Dense Visual Odometry for RGB-D Cameras Using Approximate Nearest Neighbour Fields

abstract: This paper presents a robust and efficient semi-dense visual odometry solution for RGB-D cameras. The core of our method is a 2D-3D ICP pipeline which estimates the pose of the sensor by registering the projection of a 3D semi-dense map of the reference frame with the 2D semi-dense region extracted in the current frame. The processing is speeded up by efficiently implemented approximate nearest neighbour fields under the Euclidean distance criterion, which permits the use of compact Gauss-Newton updates in the optimization. The registration is formulated as a maximum a posterior problem to deal with outliers and sensor noises, and consequently the equivalent weighted least squares problem is solved by iteratively reweighted least squares method. A variety of robust weight functions are tested and the optimum is determined based on the characteristics of the sensor model. Extensive evaluation on publicly available RGB-D datasets shows that the proposed method predominantly outperforms existing state-of-the-art methods.


title: Dynamics and stability of surface waves with bulk-soluble surfactants

abstract: In this paper we study the dynamics of a layer of incompressible viscous fluid bounded below by a rigid boundary and above by a free boundary, in the presence of a uniform gravitational field. We assume that a mass of surfactant is present both at the free surface and in the bulk of fluid, and that conversion from one species to the other is possible. The surfactants couple to the fluid dynamics through the coefficient of surface tension, which depends on the the surface density of surfactants. Gradients in this concentration give rise to Marangoni stress on the free surface. In turn, the fluids advect the surfactants and distort their concentration through geometric distortions of the free surface. We model the surfactants in a way that allows absorption and desorption of surfactant between the surface and bulk. We prove that small perturbations of the equilibrium solutions give rise to global-in-time solutions that decay to equilibrium at an exponential rate. This establishes the asymptotic stability of the equilibrium solutions.


title: Monocular LSD-SLAM Integration within AR System

abstract: In this paper, we cover the process of integrating Large-Scale Direct Simultaneous Localization and Mapping (LSD-SLAM) algorithm into our existing AR stereo engine, developed for our modified "Augmented Reality Oculus Rift". With that, we are able to track one of our realworld cameras which are mounted on the rift, within a complete unknown environment. This makes it possible to achieve a constant and full augmentation, synchronizing our 3D movement (x, y, z) in both worlds, the real world and the virtual world. The development for the basic AR setup using the Oculus Rift DK1 and two fisheye cameras is fully documented in our previous paper. After an introduction to image-based registration, we detail the LSD-SLAM algorithm and document our code implementing our integration. The AR stereo engine with Oculus Rift support can be accessed via the GIT repository https://github.com/MaXvanHeLL/ARift.git and the modified LSD-SLAM project used for the integration is available here https://github.com/MaXvanHeLL/LSD-SLAM.git.


title: Thermodynamical Effects and High Resolution Methods for Compressible Fluid Flows

abstract: One of the fundamental differences of compressible fluid flows from incompressible fluid flows is the involvement of thermodynamics. This difference should be manifested in the design of numerical methods and seems often be neglected in addition that the entropy inequality, as a conceptual derivative, is taken into account to reflect irreversible processes and verified for some first order schemes. In this paper, we refine the GRP solver to illustrate how the thermodynamical variation is integrated into the design of high resolution methods for compressible fluid flows and demonstrate numerically the importance of thermodynamic effect in the resolution of strong waves. As a by-product, we show that the GRP solver works for generic equations of state, and is independent of technical arguments.


title: Low Cost Monitoring and Intruders Detection using Wireless Video Sensor Networks

abstract: There is a growing interest in the use of video sensor networks in surveillance applications in order to detect intruders with low cost. The essential concern of such networks is whether or not a specified target can pass or intrude the monitored region without being detected. This concern forms a serious challenge to wireless video sensor networks of weak computation and battery power. In this paper, our aim is to prolong the whole network lifetime while fulfilling the surveillance application needs. We present a novel scheduling algorithm where only a subset of video nodes contribute significantly to detect intruders and prevent malicious attacker to predict the behavior of the network prior to intrusion. Our approach is chaos-based, where every node based on its last detection, a hash value and some pseudo-random numbers easily computes a decision function to go to sleep or active mode. We validate the efficiency of our approach through theoretical analysis and demonstrate the benefits of our scheduling algorithm by simulations. Results show that in addition of being able to increase the whole network lifetime and to present comparable results against random attacks (low stealth time), our scheme is also able to withstand malicious attacks due to its fully unpredictable behavior.


title: Propagation of bursting oscillations in coupled Hodgkin-Huxley Reaction-Diffusion systems

abstract: We consider networks of reaction-diffusion systems of Hodgkin-Huxley type. We give a general mathematical framework, in which we prove existence and unicity of solutions as well as existence of invariant regions and of the attractor. Then, we illustrate some relevant numerical examples and exhibit bifurcation phenomena and propagation of bursting oscillations through one and two coupled systems.


title: Non-local field-like spin-orbit torques in Rashba systems: an ab-initio study of Ag$_{2}$Bi/Ag/Fe film

abstract: We investigate from first principles the field-like spin-orbit torques (SOTs) in a Ag$_{2}$Bi-terminated Ag(111) film grown on ferromagnetic Fe(110). We find that a large part of the SOT arises from the spin-orbit interaction (SOI) in the Ag$_{2}$Bi layer far away from the Fe layers. These results clearly hint at a long range spin transfer in the direction perpendicular to the film that does not originate in the spin Hall effect. In order to bring evidence of the non-local character of the computed SOT, we show that the torque acting on the Fe layers can be engineered by the introduction of Bi vacancies in the Ag$_{2}$Bi layer. Overall, we find a drastic dependence of the SOT on the disorder type, which we explain by a complex interplay of different contributions to the SOT in the Brillouin zone.


abstract: We present Deep Generalized Canonical Correlation Analysis (DGCCA) -- a method for learning nonlinear transformations of arbitrarily many views of data, such that the resulting transformations are maximally informative of each other. While methods for nonlinear two-view representation learning (Deep CCA, (Andrew et al., 2013)) and linear many-view representation learning (Generalized CCA (Horst, 1961)) exist, DGCCA is the first CCA-style multiview representation learning technique that combines the flexibility of nonlinear (deep) representation learning with the statistical power of incorporating information from many independent sources, or views. We present the DGCCA formulation as well as an efficient stochastic optimization algorithm for solving it. We learn DGCCA representations on two distinct datasets for three downstream tasks: phonetic transcription from acoustic and articulatory measurements, and recommending hashtags and friends on a dataset of Twitter users. We find that DGCCA representations soundly beat existing methods at phonetic transcription and hashtag recommendation, and in general perform no worse than standard linear many-view techniques.


title: Recovering the HII region size statistics from 21-cm tomography

abstract: We introduce a novel technique, called "granulometry", to characterize and recover the mean size and the size distribution of HII regions from 21-cm tomography. The technique is easy to implement, but places the previously not very well defined concept of morphology on a firm mathematical foundation. The size distribution of the cold spots in 21-cm tomography can be used as a direct tracer of the underlying probability distribution of HII region sizes. We explore the capability of the method using large-scale reionization simulations and mock observational data cubes while considering capabilities of SKA1-low and a future extension to SKA2. We show that the technique allows the recovery of the HII region size distribution with a moderate signal-to-noise ratio from wide-field imaging ($\rm SNR\lesssim3$), for which the statistical uncertainty is sample variance dominated. We address the observational requirements on the angular resolution, the field-of-view, and the thermal noise limit for a successful measurement. To achieve a full scientific return from 21-cm tomography and to exploit a synergy with 21-cm power spectra, we suggest an observing strategy using wide-field imaging (several tens of square degrees) by an interferometric mosaicking/multi-beam observation with additional intermediate baselines (~2-4 km).


title: Topological protection from exceptional points in Weyl and nodal line semimetals

abstract: We investigate the topological protection of surface states in Weyl and nodal-line semimetals by characterizing them as evanescent states when the band structure is extended to complex momenta. We find in this way a sequence of exceptional points -that is, branch points with zero energy in the complex spectrum- allowing us to identify the set of surface states with complex momentum signaling the decay into the 3D semimetal. From this point of view, Weyl and nodal-line semimetals can be classified in two types depending on the way surface states decay. Type A semimetals have surface states with smaller penetration length and oscillating decay while type B semimetals have longer simple exponential decays. The difference between both types reflects in the way the branch cuts in the spectrum accommodate in the complex plane. The stability of the surface states stems in this approach from the complex structure that develops around the exceptional points, with a topological protection which is based on the fact that the branch cuts cannot be closed by small perturbations. We check this property when nodal-line semimetals are placed under circularly polarized light, where we observe that the exceptional points survive the effect of such a perturbation, though appropriate boundary conditions for zero-energy surface states cannot be satisfied in general due to the breakdown of time-reversal invariance by the radiation field.

categories: cond-mat.mes-hall math-ph math.MP quant-ph


title: Stable high-pressure phases in the H-S system determined by chemically reacting hydrogen and sulfur

abstract: Synchrotron X-ray diffraction and Raman spectroscopy have been used to study chemical reactions of molecular hydrogen (H2) with sulfur (S) at high pressures. We find theoretically predicted Cccm and Im-3m H3S to be the reaction products at 50 and 140 GPa, respectively. Im-3m H3S is a stable crystalline phase above 140 GPa and it transforms to R3m H3S on pressure release below 140 GPa. The latter phase is (meta)stable down to at least 70 GPa where it transforms to Cccm H3S upon annealing (T<1300 K) to overcome the kinetic hindrance. Cccm H3S has an extended structure with symmetric hydrogen bonds at 50 GPa and upon decompression it experiences a transformation to a molecular mixed H2S-H2 structure below 40 GPa without any apparent change in the crystal symmetry.


title: The nonlinear Schr\"odinger Equation driven by jump processes

abstract: The main result of the paper is the existence of a solution of the nonlinear Schr\"odinger equation with a \levy noise with infinite activity. To be more precise, let $A=\Delta$ be the Laplace operator with $D(A)=\{ u\in L ^2 (\mathbb{R} ^d): \Delta u \in L ^2 (\mathbb{R} ^d)\}$. Let $Z\hookrightarrow L ^2(\mathbb{R} ^d)$ be a function space and $\eta$ be a Poisson random measure on $Z$, let $g:\mathbb{R}\to\mathbb{C}$ and $h:\mathbb{R}\to\mathbb{C}$ be some given functions, satisfying certain conditions specified later. Let $\alpha\ge 1$ and $\lambda\ge 0$. We are interested in the solution of the following equation % $$ i \, d u(t,x) - \Delta u(t,x)\,dt +\lambda |u(t,x)|^{\alpha-1} u(t,x) \, dt $$ $$= \int_Z u(t,x)\, g(z(x))\,\tilde \eta (dz,dt)+\int_Z u(t,x)\, h (z(x))\, \gamma (dz, dt), $$ $$ u(0)= u_0. $$ First we consider the case, where the \levy process is a compound Poisson process. With the help of this result we can tackle the general case, and show that the equation above has a solution.


title: Flashes of Hidden Worlds at Colliders

abstract: (This is a general physics level overview article about hidden sectors, and how they motivate searches for long-lived particles. Intended for publication in Physics Today.) Searches for new physics at the Large Hadron Collider have so far come up empty, but we just might not be looking in the right place. Spectacular bursts of particles appearing seemingly out of nowhere could shed light on some of nature's most profound mysteries.


title: Measuring Majorana non-locality and spin structure with a quantum dot

abstract: Robust zero bias transport anomalies in semiconducting nanowires with proximity-induced superconductivity have been convincingly demonstrated in various experiments. While these are compatible with the existence of Majorana zero modes at the ends of the nanowire, a direct proof of their non-locality and topological protection is now needed. Here we show that a quantum dot at the end of the nanowire may be used as a powerful spectroscopic tool to quantify the degree of Majorana non-locality through a local transport measurement. Moreover, the spin polarization of dot sub-gap states at singlet-doublet transitions in the Coulomb blockade regime allows the dot to directly probe the spin structure of the Majorana wave function, and indirectly measure the spin-orbit coupling of the nanowire.


title: Deep Kernelized Autoencoders

abstract: In this paper we introduce the deep kernelized autoencoder, a neural network model that allows an explicit approximation of (i) the mapping from an input space to an arbitrary, user-specified kernel space and (ii) the back-projection from such a kernel space to input space. The proposed method is based on traditional autoencoders and is trained through a new unsupervised loss function. During training, we optimize both the reconstruction accuracy of input samples and the alignment between a kernel matrix given as prior and the inner products of the hidden representations computed by the autoencoder. Kernel alignment provides control over the hidden representation learned by the autoencoder. Experiments have been performed to evaluate both reconstruction and kernel alignment performance. Additionally, we applied our method to emulate kPCA on a denoising task obtaining promising results.


title: Control of multidimensional systems on complex network

abstract: Multidimensional systems coupled via complex networks are widespread in nature and thus frequently invoked for a large plethora of interesting applications. From ecology to physics, individual entities in mutual interactions are grouped in families, homogeneous in kind. These latter interact selectively, through a sequence of self-consistently regulated steps, whose deeply rooted architecture is stored in the assigned matrix of connections. The asymptotic equilibrium eventually attained by the system, and its associated stability, can be assessed by employing standard nonlinear dynamics tools. For many practical applications, it is however important to externally drive the system towards a desired equilibrium, which is resilient, hence stable, to external perturbations. To this end we here consider a system made up of $N$ interacting populations which evolve according to general rate equations, bearing attributes of universality. One species is added to the pool of interacting families and used as a dynamical controller to induce novel stable equilibria. Use can be made of the root locus method to shape the needed control, in terms of intrinsic reactivity and adopted protocol of injection. The proposed method is tested on both synthetic and real data, thus enabling to demonstrate its robustness and versatility.

categories: physics.soc-ph cond-mat.stat-mech nlin.AO physics.bio-ph


title: A Complete Quantitative Deduction System for the Bisimilarity Distance on Markov Chains

abstract: In this paper we propose a complete axiomatization of the bisimilarity distance of Desharnais et al. for the class of finite labelled Markov chains. Our axiomatization is given in the style of a quantitative extension of equational logic recently proposed by Mardare, Panangaden, and Plotkin (LICS 2016) that uses equality relations $t \equiv_\varepsilon s$ indexed by rationals, expressing that `$t$ is approximately equal to $s$ up to an error $\varepsilon$'. Notably, our quantitative deduction system extends in a natural way the equational system for probabilistic bisimilarity given by Stark and Smolka by introducing an axiom for dealing with the Kantorovich distance between probability distributions. The axiomatization is then used to propose a metric extension of a Kleene's style representation theorem for finite labelled Markov chains, that was proposed (in a more general coalgebraic fashion) by Silva et al. (Inf. Comput. 2011).


title: Representing the deformation $\infty$-groupoid

abstract: The goal of the present paper is to introduce a smaller, but equivalent version of the Deligne-Hinich-Getzler $\infty$-groupoid associated to a homotopy Lie algebra. In the case of differential graded Lie algebras, we represent it by a universal cosimplicial object.


title: Learning detectors of malicious web requests for intrusion detection in network traffic

abstract: This paper proposes a generic classification system designed to detect security threats based on the behavior of malware samples. The system relies on statistical features computed from proxy log fields to train detectors using a database of malware samples. The behavior detectors serve as basic reusable building blocks of the multi-level detection architecture. The detectors identify malicious communication exploiting encrypted URL strings and domains generated by a Domain Generation Algorithm (DGA) which are frequently used in Command and Control (C&C), phishing, and click fraud. Surprisingly, very precise detectors can be built given only a limited amount of information extracted from a single proxy log. This way, the computational requirements of the detectors are kept low which allows for deployment on a wide range of security devices and without depending on traffic context such as DNS logs, Whois records, webpage content, etc. Results on several weeks of live traffic from 100+ companies having 350k+ hosts show correct detection with a precision exceeding 95% of malicious flows, 95% of malicious URLs and 90% of infected hosts. In addition, a comparison with a signature and rule-based solution shows that our system is able to detect significant amount of new threats.


title: The theory of parametrically amplified electron-phonon superconductivity

abstract: The ultrafast optical manipulation of ordered phases in strongly correlated materials is a topic of significant theoretical, experimental, and technological interest. Inspired by a recent experiment on light-induced superconductivity in fullerenes [Mitrano et al., Nature 530, 2016], we develop a comprehensive theory of light-induced superconductivity in driven electron-phonon systems with lattice nonlinearities. In analogy with the operation of parametric amplifiers, we show how the interplay between the external drive and lattice nonlinearities lead to significantly enhanced effective electron-phonon couplings. We provide a detailed and unbiased study of the nonequilibrium dynamics of the driven system using the real-time Green's function technique. To this end, we develop a Floquet generalization of the Migdal-Eliashberg theory and derive a numerically tractable set of quantum Floquet-Boltzmann kinetic equations for the coupled electron-phonon system. We study the role of parametric phonon generation and electronic heating in destroying the transient superconducting state. Finally, we predict the transient formation of electronic Floquet bands in time- and angle-resolved photo-emission spectroscopy experiments as a consequence of the proposed mechanism.