title: Investigation of Dark Matter in the 3-2-3-1 Model

abstract: We prove that the $SU(3)_C\otimes SU(2)_L \otimes SU(3)_R\otimes U(1)_X$ (3-2-3-1) gauge model always contains a matter parity $W_P=(-1)^{3(B-L)+2s}$ as conserved residual gauge symmetry, where $B-L=2(\beta T_{8R}+X)$ is a $SU(3)_R\otimes U(1)_X$ charge. Due to the non-Abelian nature of $B-L$, the $W$-odd and $W$-even fields are actually unified in gauge multiplets. We investigate two viable versions for dark matter according to $\beta=\pm1/\sqrt{3}$, where the dark matter candidates can be fermion, scalar, or vector fields. We figure out the parameter spaces in the allowed regions of the relic density and direct detection cross-sections. Additionally, we examine the neutrino masses induced by the seesaw mechanism along with associated lepton flavor violation processes. The new gauge boson searches at the LEPII and LHC are discussed.


title: Survivable Network Design for Group Connectivity in Low-Treewidth Graphs

abstract: In the Group Steiner Tree problem (GST), we are given a (vertex or edge)-weighted graph $G=(V,E)$ on $n$ vertices, a root vertex $r$ and a collection of groups $\{S_i\}_{i\in[h]}: S_i\subseteq V(G)$. The goal is to find a min-cost subgraph $H$ that connects the root to every group. We consider a fault-tolerant variant of GST, which we call Restricted (Rooted) Group SNDP. In this setting, each group $S_i$ has a demand $k_i\in[k],k\in\mathbb N$, and we wish to find a min-cost $H\subseteq G$ such that, for each group $S_i$, there is a vertex in $S_i$ connected to the root via $k_i$ (vertex or edge) disjoint paths. While GST admits $O(\log^2 n\log h)$ approximation, its high connectivity variants are Label-Cover hard, and for the vertex-weighted version, the hardness holds even when $k=2$. Previously, positive results were known only for the edge-weighted version when $k=2$ [Gupta et al., SODA 2010; Khandekar et al., Theor. Comput. Sci., 2012] and for a relaxed variant where the disjoint paths may end at different vertices in a group [Chalermsook et al., SODA 2015]. Our main result is an $O(\log n\log h)$ approximation for Restricted Group SNDP that runs in time $n^{f(k, w)}$, where $w$ is the treewidth of $G$. This nearly matches the lower bound when $k$ and $w$ are constant. The key to achieving this result is a non-trivial extension of the framework in [Chalermsook et al., SODA 2017], which embeds all feasible solutions to the problem into a dynamic program (DP) table. However, finding the optimal solution in the DP table remains intractable. We formulate a linear program relaxation for the DP and obtain an approximate solution via randomized rounding. This framework also allows us to systematically construct DP tables for high-connectivity problems. As a result, we present new exact algorithms for several variants of survivable network design problems in low-treewidth graphs.


title: Measurement of the full excitation spectrum of the 7Li(p,{\gamma}){\alpha}{\alpha} reaction at 441 keV

abstract: A current challenge for ab initio calculations is systems that contain large continuum contributions such as 8Be. We report on new measurements of radiative decay widths in this nucleus that test recent Green's function Monte Carlo calculations. Traditionally, {\gamma} ray detectors have been utilized to measure the high energy photons from the 7Li(p, {\gamma}){\alpha}{\alpha} reaction. However, due to the complicated response function of these detectors it has not yet been possible to extract the full {\gamma} ray spectrum from this reaction. Here we present an alternative measurement using large area Silicon detectors to detect the two {\alpha} particles, which provides a practically background free spectrum and retains good energy resolution. The resulting spectrum is analyzed using a many-level multi channel R-matrix parametrization. Improved values for the radiative widths are extracted from the R-matrix fit. We find evidence for significant non-resonant continuum contributions and tentative evidence for a broad 0+ resonance at 12 MeV.


title: Massive MIMO for Ultra-reliable Communications with Constellations for Dual Coherent-noncoherent Detection

abstract: The stringent requirements of ultra-reliable low-latency communications (URLLC) require rethinking of the physical layer transmission techniques. Massive antenna arrays are seen as an enabler of the emerging $5^\text{th}$ generation systems, due to increases in spectral efficiency and degrees of freedom for transmissions, which can greatly improve reliability under demanding latency requirements. Massive array coherent processing relies on accurate channel state information (CSI) in order to achieve high reliability. In this paper, we investigate the impact of imperfect CSI in a single-input multiple-output (SIMO) system on the coherent receiver. An amplitude-phase keying (APK) symbol constellation is proposed, where each two symmetric symbols reside on distinct power levels. The symbols are demodulated using a dual-stage non-coherent and coherent detection strategy, in order to improve symbol reliability. By means of analysis and simulation, we find an adequate scaling of the constellation and show that for high signal-to-noise ratio (SNR) and inaccurate CSI regime, the proposed scheme enhances receiver performance.


title: The Underlying Mechanisms of Time Dilation and Doppler Effect in Curved Space-Time

abstract: In this paper, we theoretically investigate the time dilation and Doppler effect in curved space-time from the perspective of quantum field theory (QFT). A Coordinate Transformation which Maintains the Period of Clocks is introduced, and such coordinate transformation is named as CTMPC throughout this paper. By analogy with the Lorentz transformation in Minkowski space-time, CTMPC is a correct transformation in curved space-times in a sense that it shows the correct relation between the time measured by the two observers, moreover, Lorentz transformation is just a special case of CTMPC applied in Minkowski space-time. We demonstrate that the Coordinate Transformation which Maintains the Local Metric (CTMLM) is one CTMPC, while the mathematical forms of physics formulas in QFT will be maintained. As applications of CTMLM, the time dilation and Doppler effect with an arbitrary time-dependent relative velocity in curved space-time are analysed. For Minkowski space-time, the time dilation and Doppler effect agree with the clock hypothesis. For curved space-time, we show that even if the emitted wave has a narrow frequency range, the Doppler effect may, in general, broaden the frequency spectrum and, at the meantime, shift the frequencies values. These new findings will deepen our understanding on the nature of space-time and the Doppler effect in curved space-time, they may also provide theoretical guidance in future astronomical observations.


title: Short Packet Structure for Ultra-Reliable Machine-type Communication: Tradeoff between Detection and Decoding

abstract: Machine-type communication requires rethinking of the structure of short packets due to the coding limitations and the significant role of the control information. In ultra-reliable low-latency communication (URLLC), it is crucial to optimally use the limited degrees of freedom (DoFs) to send data and control information. We consider a URLLC model for short packet transmission with acknowledgement (ACK). We compare the detection/decoding performance of two short packet structures: (1) time-multiplexed detection sequence and data; and (2) structure in which both packet detection and data decoding use all DoFs. Specifically, as an instance of the second structure we use superimposed sequences for detection and data. We derive the probabilities of false alarm and misdetection for an AWGN channel and numerically minimize the packet error probability (PER), showing that for delay-constrained data and ACK exchange, there is a tradeoff between the resources spent for detection and decoding. We show that the optimal PER for the superimposed structure is achieved for higher detection overhead. For this reason, the PER is also higher than in the preamble case. However, the superimposed structure is advantageous due to its flexibility to achieve optimal operation without the need to use multiple codebooks.


title: A Neurorobotic Experiment for Crossmodal Conflict Resolution in Complex Environments

abstract: Crossmodal conflict resolution is crucial for robot sensorimotor coupling through the interaction with the environment, yielding swift and robust behaviour also in noisy conditions. In this paper, we propose a neurorobotic experiment in which an iCub robot exhibits human-like responses in a complex crossmodal environment. To better understand how humans deal with multisensory conflicts, we conducted a behavioural study exposing 33 subjects to congruent and incongruent dynamic audio-visual cues. In contrast to previous studies using simplified stimuli, we designed a scenario with four animated avatars and observed that the magnitude and extension of the visual bias are related to the semantics embedded in the scene, i.e., visual cues that are congruent with environmental statistics (moving lips and vocalization) induce the strongest bias. We implement a deep learning model that processes stereophonic sound, facial features, and body motion to trigger a discrete behavioural response. After training the model, we exposed the iCub to the same experimental conditions as the human subjects, showing that the robot can replicate similar responses in real time. Our interdisciplinary work provides important insights into how crossmodal conflict resolution can be modelled in robots and introduces future research directions for the efficient combination of sensory observations with internally generated knowledge and expectations.


title: Solving determinantal systems using homotopy techniques

abstract: Let $\K$ be a field of characteristic zero and $\Kbar$ be an algebraic closure of $\K$. Consider a sequence of polynomials$G=(g\_1,\dots,g\_s)$ in $\K[X\_1,\dots,X\_n]$, a polynomial matrix $\F=[f\_{i,j}] \in \K[X\_1,\dots,X\_n]^{p \times q}$, with $p \leq q$,and the algebraic set $V\_p(F, G)$ of points in $\KKbar$ at which all polynomials in $\G$ and all $p$-minors of $\F$vanish. Such polynomial systems appear naturally in e.g. polynomial optimization, computational geometry.We provide bounds on the number of isolated points in $V\_p(F, G)$ depending on the maxima of the degrees in rows (resp. columns) of $\F$. Next, we design homotopy algorithms for computing those points. These algorithms take advantage of the determinantal structure of the system defining $V\_p(F, G)$. In particular, the algorithms run in time that is polynomial in the bound on the number of isolated points.


title: Tensor Decomposition for Compressing Recurrent Neural Network

abstract: In the machine learning fields, Recurrent Neural Network (RNN) has become a popular architecture for sequential data modeling. However, behind the impressive performance, RNNs require a large number of parameters for both training and inference. In this paper, we are trying to reduce the number of parameters and maintain the expressive power from RNN simultaneously. We utilize several tensor decompositions method including CANDECOMP/PARAFAC (CP), Tucker decomposition and Tensor Train (TT) to re-parameterize the Gated Recurrent Unit (GRU) RNN. We evaluate all tensor-based RNNs performance on sequence modeling tasks with a various number of parameters. Based on our experiment results, TT-GRU achieved the best results in a various number of parameters compared to other decomposition methods.


