title: Skin lesion detection based on an ensemble of deep convolutional neural network

abstract: Skin cancer is a major public health problem, with over 5 million newly diagnosed cases in the United States each year. Melanoma is the deadliest form of skin cancer, responsible for over 9,000 deaths each year. In this paper, we propose an ensemble of deep convolutional neural networks to classify dermoscopy images into three classes. To achieve the highest classification accuracy, we fuse the outputs of the softmax layers of four different neural architectures. For aggregation, we consider the individual accuracies of the networks weighted by the confidence values provided by their final softmax layers. This fusion-based approach outperformed all the individual neural networks regarding classification accuracy.


title: Origin of quantum oscillations in doped cuprates

abstract: It is proposed that Fermi-surface reconstruction in electron-doped Ln_{2-x}Ce_xCuO_4 (Ln = Pr, Nd) and in hole-doped YBa_2Cu_3O_{6+y} and YBa_2Cu_4O_8 occurs when the Fermi arcs extend into the second Brillouin zone (BZ). The criterion employs the axial component of the Fermi-arc tips, \hat{q} > 0.5, depending on both the position of the Fermi arc's center \dot{Q} and the incommensurabity \delta_c of unidirectional (striped) charge-density waves (CDWs). Qualitatively, the concave end-pieces of the Fermi arcs, terminated by Bragg-reflection mirrors due to the CDWs and severed at the boundary of the first BZ by lattice Bragg reflection, are assumed to join and relax to convex loops. Those entities may correspond to the electron pockets attributed to the quantum oscillations observed in these compounds. The criterion also explains why no quantum oscillations are found in the simple hole-doped lanthanum cuprates, La_{2-x}Ae_xCuO_4 (Ae = Sr, Ba), and in the bismuth cuprates Bi_2Sr_{2-x}La_xCuO_{6+y} and Bi_2Sr_2CaCu_2O_{8+y}. The possibility of quantum oscillations in hole-doped, partly substituted La_{2-y-x}Ln_ySr_xCuO_4 (Ln = Nd, Eu; y = 0.4, 0.2) in the high-end doping interval of their pseudogap phase, 0.182 < x < 0.235, is raised. A geometric modification of Bragg-reflection mirrors applies to HgBa_2CuO_{4+y} where CDWs are bidirectional (checkerboard-like).


title: Least action principle for Lorentz force in dilaton-Maxwell electrodynamics

abstract: The least action principle is established for the dynamics of a test particle in a dilaton-Maxwell background. These dynamics and background are invariant under the action of the dilatation transformation; explicit form of the corresponding generalization of the Lorentz force is established for the considered model. On a stationary background, we have found the integral of motion of the energy type. This integral is used to resolve the radial dynamics of test particles in a spherically symmetric electrostatic background.


title: BPS states of Fourfolds as candidates for Kaluza-Klein modes

abstract: Within the framework of the cosmological theory of the Big Bang, F-theory that unifies all four types of fundamental interactions is represented. Among the most exciting predictions of physics beyond the Standard Model is the assumption of the space of extra dimensions that solves the hierarchy problem. With the presence of this extra dimensions are connected the searches for Kaluza-Klein partners of gravitons, gauge bosons and microscopic black holes at the LHC. In the framework of F-theory, Calabi-Yau fourfold is considered as a space of extra dimensions. We study the duality between the F-theory compactified on the K3-surface and $E_8\times E_8$ heterotic string compactified on the torus, $T^2$. The set of BPS states corresponding to the Calabi-Yau fourfolds, which has either an elliptic curve or a K3-fibration as a layer, is studied in the aspect of correspondence to the KK modes of the M-theory on $R^8\times T^2 \times S^1 / Z_2$. The singularities of the moduli space of the Calabi-Yau fourfold make it possible to observe massive KK modes, the masses of which are obtained from the M-theory of supergravity. The result is of interest for a theoretical understanding of the KK modes, the experimental searches for which are carried out at the LHC.


title: Coexistence of Radar and Communication Systems in CBRS Bands Through Downlink Power Control

abstract: Citizen Broadband Radio Service band (3550 - 3700 GHz) is seen as one of the key frequency bands to enable improvements in performance of wireless broadband and cellular systems. A careful study of interference caused by a secondary cellular communication system coexisting with an incumbent naval radar is required to establish a pragmatic protection distance, which not only protects the incumbent from harmful interference but also increases the spectrum access opportunity for the secondary system. In this context, this paper investigates the co-channel and adjacent channel coexistence of a ship-borne naval radar and a wide-area cellular communication system and presents the analysis of interference caused by downlink transmission in the cellular system on the naval radar for different values of radar protection distance. The results of such analysis suggest that maintaining a protection distance of 30 km from the radar will ensure the required INR protection criterion of -6 dB at the radar receiver with > 0.9 probability, even when the secondary network operates in the same channel as the radar. Novel power control algorithms to assign operating powers to the coexisting cellular devices are also proposed to further reduce the protection distance from radar while still meeting the radar INR protection requirement.


title: Dynamical dissociation of quarkonia by wave function decoherence

abstract: We investigate the real-time evolution of quarkonium bound states in a quark-gluon plasma in one dimension using an improved QCD based stochastic potential model. This model describes the quarkonium dynamics in terms of a Schr\"odinger equation with an in-medium potential and two noise terms encoding the residual interactions between the heavy quarks and the medium. The probabilities of bound states in a static medium and in a boost-invariantly expanding quark-gluon plasma are discussed. We draw two conclusions from our results: One is that the outcome of the stochastic potential model is qualitatively consistent with the experimental data in relativistic heavy-ion collisions. The other is that the noise plays an important role in order to describe quarkonium dynamics in medium, in particular it causes decoherence of the quarkonium wave function. The effectiveness of decoherence is controlled by a new length scale $l_{\rm corr}$. It represents the noise correlation length and its effect has not been included in existing phenomenological studies.


title: Frequency Switching for Simultaneous Wireless Information and Power Transfer

abstract: A new frequency switching receiver structure is proposed for simultaneous wireless information and power transfer in multi-carrier communication systems. Each subcarrier is switched to either the energy harvesting unit or the information decoding unit, according to the optimal subcarrier allocation. To implement the system, one-bit feedback is required for each subcarrier. Two optimization problems are defined, converted to binary knapsack problems, and solved using dynamic programming approaches. Upper bounds are obtained using continuous relaxations. Power allocation is integrated to further increase the performance. Numerical studies show that the proposed frequency switching based model is better than existing models in a wide range of parameters.


title: Special tilting modules for algebras with positive dominant dimension

abstract: We study certain special tilting and cotilting modules for an algebra with positive dominant dimension, each of which is generated or cogenerated (and usually both) by projective-injectives. These modules have various interesting properties, for example that their endomorphism algebras always have global dimension at most that of the original algebra. We characterise minimal d-Auslander-Gorenstein algebras and d-Auslander algebras via the property that these special tilting and cotilting modules coincide. By the Morita-Tachikawa correspondence, any algebra of dominant dimension at least 2 may be expressed (essentially uniquely) as the endomorphism algebra of a generator-cogenerator for another algebra, and we also study our special tilting and cotilting modules from this point of view, via the theory of recollements and intermediate extension functors.


title: Diagrammatic Monte-Carlo for weak-coupling expansion of non-Abelian lattice field theories: large-N U(N)xU(N) principal chiral model

abstract: We develop numerical tools for Diagrammatic Monte-Carlo simulations of non-Abelian lattice field theories in the t'Hooft large-N limit based on the weak-coupling expansion. First we note that the path integral measure of such theories contributes a bare mass term in the effective action which is proportional to the bare coupling constant. This mass term renders the perturbative expansion infrared-finite and allows to study it directly in the large-N and infinite-volume limits using the Diagrammatic Monte-Carlo approach. On the exactly solvable example of a large-N O(N) sigma model in D=2 dimensions we show that this infrared-finite weak-coupling expansion contains, in addition to powers of bare coupling, also powers of its logarithm, reminiscent of re-summed perturbation theory in thermal field theory and resurgent trans-series without exponential terms. We numerically demonstrate the convergence of these double series to the manifestly non-perturbative dynamical mass gap. We then develop a Diagrammatic Monte-Carlo algorithm for sampling planar diagrams in the large-N matrix field theory, and apply it to study this infrared-finite weak-coupling expansion for large-N U(N)xU(N) nonlinear sigma model (principal chiral model) in D=2. We sample up to 12 leading orders of the weak-coupling expansion, which is the practical limit set by the increasingly strong sign problem at high orders. Comparing Diagrammatic Monte-Carlo with conventional Monte-Carlo simulations extrapolated to infinite N, we find a good agreement for the energy density as well as for the critical temperature of the "deconfinement" transition. Finally, we comment on the applicability of our approach to planar QCD at zero and finite density.


title: Statistical Physics of Balance Theory

abstract: Triadic relationships are accepted to play a key role in the dynamics of social and political networks. Building on insights gleaned from balance theory in social network studies and from Boltzmann-Gibbs statistical physics, we propose a model to quantitatively capture the dynamics of the four types of triadic relationships in a network. Central to our model are the triads' incidence rates and the idea that those can be modeled by assigning a specific triadic energy to each type of triadic relation. We emphasize the role of the degeneracy of the different triads and how it impacts the degree of frustration in the political network. In order to account for a persistent form of disorder in the formation of the triadic relationships, we introduce the systemic variable temperature. In order to learn about the dynamics and motives, we propose a generic Hamiltonian with three terms to model the triadic energies. One term is connected with a three-body interaction that captures balance theory. The other terms take into account the impact of heterogeneity and of negative edges in the triads. The validity of our model is tested on four datasets including the time series of triadic relationships for the standings between two classes of alliances in a massively multiplayer online game (MMOG). We also analyze real-world data for the relationships between the "agents" involved in the Syrian civil war, and in the relations between countries during the Cold War era. We find emerging properties in the triadic relationships in a political network, for example reflecting itself in a persistent hierarchy between the four triadic energies, and in the consistency of the extracted parameters from comparing the model Hamiltonian to the data.