title: Distance entropy cartography characterises centrality in complex networks

abstract: We introduce distance entropy as a measure of homogeneity in the distribution of path lengths between a given node and its neighbours in a complex network. Distance entropy defines a new centrality measure whose properties are investigated for a variety of synthetic network models. By coupling distance entropy information with closeness centrality, we introduce a network cartography which allows one to reduce the degeneracy of ranking based on closeness alone. We apply this methodology to the empirical multiplex lexical network encoding the linguistic relationships known to English speaking toddlers. We show that the distance entropy cartography better predicts how children learn words compared to closeness centrality. Our results highlight the importance of distance entropy for gaining insights from distance patterns in complex networks.

categories: physics.soc-ph cond-mat.stat-mech cs.CL cs.SI physics.data-an


title: Reaction-diffusion on random spatial networks with scale-free jumping rates via effective medium theory

abstract: We study epidemic processes using a metapopulation approach on the line featuring random transport rates between arbitrarily distant sites. An average transport network is found using a recently developed variant of the effective medium approximation (EMA) that is capable of dealing with these long-range connections. Using a Feynman-Kac argument in the effective medium, we derive an estimate on the size of the infected domain, and reproduce the known result of its exponential growth in time. We hereby demonstrate the applicability of long-range EMA to dynamical processes on networks more intricate than simple diffusion.


title: Towards Open Boundary Molecular Dynamics Simulation of Ionic Liquids

abstract: We extend the use of the adaptive resolution method (AdResS) in its Grand Canonical-like version (GC-AdResS) to the molecular dynamics simulation of 1,3-dimethylimidazolium chloride. We show that the partitioning of the total system in a subsystem of interest with atomistic details and a reservoir of coarse-grained particles leads to satisfactory results. The challenging aspect of this study, compared to previous AdResS simulations, is the presence of charged particles and the necessity of addressing the question about the minimal physical input needed to model the coarse-grained particles in the reservoir. We propose two different approaches and show that in both cases they are sufficient to capture the decisive physical characteristics that allow a valid system-reservoir coupling. The technical satisfactory result paves the way for multiscale analysis of ionic liquids and for truly open boundary molecular simulations.


title: Multigap Superconductivity in a charge density wave superconductor LaPt2Si2

abstract: The superconducting gap structure of a charge density wave (CDW) superconductor LaPt$_2$Si$_2$ ($T_c$ = 1.6~K) having a quasi two dimensional crystal structure has been investigated using muon spin rotation/relaxation ($\mu$SR) measurements carried out in transverse field (TF), zero field (ZF) and longitudinal field (LF) geometries. Rigorous analysis of TF-$\mu$SR spectra in the superconducting state corroborates that the temperature dependence of the effective penetration depth, $\lambda_L$, derived from muon spin depolarization, fits to an isotropic $s+s-$wave model suggesting that the Fermi surface contains two gaps of different magnitude rather than an isotropic gap expected for a conventional $s-$wave superconductor. On the other hand, ZF $\mu$SR data do not show any significant change in muon spin relaxation rate above and below the superconducting transition indicating the fact that time-reversal symmetry is preserved in the system.


title: Independence and totalness of subspaces in phase space methods

abstract: The concepts of independence and totalness of subspaces are introduced in the context of quasi-probability distributions in phase space, for quantum systems with finite-dimensional Hilbert space. It is shown that due to the non-distributivity of the lattice of subspaces, there are various levels of independence, from pairwise independence up to (full) independence. Pairwise totalness, totalness and other intermediate concepts are also introduced, which roughly express that the subspaces overlap strongly among themselves, and they cover the full Hilbert space. A duality between independence and totalness, that involves orthocomplementation (logical NOT operation), is discussed. Another approach to independence is also studied, using Rota's formalism on independent partitions of the Hilbert space. This is used to define informational independence, which is proved to be equivalent to independence. As an application, the pentagram (used in discussions on contextuality) is analyzed using these concepts.


title: Inertial drag on a sphere settling in a stratified fluid

abstract: We compute the drag force on a sphere settling slowly in a quiescent, linearly stratified fluid. Stratification can significantly enhance the drag experienced by the settling particle. The magnitude of this effect depends on whether fluid-density transport around the settling particle is due to diffusion, to advection by the disturbance flow caused by the particle, or due to both. It therefore matters how efficiently the fluid disturbance is convected away from the particle by fluid-inertial terms. When these terms dominate, the Oseen drag force must be recovered. We compute by perturbation theory how the Oseen drag is modified by diffusion and stratification. Our results are in good agreement with recent direct-numerical simulation studies of the problem at small Reynolds numbers and large (but not too large) Froude numbers.


title: WACA: Wearable-Assisted Continuous Authentication

abstract: One-time login process in conventional authentication systems does not guarantee that the identified user is the actual user throughout the session. However, it is necessary to re-verify the user identity periodically throughout a login session without reducing the user convenience. Continuous authentication can address this issue. However, existing methods are either not reliable or not usable. In this paper, we introduce a usable and reliable method called Wearable Assisted Continuous Authentication (WACA). WACA relies on the sensor based keystroke dynamics, where the authentication data is acquired through the built in sensors of a wearable (e.g., smartwatch) while the user is typing. We implemented the WACA framework and evaluated its performance on real devices with real users. The empirical evaluation of WACA reveals that WACA is feasible and its error rate is as low as 1 percent with 30 seconds of processing time and 2 3% for 20 seconds. The computational overhead is minimal. Furthermore, we tested WACA against different attack scenarios. WACA is capable of identifying insider threats with very high accuracy (99.2%) and also robust against powerful adversaries such as imitation and statistical attackers.


title: On the Sublinear Convergence of Randomly Perturbed Alternating Gradient Descent to Second Order Stationary Solutions

abstract: The alternating gradient descent (AGD) is a simple but popular algorithm which has been applied to problems in optimization, machine learning, data ming, and signal processing, etc. The algorithm updates two blocks of variables in an alternating manner, in which a gradient step is taken on one block, while keeping the remaining block fixed. When the objective function is nonconvex, it is well-known the AGD converges to the first-order stationary solution with a global sublinear rate. In this paper, we show that a variant of AGD-type algorithms will not be trapped by "bad" stationary solutions such as saddle points and local maximum points. In particular, we consider a smooth unconstrained optimization problem, and propose a perturbed AGD (PA-GD) which converges (with high probability) to the set of second-order stationary solutions (SS2) with a global sublinear rate. To the best of our knowledge, this is the first alternating type algorithm which takes $\mathcal{O}(\text{polylog}(d)/\epsilon^{7/3})$ iterations to achieve SS2 with high probability [where polylog$(d)$ is polynomial of the logarithm of dimension $d$ of the problem].


title: Convolutional Neural Networks with Alternately Updated Clique

abstract: Improving information flow in deep networks helps to ease the training difficulties and utilize parameters more efficiently. Here we propose a new convolutional neural network architecture with alternately updated clique (CliqueNet). In contrast to prior networks, there are both forward and backward connections between any two layers in the same block. The layers are constructed as a loop and are updated alternately. The CliqueNet has some unique properties. For each layer, it is both the input and output of any other layer in the same block, so that the information flow among layers is maximized. During propagation, the newly updated layers are concatenated to re-update previously updated layer, and parameters are reused for multiple times. This recurrent feedback structure is able to bring higher level visual information back to refine low-level filters and achieve spatial attention. We analyze the features generated at different stages and observe that using refined features leads to a better result. We adopt a multi-scale feature strategy that effectively avoids the progressive growth of parameters. Experiments on image recognition datasets including CIFAR-10, CIFAR-100, SVHN and ImageNet show that our proposed models achieve the state-of-the-art performance with fewer parameters.


title: Limits on Inferring the Past

abstract: Here we define and study the properties of retrodictive inference. We derive equations relating retrodiction entropy and thermodynamic entropy, and as a special case, show that under equilibrium conditions, the two are identical. We demonstrate relations involving the KL-divergence and retrodiction probability, and bound the time rate of change of retrodiction entropy. As a specific case, we invert various Langevin processes, inferring the initial condition of \(N\) particles given their final positions at some later time. We evaluate the retrodiction entropy for Langevin dynamics exactly for special cases, and find that one's ability to infer the initial state of a system can exhibit two possible qualitative behaviors depending on the potential energy landscape, either decreasing indefinitely, or asymptotically approaching a fixed value. We also study how well we can retrodict points that evolve based on the logistic map. We find singular changes in the retrodictivity near bifurcations. Counterintuitively, the transition to chaos is accompanied by maximal retrodictability.


title: Late metal-silicate separation on the IAB parent asteroid: Constraints from combined W and Pt isotopes and thermal modelling

abstract: The short-lived $^{182}$Hf-$^{182}$W decay system is a powerful chronometer for constraining the timing of metal-silicate separation and core formation in planetesimals and planets. Neutron capture effects on W isotopes, however, significantly hamper the application of this tool. In order to correct for neutron capture effects, Pt isotopes have emerged as a reliable in-situ neutron dosimeter. This study applies this method to IAB iron meteorites, in order to constrain the timing of metal segregation on the IAB parent body. The $\epsilon^{182}$W values obtained for the IAB iron meteorites range from -3.61 $\pm$ 0.10 to -2.73 $\pm$ 0.09. Correlating $\epsilon^{\mathrm{i}}$Pt with $^{182}$W data yields a pre-neutron capture $^{182}$W of -2.90 $\pm$ 0.06. This corresponds to a metal-silicate separation age of 6.0 $\pm$ 0.8 Ma after CAI for the IAB parent body, and is interpreted to represent a body-wide melting event. Later, between 10 and 14 Ma after CAI, an impact led to a catastrophic break-up and subsequent reassembly of the parent body. Thermal models of the interior evolution that are consistent with these estimates suggest that the IAB parent body underwent metal-silicate separation as a result of internal heating by short-lived radionuclides and accreted at around 1.4 $\pm$ 0.1 Ma after CAIs with a radius of greater than 60 km.


title: Effect of Mars Atmospheric Loss on Snow Melt Potential in a 3.5-Gyr Mars Climate Evolution Model

abstract: Post-Noachian Martian paleochannels indicate the existence of liquid water on the surface of Mars after about 3.5 Gya (Irwin et al., 2015; Palucis et al., 2016). In order to explore the effects of variations in CO$_{2}$ partial pressure and obliquity on the possibility of surface water, we created a zero-dimensional surface energy balance model. We combine this model with physically consistent orbital histories to track conditions over the last 3.5 Gyr of Martian history. We find that melting is allowed for atmospheric pressures corresponding to exponential loss rates of $dP/dt \propto t^{-3.73}$ or faster, but this rate is within $0.5 \sigma$ of the rate calculated from initial measurements made by the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission, if we assume all the escaping oxygen measured by MAVEN comes from atmospheric CO$_{2}$ (Lillis et al., 2017; Tu et al., 2015). Melting at this loss rate matches selected key geologic constraints on the formation of Hesperian river networks, assuming optimal melt conditions during the warmest part of each Mars year (Irwin et al., 2015; Stopar et al., 2006; Kite et al., 2017a,b). The atmospheric pressure has a larger effect on the surface energy than changes in Mars's mean obliquity. These results show that initial measurements of atmosphere loss by MAVEN are consistent with atmospheric loss being the dominant process that switched Mars from a melt-permitting to a melt-absent climate (Jakosky et al., 2017), but non-CO$_{2}$ warming will be required if $<2$ Gya paleochannels are confirmed, or if most of the escaping oxygen measured by MAVEN comes from H$_{2}$O.


title: $\rho J/\Psi$ Scattering in an Improved Many Body Potential

abstract: We calculate the cross-sections for the processes $\rho J/\psi\rightarrow D^0\bar{D}^0$, $\rho J/\psi\rightarrow D^0\bar{D}^{0*}$ ($D^{0*}\bar{D}^{0}$) and $\rho J/\psi\rightarrow D^{0*}\bar{D}^{0*}$ using a QCD-motivated many-body overlap factor to modify the usual sum of two-body interaction model. The realistic Cornell potential has been used for pairwise interaction in the four quark Hamiltonian and noted to give lesser cross-sections as compared to the quadratic potential. The Resonating group method is employed along with the Born approximation which decouples its integral equations. It is pointed out that the additional QCD effect (a gluonic field ovelap factor) result in a significant suppression in the cross sections as compared to the more popular sum of two-body interaction.


title: Analysis of imperfections in the coherent optical excitation of single atoms to Rydberg states

abstract: We study experimentally various physical limitations and technical imperfections that lead to damping and finite contrast of optically-driven Rabi oscillations between ground and Rydberg states of a single atom. Finite contrast is due to preparation and detection errors and we show how to model and measure them accurately. Part of these errors originates from the finite lifetime of Rydberg states and we observe its $n^3$-scaling with the principal quantum number $n$. To explain the damping of Rabi oscillations, we use simple numerical models, taking into account independently measured experimental imperfections, and show that the observed damping actually results from the accumulation of several small effects, each at the level of a few percents. We discuss prospects for improving the coherence of ground-Rydberg Rabi oscillations in view of applications in quantum simulation and quantum information processing with arrays of single Rydberg atoms.


title: Exponentiable Higher Toposes

abstract: We characterise the class of exponentiable $\infty$-toposes: $\mathcal X$ is exponentiable if and only if $\mathcal S\mathrm{h}(\mathcal X)$ is a continuous $\infty$-category. The heart of the proof is the description of the $\infty$-category of $\mathcal C$-valued sheaves on $\mathcal X$ as an $\infty$-category of functors that satisfy finite limits conditions as well as filtered colimits conditions (instead of limits conditions purely); we call such functors $\omega$-continuous sheaves. As an application, we show that when $\mathcal X$ is exponentiable, its $\infty$-category of stable sheaves $\mathcal S\mathrm{h}(\mathcal X, \mathrm{Sp})$ is a dualisable object in the $\infty$-category of presentable stable $\infty$-categories.


title: Fine-grained wound tissue analysis using deep neural network

abstract: Tissue assessment for chronic wounds is the basis of wound grading and selection of treatment approaches. While several image processing approaches have been proposed for automatic wound tissue analysis, there has been a shortcoming in these approaches for clinical practices. In particular, seemingly, all previous approaches have assumed only 3 tissue types in the chronic wounds, while these wounds commonly exhibit 7 distinct tissue types that presence of each one changes the treatment procedure. In this paper, for the first time, we investigate the classification of 7 wound issue types. We work with wound professionals to build a new database of 7 types of wound tissue. We propose to use pre-trained deep neural networks for feature extraction and classification at the patch-level. We perform experiments to demonstrate that our approach outperforms other state-of-the-art. We will make our database publicly available to facilitate research in wound assessment.


title: A few remarks on invariable generation in infinite groups

abstract: A subset $S$ of a group $G$ invariably generates $G$ if $G$ is generated by $\{ s^g(s) | s\in S\} $ for any choice of $g(s)\in G, s\in S$. In case $G$ is topological one defines similarly the notion of topological invariable generation. A topological group $G$ is said to be $\mathcal{TIG}$ if it is topologically invariably generated by some subset $S\subseteq G$. In this paper we study the problem of (topological) invariable generation for linear groups and for automorphism groups of trees. Our main results show that the Lie group $\mathrm{SL}_{2}(\mathbb{R})$ and the automorphism group of a regular tree are $\mathcal{TIG}$, and that the groups $PSL_m(K) ,m\geq 2$ are not $\mathcal{IG}$ for certain countable fields of infinite transcendence degree over the prime field.


title: Stresses in curved nematic membranes

abstract: Ordering configurations of a director field on a curved membrane induce stress. In this work, we present a theoretical framework to calculate the stress tensor and the torque as a consequence of the nematic ordering; we use the variational principle and invariance of the energy under Euclidean motions. Euler-Lagrange equations of the membrane as well as the corresponding boundary conditions also appear as natural results. The stress tensor found includes attraction-repulsion forces between defects; likewise, defects are attracted to patches with the same sign in gaussian curvature. These forces are mediated by the Green function of Laplace-Beltrami operator of the surface. In addition, we find non-isotropic forces that involve derivatives of the Green function and the gaussian curvature, even in the normal direction to the membrane. We examine the case of axial membranes to analyze the spherical one. For spherical vesicles we find the modified Young-Laplace law as a consequence of the nematic texture. In the case of spherical cap with defect at the north pole, we find that the force is repulsive respect to the north pole, indicating that it is an unstable equilibrium point.


title: Long lifetimes of ultra-hot particles in interacting Fermi systems

abstract: The energy dependence of the relaxation rate of hot electrons due to interaction with the Fermi sea is studied. We consider 2D and 3D systems, quasi-1D quantum wires with multiple transverse bands, as well as single-channel 1D wires. Our analysis includes both spinful and spin-polarized setups, with short-range and Coulomb interactions. We show that, quite generally, the relaxation rate is a non-monotonic function of the electron energy and decays as a power-law at high energies. In other words, ultra-hot electrons regain their coherence with increasing energy. Such a behavior was observed in a recent experiment on multi-band quantum wires, J. Reiner et al, Phys. Rev. X {\bf 7}, 021016 (2017).


title: Error analysis of Nitsche's mortar method

abstract: Optimal a priori and a posteriori error estimates are derived for Nitsche's mortar finite elements. The analysis is based on the equivalence of the Nitsche's method and the stabilised mixed method. The Nitsche's method is defined so that it is robust with respect to large jumps in the material and mesh parameters over the interface. Numerical results demonstrate the robustness of the a posteriori estimators.


title: Capacitively Driven Global Interconnect with Magnetoelectric Switching Based Receiver for Higher Energy Efficiency

abstract: We propose capacitively driven low-swing global interconnect circuit using a receiver that utilizes magnetoelectric (ME) effect induced magnetization switching to reduce the energy consumption. Capacitively driven wire has recently been shown to be effective in improving the performance of global interconnects. Such techniques can reduce the signal swing in the interconnect by using a capacitive divider network and does not require an additional voltage supply. However, the large reduction in signal swing makes it necessary to use differential signaling and amplification for successful regeneration at the receiver, which add area and static power. ME effect induced magnetization reversal has recently been proposed which shows the possibility of using a low voltage to switch a nanomagnet adjacent to a multi-ferroic oxide. Here, we propose an ME effect based receiver that uses the low voltage at the receiving end of the global wire to switch a nanomagnet. The nanomagnet is also used as the free layer of a magnetic tunnel junction (MTJ), the resistance of which is tuned through the ME effect. This change in MTJ resistance is converted to full swing binary signals by using simple digital components. This process allows capacitive low swing interconnection without differential signaling or amplification, which leads to significant energy efficiency. Our simulation results indicate that for 5-10 mm long global wires in IBM 45 nm technology, capacitive ME design consumes 3x lower energy compared to full-swing CMOS design and 2x lower energy compared to differential amplifier based low-swing capacitive CMOS design.


title: Talking about Probability, Inference and Decisions. Part 1: The Witches of Bayes

abstract: In October 2017 the Italian National Institute of Statistics (ISTAT), Italy's body for official statistics, has published the book of fairy tales Le streghe di Bayes (The witches of Bayes) written by ISTAT staff members with the commendable aim of introducing statistical and probabilistic reasoning to children. In this paper the fairy tale which gives the name to the book is analyzed in a dialog between three teachers with different background and expertise. The outcomes are definitively discouraging, especially when the story is compared to the appendix of the book, in which the teaching power of every story is indeed explained (as a matter of fact, without the appendix the fairy tale of the witches seemed to be written with the purpose of make the 'Bayesians', meant as the villagers from 'Bayes', ridiculous). In fact the fairy tale of the witches does not contain any Bayesian reasoning, the suggested decision strategy is simply wrong and the story does not even seem to be easily modifiable (besides the trivial correction of the decision strategy) in order to make it usable as a teaching tool. As it happens in real dialogues, besides the fairy tale in question, the dialogue touches several issues somehow related to the story and concerning probability, inference, prediction and decision making. The present paper is an indirect response to the invitation by the ISBA bulletin to comment on the fairy tale.