title: Massive Star Cluster Formation and Destruction in Luminous Infrared Galaxies in GOALS

abstract: We present the results of a {\it Hubble Space Telescope} ACS/HRC FUV, ACS/WFC optical study into the cluster populations of a sample of 22 Luminous Infrared Galaxies in the Great Observatories All-Sky LIRG Survey. Through integrated broadband photometry we have derived ages and masses for a total of 484 star clusters contained within these systems. This allows us to examine the properties of star clusters found in the extreme environments of LIRGs relative to lower luminosity star-forming galaxies in the local Universe. We find that by adopting a Bruzual \& Charlot simple stellar population (SSP) model and Salpeter initial mass function, the age distribution of clusters declines as $dN/d\tau = \tau^{-0.9 +/- 0.3}$, consistent with the age distribution derived for the Antennae Galaxies, and interpreted as evidence for rapid cluster disruption occuring in the strong tidal fields of merging galaxies. The large number of $10^{6} M_{\odot}$ young clusters identified in the sample also suggests that LIRGs are capable of producing more high-mass clusters than what is observed to date in any lower luminosity star-forming galaxy in the local Universe. The observed cluster mass distribution of $dN/dM = M^{-1.95 +/- 0.11}$ is consistent with the canonical -2 power law used to describe the underlying initial cluster mass function (ICMF) for a wide range of galactic environments. We interpret this as evidence against mass-dependent cluster disruption, which would flatten the observed CMF relative to the underlying ICMF distribution.


title: Ab initio study of magnetocrystalline anisotropy, magnetostriction, and Fermi surface of L10 FeNi (tetrataenite)