title: How long, O Bayesian network, will I sample thee? A program analysis perspective on expected sampling times

abstract: Bayesian networks (BNs) are probabilistic graphical models for describing complex joint probability distributions. The main problem for BNs is inference: Determine the probability of an event given observed evidence. Since exact inference is often infeasible for large BNs, popular approximate inference methods rely on sampling. We study the problem of determining the expected time to obtain a single valid sample from a BN. To this end, we translate the BN together with observations into a probabilistic program. We provide proof rules that yield the exact expected runtime of this program in a fully automated fashion. We implemented our approach and successfully analyzed various real-world BNs taken from the Bayesian network repository.


title: The phase structure of Causal Dynamical Triangulations with toroidal spatial topology

abstract: We investigate the impact of topology on the phase structure of four-dimensional Causal Dynamical Triangulations (CDT). Using numerical Monte Carlo simulations we study CDT with toroidal spatial topology. We confirm existence of all four distinct phases of quantum geometry earlier observed in CDT with spherical spatial topology. We plot the toroidal CDT phase diagram and find that it looks very similar to the case of the spherical spatial topology.


title: Revealing a mode interplay that controls second harmonic radiation in gold nanoantennas

abstract: In this work, we investigate the generation of second harmonic light by gold nanorods and demonstrate that the collected nonlinear intensity depends upon a phase interplay between different modes available in the nanostructure. By recording the backward and forward emitted second harmonic signals from nanorods with various lengths, we find that the maximum nonlinear signal emitted in the forward and backward directions is not obtained for the same nanorod length. We confirm the experimental results with the help of full-wave computations done with a surface integral equation method. These observations are explained by the multipolar nature of the second harmonic emission, which emphasizes the role played by the relative phase between the second harmonic modes. Our findings are of a particular importance for the design of plasmonic nanostructures with controllable nonlinear emission and nonlinear plasmonic sensors as well as for the coherent control of harmonic generations in plasmonic nanostructures.


title: Black hole information paradox without Hawking radiation

abstract: By entangling soft massless particles one can create an arbitrarily large amount of entanglement entropy that carries an arbitrarily small amount of energy. Dropping this entropy into the black hole (b.h.) one can increase the b.h. entropy by an amount that violates Bekenstein bound or any other reasonable bound, leading to a version of b.h. information paradox that does not involve Hawking radiation. Among many proposed solutions of the standard b.h. information paradox with Hawking radiation, only a few can also resolve this version without the Hawking radiation. The assumption that both versions should be resolved in the same way significantly helps to reduce the space of possible resolutions.


title: A Simple Method to improve Initialization Robustness for Active Contours driven by Local Region Fitting Energy

abstract: Active contour models based on local region fitting energy can segment images with intensity inhomogeneity effectively, but their segmentation results are easy to error if the initial contour is inappropriate. In this paper, we present a simple and universal method of improving the robustness of initial contour for these local fitting-based models. The core idea of proposed method is exchanging the fitting values on the two sides of contour, so that the fitting values inside the contour are always larger (or smaller) than the values outside the contour in the process of curve evolution. In this way, the whole curve will evolve along the inner (or outer) boundaries of object, and less likely to be stuck in the object or background. Experimental results have proved that using the proposed method can enhance the robustness of initial contour and meanwhile keep the original advantages in the local fitting-based models.


title: Sensor and Sink Placement, Scheduling and Routing Algorithms for Connected Coverage of Wireless Sensor Networks

abstract: A sensor is a small electronic device which has the ability to sense, compute and communicate either with other sensors or directly with a base station (sink). In a wireless sensor network (WSN), the sensors monitor a region and transmit the collected data packets through routes to the sinks. In this study, we propose a mixed--integer linear programming (MILP) model to maximize the number of time periods that a WSN carries out the desired tasks with limited energy and budget. Our sink and sensor placement, scheduling, routing with connected coverage ($SPSRC$) model is the first in the literature that combines the decisions for the locations of sinks and sensors, activity schedules of the deployed sensors, and data flow routes from each active sensor to its assigned sink for connected coverage of the network over a finite planning horizon. The problem is NP--hard and difficult to solve even for small instances. Assuming that the sink locations are known, we develop heuristics which construct a feasible solution of the problem by gradually satisfying the constraints. Then, we introduce search heuristics to determine the locations of the sinks to maximize the network lifetime. Computational experiments reveal that our heuristic methods can find near optimal solutions in an acceptable amount of time compared to the commercial solver CPLEX 12.7.0.


title: On the norm attainment set of a bounded linear operator and semi-inner-products in normed spaces

abstract: We obtain a complete characterization of the norm attainment set of a bounded linear functional on a normed space, in terms of a semi-inner-product defined on the space. Motivated by this result, we further apply the concept of semi-inner-products to obtain a complete characterization of the norm attainment set of a bounded linear operator between any two real normed spaces. In particular, this answers an open question raised recently in [Sain, D., \textit{On the norm attainment set of a bounded linear operator}, J. Math. Anal. Appl., \textbf{457} (2018), 67-76.]. Our results illustrate the applicability of semi-inner-products towards a better understanding of the geometry of normed spaces.


title: Precision medicine as a control problem: Using simulation and deep reinforcement learning to discover adaptive, personalized multi-cytokine therapy for sepsis

abstract: Sepsis is a life-threatening condition affecting one million people per year in the US in which dysregulation of the body's own immune system causes damage to its tissues, resulting in a 28 - 50% mortality rate. Clinical trials for sepsis treatment over the last 20 years have failed to produce a single currently FDA approved drug treatment. In this study, we attempt to discover an effective cytokine mediation treatment strategy for sepsis using a previously developed agent-based model that simulates the innate immune response to infection: the Innate Immune Response agent-based model (IIRABM). Previous attempts at reducing mortality with multi-cytokine mediation using the IIRABM have failed to reduce mortality across all patient parameterizations and motivated us to investigate whether adaptive, personalized multi-cytokine mediation can control the trajectory of sepsis and lower patient mortality. We used the IIRABM to compute a treatment policy in which systemic patient measurements are used in a feedback loop to inform future treatment. Using deep reinforcement learning, we identified a policy that achieves 0% mortality on the patient parameterization on which it was trained. More importantly, this policy also achieves 0.8% mortality over 500 randomly selected patient parameterizations with baseline mortalities ranging from 1 - 99% (with an average of 49%) spanning the entire clinically plausible parameter space of the IIRABM. These results suggest that adaptive, personalized multi-cytokine mediation therapy could be a promising approach for treating sepsis. We hope that this work motivates researchers to consider such an approach as part of future clinical trials. To the best of our knowledge, this work is the first to consider adaptive, personalized multi-cytokine mediation therapy for sepsis, and is the first to exploit deep reinforcement learning on a biological simulation.


title: A Primer on Resurgent Transseries and Their Asymptotics

abstract: The computation of observables in general interacting theories, be them quantum mechanical, field, gauge or string theories, is a non-trivial problem which in many cases can only be addressed by resorting to perturbative methods. In most physically interesting problems these perturbative expansions result in asymptotic series with zero radius of convergence. These asymptotic series then require the use of resurgence and transseries in order for the associated observables to become nonperturbatively well-defined. Resurgence encodes the complete large-order asymptotic behaviour of the coefficients from a perturbative expansion, generically in terms of (multi) instanton sectors and for each problem in terms of its Stokes constants. Some observables arise from linear problems, and have a finite number of instanton sectors and associated Stokes constants; some other observables arise from nonlinear problems, and have an infinite number of instanton sectors and Stokes constants. By means of two very explicit examples, and with emphasis on a pedagogical style of presentation, this work aims at serving as a primer on the aforementioned resurgent, large-order asymptotics of general perturbative expansions. This includes discussions of transseries, Stokes phenomena, generalized steepest-descent methods, Borel transforms, nonlinear resonance, and alien calculus. Furthermore, resurgent properties of transseries---usually described mathematically via alien calculus---are recast in equivalent physical languages: either a "statistical mechanical" language, as motions in chains and lattices; or a "conformal field theoretical" language, with underlying Virasoro-like algebraic structures.

categories: hep-th hep-lat hep-ph math-ph math.MP quant-ph


title: Three-particle $N\pi\pi$ state contribution to the nucleon two-point function in lattice QCD

abstract: The three-particle $N\pi\pi$-state contribution to the QCD two-point function of standard nucleon interpolating fields is computed to leading order in chiral perturbation theory. Using the experimental values for two low-energy coefficients the impact of this contribution on lattice QCD calculations of the nucleon mass is estimated. The impact is found to be at the per mille level at most and negligible in practice.


title: Modeling our survival in a zombie apocalypse

abstract: In this work we applied several concepts on the modeling of complex systems in an attempt to save mankind in the case of a zombie apocalypse. We developed a dynamical system in order to model a zombie outbreak in which we, regular humans, are aided by military personnel in the fight against the zombies. Our analysis has shown that the initial amount of military personnel play a key rule on our survival, even when the zombies are extremely aggressive and in large advantage.This document is a pre-print version of the paper de Mendon\c{c}a, J.P.A., Teixeira, L.M.V., Sato, F. et al. Math Intelligencer (2019). https://doi.org/10.1007/s00283-019-09893-9.


title: On the Low-Complexity, Hardware-Friendly Tridiagonal Matrix Inversion for Correlated Massive MIMO Systems

abstract: In massive MIMO (M-MIMO) systems, one of the key challenges in the implementation is the large-scale matrix inversion operation, as widely used in channel estimation, equalization, detection, and decoding procedures. Traditionally, to handle this complexity issue, several low-complexity matrix inversion approximation methods have been proposed, including the classic Cholesky decomposition and the Neumann series expansion (NSE). However, the conventional approaches failed to exploit neither the special structure of channel matrices nor the critical issues in the hardware implementation, which results in poorer throughput performance and longer processing delay. In this paper, by targeting at the correlated M-MIMO systems, we propose a modified NSE based on tridiagonal matrix inversion approximation (TMA) to accommodate the complexity as well as the performance issue in the conventional hardware implementation, and analyze the corresponding approximation errors. Meanwhile, we investigate the VLSI implementation for the proposed detection algorithm based on a Xilinx Virtex-7 XC7VX690T FPGA platform. It is shown that for correlated massive MIMO systems, it can achieve near-MMSE performance and $630$ Mb/s throughput. Compared with other benchmark systems, the proposed pipelined TMA detector can get high throughput-to-hardware ratio. Finally, we also propose a fast iteration structure for further research.


title: Changes in the metallicity of gas giant planets due to pebble accretion

abstract: We run numerical simulations to study the accretion of gas and dust grains onto gas giant planets embedded into massive protoplanetary discs. The outcome is found to depend on the disc cooling rate, planet mass, grain size and irradiative feedback from the planet. If radiative cooling is efficient, planets accrete both gas and pebbles rapidly, open a gap and usually become massive brown dwarfs. In the inefficient cooling case, gas is too hot to accrete onto the planet but pebble accretion continues and the planets migrate inward rapidly. Radiative feedback from the planet tends to suppress gas accretion. Our simulations predict that metal enrichment of planets by dust grain accretion inversely correlates with the final planet mass, in accordance with the observed trend in the inferred bulk composition of Solar System and exosolar giant planets. To account for observations, however, as much as ~30-50% of the dust mass should be in the form of large grains.


title: Identifying Sources and Sinks in the Presence of Multiple Agents with Gaussian Process Vector Calculus

abstract: In systems of multiple agents, identifying the cause of observed agent dynamics is challenging. Often, these agents operate in diverse, non-stationary environments, where models rely on hand-crafted environment-specific features to infer influential regions in the system's surroundings. To overcome the limitations of these inflexible models, we present GP-LAPLACE, a technique for locating sources and sinks from trajectories in time-varying fields. Using Gaussian processes, we jointly infer a spatio-temporal vector field, as well as canonical vector calculus operations on that field. Notably, we do this from only agent trajectories without requiring knowledge of the environment, and also obtain a metric for denoting the significance of inferred causal features in the environment by exploiting our probabilistic method. To evaluate our approach, we apply it to both synthetic and real-world GPS data, demonstrating the applicability of our technique in the presence of multiple agents, as well as its superiority over existing methods.


title: Oriented Convex Containers of Polygons

abstract: We consider the optimal containment of polygonal regions within convex containers with the special property of 'orientedness' - an oriented region enables us to choose a preferred direction on the plane (this direction is not necessarily an axis of symmetry) - and derive preliminary results.


title: Quantum cognition goes beyond-quantum: modeling the collective participant in psychological measurements

abstract: In psychological measurements, two levels should be distinguished: the 'individual level', relative to the different participants in a given cognitive situation, and the 'collective level', relative to the overall statistics of their outcomes, which we propose to associate with a notion of 'collective participant'. When the distinction between these two levels is properly formalized, it reveals why the modeling of the collective participant generally requires beyond-quantum - non-Bornian - probabilistic models, when sequential measurements at the individual level are considered, and this though a pure quantum description remains valid for single measurement situations.


title: Classical electrodynamics and gauge symmetry of the X-boson

abstract: In this paper we have obtained several new results concerning the X-boson, which is being considered recently as one of the main candidate of the dark matter particle content. The classical electrodynamics for the X-boson model was explored to understand its propagation in space-time. The field equations (Maxwell-type) and the corresponding wave equations were obtained. They indicate the dispersion relations of both a massive and massless particles. These results can be interpreted as a photon and the X-boson (with a mass of 17 MeV), respectively. The gauge symmetry and gauge transformations were discussed. A full model diagonalization was introduced to obtain a Maxwell sector added up to a Proca sector. After that, we have obtained the retarded Green functions which yield the inhomogeneous solutions of the wave equations for the X-boson fields. We have applied this solution to the charge point motion, and we transformed the potential solution to a quadrature problem. In particular, as an example, we have analyzed the point charge in a straight line uniform motion, where the corresponding electromagnetic charge were calculated.


title: Thermodynamics of Ricci-Gauss-Bonnet Dark Energy

abstract: We investigate the validity of generalized second law of thermodynamics of a physical system comprising of newly proposed dark energy model called Ricci Gauss-Bonnet and cold dark matter enveloped by apparent horizon and event horizon in flat Friedmann-Robertson-Walker (FRW) universe. For this purpose, Bekenstein entropy, Renyi, logarithmic and power law entropic corrections are used. It is found that this law exhibits the validity on both apparent and event horizons except for the case of logarithmic entropic correction at apparent horizon. Also, we check the thermodynamical equilibrium condition for all cases of entropy and found its vitality in all cases of entropy.


title: Reconstruction of brain networks involved in magnetophosphene perception using dense electroencephalography

abstract: Characterizing functional brain networks in humans during magnetophosphene perception. Dense electroencephalography (EEG, 128 channels) was performed in N=3 volunteers during high-level (50 mT) magnetic field (MF) exposure. Functional brain networks were reconstructed, at the cortical level from scalp recordings, using the EEG source connectivity method. Magnetophosphene perception appears to consistently activate the right inferior occipito-temporal pathway. This study provides the very first neuroimaging results characterizing magnetophosphene perception in humans. The use of dense-EEG source connectivity is a promising approach in the field of bioelectromagnetics.


title: Influences of the coordinate dependent noncommutative space on charged and spin currents

abstract: We study the charged and spin currents on a coordinate dependent noncommutative space. Starting from the noncommutative extended relativistic equation of motion, the non-relativistic approximation is obtained by using the Foldy-Wouthuysen transformation, and then the charged and spin currents are derived by using the extended Drude model. We find that the charged current is twisted by modifying the off-diagonal elements of the Hall conductivity, however, the spin current is not affected up to leading order of the noncommutative parameter.


title: Symmetric indefinite triangular factorization revealing the rank profile matrix

abstract: We present a novel recursive algorithm for reducing a symmetric matrix to a triangular factorization which reveals the rank profile matrix. That is, the algorithm computes a factorization $\mathbf{P}^T\mathbf{A}\mathbf{P} = \mathbf{L}\mathbf{D}\mathbf{L}^T$ where $\mathbf{P}$ is a permutation matrix, $\mathbf{L}$ is lower triangular with a unit diagonal and $\mathbf{D}$ is symmetric block diagonal with $1{\times}1$ and $2{\times}2$ antidiagonal blocks. The novel algorithm requires $O(n^2r^{\omega-2})$ arithmetic operations. Furthermore, experimental results demonstrate that our algorithm can even be slightly more than twice as fast as the state of the art unsymmetric Gaussian elimination in most cases, that is it achieves approximately the same computational speed. By adapting the pivoting strategy developed in the unsymmetric case, we show how to recover the rank profile matrix from the permutation matrix and the support of the block-diagonal matrix. There is an obstruction in characteristic $2$ for revealing the rank profile matrix which requires to relax the shape of the block diagonal by allowing the 2-dimensional blocks to have a non-zero bottom-right coefficient. This relaxed decomposition can then be transformed into a standard $\mathbf{P}\mathbf{L}\mathbf{D}\mathbf{L}^T\mathbf{P}^T$ decomposition at a negligible cost.


title: Dissipative operators and operator Lipschitz functions

abstract: The purpose of this paper is to obtain an integral representation for the difference $f(L_1)-f(L_2)$ of functions of maximal dissipative operators. This representation in terms of double operator integrals will allow us to establish Lipschitz type estimates for functions of maximal dissipative operators. We also consider a similar problem for quasicommutators, i.e., operators of the form $f(L_1)R-Rf(L_2)$.


title: Parallel Transport of Higher Flat Gerbes as an Extended Homotopy Quantum Field Theory

abstract: We prove that the parallel transport of a flat $n-1$-gerbe on any given target space gives rise to an $n$-dimensional extended homotopy quantum field theory. In case the target space is the classifying space of a finite group, we provide explicit formulae for this homotopy quantum field theory in terms of transgression. Moreover, we use the geometric theory of orbifolds to give a dimension-independent version of twisted and equivariant Dijkgraaf-Witten models. Finally, we introduce twisted equivariant Dijkgraaf-Witten theories giving us in the 3-2-1-dimensional case a new class of equivariant modular tensor categories which can be understood as twisted versions of the equivariant modular categories constructed by Maier, Nikolaus and Schweigert.


title: Room-temperature lasing action in GaN quantum wells in the infrared 1.5 micron region

abstract: Large-scale optoelectronics integration is strongly limited by the lack of efficient light sources, which could be integrated with the silicon complementary metal-oxide-semiconductor (CMOS) technology. Persistent efforts continue to achieve efficient light emission from silicon in the extending the silicon technology into fully integrated optoelectronic circuits. Here, we report the realization of room-temperature stimulated emission in the technologically crucial 1.5 micron wavelength range from Er-doped GaN multiple-quantum wells on silicon and sapphire. Employing the well-acknowledged variable stripe technique, we have demonstrated an optical gain up to 170 cm-1 in the multiple-quantum well structures. The observation of the stimulated emission is accompanied by the characteristic threshold behavior of emission intensity as a function of pump fluence, spectral linewidth narrowing and excitation length. The demonstration of room-temperature lasing at the minimum loss window of optical fibers and in the eye-safe wavelength region of 1.5 micron are highly sought-after for use in many applications including defense, industrial processing, communication, medicine, spectroscopy and imaging. As the synthesis of Er-doped GaN epitaxial layers on silicon and sapphire has been successfully demonstrated, the results laid the foundation for achieving hybrid GaN-Si lasers providing a new pathway towards full photonic integration for silicon optoelectronics.


title: The density of expected persistence diagrams and its kernel based estimation

abstract: Persistence diagrams play a fundamental role in Topological Data Analysis where they are used as topological descriptors of filtrations built on top of data. They consist in discrete multisets of points in the plane $\mathbb{R}^2$ that can equivalently be seen as discrete measures in $\mathbb{R}^2$. When the data come as a random point cloud, these discrete measures become random measures whose expectation is studied in this paper. First, we show that for a wide class of filtrations, including the \v{C}ech and Rips-Vietoris filtrations, the expected persistence diagram, that is a deterministic measure on $\mathbb{R}^2$ , has a density with respect to the Lebesgue measure. Second, building on the previous result we show that the persistence surface recently introduced in [Adams & al., Persistence images: a stable vector representation of persistent homology] can be seen as a kernel estimator of this density. We propose a cross-validation scheme for selecting an optimal bandwidth, which is proven to be a consistent procedure to estimate the density.