abstract: The ordered L1$_0$ FeNi phase (tetrataenite) is recently considered as a promising candidate for the rare-earth free permanent magnets applications. In this work we calculate several characteristics of the L1$_0$ FeNi, where most of the results come form the fully relativistic full potential FPLO method with the generalized gradient approximation (GGA). A special attention deserves the summary of the magnetocrystalline anisotropy energies (MAE's), the full potential calculations of the anisotropy constant $K_3$, and the combined analysis of the Fermi surface and three-dimensional $\mathbf{k}$-resolved MAE. Other calculated parameters presented in this article are the magnetic moments $m_{s}$ and $m_{l}$, magnetostrictive coefficient $\lambda_{001}$, bulk modulus B$_0$, and lattice parameters. The MAE's summary shows rather big discrepancies between the experimental MAE's from literature and also between the calculated MAE's.


title: Bayesian Joint Topic Modelling for Weakly Supervised Object Localisation

abstract: We address the problem of localisation of objects as bounding boxes in images with weak labels. This weakly supervised object localisation problem has been tackled in the past using discriminative models where each object class is localised independently from other classes. We propose a novel framework based on Bayesian joint topic modelling. Our framework has three distinctive advantages over previous works: (1) All object classes and image backgrounds are modelled jointly together in a single generative model so that "explaining away" inference can resolve ambiguity and lead to better learning and localisation. (2) The Bayesian formulation of the model enables easy integration of prior knowledge about object appearance to compensate for limited supervision. (3) Our model can be learned with a mixture of weakly labelled and unlabelled data, allowing the large volume of unlabelled images on the Internet to be exploited for learning. Extensive experiments on the challenging VOC dataset demonstrate that our approach outperforms the state-of-the-art competitors.


title: Small-scale filament eruptions as the driver of solar coronal hole X-ray jets

abstract: Solar X-ray jets are evidently made by a burst of reconnection of closed magnetic field in a jet's base with ambient "open" field (1,2). In the widely-accepted version of the "emerging-flux" model, that reconnection occurs at a current sheet between the open field and emerging closed field and also makes a compact hot brightening that is usually observed at the edge of the jet's base (1,3). Here we report on high-resolution X-ray and EUV observations of 20 randomly-selected X-ray jets in polar coronal holes. In each jet, contrary to the emerging-flux model, a miniature version of the filament eruptions that initiate coronal mass ejections (CMEs) (4-7) drives the jet-producing reconnection, and the compact hot brightening is made by internal reconnection of the legs of the minifilament-carrying erupting closed field, analogous to solar flares of larger-scale eruptions. Previous observations have found that some jets are driven by base-field eruptions (8-10,12), but only one such study, of only one jet, provisionally questioned the emerging-flux model (13). Our observations support the view that solar filament eruptions are made by a fundamental explosive magnetic process that occurs on a vast range of scales, from the biggest CME/flare eruptions down to X-ray jets, and perhaps down to even smaller jets that are candidates for powering coronal heating (10,14,15). A picture similar to that suggested by our observations was drawn before, inferred from different observations and based on a different origin of the erupting minifilament flux rope (11) (see Methods).


title: Measurement of $b$-hadron pair production with the ATLAS detector in proton-proton collisions at $\sqrt{s}=8$ TeV

abstract: A measurement of $b$-hadron pair production is presented, based on a data set corresponding to an integrated luminosity of 11.4 fb$^{-1}$ of proton--proton collisions recorded at $\sqrt{s}=8$ TeV with the ATLAS detector at the LHC. Events are selected in which a $b$-hadron is reconstructed in a decay channel containing $J/\psi \rightarrow \mu\mu$, and a second $b$-hadron is reconstructed in a decay channel containing a muon. Results are presented in a fiducial volume defined by kinematic requirements on three muons based on those used in the analysis. The fiducial cross section is measured to be $17.7 \pm 0.1 ($stat.$) \pm 2.0 ($syst.$)$ nb. A number of normalised differential cross sections are also measured, and compared to predictions from the Pythia8, Herwig++, MadGraph5\_aMC@NLO+Pythia8 and Sherpa event generators, providing new constraints on heavy flavour production.


title: The effect of finite-conductvity Hartmann walls on the linear stability of Hunt's flow

abstract: We analyse numerically the linear stability of the fully developed liquid metal flow in a square duct with insulating side walls and thin electrically conducting horizontal walls with the wall conductance ratio $c=0.01\cdots 1$ subject to a vertical magnetic field with the Hartmann numbers up to $Ha=10^{4}.$ In a sufficiently strong magnetic field, the flow consists of two jets at the side walls walls and a near-stagnant core with the relative velocity $\sim(cHa)^{-1}.$ We find that for $Ha\gtrsim300,$ the effect of wall conductivity on the stability of the flow is mainly determined by the effective Hartmann wall conductance ratio $cHa.$ For $c\ll 1,$ the increase of the magnetic field or that of the wall conductivity has a destabilizing effect on the flow. Maximal destabilization of the flow occurs at $Ha\approx30/c.$ In a stronger magnetic field with $cHa\gtrsim 30,$ the destabilizing effect vanishes and the asymptotic results of Priede et al. [J. Fluid Mech. 649, 115, 2010] for the ideal Hunt's flow with perfectly conducting Hartmann walls are recovered.


title: Optimal frame designs for multitasking devices with weight restrictions

abstract: Let $\mathbf d=(d_j)_{j\in\mathbb I_m}\in\mathbb N^m$ be a finite sequence (of dimensions) and $\alpha=(\alpha_i)_{i\in\mathbb I_n}$ be a sequence of positive numbers (of weights), where $\mathbb I_k=\{1,\ldots,k\}$ for $k\in\mathbb N$. We introduce the $(\alpha\, , \,\mathbf d)$-designs i.e., $m$-tuples $\Phi=(\mathcal F_j)_{j\in\mathbb I_m}$ such that $\mathcal F_j=\{f_{ij}\}_{i\in\mathbb I_n}$ is a finite sequence in $\mathbb C^{d_j}$, $j\in\mathbb I_m$, and such that the sequence of non-negative numbers $(\|f_{ij}\|^2)_{j\in\mathbb I_m}$ forms a partition of $\alpha_i$, $i\in\mathbb I_n$. We characterize the existence of $(\alpha\, , \, \mathbf d)$-designs with prescribed properties in terms of majorization relations. We show, by means of a finite-step algorithm, that there exist $(\alpha\, , \, \mathbf d)$-designs $\Phi^{\rm op}=(\mathcal F_j^{\rm op})_{j\in\mathbb I_m}$ that are universally optimal; that is, for every convex function $\varphi:[0,\infty)\rightarrow [0,\infty)$ then $\Phi^{\rm op}$ minimizes the joint convex potential induced by $\varphi$ among $(\alpha\, , \, \mathbf d)$-designs, namely $$ \sum_{j\in\mathbb I_m}\text{P}_\varphi(\mathcal F_j^{\rm op})\leq \sum_{j\in \mathbb I_m}\text{P}_\varphi(\mathcal F_j) $$ for every $(\alpha\, , \, \mathbf d)$-design $\Phi=(\mathcal F_j)_{j\in\mathbb I_m}$, where $\text{P}_\varphi(\mathcal F)=tr(\varphi(S_{\mathcal F}))$; in particular, $\Phi^{\rm op}$ minimizes both the joint frame potential and the joint mean square error among $(\alpha\, , \, \mathbf d)$-designs. We show that in this case $\mathcal F_j^{\rm op}$ is a frame for $\mathbb C^{d_j}$, for $j\in\mathbb I_m$. This corresponds to the existence of optimal encoding-decoding schemes for multitasking devices with energy restrictions.


title: Real time observation of granular rock analogue material deformation and failure using nonlinear laser interferometry

abstract: A better understanding and anticipation of natural processes such as landsliding or seismic fault activity requires detailed theoretical and experimental analysis of rock mechanics and geomaterial dynamics. These last decades, considerable progress has been made towards understanding deformation and fracture process in laboratory experiment on granular rock materials, as the well-known shear banding experiment. One of the reasons for this progress is the continuous improvement in the instrumental techniques of observation. But the lack of real time methods does not allow the detection of indicators of the upcoming fracture process and thus to anticipate the phenomenon. Here, we have performed uniaxial compression experiments to analyse the response of a granular rock material sample to different shocks. We use a novel interferometric laser sensor based on the nonlinear self-mixing interferometry technique to observe in real time the deformations of the sample and assess its usefulness as a diagnostic tool for the analysis of geomaterial dynamics. Due to the high spatial and temporal resolution of this approach, we observe both vibrations processes in response to a dynamic loading and the onset of failure. The latter is preceded by a continuous variation of vibration period of the material. After several shocks, the material response is no longer reversible and we detect a progressive accumulation of irreversible deformation leading to the fracture process. We demonstrate that material failure is anticipated by the critical slowing down of the surface vibrational motion, which may therefore be envisioned as an early warning signal or predictor to the macroscopic failure of the sample. The nonlinear self-mixing interferometry technique is readily extensible to fault propagation measurements. As such, it opens a new window of observation for the study of geomaterial deformation and failure.


title: Banach Algebras Associated to Metric Operator Fields

abstract: Motivated by noncommutative geometry and quantum physics, the concept of `metric operator field' is introduced. Roughly speaking, a metric operator field is a vector field on a set with values in self tensor product of a bundle of C*-algebras, satisfying properties similar to an ordinary metric (distance function). It is proved that to any such object there naturally correspond a Banach *-algebra that we call Lipschitz algebra, a class of probabilistic metrics, and (under some conditions) a (nontrivial) continuous field of C*-algebras in the sense of Dixmier. It is proved that for metric operator fields with values in von Neumann algebras the associated Lipschitz algebras are dual Banach spaces, and under some conditions, they are not amenable Banach algebras. Some examples and constructions are considered. We also discuss very briefly a possible application to quantum gravity.


title: Multi-mode active control of friction, dynamic ratchets and actuators

abstract: Active control of friction by ultrasonic vibration is a well-known effect with numerous technical applications ranging from press forming to micromechanical actuators. Reduction of friction is observed with vibration applied in any of the three possible directions (normal to the contact plane, in the direction of motion and in-plane transverse). In this work, we consider the multi-mode active control of sliding friction, where phase-shifted oscillations in two or more directions act at the same time. Our analysis is based on a macroscopic contact-mechanical model that was recently shown to be well-suited for describing dynamic frictional processes. For simplicity, we limit our analysis to a constant, load-independent normal and tangential stiffness and two superimposed phase-shifted harmonic oscillations, one of them being normal to the plane and the other in the direction of motion. As in previous works utilizing the present model, we assume a constant local coefficient of friction, with reduction of the observed force of friction arising entirely from the macroscopic dynamics of the system. Our numerical simulations show that the resulting law of friction is determined by just three dimensionless parameters. Depending on the values of these parameters, three qualitatively different types of behavior are observed: (a) symmetric velocity-dependence of the coefficient of friction (same for positive and negative velocities), (b) asymmetric dependence with respect to the sign of the velocity, but with zero force at zero velocity, and (c) asymmetric dependence with non-zero force at zero velocity. The latter two cases can be interpreted as a "dynamic ratchet" (b) and an actuator (c).


title: Limits on spin-dependent WIMP-nucleon cross section obtained from the complete LUX exposure

abstract: We present experimental constraints on the spin-dependent WIMP-nucleon elastic cross sections from the total 129.5 kg-year exposure acquired by the Large Underground Xenon experiment (LUX), operating at the Sanford Underground Research Facility in Lead, South Dakota (USA). A profile likelihood ratio analysis allows 90% CL upper limits to be set on the WIMP-neutron (WIMP-proton) cross section of $\sigma_n$ = 1.6$\times 10^{-41}$ cm$^{2}$ ($\sigma_p$ = 5$\times 10^{-40}$ cm$^{2}$) at 35 GeV$c^{-2}$, almost a sixfold improvement over the previous LUX spin-dependent results. The spin-dependent WIMP-neutron limit is the most sensitive constraint to date.


title: A note on the uniqueness of models in social abstract argumentation

abstract: Social abstract argumentation is a principled way to assign values to conflicting (weighted) arguments. In this note we discuss the important property of the uniqueness of the model.


title: Feasibility of measuring the magnetic dipole moments of the charm baryons at the LHC using bent crystals

abstract: In this paper we revisit the idea of measuring the magnetic dipole moments of the charm baryons and, in particular, of charmed Lambda by studying the spin precession induced by the strong effective magnetic field inside the channels of a bent crystal. We present a detailed sensitivity study showing the feasibility of such an experiment at the LHC in the coming years.


title: Investigating the Role of Coherence Effects on Jet Quenching in Pb-Pb Collisions at $\sqrt{s_{NN}} =2.76$ TeV using Jet Substructure

abstract: We report measurements of two jet shapes, the ratio of 2-Subjettiness to 1-Subjettiness ($\it{\tau_{2}}/\it{\tau_{1}}$) and the opening angle between the two axes of the 2-Subjettiness jet shape, which is obtained by reclustering the jet with the exclusive-$\it{k}_{T}$ algorithm~\cite{kT} and undoing the final clustering step. The aim of this measurement is to explore a possible change in the rate of 2-pronged objects in Pb-Pb compared to pp due to colour coherence. Coherence effects~\cite{ColourCoherence} relate to the ability of the medium to resolve a jet's substructure, which has an impact on the energy loss magnitude and mechanism of the traversing jet. In both collision systems charged jets are found with the anti-$\it{k}_{T}$ algorithm~\cite{Anti-kT}, a resolution parameter of $\it{R}=$0.4\ and a constituent cut off of $0.15$ GeV. This analysis uses hadron-jet coincidence techniques in Pb-Pb collisions to reject the combinatorial background and corrects further for background effects by employing various jet shape subtraction techniques and two dimensional unfolding. Measurements of the Nsubjettiness for jet momenta of $40-60$ GeV/$\it{c}$ in Pb-Pb collisions at $\sqrt{s_{NN}} =2.76$ TeV and pp collisions at $\sqrt{s} =7$ TeV will be presented and compared to PYTHIA simulations.


title: Time evolution of the spectral break in the high-energy extra component of GRB 090926A

abstract: We report on a new time-resolved analysis of the GRB 090926A broad-band spectrum during its prompt phase, and on its interpretation in the framework of prompt emission models. We characterize the high-energy extra component of GRB 090926A using Pass 8 data from the Fermi Large Area Telescope (LAT), which offer a greater sensitivity than any data set used in previous studies of this burst, particularly in the 30-100 MeV energy band. Then, we combine the LAT data with the Fermi Gamma-ray Burst Monitor (GBM) in joint spectral fits in order to characterize the time evolution of the broad-band spectrum from keV to GeV energies. Careful attention is paid to the LAT systematic effects. Finally, we perform a temporal analysis of the light curves and we compute the variability time scales during and after the bright spike at ~10 s post-trigger. Our analysis confirms and better constrains the spectral break in the high-energy extra component, which has been previously reported during the bright spike. Furthermore, it reveals that the spectral attenuation persists at later times, with an increase of the break characteristic energy up to the GeV domain until the end of the prompt phase. We interpret the high-energy spectral break as caused by photon opacity to pair creation. Requiring that all emissions are produced above the photosphere of GRB 090926A, we compute the bulk Lorentz factor of the outflow, $\Gamma$. The latter decreases from 230 during the spike to 100 at the end of the prompt emission. Assuming, instead, that the spectral break reflects the natural curvature of the inverse Compton spectrum, lower limits corresponding to larger values of $\Gamma$ are also derived. Combined with the extreme temporal variability of GRB 090926A, these Lorentz factors lead to emission radii $R\sim10^{14}$ cm which are consistent with an internal origin of both the keV-MeV and GeV prompt emissions.


title: Optimal Computation of Overabundant Words

abstract: The observed frequency of the longest proper prefix, the longest proper suffix, and the longest infix of a word $w$ in a given sequence $x$ can be used for classifying $w$ as avoided or overabundant. The definitions used for the expectation and deviation of $w$ in this statistical model were described and biologically justified by Brendel et al. (J Biomol Struct Dyn 1986). We have very recently introduced a time-optimal algorithm for computing all avoided words of a given sequence over an integer alphabet (Algorithms Mol Biol 2017). In this article, we extend this study by presenting an $\mathcal{O}(n)$-time and $\mathcal{O}(n)$-space algorithm for computing all overabundant words in a sequence $x$ of length $n$ over an integer alphabet. Our main result is based on a new non-trivial combinatorial property of the suffix tree $\mathcal{T}$ of $x$: the number of distinct factors of $x$ whose longest infix is the label of an explicit node of $\mathcal{T}$ is no more than $3n-4$. We further show that the presented algorithm is time-optimal by proving that $\mathcal{O}(n)$ is a tight upper bound for the number of overabundant words. Finally, we present experimental results, using both synthetic and real data, which justify the effectiveness and efficiency of our approach in practical terms.


title: Cell Tracking via Proposal Generation and Selection

abstract: Microscopy imaging plays a vital role in understanding many biological processes in development and disease. The recent advances in automation of microscopes and development of methods and markers for live cell imaging has led to rapid growth in the amount of image data being captured. To efficiently and reliably extract useful insights from these captured sequences, automated cell tracking is essential. This is a challenging problem due to large variation in the appearance and shapes of cells depending on many factors including imaging methodology, biological characteristics of cells, cell matrix composition, labeling methodology, etc. Often cell tracking methods require a sequence-specific segmentation method and manual tuning of many tracking parameters, which limits their applicability to sequences other than those they are designed for. In this paper, we propose 1) a deep learning based cell proposal method, which proposes candidates for cells along with their scores, and 2) a cell tracking method, which links proposals in adjacent frames in a graphical model using edges representing different cellular events and poses joint cell detection and tracking as the selection of a subset of cell and edge proposals. Our method is completely automated and given enough training data can be applied to a wide variety of microscopy sequences. We evaluate our method on multiple fluorescence and phase contrast microscopy sequences containing cells of various shapes and appearances from ISBI cell tracking challenge, and show that our method outperforms existing cell tracking methods. Code is available at: https://github.com/SaadUllahAkram/CellTracker


title: Generative Adversarial Trainer: Defense to Adversarial Perturbations with GAN

abstract: We propose a novel technique to make neural network robust to adversarial examples using a generative adversarial network. We alternately train both classifier and generator networks. The generator network generates an adversarial perturbation that can easily fool the classifier network by using a gradient of each image. Simultaneously, the classifier network is trained to classify correctly both original and adversarial images generated by the generator. These procedures help the classifier network to become more robust to adversarial perturbations. Furthermore, our adversarial training framework efficiently reduces overfitting and outperforms other regularization methods such as Dropout. We applied our method to supervised learning for CIFAR datasets, and experimantal results show that our method significantly lowers the generalization error of the network. To the best of our knowledge, this is the first method which uses GAN to improve supervised learning.


title: The Mira-Titan Universe II: Matter Power Spectrum Emulation

abstract: We introduce a new cosmic emulator for the matter power spectrum covering eight cosmological parameters. Targeted at optical surveys, the emulator provides accurate predictions out to a wavenumber k~5/Mpc and redshift z<=2. Besides covering the standard set of LCDM parameters, massive neutrinos and a dynamical dark energy of state are included. The emulator is built on a sample set of 36 cosmological models, carefully chosen to provide accurate predictions over the wide and large parameter space. For each model, we have performed a high-resolution simulation, augmented with sixteen medium-resolution simulations and TimeRG perturbation theory results to provide accurate coverage of a wide k-range; the dataset generated as part of this project is more than 1.2Pbyte. With the current set of simulated models, we achieve an accuracy of approximately 4%. Because the sampling approach used here has established convergence and error-control properties, follow-on results with more than a hundred cosmological models will soon achieve ~1% accuracy. We compare our approach with other prediction schemes that are based on halo model ideas and remapping approaches. The new emulator code is publicly available.


title: Logical Parsing from Natural Language Based on a Neural Translation Model

abstract: Semantic parsing has emerged as a significant and powerful paradigm for natural language interface and question answering systems. Traditional methods of building a semantic parser rely on high-quality lexicons, hand-crafted grammars and linguistic features which are limited by applied domain or representation. In this paper, we propose a general approach to learn from denotations based on Seq2Seq model augmented with attention mechanism. We encode input sequence into vectors and use dynamic programming to infer candidate logical forms. We utilize the fact that similar utterances should have similar logical forms to help reduce the searching space. Under our learning policy, the Seq2Seq model can learn mappings gradually with noises. Curriculum learning is adopted to make the learning smoother. We test our method on the arithmetic domain which shows our model can successfully infer the correct logical forms and learn the word meanings, compositionality and operation orders simultaneously.


title: Constructing current singularity in a 3D line-tied plasma

abstract: We revisit Parker's conjecture of current singularity formation in 3D line-tied plasmas using a recently developed numerical method, variational integration for ideal magnetohydrodynamics in Lagrangian labeling. With the frozen-in equation built-in, the method is free of artificial reconnection, and hence it is arguably an optimal tool for studying current singularity formation. Using this method, the formation of current singularity has previously been confirmed in the Hahm--Kulsrud--Taylor problem in 2D. In this paper, we extend this problem to 3D line-tied geometry. The linear solution, which is singular in 2D, is found to be smooth for arbitrary system length. However, with finite amplitude, the linear solution can become pathological when the system is sufficiently long. The nonlinear solutions turn out to be smooth for short systems. Nonetheless, the scaling of peak current density versus system length suggests that the nonlinear solution may become singular at finite length. With the results in hand, we can neither confirm nor rule out this possibility conclusively, since we cannot obtain solutions with system length near the extrapolated critical value.


title: Increasing the Discovery Power and Confidence Levels of Disease Association Studies: A Survey

abstract: The majority of common diseases are influenced by multiple genetic and environmental factors such as Cancer. Even though uncovering the main causes of disease is deemed difficult due to the complexity of gene-gene and gene-environment interactions, major research efforts aim at identifying disease risk factors, especially genetic ones. Over the past decade, disease association studies have been used to uncover the susceptibility, aetiology and mechanisms of action pertaining to common diseases. In disease association studies, genetic data is analyzed in order to reveal the relationship between different types of variants, and a disease of interest. The ultimate goal of association studies is to facilitate susceptibility testing for disease prediction, early diagnosis and enhanced prognosis . Susceptibility testing and disease prediction are particularly important for diseases that can be prevented by diet, drugs or change in lifestyle. The discovered associations assist in understanding the molecular mechanisms influenced by the reported variants, and in identifying important risk factors. Current association studies suffer from several shortcomings. This report surveys the literature that addresses the shortcomings of current methods the identify genetic disease associations. In addition, it reviews the suggested solutions that either enhance some aspect of the methodologies, or complement them.


title: Asynchronous Announcements

abstract: We propose a multi-agent epistemic logic of asynchronous announcements, where truthful announcements are publicly sent but individually received by agents, and in the order in which they were sent. Additional to epistemic modalities the logic contains dynamic modalities for making announcements and for receiving them. What an agent believes is a function of her initial uncertainty and of the announcements she has received. Beliefs need not be truthful, because announcements already made may not yet have been received. As announcements are true when sent, certain message sequences can be ruled out, just like inconsistent cuts in distributed computing. We provide a complete axiomatization for this \emph{asynchronous announcement logic} (AA). It is a reduction system that also demonstrates that any formula in $AA$ is equivalent to one without dynamic modalities, just as for public announcement logic. A detailed example modelling message exchanging processes in distributed computing in $AA$ closes our investigation.


title: Simple Witt modules that are finitely generated over the cartan subalgebra

abstract: Let $d\ge1$ be an integer, $W_d$ and $\mathcal{K}_d$ be the Witt algebra and the weyl algebra over the Laurent polynomial algebra $A_d=\mathbb{C} [x_1^{\pm1}, x_2^{\pm1}, ..., x_d^{\pm1}]$, respectively. For any $\mathfrak{gl}_d$-module $M$ and any admissible module $P$ over the extended Witt algebra $\widetilde W_d$, we define a $W_d$-module structure on the tensor product $P\otimes M$. We prove in this paper that any simple $W_d$-module that is finitely generated over the cartan subalgebra is a quotient module of the $W_d$-module $P \otimes M$ for a finite dimensional simple $\mathfrak{gl}_d$-module $M$ and a simple $\mathcal{K}_d$-module $P$ that are finitely generated over the cartan subalgebra. We also characterize all simple $\mathcal{K}_d$-modules and all simple admissible $\widetilde W_d$-modules that are finitely generated over the cartan subalgebra.


title: That is not dead which can eternal lie: the aestivation hypothesis for resolving Fermi's paradox

abstract: If a civilization wants to maximize computation it appears rational to aestivate until the far future in order to exploit the low temperature environment: this can produce a $10^{30}$ multiplier of achievable computation. We hence suggest the "aestivation hypothesis": the reason we are not observing manifestations of alien civilizations is that they are currently (mostly) inactive, patiently waiting for future cosmic eras. This paper analyzes the assumptions going into the hypothesis and how physical law and observational evidence constrain the motivations of aliens compatible with the hypothesis.


title: Chimera states in multi-strain epidemic models with temporary immunity

abstract: We investigate a time-delayed epidemic model for multi-strain diseases with temporary immunity. In the absence of cross-immunity between strains, dynamics of each individual strain exhibits emergence and anni- hilation of limit cycles due to a Hopf bifurcation of the endemic equilibrium, and a saddle-node bifurcation of limit cycles depending on the time delay associated with duration of temporary immunity. Effects of all-to-all and non-local coupling topologies are systematically investigated by means of numerical simulations, and they suggest that cross-immunity is able to induce a diverse range of complex dynamical behaviors and synchro- nization patterns, including discrete traveling waves, solitary states, and amplitude chimeras. Interestingly, chimera states are observed for narrower cross-immunity kernels, which can have profound implications for understanding the dynamics of multi-strain diseases.


title: Machine Learning Techniques for Mortality Modeling

abstract: Various stochastic models have been proposed to estimate mortality rates. In this paper we illustrate how machine learning techniques allow us to analyze the quality of such mortality models. In addition, we present how these techniques can be used for differentiating the different causes of death in mortality modeling.


title: Topological spin dynamics in cubic FeGe near room temperature

abstract: Understanding spin-wave dynamics in chiral magnets is a key step for the development of high-speed, spin-wave based spintronic devices that take advantage of chiral and topological spin textures for their operation. Here we present an experimental and theoretical study of spin-wave dynamics in a cubic B20 FeGe single crystal. Using the combination of waveguide microwave absorption spectroscopy (MAS), micromagnetic simulations, and analytical theory, we identify the resonance dynamics in all magnetic phases (field polarized, conical, helical, and skyrmion phases). Because the resonance frequencies of specific chiral spin textures are unique, quantitative agreement between our theoretical predictions and experimental findings for all resonance frequencies and spin wave modes enables us to unambiguously identify chiral magnetic phases and to demonstrate that MAS is a powerful tool to efficiently extract a magnetic phase diagram. These results provide a new tool to accelerate the integration of chiral magnetic materials into spintronic devices.


title: The evolving system of Trypillian settlements

abstract: Archaeological settlement systems are usually analysed in terms of the relation between the rank of a settlement and its size (Zipf's law). We argue that this approach is unreliable, and can be misleading, in application to archaeological data where the recovery rate of settlements can be low and their size estimates are often approximate at best. An alternative framework for the settlement data interpretation, better connected with the theoretical concepts of the stochastic evolution of settlement systems, is applied to the evolving system of settlements of the Late Neolithic-Bronze Age Trypillia cultural complex (5,400-2,800 BC) in modern Ukraine. The stochastic evolution model provides a consistent and accurate explanation of the frequency of occurrence of Trypillian settlement areas in the range from 0.05 to 500 ha. Thus validated, the model leads to reliable estimates of the typical size of a newly formed settlement as well as the growth rates of the total number of settlements and their areas. The parameters of the settlement system thus revealed are consistent with palaeoeconomy reconstructions for the Trypillia area.


title: Nucleon axial form factors using lattice QCD simulations with a physical value of the pion mass

abstract: We present results on the nucleon axial and induced pseudo-scalar form factors using an ensemble of two degenerate twisted mass clover-improved fermions generated at the physical value of the pion mass. We evaluate the isovector and the isoscalar, as well as, the strange and the charm axial form factors. The disconnected contributions are evaluated using recently developed methods that include deflation of the lower eigenstates, allowing us to extract the isoscalar, strange and charm axial form factors. We find that the disconnected quark loop contributions are non-zero and particularly large for the induced pseudo-scalar form factor.


title: The isoperimetric problem in the 2-dimensional Finsler space forms with k = 0. I

abstract: In this paper, the isoperimetric problem in the 2-dimensional Finsler space form $(F_B,B^2(1))$ with k = 0 by using the Busemann-Hausdorff area is investigated. We prove that the circle centered the origin achieves the local maximum area of the isoperimetric problem.


title: A covariant Lagrangian for stable nonsingular bounce

abstract: The nonsingular bounce models usually suffer from the ghost or gradient instabilities, as has been proved recently. In this paper, we propose a covariant effective theory for stable nonsingular bounce, which has the quadratic order of the second order derivative of the field $\phi$ but the background set only by $P(\phi,X)$. With it, we explicitly construct a fully stable nonsingular bounce model for the ekpyrotic scenario.


title: The effect of the environment on the structure, morphology and star-formation history of intermediate-redshift galaxies

abstract: With the aim of understanding the effect of the environment on the star formation history and morphological transformation of galaxies, we present a detailed analysis of the colour, morphology and internal structure of cluster and field galaxies at $0.4 \le z \le 0.8$. We use {\em HST} data for over 500 galaxies from the ESO Distant Cluster Survey (EDisCS) to quantify how the galaxies' light distribution deviate from symmetric smooth profiles. We visually inspect the galaxies' images to identify the likely causes for such deviations. We find that the residual flux fraction ($RFF$), which measures the fractional contribution to the galaxy light of the residuals left after subtracting a symmetric and smooth model, is very sensitive to the degree of structural disturbance but not the causes of such disturbance. On the other hand, the asymmetry of these residuals ($A_{\rm res}$) is more sensitive to the causes of the disturbance, with merging galaxies having the highest values of $A_{\rm res}$. Using these quantitative parameters we find that, at a fixed morphology, cluster and field galaxies show statistically similar degrees of disturbance. However, there is a higher fraction of symmetric and passive spirals in the cluster than in the field. These galaxies have smoother light distributions than their star-forming counterparts. We also find that while almost all field and cluster S0s appear undisturbed, there is a relatively small population of star-forming S0s in clusters but not in the field. These findings are consistent with relatively gentle environmental processes acting on galaxies infalling onto clusters.


title: The K-Pg event as a key to bat evolution

abstract: Bats are unique mammals. This note discusses some questions regarding bat evolution including why they are nocturnal and why they can echolocate. It is hypothesized that echolocation was necessary for bats to survive the period of limited visibility that followed the Cretaceous-Paleogene (K-Pg) extinction event.


title: A broadband self-consistent modelling of the X-ray spectrum of 4U 1626--67

abstract: The accretion-powered X-ray pulsar 4U 1626--67 is one of the few highly magnetized pulsars that accretes through Roche-lobe overflow from a low-mass companion. The characteristics of its broadband spectrum are similar to those of X-ray pulsars hosted in a high-mass X-ray binary systems, with a broad resonant cyclotron scattering feature (CRSF) at $\sim$37 keV. In this work, we examine the pulse-resolved and the pulse-averaged broadband spectrum using data from NuSTAR and Swift. We use the Becker & Wolff model of bulk+thermal Comptonization to infer key physical parameters of the accretion column flow and a broadband model for the disk reflected spectrum. In the softer X-ray band, we need to add a soft black-body component with $kT_{\rm bb}$ $\sim$0.5 keV, whose characteristics indicate a possible origin from the neutron star surface. Residuals suggest that the shape of the cyclotron line could be more satisfactorily fitted using a narrow core and broader wings and, at higher energies, a second harmonic could be present at $\sim$61 keV.


title: A complex network analysis of ethnic conflicts and human rights violations

abstract: News reports in media contain records of a wide range of socio-economic and political events in time. Using a publicly available, large digital database of news records, and aggregating them over time, we study the network of ethnic conflicts and human rights violations. Complex network analyses of the events and the involved actors provide important insights on the engaging actors, groups, establishments and sometimes nations, pointing at their long range effect over space and time. We find power law decays in distributions of actor mentions, co-actor mentions and degrees and dominance of influential actors and groups. Most influential actors or groups form a giant connected component which grows in time, and is expected to encompass all actors globally in the long run. We demonstrate how targeted removal of actors may help stop spreading unruly events. We study the cause-effect relation between types of events, and our quantitative analysis confirm that ethnic conflicts lead to human rights violations, while it does not support the converse.


title: Competition between Chaotic and Non-Chaotic Phases in a Quadratically Coupled Sachdev-Ye-Kitaev Model

abstract: The Sachdev-Ye-Kitaev (SYK) model is a concrete solvable model to study non-Fermi liquid properties, holographic duality and maximally chaotic behavior. In this work, we consider a generalization of the SYK model that contains two SYK models with different number of Majorana modes coupled by quadratic terms. This model is also solvable, and the solution shows a zero-temperature quantum phase transition between two non-Fermi liquid chaotic phases. This phase transition is driven by tuning the ratio of two mode numbers, and a Fermi liquid non-chaotic phase sits at the critical point with equal mode number. At finite temperature, the Fermi liquid phase expands to a finite regime. More intriguingly, a different non-Fermi liquid phase emerges at finite temperature. We characterize the phase diagram in term of the spectral function, the Lyapunov exponent and the entropy. Our results illustrate a concrete example of quantum phase transition and critical regime between two non-Fermi liquid phases.


title: Inferring gene regulatory networks from single-cell data: a mechanistic approach

abstract: The recent development of single-cell transcriptomics has enabled gene expression to be measured in individual cells instead of being population-averaged. Despite this considerable precision improvement, inferring regulatory networks remains challenging because stochasticity now proves to play a fundamental role in gene expression. In particular, mRNA synthesis is now acknowledged to occur in a highly bursty manner. We propose to view the inference problem as a fitting procedure for a mechanistic gene network model that is inherently stochastic and takes not only protein, but also mRNA levels into account. We first explain how to build and simulate this network model based upon the coupling of genes that are described as piecewise-deterministic Markov processes. Our model is modular and can be used to implement various biochemical hypotheses including causal interactions between genes. However, a naive fitting procedure would be intractable. By performing a relevant approximation of the stationary distribution, we derive a tractable procedure that corresponds to a statistical hidden Markov model with interpretable parameters. This approximation turns out to be extremely close to the theoretical distribution in the case of a simple toggle-switch, and we show that it can indeed fit real single-cell data. As a first step toward inference, our approach was applied to a number of simple two-gene networks simulated in silico from the mechanistic model and satisfactorily recovered the original networks. Our results demonstrate that functional interactions between genes can be inferred from the distribution of a mechanistic, dynamical stochastic model that is able to describe gene expression in individual cells. This approach seems promising in relation to the current explosion of single-cell expression data.


title: Six-vertex model and non-linear differential equations I. Spectral problem

abstract: In this work we relate the spectral problem of the toroidal six-vertex model's transfer matrix with the theory of integrable non-linear differential equations. More precisely, we establish an analogy between the Classical Inverse Scattering Method and previously proposed functional equations originating from the Yang-Baxter algebra. The latter equations are then regarded as an Auxiliary Linear Problem allowing us to show that the six-vertex model's spectrum solves Riccati-type non-linear differential equations. Generating functions of conserved quantities are expressed in terms of determinants and we also discuss a relation between our Riccati equations and a stationary Schr\"odinger equation.


title: On the solutions of the critical Lane-Emden equation in higher space dimensions

abstract: In this paper we study solutions of the critical Lane-Emden equation in higher space dimensions. We show that after certain transformations the general solution can be written in terms of elliptic functions. We restrict ourselves to real solutions which can be used in physical applications.


title: Accelerating quantum instanton calculations of the kinetic isotope effects

abstract: Path integral implementation of the quantum instanton approximation currently belongs among the most accurate methods for computing quantum rate constants and kinetic isotope effects, but its use has been limited due to the rather high computational cost. Here we demonstrate that the efficiency of quantum instanton calculations of the kinetic isotope effects can be increased by orders of magnitude by combining two approaches: The convergence to the quantum limit is accelerated by employing high-order path integral factorizations of the Boltzmann operator, while the statistical convergence is improved by implementing virial estimators for relevant quantities. After deriving several new virial estimators for the high-order factorization and evaluating the resulting increase in efficiency, using $\mathrm{\cdot H_{\alpha}+H_{\beta}H_{\gamma}\rightarrow H_{\alpha}H_{\beta}+\cdot H_{\gamma }}$ reaction as an example, we apply the proposed method to obtain several kinetic isotope effects on $\mathrm{CH_{4}+\cdot H\rightleftharpoons\cdot CH_{3}+H_{2}}$ forward and backward reactions.


title: Accelerating equilibrium isotope effect calculations: I. Stochastic thermodynamic integration with respect to mass

abstract: Accurate path integral Monte Carlo or molecular dynamics calculations of isotope effects have until recently been expensive because of the necessity to reduce three types of errors present in such calculations: statistical errors due to sampling, path integral discretization errors, and thermodynamic integration errors. While the statistical errors can be reduced with virial estimators and path integral discretization errors with high-order factorization of the Boltzmann operator, here we propose a method for accelerating isotope effect calculations by eliminating the integration error. We show that the integration error can be removed entirely by changing particle masses stochastically during the calculation and by using a piecewise linear umbrella biasing potential. Moreover, we demonstrate numerically that this approach does not increase the statistical error. The resulting acceleration of isotope effect calculations is demonstrated on a model harmonic system and on deuterated species of methane.


title: Nonconvex Generalization of Alternating Direction Method of Multipliers for Nonlinear Equality Constrained Problems

abstract: The classic Alternating Direction Method of Multipliers (ADMM) is a popular framework to solve linear-equality constrained problems. In this paper, we extend the ADMM naturally to nonlinear equality-constrained problems, called neADMM. The difficulty of neADMM is to solve nonconvex subproblems. We provide globally optimal solutions to them in two important applications. Experiments on synthetic and real-world datasets demonstrate excellent performance and scalability of our proposed neADMM over existing state-of-the-start methods.


title: Classifying Radio Galaxies with Convolutional Neural Network

abstract: We present the application of deep machine learning technique to classify radio images of extended sources on a morphological basis using convolutional neural networks. In this study, we have taken the case of Fanaroff-Riley (FR) class of radio galaxies as well as radio galaxies with bent-tailed morphology. We have used archival data from the Very Large Array (VLA) - Faint Images of the Radio Sky at Twenty Centimeters (FIRST) survey and existing visually classified samples available in literature to train a neural network for morphological classification of these categories of radio sources. Our training sample size for each of these categories is approximately 200 sources, which has been augmented by rotated versions of the same. Our study shows that convolutional neural networks can classify images of the FRI and FRII and bent-tailed radio galaxies with high accuracy (maximum precision at 95%) using well-defined samples and fusion classifier, which combines the results of binary classifications, while allowing for a mechanism to find sources with unusual morphologies. The individual precision is highest for bent-tailed radio galaxies at 95% and is 91% and 75% for the FRI and FRII classes, respectively, whereas the recall is highest for FRI and FRIIs at 91% each, while bent-tailed class has a recall of 79%. These results show that our results are comparable to that of manual classification while being much faster. Finally, we discuss the computational and data-related challenges associated with morphological classification of radio galaxies with convolutional neural networks.


title: A Distributed Learning Dynamics in Social Groups

abstract: We study a distributed learning process observed in human groups and other social animals. This learning process appears in settings in which each individual in a group is trying to decide over time, in a distributed manner, which option to select among a shared set of options. Specifically, we consider a stochastic dynamics in a group in which every individual selects an option in the following two-step process: (1) select a random individual and observe the option that individual chose in the previous time step, and (2) adopt that option if its stochastic quality was good at that time step. Various instantiations of such distributed learning appear in nature, and have also been studied in the social science literature. From the perspective of an individual, an attractive feature of this learning process is that it is a simple heuristic that requires extremely limited computational capacities. But what does it mean for the group -- could such a simple, distributed and essentially memoryless process lead the group as a whole to perform optimally? We show that the answer to this question is yes -- this distributed learning is highly effective at identifying the best option and is close to optimal for the group overall. Our analysis also gives quantitative bounds that show fast convergence of these stochastic dynamics. Prior to our work the only theoretical work related to such learning dynamics has been either in deterministic special cases or in the asymptotic setting. Finally, we observe that our infinite population dynamics is a stochastic variant of the classic multiplicative weights update (MWU) method. Consequently, we arrive at the following interesting converse: the learning dynamics on a finite population considered here can be viewed as a novel distributed and low-memory implementation of the classic MWU method.


title: 3D Placement of an Unmanned Aerial Vehicle Base Station (UAV-BS) for Energy-Efficient Maximal Coverage

abstract: Unmanned Aerial Vehicle mounted base stations (UAV-BSs) can provide wireless services in a variety of scenarios. In this letter, we propose an optimal placement algorithm for UAV-BSs that maximizes the number of covered users using the minimum transmit power. We decouple the UAV-BS deployment problem in the vertical and horizontal dimensions without any loss of optimality. Furthermore, we model the UAV-BS deployment in the horizontal dimension as a circle placement problem and a smallest enclosing circle problem. Simulations are conducted to evaluate the performance of the proposed method for different spatial distributions of the users.


title: Low spin wave damping in the insulating chiral magnet Cu$_{2}$OSeO$_{3}$

abstract: Chiral magnets with topologically nontrivial spin order such as Skyrmions have generated enormous interest in both fundamental and applied sciences. We report broadband microwave spectroscopy performed on the insulating chiral ferrimagnet Cu$_{2}$OSeO$_{3}$. For the damping of magnetization dynamics we find a remarkably small Gilbert damping parameter of about $1\times10^{-4}$ at 5 K. This value is only a factor of 4 larger than the one reported for the best insulating ferrimagnet yttrium iron garnet. We detect a series of sharp resonances and attribute them to confined spin waves in the mm-sized samples. Considering the small damping, insulating chiral magnets turn out to be promising candidates when exploring non-collinear spin structures for high frequency applications.


title: Explicitly correlated Helium wave function and hyperspherical coordinates

abstract: Wave functions of a new functional kind have been proposed for Helium-like atoms in this work . These functions explicitly depend on interelectronic and hyperspherical coordinates. The best ground state energy for the Helium atom $ -2.903724376677 a.u.$ has been calculated with variational method with basis set of simple functions with a single exponential parameter. To the author's knowledge, this is the best result with use of hyperspherical coordinates so far. Comparable result has been obtained for the hydrogen anion. For Helium atom, our best wave functions matched the Kato cusp conditions within the accuracy below $6.10^{-4} $. An important feature of proposed wave functions is the inclusion of negative powers of $R=\sqrt(r^{2}_{1}+r^{2}_{2})$ in combination with positive powers of $r_{12}$ into the wave function. We showed that this is necessary condition for proposed wave function to be a formal solution of Schr\"odinger equation.

categories: physics.chem-ph physics.atom-ph physics.comp-ph quant-ph


title: A notion of minor-based matroid connectivity

abstract: For a matroid $N$, a matroid $M$ is $N$-connected if every two elements of $M$ are in an $N$-minor together. Thus a matroid is connected if and only if it is $U_{1,2}$-connected. This paper proves that $U_{1,2}$ is the only connected matroid $N$ such that if $M$ is $N$-connected with $|E(M)| > |E(N)|$, then $M \backslash e$ or $M / e$ is $N$-connected for all elements $e$. Moreover, we show that $U_{1,2}$ and $M(\mathcal{W}_2)$ are the only connected matroids $N$ such that, whenever a matroid has an $N$-minor using $\{e,f\}$ and an $N$-minor using $\{f,g\}$, it also has an $N$-minor using $\{e,g\}$. Finally, we show that $M$ is $U_{0,1} \oplus U_{1,1}$-connected if and only if every clonal class of $M$ is trivial.


title: Learning Deep Networks from Noisy Labels with Dropout Regularization

abstract: Large datasets often have unreliable labels-such as those obtained from Amazon's Mechanical Turk or social media platforms-and classifiers trained on mislabeled datasets often exhibit poor performance. We present a simple, effective technique for accounting for label noise when training deep neural networks. We augment a standard deep network with a softmax layer that models the label noise statistics. Then, we train the deep network and noise model jointly via end-to-end stochastic gradient descent on the (perhaps mislabeled) dataset. The augmented model is overdetermined, so in order to encourage the learning of a non-trivial noise model, we apply dropout regularization to the weights of the noise model during training. Numerical experiments on noisy versions of the CIFAR-10 and MNIST datasets show that the proposed dropout technique outperforms state-of-the-art methods.


title: Compressive Estimation of a Stochastic Process with Unknown Autocorrelation Function

abstract: In this paper, we study the prediction of a circularly symmetric zero-mean stationary Gaussian process from a window of observations consisting of finitely many samples. This is a prevalent problem in a wide range of applications in communication theory and signal processing. Due to stationarity, when the autocorrelation function or equivalently the power spectral density (PSD) of the process is available, the Minimum Mean Squared Error (MMSE) predictor is readily obtained. In particular, it is given by a linear operator that depends on autocorrelation of the process as well as the noise power in the observed samples. The prediction becomes, however, quite challenging when the PSD of the process is unknown. In this paper, we propose a blind predictor that does not require the a priori knowledge of the PSD of the process and compare its performance with that of an MMSE predictor that has a full knowledge of the PSD. To design such a blind predictor, we use the random spectral representation of a stationary Gaussian process. We apply the well-known atomic-norm minimization technique to the observed samples to obtain a discrete quantization of the underlying random spectrum, which we use to predict the process. Our simulation results show that this estimator has a good performance comparable with that of the MMSE estimator.


title: Creation of a Bose-condensed gas of rubidium 87 by laser cooling

abstract: We demonstrate direct laser cooling of a gas of rubidium 87 atoms to quantum degeneracy. The method does not involve evaporative cooling, is fast, and induces little atom loss. The atoms are trapped in a two-dimensional optical lattice that enables cycles of cloud compression to increase the density, followed by degenerate Raman sideband cooling to decrease the temperature. Light-induced loss at high atomic density is substantially reduced by using far red detuned optical pumping light. Starting with 2000 atoms, we prepare 1400 atoms in 300 ms at quantum degeneracy, as confirmed by the appearance of a bimodal velocity distribution as the system crosses over from a classical gas to a Bose-condensed, interacting one-dimensional gas with a macroscopic population of the quantum ground state. The method should be broadly applicable to many bosonic and fermionic species, and to systems where evaporative cooling is not possible.


title: Adjustments to Computer Models via Projected Kernel Calibration

abstract: Identification of model parameters in computer simulations is an important topic in computer experiments. We propose a new method, called the projected kernel calibration method, to estimate these model parameters. The proposed method is proven to be asymptotic normal and semi-parametric efficient. As a frequentist method, the proposed method is as efficient as the $L_2$ calibration method proposed by Tuo and Wu [Ann. Statist. 43 (2015) 2331-2352]. On the other hand, the proposed method has a natural Bayesian version, which the $L_2$ method does not have. This Bayesian version allows users to calculate the credible region of the calibration parameters without using a large sample approximation. We also show that, the inconsistency problem of the calibration method proposed by Kennedy and O'Hagan [J. R. Stat. Soc. Ser. B. Stat. Methodol. 63 (2001) 425-464] can be rectified by a simple modification of the kernel matrix.


title: Wright meets Markowitz: How standard portfolio theory changes when assets are technologies following experience curves

abstract: We consider how to optimally allocate investments in a portfolio of competing technologies using the standard mean-variance framework of portfolio theory. We assume that technologies follow the empirically observed relationship known as Wright's law, also called a "learning curve" or "experience curve", which postulates that costs drop as cumulative production increases. This introduces a positive feedback between cost and investment that complicates the portfolio problem, leading to multiple local optima, and causing a trade-off between concentrating investments in one project to spur rapid progress vs. diversifying over many projects to hedge against failure. We study the two-technology case and characterize the optimal diversification in terms of progress rates, variability, initial costs, initial experience, risk aversion, discount rate and total demand. The efficient frontier framework is used to visualize technology portfolios and show how feedback results in nonlinear distortions of the feasible set. For the two-period case, in which learning and uncertainty interact with discounting, we compare different scenarios and find that the discount rate plays a critical role.


title: Expected seismicity and the seismic noise environment of Europa

abstract: Seismic data will be a vital geophysical constraint on internal structure of Europa if we land instruments on the surface. Quantifying expected seismic activity on Europa both in terms of large, recognizable signals and ambient background noise is important for understanding dynamics of the moon, as well as interpretation of potential future data. Seismic energy sources will likely include cracking in the ice shell and turbulent motion in the oceans. We define a range of models of seismic activity in Europa's ice shell by assuming each model follows a Gutenberg-Richter relationship with varying parameters. A range of cumulative seismic moment release between $10^{16}$ and $10^{18}$ Nm/yr is defined by scaling tidal dissipation energy to tectonic events on the Earth's moon. Random catalogs are generated and used to create synthetic continuous noise records through numerical wave propagation in thermodynamically self-consistent models of the interior structure of Europa. Spectral characteristics of the noise are calculated by determining probabilistic power spectral densities of the synthetic records. While the range of seismicity models predicts noise levels that vary by 80 dB, we show that most noise estimates are below the self-noise floor of high-frequency geophones, but may be recorded by more sensitive instruments. The largest expected signals exceed background noise by $\sim$50 dB. Noise records may allow for constraints on interior structure through autocorrelation. Models of seismic noise generated by pressure variations at the base of the ice shell due to turbulent motions in the subsurface ocean may also generate observable seismic noise.


title: Characterizing Lie groups by controlling their zero-dimensional subgroups

abstract: We provide characterizations of Lie groups as compact-like groups in which all closed zero-dimensional metric (compact) subgroups are discrete. The "compact-like" properties we consider include (local) compactness, (local) omega-boundedness, (local) countable compactness, (local) precompactness, (local) minimality and sequential completeness. Below is a sample of our characterizations: (i) A topological group is a Lie group if and only if it is locally compact and has no infinite compact metric zero-dimensional subgroups. (ii) An abelian topological group G is a Lie group if and only if G is locally minimal, locally precompact and all closed metric zero-dimensional subgroups of G are discrete. (iii) An abelian topological group is a compact Lie group if and only if it is minimal and has no infinite closed metric zero-dimensional subgroups. (iv) An infinite topological group is a compact Lie group if and only if it is sequentially complete, precompact, locally minimal, contains a non-empty open connected subset and all its compact metric zero-dimensional subgroups are finite.


title: PT -symmetric rational Calogero model with balanced loss and gain

abstract: A two body rational Calogero model with balanced loss and gain is investigated. The system yields a Hamiltonian which is symmetric under the combined operation of parity (P) and time reversal (T ) symmetry. It is shown that the system is integrable and exact, stable classical solutions are obtained for particular ranges of the parameters. The corresponding quantum system admits bound state solutions for exactly the same ranges of the parameters for which the classical solutions are stable. The eigen spectra of the system is presented with a discussion on the normalization of the wave functions in proper Stokes wedges. Finally, the Calogero model with balanced loss and gain is studied classically, when the pair-wise harmonic interaction term is replaced by a common confining harmonic potential. The system admits stable solutions for particular ranges of the parameters. However, the integrability and/or exact solvability of the system is obscure due to the presence of the loss and gain terms. The perturbative solutions are obtained and are compared with the numerical results.


title: Rapid Mixing of Local Graph Dynamics

abstract: Graph dynamics arise naturally in many contexts. For instance in peer-to-peer networks, a participating peer may replace an existing connection with one neighbour by a new connection with a neighbour's neighbour. Several such local rewiring rules have been proposed to ensure that peer-to-peer networks achieve good connectivity properties (e.g. high expansion) in equilibrium. However it has remained an open question whether there existed such rules that also led to fast convergence to equilibrium. In this work we provide an affirmative answer: We exhibit a local rewiring rule that converges to equilibrium after each participating node has undergone only a number of rewirings that is poly-logarithmic in the system size. The proof involves consideration of the whole isoperimetric profile of the graph, and may be of independent interest.


title: Deep Projective 3D Semantic Segmentation

abstract: Semantic segmentation of 3D point clouds is a challenging problem with numerous real-world applications. While deep learning has revolutionized the field of image semantic segmentation, its impact on point cloud data has been limited so far. Recent attempts, based on 3D deep learning approaches (3D-CNNs), have achieved below-expected results. Such methods require voxelizations of the underlying point cloud data, leading to decreased spatial resolution and increased memory consumption. Additionally, 3D-CNNs greatly suffer from the limited availability of annotated datasets. In this paper, we propose an alternative framework that avoids the limitations of 3D-CNNs. Instead of directly solving the problem in 3D, we first project the point cloud onto a set of synthetic 2D-images. These images are then used as input to a 2D-CNN, designed for semantic segmentation. Finally, the obtained prediction scores are re-projected to the point cloud to obtain the segmentation results. We further investigate the impact of multiple modalities, such as color, depth and surface normals, in a multi-stream network architecture. Experiments are performed on the recent Semantic3D dataset. Our approach sets a new state-of-the-art by achieving a relative gain of 7.9 %, compared to the previous best approach.


title: Incentive Mechanism Design for Cache-Assisted D2D Communications: A Mobility-Aware Approach

abstract: Caching popular contents at mobile devices, assisted by device-to-device (D2D) communications, is considered as a promising technique for mobile content delivery. It can effectively reduce backhaul traffic and service cost, as well as improving the spectrum efficiency. However, due to the selfishness of mobile users, incentive mechanisms will be needed to motivate device caching. In this paper, we investigate incentive mechanism design in cache-assisted D2D networks, taking advantage of the user mobility information. An inter-contact model is adopted to capture the average time between two consecutive contacts of each device pair. A Stackelberg game is formulated, where each user plays as a follower aiming at maximizing its own utility and the mobile network operator (MNO) plays as a leader aiming at minimizing the cost. Assuming that user responses can be predicted by the MNO, a cost minimization problem is formulated. Since this problem is NP-hard, we reformulate it as a non-negative submodular maximization problem and develop $(\frac{1}{4+\epsilon})$-approximation local search algorithm to solve it. In the simulation, we demonstrate that the local search algorithm provides near optimal performance. By comparing with other caching strategies, we validate the effectiveness of the proposed incentive-based mobility-aware caching strategy.


title: Analysis of Channel-Based User Authentication by Key-Less and Key-Based Approaches

abstract: User authentication (UA) supports the receiver in deciding whether a message comes from the claimed transmitter or from an impersonating attacker. In cryptographic approaches messages are signed with either an asymmetric or symmetric key, and a source of randomness is required to generate the key. In physical layer authentication (PLA) instead the receiver checks if received messages presumably coming from the same source undergo the same channel. We compare these solutions by considering the physical-layer channel features as randomness source for generating the key, thus allowing an immediate comparison with PLA (that already uses these features). For the symmetric-key approach we use secret key agreement, while for asymmetric-key the channel is used as entropy source at the transmitter. We focus on the asymptotic case of an infinite number of independent and identically distributed channel realizations, showing the correctness of all schemes and analyzing the secure authentication rate, that dictates the rate at which the probability that UA security is broken goes to zero as the number of used channel resources (to generate the key or for PLA) goes to infinity. Both passive and active attacks are considered and by numerical results we compare the various systems.


title: Rapid X-ray variability properties during the unusual very hard state in neutron-star low-mass X-ray binaries

abstract: Here we study the rapid X-ray variability (using XMM-Newton observations) of three neutron-star low-mass X-ray binaries (1RXS J180408.9-342058, EXO 1745-248, and IGR J18245-2452) during their recently proposed very hard spectral state (Parikh et al. 2017). All our systems exhibit a strong to very strong noise component in their power density spectra (rms amplitudes ranging from 34% to 102%) with very low characteristic frequencies (as low as 0.01 Hz). These properties are more extreme than what is commonly observed in the canonical hard state of neutron-star low-mass X-ray binaries observed at X-ray luminosities similar to those we observe from our sources. This suggests that indeed the very hard state is a distinct spectral-timing state from the hard state, although we argue that the variability behaviour of IGR J18245-2452 is very extreme and possibly this source was in a very unusual state. We also compare our results with the rapid X-ray variability of the accreting millisecond X-ray pulsars IGR J00291+5934 and Swift J0911.9-6452 (also using XMM-Newton data) for which previously similar variability phenomena were observed. Although their energy spectra (as observed using the Swift X-ray telescope) were not necessarily as hard (i.e., for Swift J0911.9-6452) as for our other three sources, we conclude that likely both sources were also in very similar state during their XMM-Newton observations. This suggest that different sources that are found in this new state might exhibit different spectral hardness and one has to study both the spectral as well as rapid variability to identify this unusual state.


title: Evolution of tripartite entangled states in a decohering environment and their experimental protection using dynamical decoupling

abstract: We embarked upon the task of experimental protection of different classes of tripartite entangled states, namely the maximally entangled GHZ and W states and the ${\rm W \bar{W}}$ state, using dynamical decoupling. The states were created on a three-qubit NMR quantum information processor and allowed to evolve in the naturally noisy NMR environment. Tripartite entanglement was monitored at each time instant during state evolution, using negativity as an entanglement measure. It was found that the W state is most robust while the GHZ-type states are most fragile against the natural decoherence present in the NMR system. The ${\rm W \bar{W}}$ state which is in the GHZ-class, yet stores entanglement in a manner akin to the W state, surprisingly turned out to be more robust than the GHZ state. The experimental data were best modeled by considering the main noise channel to be an uncorrelated phase damping channel acting independently on each qubit, alongwith a generalized amplitude damping channel. Using dynamical decoupling, we were able to achieve a significant protection of entanglement for GHZ states. There was a marginal improvement in the state fidelity for the W state (which is already robust against natural system decoherence), while the ${\rm W \bar{W}}$ state showed a significant improvement in fidelity and protection against decoherence.


title: Comparing topologies on linearly recursive sequences

abstract: The space of linearly recursive sequences of complex numbers admits two distinguished topologies. Namely, the adic topology induced by the ideal of those sequences whose first term is $0$ and the topology induced from the Krull topology on the space of complex power series via a suitable embedding. We show that these topologies are not equivalent.


title: Two-component domain decomposition scheme with overlapping subdomains for parabolic equations

abstract: An iteration-free method of domain decomposition is considered for approximate solving a boundary value problem for a second-order parabolic equation. A standard approach to constructing domain decomposition schemes is based on a partition of unity for the domain under the consideration. Here a new general approach is proposed for constructing domain decomposition schemes with overlapping subdomains based on indicator functions of subdomains. The basic peculiarity of this method is connected with a representation of the problem operator as the sum of two operators, which are constructed for two separate subdomains with the subtraction of the operator that is associated with the intersection of the subdomains. There is developed a two-component factorized scheme, which can be treated as a generalization of the standard Alternating Direction Implicit (ADI) schemes to the case of a special three-component splitting. There are obtained conditions for the unconditional stability of regionally additive schemes constructed using indicator functions of subdomains. Numerical results are presented for a model two-dimensional problem.


title: A conjectural Peterson isomorphism in K-theory

abstract: We state a precise conjectural isomorphism between localizations of the equivariant quantum K-theory ring of a flag variety and the equivariant K-homology ring of the affine Grassmannian, in particular relating their Schubert bases and structure constants. This generalizes Peterson's isomorphism in (co)homology. We prove a formula for the Pontryagin structure constants in the K-homology ring, and we use it to check our conjecture in few situations.


title: Atomically flat single terminated oxide substrate surfaces

abstract: To achieve high quality epitaxial thin films and heterostructures of transition metal oxides with atomically controlled interfaces, one critical requirement is the use of atomically flat single terminated oxide substrates since the atomic arrangements and the reaction chemistry of the topmost surface layer of substrates determine the growth and consequent properties of the overlying films. Achieving the atomically flat and chemically single terminated surface state of commercially available substrates, however, requires judicious efforts because the surface of as-received substrates is of chemically mixed nature and also often polar. In this review, we summarize the surface treatment procedures to accomplish atomically flat surfaces with single terminating layer for various metal oxide substrates. We particularly focus on the substrates with lattice constant ranging from 4.00 to 3.70 angstrom, as the lattice constant of most perovskite materials falls into this range. For materials outside the range, one can utilize the substrates to induce compressive or tensile strain on the films and explore new states not available in bulk. The substrates covered in this review, which have been chosen with commercial availability and, most importantly, experimental practicality as a criterion, are KTaO3, REScO3 (RE = Rare-earth elements), SrTiO3, La0.18Sr0.82Al0.59Ta0.41O3 (LSAT), NdGaO3, LaAlO3, SrLaAlO4, and YAlO3. Analyzing all the established procedures, we conclude that atomically flat surfaces with selective A- or B-site single termination would be obtained for most commercially available oxide substrates. We further note that this topmost surface layer selectivity would provide an additional degree of freedom in searching for unforeseen emergent phenomena and functional applications in epitaxial oxide thin films and heterostructures with atomically controlled interfaces.


title: Optimal properties of the canonical tight probabilistic frame

abstract: A probabilistic frame is a Borel probability measure with finite second moment whose support spans $\mathbb{R}^d$. A Parseval probabilistic frame is one for which the associated matrix of the second moments is the identity matrix in $\mathbb{R}^d$. Each probabilistic frame is canonically associated to a Parseval probabilistic frame. In this paper, we show that this canonical Parseval probabilistic frame is the closest Parseval probabilistic frame to a given probabilistic frame in the $2-$Wasserstein distance. Our proof is based on two main ingredients. On the one hand, we show that a probabilistic frame can be approximated in the $2-$Wasserstein metric with (compactly supported) finite frames whose bounds can be controlled. On the other hand, we establish some fine continuity properties of the function that maps a probabilistic frame to its canonical Parseval probabilistic frame. Our results generalize similar ones for finite frames and their associated Parseval frames.


title: A higher-order Skyrme model

abstract: We propose a higher-order Skyrme model with derivative terms of eighth, tenth and twelfth order. Our construction yields simple and easy-to-interpret higher-order Lagrangians. We first show that a Skyrmion with higher-order terms proposed by Marleau has an instability in the form of a baby-Skyrmion string, while the static energies of our construction are positive definite, implying stability against time-independent perturbations. However, we also find that the Hamiltonians of our construction possess two kinds of dynamical instabilities, which may indicate the instability with respect to time-dependent perturbations. Different from the well-known Ostrogradsky instability, the instabilities that we find are intrinsically of nonlinear nature and also due to the fact that even powers of the inverse metric gives a ghost-like higher-order kinetic-like term. The vacuum state is, however, stable. Finally, we show that at sufficiently low energies, our Hamiltonians in the simplest cases, are stable against time-dependent perturbations.


title: Frequentist Consistency of Variational Bayes

abstract: A key challenge for modern Bayesian statistics is how to perform scalable inference of posterior distributions. To address this challenge, variational Bayes (VB) methods have emerged as a popular alternative to the classical Markov chain Monte Carlo (MCMC) methods. VB methods tend to be faster while achieving comparable predictive performance. However, there are few theoretical results around VB. In this paper, we establish frequentist consistency and asymptotic normality of VB methods. Specifically, we connect VB methods to point estimates based on variational approximations, called frequentist variational approximations, and we use the connection to prove a variational Bernstein-von Mises theorem. The theorem leverages the theoretical characterizations of frequentist variational approximations to understand asymptotic properties of VB. In summary, we prove that (1) the VB posterior converges to the Kullback-Leibler (KL) minimizer of a normal distribution, centered at the truth and (2) the corresponding variational expectation of the parameter is consistent and asymptotically normal. As applications of the theorem, we derive asymptotic properties of VB posteriors in Bayesian mixture models, Bayesian generalized linear mixed models, and Bayesian stochastic block models. We conduct a simulation study to illustrate these theoretical results.


title: Test of the local form of higher-spin equations via AdS/CFT

abstract: The local form of higher-spin equations found recently to the second order [1] is shown to properly reproduce the anticipated $AdS/CFT$ correlators for appropriate boundary conditions. It is argued that consistent $AdS/CFT$ holography for the parity-broken boundary models needs a nontrivial modification of the bosonic truncation of the original higher-spin theory with the doubled number of fields, as well as a nonlinear deformation of the boundary conditions in the higher orders.


title: Improved Computation of Involutive Bases

abstract: In this paper, we describe improved algorithms to compute Janet and Pommaret bases. To this end, based on the method proposed by Moller et al., we present a more efficient variant of Gerdt's algorithm (than the algorithm presented by Gerdt-Hashemi-M.Alizadeh) to compute minimal involutive bases. Further, by using the involutive version of Hilbert driven technique, along with the new variant of Gerdt's algorithm, we modify the algorithm, given by Seiler, to compute a linear change of coordinates for a given homogeneous ideal so that the new ideal (after performing this change) possesses a finite Pommaret basis. All the proposed algorithms have been implemented in Maple and their efficiency is discussed via a set of benchmark polynomials.


title: The Interactive Sum Choice Number of Trees

abstract: We study the interactive sum choice number, a game coloring parameter introduced by Bonamy and Meeks, and obtain a recursive formula for the interactive sum choice number of forests. This formula coincides with a formula for the slow coloring cost of forests, a parameter introduced by Mahoney, Puleo, and West, and shows that these parameters are equal on forests. This answers a question of Bonamy and Meeks.


title: Alteration of cerebrovascular haemodynamic patterns due to atrial fibrillation: an in silico investigation

abstract: There has recently been growing evidence that atrial fibrillation (AF), the most common cardiac arrhythmia, is independently associated with the risk of dementia. This represents a very recent frontier with high social impact for the number of individuals involved and for the expected increase in AF incidence in the next 40 years. Although a number of potential haemodynamic processes, such as microembolisms, altered cerebral blood flow, hypoperfusion and microbleeds, arise as connecting links between the two pathologies, the causal mechanisms are far from clear. An in silico approach is proposed that combines in sequence two lumped-parameter schemes, for the cardiovascular system and the cerebral circulation. The systemic arterial pressure is obtained from the cardiovascular system and used as the input for the cerebral circulation, with the aim of studying the role of AF on the cerebral haemodynamics with respect to normal sinus rhythm (NSR), over a 5000 beat recording. In particular, the alteration of the haemodynamic (pressure and flowrate) patterns in the microcirculation during AF is analysed by means of different statistical tools, from correlation coefficients to autocorrelation functions, crossing times, extreme values analysis and multivariate linear regression models. A remarkable signal alteration, such as a reduction in signal correlation (NSR, about 3 s; AF, less than 1 s) and increased probability (up to three to four times higher in AF than in NSR) of extreme value events, emerges for the peripheral brain circulation. The described scenario offers a number of plausible cause-effect mechanisms that might explain the occurrence of critical events and the haemodynamic links relating to AF and dementia.


title: Four-Loop Nonplanar Cusp Anomalous Dimension in N=4 Supersymmetric Yang-Mills Theory

abstract: The light-like cusp anomalous dimension is a universal function that controls infrared divergences in quite general quantum field theories. In the maximally supersymmetric Yang-Mills theory this function is fixed fully by integrability to the three-loop order. At four loops a non-planar correction appears which we obtain for the first time from a numerical computation of the Sudakov form factor. Key ingredients are widely applicable methods to control the number-theoretic aspects of the appearing integrals. Our result shows explicitly that quadratic Casimir scaling breaks down at four loops.