title: A High GOPs/Slice Time Series Classifier for Portable and Embedded Biomedical Applications

abstract: Nowadays a diverse range of physiological data can be captured continuously for various applications in particular wellbeing and healthcare. Such data require efficient methods for classification and analysis. Deep learning algorithms have shown remarkable potential regarding such analyses, however, the use of these algorithms on low-power wearable devices is challenged by resource constraints such as area and power consumption. Most of the available on-chip deep learning processors contain complex and dense hardware architectures in order to achieve the highest possible throughput. Such a trend in hardware design may not be efficient in applications where on-node computation is required and the focus is more on the area and power efficiency as in the case of portable and embedded biomedical devices. This paper presents an efficient time-series classifier capable of automatically detecting effective features and classifying the input signals in real-time. In the proposed classifier, throughput is traded off with hardware complexity and cost using resource sharing techniques. A Convolutional Neural Network (CNN) is employed to extract input features and then a Long-Short-Term-Memory (LSTM) architecture with ternary weight precision classifies the input signals according to the extracted features. Hardware implementation on a Xilinx FPGA confirm that the proposed hardware can accurately classify multiple complex biomedical time series data with low area and power consumption and outperform all previously presented state-of-the-art records. Most notably, our classifier reaches 1.3$\times$ higher GOPs/Slice than similar state of the art FPGA-based accelerators.


title: Phase Transition for the Contact Process in a Random Environment on Zd*Z+

abstract: We review the results in Chen & Yao(2009,2012) which concern the contact process in a static random environment on the half space Z^d*Z^+ and make some addition to them. Furthermore, we explain why our methods cannot apply to the whole space case and compare our results with some related works.


title: Characteristics of type III radio bursts and solar S bursts

abstract: The Sun is an active source of radio emission which is often associated with the acceleration of electrons arising from processes such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), numerous solar S bursts (where S stands for short) and storms of Type III radio bursts have been observed, that are not directly relates to flares and CMEs. Here, we expand our understanding on the spectral characteristic of these two different types of radio bursts based on observations from the Low Frequency Array (LOFAR). On 9 July 2013, over 3000 solar S bursts accompanied by over 800 Type III radio bursts were observed over a time period of ~8 hours. The characteristics of Type III radio bursts are consistent to previous studies, while S bursts show narrow bandwidths, durations and drift rates of about 1/2 the drift rate of Type III bursts. Type III bursts and solar S bursts occur in a region in the corona where plasma emission is the dominant emission mechanism as determined by data constrained density and magnetic field models.


title: Time Varying Channel Tracking with Spatial and Temporal BEM for Massive MIMO Systems

abstract: In this paper, we propose a channel tracking method for massive multi-input and multi-output systems under both time-varying and spatial-varying circumstance. Exploiting the characteristics of massive antenna array, a spatial-temporal basis expansion model is designed to reduce the effective dimensions of up-link and down-link channel, which decomposes channel state information into the time-varying spatial information and gain information. We firstly model the users movements as a one-order unknown Markov process, which is blindly learned by the expectation and maximization (EM) approach. Then, the up-link time varying spatial information can be blindly tracked by Taylor series expansion of the steering vector, while the rest up-link channel gain information can be trained by only a few pilot symbols. Due to angle reciprocity (spatial reciprocity), the spatial information of the down-link channel can be immediately obtained from the up-link counterpart, which greatly reduces the complexity of down-link channel tracking. Various numerical results are provided to demonstrate the effectiveness of the proposed method.


title: Inflation with explicit parametric connection between General Relativity and Scalar-Tensor gravity

abstract: In this paper we consider the cosmological inflation within scalar-tensor gravity and compare it with standard inflation based on general relativity. The difference is determined by the value of the parameter $\Delta$. This approach is associated with using the special ansatz which links a function that defines a type of gravity with a scale factor of the universe.


title: DiGrad: Multi-Task Reinforcement Learning with Shared Actions

abstract: Most reinforcement learning algorithms are inefficient for learning multiple tasks in complex robotic systems, where different tasks share a set of actions. In such environments a compound policy may be learnt with shared neural network parameters, which performs multiple tasks concurrently. However such compound policy may get biased towards a task or the gradients from different tasks negate each other, making the learning unstable and sometimes less data efficient. In this paper, we propose a new approach for simultaneous training of multiple tasks sharing a set of common actions in continuous action spaces, which we call as DiGrad (Differential Policy Gradient). The proposed framework is based on differential policy gradients and can accommodate multi-task learning in a single actor-critic network. We also propose a simple heuristic in the differential policy gradient update to further improve the learning. The proposed architecture was tested on 8 link planar manipulator and 27 degrees of freedom(DoF) Humanoid for learning multi-goal reachability tasks for 3 and 2 end effectors respectively. We show that our approach supports efficient multi-task learning in complex robotic systems, outperforming related methods in continuous action spaces.


title: Nonlinear Schr\"odinger equation, differentiation by parts and modulation spaces

abstract: We show the existence of weak solutions in the extended sense of the Cauchy problem for the cubic nonlinear Schr\"odinger equation in the modulation space $M_{p,q}^{s}(\mathbb R)$ where $1\leq q\leq2$, $2\leq p<\frac{10q'}{q'+6}$ and $s\geq0$. Moreover, for either $1\leq q\leq\frac32, s\geq0$ and $2\leq p\leq 3$ or $\frac32<q\leq\frac{18}{11}, s>\frac23-\frac1{q}$ and $2\leq p\leq 3$ or $\frac{18}{11}<q\leq2, s>\frac23-\frac1{q}$ and $2\leq p<\frac{10q'}{q'+6}$ we show that the Cauchy problem is unconditionally wellposed in $M_{p,q}^{s}(\mathbb R).$ This improves \cite{NP}, where the case $p=2$ was considered and the differentiation by parts technique was introduced to a problem with continuous Fourier variable. Here the same technique is used, but more delicate estimates are necessary for $p\neq2$.


title: Leakage and Protocol Composition in a Game-Theoretic Perspective

abstract: In the inference attacks studied in Quantitative Information Flow (QIF), the adversary typically tries to interfere with the system in the attempt to increase its leakage of secret information. The defender, on the other hand, typically tries to decrease leakage by introducing some controlled noise. This noise introduction can be modeled as a type of protocol composition, i.e., a probabilistic choice among different protocols, and its effect on the amount of leakage depends heavily on whether or not this choice is visible to the adversary. In this work we consider operators for modeling visible and invisible choice in protocol composition, and we study their algebraic properties. We then formalize the interplay between defender and adversary in a game-theoretic framework adapted to the specific issues of QIF, where the payoff is information leakage. We consider various kinds of leakage games, depending on whether players act simultaneously or sequentially, and on whether or not the choices of the defender are visible to the adversary. Finally, we establish a hierarchy of these games in terms of their information leakage, and provide methods for finding optimal strategies (at the points of equilibrium) for both attacker and defender in the various cases. The full version of this paper can be found in arXiv:1803.10042


title: Generalized boosts with shell structure of the parameter space

abstract: A modification of boost transformation in arbitrary pseudo-Euclidean space is suggested, which in the case of the Minkowski space admits the existence of inertial reference frames moving with velocities taking values in a certain bounded interval. The velocity space may be partitioned by hypersurfaces \beta^2=\beta^2_k=const, k=1,2,3,..., into a finite or countable number of domains (shells), each of which has own class of inertial "reference frames" and the velocity composition law. These shells are in one-to-one correspondence. A set of mappings of shells to each other forms the group, isomorphic to permutation group in the case of finite number of shells, or the group of integers in the case of countable number of shells in the velocity space


title: Quantitative Separation Logic - A Logic for Reasoning about Probabilistic Programs

abstract: We present quantitative separation logic ($\mathsf{QSL}$). In contrast to classical separation logic, $\mathsf{QSL}$ employs quantities which evaluate to real numbers instead of predicates which evaluate to Boolean values. The connectives of classical separation logic, separating conjunction and separating implication, are lifted from predicates to quantities. This extension is conservative: Both connectives are backward compatible to their classical analogs and obey the same laws, e.g. modus ponens, adjointness, etc. Furthermore, we develop a weakest precondition calculus for quantitative reasoning about probabilistic pointer programs in $\mathsf{QSL}$. This calculus is a conservative extension of both Reynolds' separation logic for heap-manipulating programs and Kozen's / McIver and Morgan's weakest preexpectations for probabilistic programs. Soundness is proven with respect to an operational semantics based on Markov decision processes. Our calculus preserves O'Hearn's frame rule, which enables local reasoning. We demonstrate that our calculus enables reasoning about quantities such as the probability of terminating with an empty heap, the probability of reaching a certain array permutation, or the expected length of a list.


title: Eigenstate thermalization hypothesis, time operator, and extremely quick relaxation of fidelity

abstract: The eigenstate thermalization hypothesis (ETH) insists that for nonintegrable systems each energy eigenstate accurately gives microcanonical expectation values for a class of observables. As a mechanism for ETH to hold, we show that the energy eigenstates are superposition of uncountably many quasi eigenstates of operationally defined "time operator", which are thermal for thermodynamic isolated quantum many-body systems and approximately orthogonal in terms of extremely short relaxation time of the fidelity. In this way, our scenario provides a theoretical explanation of ETH.


title: Turing-Hopf bifurcation and spatiotemporal patterns in a ratio-dependent diffusive Holling-Tanner system with time delay

abstract: The Turing-Hopf type spatiotemporal patterns in a diffusive Holling-Tanner model with discrete time delay is considered. A global Turing bifurcation theorem for $\tau=0$ and a local Turing bifurcation theorem for $\tau>0$ are given by the method of eigenvalue analysis and prior estimation. Further considering the degenerated situation, the existence of Bogdanov-Takens bifurcation and Turing-Hopf bifurcation are obtained. The normal form method is used to study the explicit dynamics near the Turing-Hopf singularity, and we establish the existence of various self-organized spatiotemporal patterns, such as two non-constant steady states (stripe patterns) coexist and two spatially inhomogeneous periodic solutions (spot patterns) coexist. Moreover, the Turing-Turing-Hopf type spatiotemporal patterns, that is a subharmonic phenomenon with two spatial wave numbers and one temporal frequency, are also found and theoretically explained, when there is another Turing bifurcation curve which is relatively closed to the studied one.


title: QCH Kahler surfaces II

abstract: In this paper we give new examples of QCH Kahler surfaces whose opposite almost Hermitian strucure is Hermitian and not locally conformally Kahler. In this way we give also a large class of examples of Hermitian surfaces with J-invariant Ricci tensor which are not l.c.k.


title: Approximate subloops and Freiman's theorem in finitely generated commutative moufang loops

abstract: Fix a parameter $K\geqslant 1$. A $K$-approximate subgroup is a finite set $A$ in a group $G$ which contains the identity, is symmetric and such that $A.A$ can covered by $K$ left translates of $A$. This article deals with the generalisation of the concept of approximate groups in the case of loops which we call approximate loops and the description of $K$-approximate subloops when the ambient loop is a finitely generated commutative moufang loop. Specifically we have a Freiman type theorem where such approximate subloops are controlled by arithmetic progressions defined in the commutative moufang loops.


title: Efficient V2V Communication Scheme for 5G MmWave Hyper-Connected CAVs

abstract: Connected and Autonomous Vehicles (CAVs) require continuous access to sensory data to perform complex high-speed maneuvers and advanced trajectory planning. High priority CAVs are particularly reliant on extended perception horizon facilitated by sensory data exchange between CAVs. Existing technologies such as the Dedicated Short Range Communications (DSRC) are ill-equipped to provide advanced cooperative perception service. This creates the need for more sophisticated technologies such as the 5G Millimetre-Waves (mmWaves). In this work, we propose a distributed Vehicle-to-Vehicle (V2V) mmWaves association scheme operating in a heterogeneous manner. Our system utilises the information exchanged within the DSRC frequency band to bootstrap the best CAV pairs formation. Using a Stable Fixtures Matching Game, we form V2V multipoint-to-multipoint links. Compared to more traditional point-to-point links, our system provides almost twice as much sensory data exchange capacity for high priority CAVs while doubling the mmWaves channel utilisation for all the vehicles in the network.


title: Low-Overhead Coordination in Sub-28 Millimeter-Wave Networks

abstract: In this paper, we present some contributions from our recent investigation. We address the open issue of interference coordination for sub-28 GHz millimeter-wave communication, by proposing fast-converging coordination algorithms, for dense multi-user multi-cell networks. We propose to optimize a lower bound on the network sum-rate, after investigating its tightness. The bound in question results in distributed optimization, requiring local information at each base station and user. We derive the optimal solution to the transmit and receive filter updates, that we dub non-homogeneous waterfilling, and show its convergence to a stationary point of the bound. We also underline a built-in mechanism to turn-off data streams with low-SINR, and allocate power to high-SNR streams. This "stream control" is a at the root of the fast-converging nature of the algorithm. Our numerical result conclude that low-overhead coordination offers large gains, for dense sub-$28$ GHz systems. These findings bear direct relevance to the ongoing discussions around 5G New Radio.


title: On the Benefits of Asymmetric Coded Cache Placement in Combination Networks with End-User Caches

abstract: This paper investigates the fundamental tradeoff between cache size and download time in the (H;r;M;N) combination network, where a server with N files is connected to H relays (without caches) and each of the K:=\binom{H}{r} users (with caches of size M files) is connected to a different subset of r relays. Existing schemes fall within two categories: either use the uncoded symmetric cache placement originally proposed for the shared-link model and design delivery phase dependent on the network topology, or effectively divide the combination network into H independent shared-link networks each serving \binom{H-1}{r-1} users; in either case, the placement phase is independent of network topology. In this paper, a novel strategy is proposed where the coded cache placement is dependent on network topology. The proposed scheme is shown to be information theoretically optimal for large cache size. In addition, when not exactly optimal, the proposed scheme can also outperform existing schemes.


title: Time-diffracting beams: On their nature, diffraction-free propagation as needles of light, and nonlinear generation

abstract: We investigate on the properties of the recently introduced time-diffracting (TD) beams in free space. They are shown to be paraxial and quasi-monochromatic realizations of localized waves, spatiotemporal localized waves travelling undistorted at arbitrary speeds. The paraxial and quasi-monochromatic regime is shown to be necessary to observe what can properly be named diffraction in time. TD beams of finite energy travelling at quasi-luminal velocities are seen to form substantially longer foci or needles of light than the so-called abruptly focusing and defocusing needle of light, or limiting TD beam of infinite speed. Exploring the properties of TD beams under Lorentz transformations and transformation by paraxial optical systems, we realize that the relativistically moving nonlinear polarization of material media induced by a strongly localized fundamental pump wave generates a TD beam at its second harmonic, whose diffraction-free behavior as a needle of light in free space can be optimized with a standard $4f$-imager system.


title: Coexistence and duality in competing species models

abstract: We discuss some stochastic spatial generalizations of the Lotka--Volterra model for competing species. The generalizations take the forms of spin systems on general discrete sets and interacting diffusions on integer lattices. Methods for proving coexistence in these generalizations and some related open questions are discussed. We use duality as the central point of view. It relates coexistence of the models to survival of their dual processes.


title: Homogenization and Scattering Analysis of Second-Harmonic Generation in Nonlinear Metasurfaces

abstract: We propose an extensive discussion on the homogenization and scattering analysis of second-order nonlinear metasurfaces. Our developments are based on the generalized sheet transition conditions (GSTCs) which are used to model the electromagnetic responses of nonlinear metasurfaces. The GSTCs are solved both in the frequency domain, assuming an undepleted pump regime, and in the time-domain, assuming dispersionless material properties but a possible depleted pump regime. Based on these two modeling approaches, we derive the general second-harmonic reflectionless and transmissionless conditions as well as the conditions of asymmetric reflection and transmission. We also discuss and clarify the concept of nonreciprocal scattering pertaining to nonlinear metasurfaces.


title: HSI-CNN: A Novel Convolution Neural Network for Hyperspectral Image

abstract: With the development of deep learning, the performance of hyperspectral image (HSI) classification has been greatly improved in recent years. The shortage of training samples has become a bottleneck for further improvement of performance. In this paper, we propose a novel convolutional neural network framework for the characteristics of hyperspectral image data, called HSI-CNN. Firstly, the spectral-spatial feature is extracted from a target pixel and its neighbors. Then, a number of one-dimensional feature maps, obtained by convolution operation on spectral-spatial features, are stacked into a two-dimensional matrix. Finally, the two-dimensional matrix considered as an image is fed into standard CNN. This is why we call it HSI-CNN. In addition, we also implements two depth network classification models, called HSI-CNN+XGBoost and HSI-CapsNet, in order to compare the performance of our framework. Experiments show that the performance of hyperspectral image classification is improved efficiently with HSI-CNN framework. We evaluate the model's performance using four popular HSI datasets, which are the Kennedy Space Center (KSC), Indian Pines (IP), Pavia University scene (PU) and Salinas scene (SA). As far as we concerned, HSI-CNN has got the state-of-art accuracy among all methods we have known on these datasets of 99.28%, 99.09%, 99.42%, 98.95% separately.


title: Caching in Combination Networks: A Novel Delivery by Leveraging the Network Topology

abstract: Maddah-Ali and Niesen (MAN) in 2014 surprisingly showed that it is possible to serve an arbitrarily large number of cache-equipped users with a constant number of transmissions by using coded caching in shared-link broadcast networks. This paper studies the tradeoff between the user's cache size and the file download time for combination networks, where users with caches communicate with the servers through intermediate relays. Motivated by the so-called separation approach, it is assumed that placement and multicast message generation are done according to the MAN original scheme and regardless of the network topology. The main contribution of this paper is the design of a novel two-phase delivery scheme that, accounting to the network topology, outperforms schemes available in the literature. The key idea is to create additional (compared to MAN) multicasting opportunities: in the first phase coded messages are sent with the goal of increasing the amount of `side information' at the users, which is then leveraged during the second phase. The download time with the novel scheme is shown to be proportional to 1=H (with H being the number or relays) and to be order optimal under the constraint of uncoded placement for some parameter regimes.


title: A two-stage ensemble Kalman filter based on multiscale model reduction for inverse problems in time fractional diffusion-wave equations

abstract: Ensemble Kalman filter (EnKF) has been widely used in state estimation and parameter estimation for the dynamic system where observational data is obtained sequentially in time. To reduce uncertainty and accelerate posterior inference, a two-stage ensemble Kalman filter is presented to improve the sequential analysis of EnKF. It is known that the final posterior ensemble may be concentrated in a small portion of the entire support of the initial prior ensemble. It will be much more efficient if we first build a new prior by some partial observations, and construct a surrogate only over the significant region of the new prior. To this end, we construct a very coarse model using generalized multiscale finite element method (GMsFEM) and generate a new prior ensemble in the first stage. GMsFEM provides a set of hierarchical multiscale basis functions supported in coarse blocks. This gives flexibility and adaptivity to choosing degree of freedoms to construct a reduce model. In the second stage, we build an initial surrogate model based on the new prior by using GMsFEM and sparse generalized polynomial chaos (gPC)-based stochastic collocation methods. To improve the initial surrogate model, we dynamically update the surrogate model, which is adapted to the sequential availability of data and the updated analysis. The two-stage EnKF can achieve a better estimation than standard EnKF, and significantly improve the efficiency to update the ensemble analysis (posterior exploration). To enhance the applicability and flexibility in Bayesian inverse problems, we extend the two-stage EnKF to non-Gaussian models and hierarchical models. In the paper, we focus on the time fractional diffusion-wave models in porous media and investigate their Bayesian inverse problems using the proposed two-stage EnKF.


title: A Novel Asymmetric Coded Placement in Combination Networks with end-user Caches

abstract: The tradeoff between the user's memory size and the worst-case download time in the $(H,r,M,N)$ combination network is studied, where a central server communicates with $K$ users through $H$ immediate relays, and each user has local cache of size $M$ files and is connected to a different subset of $r$ relays. The main contribution of this paper is the design of a coded caching scheme with asymmetric coded placement by leveraging coordination among the relays, which was not exploited in past work. Mathematical analysis and numerical results show that the proposed schemes outperform existing schemes.


title: Importance sampling with imperfect cloning for the computation of generalized Lyapunov exponents

abstract: We revisit the numerical calculation of generalized Lyapunov exponents, $L$($q$), in deterministic dynamical systems. The standard method consists of adding noise to the dynamics in order to use importance sampling algorithms. Then $L$($q$) is obtained by taking the limit noise-amplitude $\to$0 after the calculation. We focus on a particular method that involves periodic cloning and pruning of a set of trajectories. However, instead of considering a noisy dynamics, we implement an imperfect (noisy) cloning. This alternative method is compared with the standard one and, when possible, with analytical results. As a workbench, we use the asymmetric tent map, the standard map, and a system of coupled symplectic maps. The general conclusion of this study is that the imperfect-cloning method performs as well as the standard one, with the advantage of preserving the deterministic dynamics.


title: Transient fractality as a mechanism for emergent irreversibility in chaotic Hamiltonian dynamics

abstract: Understanding irreversibility in macrophysics from reversible microphysics has been the holy grail in statistical physics ever since the mid-19th century. Here the central question concerns the arrow of time, which boils down to deriving macroscopic emergent irreversibility from microscopic reversible equations of motion. As suggested by Boltzmann, this irreversibility amounts to improbability (rather than impossibility) of the second-law-violating events. Later studies suggest that this improbability arises from a fractal attractor which is dynamically generated in phase space in reversible dissipative systems. However, the same mechanism seems inapplicable to reversible conservative systems, since a zero-volume fractal attractor is incompatible with the nonzero phase-space volume, which is a constant of motion due to the Liouville theorem. Here we demonstrate that in a Hamiltonian system the fractal scaling emerges transiently over an intermediate length scale. Notably, this transient fractality is unveiled by invoking the Loschmidt demon with an imperfect accuracy. Moreover, we show that irreversibility from the fractality can be evaluated by means of information theory and the fluctuation theorem. The fractality provides a unified understanding of emergent irreversibility over an intermediate time scale regardless of whether the underlying reversible dynamics is dissipative or conservative.


title: Accelerator modes and anomalous diffusion in 3D volume-preserving maps

abstract: Angle-action maps that are periodic in the action direction can have accelerator modes: orbits that are periodic when projected onto the torus, but that lift to unbounded orbits in an action variable. In this paper we construct a volume-preserving family of maps, with two angles and one action, that have accelerator modes created at Hopf-one (or saddle-center-Hopf) bifurcations. Near such a bifurcation we show that there is often a bubble of invariant tori. Computations of chaotic orbits near such a bubble show that the trapping times have an algebraic decay similar to that seen around stability islands in area-preserving maps. As in the 2D case, this gives rise to anomalous diffusive properties of the action in our 3D map.


title: Optical helicity and Hertz vectors

abstract: We study the conserved quantity associated with the dual symmetry of the Maxwell equations, called the optical helicity, by means of transverse Hertz vectors. In the presence of charges, its evolution yields the integral of $\bm{E}\cdot\bm{B}$ which is the anomaly term for chiral fermions. We also discuss the helicity change in condensed matter systems where topological magnetoelectric effect emerges. An alternative expression of the optical helicity is also found. Lastly, a dual symmetric Hertz Lagrangian is constructed and its conserved charge is derived.


title: On quadratic curves over finite fields

abstract: The geometry of algebraic curves over finite fields is a rich area of research. In previous work, the authors investigated a particular aspect of the geometry over finite fields of the classical unit circle, namely how the number of solutions of the circle equation depends on the characteristic $p$ and the degree $n\geq 1$ of the finite field $\mathbb{F}_{p^n}$. In this paper, we make a similar study of the geometry over finite fields of the quadratic curves defined by the quadratic equations in two variables for the classical conic sections. In particular the quadratic equation with mixed term is interesting, and our results display a rich variety of possibilities for the number of solutions to this equation over a finite field.


title: Goal-oriented adaptive surrogate construction for stochastic inversion

abstract: Stochastic inverse problems are generally solved by some form of finite sampling of a space of uncertain parameters. For computationally expensive models, surrogate response surfaces are often employed to increase the number of samples used in approximating the solution. The result is generally a trade off in errors where the stochastic error is reduced at the cost of an increase in deterministic/discretization errors in the evaluation of the surrogate. Such stochastic errors pollute predictions based on the stochastic inverse. In this work, we formulate a method for adaptively creating a special class of surrogate response surfaces with this stochastic error in mind. Adjoint techniques are used to enhance the local approximation properties of the surrogate allowing the construction of a higher-level enhanced surrogate. Using these two levels of surrogates, appropriately derived local error indicators are computed and used to guide refinement of both levels of the surrogates. Three types of refinement strategies are presented and combined in an iterative adaptive surrogate construction algorithm. Numerical examples, including a complex vibroacoustics application, demonstrate how this adaptive strategy allows for accurate predictions under uncertainty for a much smaller computational cost than uniform refinement.


title: An Approximate Pareto Set for Minimizing the Maximum Lateness and Makespan on Parallel Machines

abstract: We consider the two-parallel machines scheduling problem, with the aim of minimizing the maximum lateness and the makespan. Formally, the problem is defined as follows. We have to schedule a set J of n jobs on two identical machines. Each job i in J has a processing time p_i and a delivery time q_i. Each machine can only perform one job at a given time. The machines are available at time t=0 and each of them can process at most one job at a given time. The problem is to find a sequence of jobs, with the objective of minimizing the maximum lateness L_max and the makespan C_max. With no loss of generality, we consider that all data are integers and that jobs are indexed in non-increasing order of their delivery times: q_1 >= q_2 >= ... >= q_n. This paper proposes an exact algorithm (based on a dynamic programming) to generate the complete Pareto Frontier in a pseudo-polynomial time. Then, we present an FPTAS (Fully Polynomial Time Approximation Scheme) to generate an approximate Pareto Frontier, based on the conversion of the dynamic programming. The proposed FPTAS is strongly polynomial. Some numerical experiments are provided in order to compare the two proposed approaches.


title: As you like it: Localization via paired comparisons

abstract: Suppose that we wish to estimate a vector $\mathbf{x}$ from a set of binary paired comparisons of the form "$\mathbf{x}$ is closer to $\mathbf{p}$ than to $\mathbf{q}$" for various choices of vectors $\mathbf{p}$ and $\mathbf{q}$. The problem of estimating $\mathbf{x}$ from this type of observation arises in a variety of contexts, including nonmetric multidimensional scaling, "unfolding," and ranking problems, often because it provides a powerful and flexible model of preference. We describe theoretical bounds for how well we can expect to estimate $\mathbf{x}$ under a randomized model for $\mathbf{p}$ and $\mathbf{q}$. We also present results for the case where the comparisons are noisy and subject to some degree of error. Additionally, we show that under a randomized model for $\mathbf{p}$ and $\mathbf{q}$, a suitable number of binary paired comparisons yield a stable embedding of the space of target vectors. Finally, we also show that we can achieve significant gains by adaptively changing the distribution for choosing $\mathbf{p}$ and $\mathbf{q}$.


title: Partial Identification of Expectations with Interval Data

abstract: A conditional expectation function (CEF) can at best be partially identified when the conditioning variable is interval censored. When the number of bins is small, existing methods often yield minimally informative bounds. We propose three innovations that make meaningful inference possible in interval data contexts. First, we prove novel nonparametric bounds for contexts where the distribution of the censored variable is known. Second, we show that a class of measures that describe the conditional mean across a fixed interval of the conditioning space can often be bounded tightly even when the CEF itself cannot. Third, we show that a constraint on CEF curvature can either tighten bounds or can substitute for the monotonicity assumption often made in interval data applications. We derive analytical bounds that use the first two innovations, and develop a numerical method to calculate bounds under the third. We show the performance of the method in simulations and then present two applications. First, we resolve a known problem in the estimation of mortality as a function of education: because individuals with high school or less are a smaller and thus more negatively selected group over time, estimates of their mortality change are likely to be biased. Our method makes it possible to hold education rank bins constant over time, revealing that current estimates of rising mortality for less educated women are biased upward in some cases by a factor of three. Second, we apply the method to the estimation of intergenerational mobility, where researchers frequently use coarsely measured education data in the many contexts where matched parent-child income data are unavailable. Conventional measures like the rank-rank correlation may be uninformative once interval censoring is taken into account; CEF interval-based measures of mobility are bounded tightly.


title: Exact Controllability of linear KP-I equation

abstract: We prove the exact controllability of linear KP-I equation if the control input is added on a vertical domain. More generally, we have obtained the least dispersion needed to insure observability for fractional linear KP I equation.


title: Multicompartment Magnetic Resonance Fingerprinting

abstract: Magnetic resonance fingerprinting (MRF) is a technique for quantitative estimation of spin-relaxation parameters from magnetic-resonance data. Most current MRF approaches assume that only one tissue is present in each voxel, which neglects the tissue's microstructure, and may lead to artifacts in the recovered parameter maps at boundaries between tissues. In this work, we propose a multicompartment MRF model that accounts for the presence of multiple tissues per voxel. The model is fit to the data by iteratively solving a sparse linear inverse problem at each voxel, in order to express the magnetization signal as a linear combination of a few fingerprints in the precomputed dictionary. Thresholding-based methods commonly used for sparse recovery and compressed sensing do not perform well in this setting due to the high local coherence of the dictionary. Instead, we solve this challenging sparse-recovery problem by applying reweighted-l1-norm regularization, implemented using an efficient interior-point method. The proposed approach is validated with simulated data at different noise levels and undersampling factors, as well as with a controlled phantom imaging experiment on a clinical magnetic-resonance system.