title: Repulsive polarons in alkaline-earth(-like) atoms across an orbital Feshbach resonance

abstract: We characterize properties of the so-called repulsive polaron across the recently discovered orbital Feshbach resonance in alkaline-earth(-like) atoms. Being a metastable quasiparticle excitation at the positive energy, the repulsive polaron is induced by the interaction between an impurity atom and a Fermi sea. By analyzing in detail the energy, the polaron residue, the effective mass, and the decay rate of the repulsive polaron, we reveal interesting features that are intimately related to the two-channel nature of the orbital Feshbach resonance. In particular, we find that the life time of the repulsive polaron is non-monotonic in the Zeeman-field detuning bewteen the two channels, and has a maximum on the BEC-side of the resonance. Further, by considering the stability of a mixture of the impurity and the majority atoms against phase separation, we show that the itin- erant ferromagnetism may exist near the orbital Feshbach resonance at appropriate densities. Our results can be readily probed experimentally, and have interesting implications for the observation of itinerant ferromagnetism near an orbital Feshbach resonance.


title: Measurement of the Drell--Yan triple-differential cross section in $pp$ collisions at $\sqrt{s} = 8$ TeV

abstract: This paper presents a measurement of the triple-differential cross section for the Drell--Yan process $Z/\gamma^*\rightarrow \ell^+\ell^-$ where $\ell$ is an electron or a muon. The measurement is performed for invariant masses of the lepton pairs, $m_{\ell\ell}$, between $46$ and $200$ GeV using a sample of $20.2$ fb$^{-1}$ of $pp$ collisions data at a centre-of-mass energy of $\sqrt{s}=8$ TeV collected by the ATLAS detector at the LHC in 2012. The data are presented in bins of invariant mass, absolute dilepton rapidity, $|y_{\ell\ell}|$, and the angular variable $\cos\theta^{*}$ between the outgoing lepton and the incoming quark in the Collins--Soper frame. The measurements are performed in the range $|y_{\ell\ell}|<2.4$ in the muon channel, and extended to $|y_{\ell\ell}|<3.6$ in the electron channel. The cross sections are used to determine the $Z$ boson forward-backward asymmetry as a function of $|y_{\ell\ell}|$ and $m_{\ell\ell}$. The measurements achieve high-precision, below the percent level in the pole region, excluding the uncertainty in the integrated luminosity, and are in agreement with predictions. These precision data are sensitive to the parton distribution functions and the effective weak mixing angle.


title: Dynamic Portfolio Optimization with Looping Contagion Risk

abstract: In this paper we consider a utility maximization problem with defaultable stocks and looping contagion risk. We assume that the default intensity of one company depends on the stock prices of itself and other companies, and the default of the company induces immediate drops in the stock prices of the surviving companies. We prove that the value function is the unique viscosity solution of the HJB equation. We also perform some numerical tests to compare and analyse the statistical distributions of the terminal wealth of log utility and power utility based on two strategies, one using the full information of intensity process and the other a proxy constant intensity process.


title: Doubly Damped Stochastic Parallel Translations and Hessian Formulas

abstract: We study the Hessian of the solutions of time-independent Schr\"odinger equations, aiming to obtain as large a class as possible of complete Riemannian manifolds for which the estimate $C(\frac 1 t +\frac {d^2}{t^2})$ holds. For this purpose we introduce the doubly damped stochastic parallel transport equation, study them and make exponential estimates on them, deduce a second order Feynman-Kac formula and obtain the desired estimates. Our aim here is to explain the intuition, the basic techniques, and the formulas which might be useful in other studies.


title: Non-s wave superconductivity in boron-doped nanodiamond films with 0-{\pi} Josephson junction array

abstract: Superconducting transport properties of granular materials are greatly influenced by the microstructure. We show that in heavily boron-doped diamond films (HBDDF) films some sharp transport features can be manipulated by applying a magnetic field and controlled finite bias current. We demonstrate the conductivity cross-over from dirty metal to the superconducting state through an insulating peak arising at a very low current or magnetic field region and particularly pronounced negative magnetoresistance with periodic oscillatory features. The current-voltage characteristics show features of the Berezinskii-Kosterlitz-Thouless (BKT) phase transitions which verifies the two-dimensional structure in HBDDF observed recently. A zero bias conductance peak can be attributed to the Andreev bound state formed at the grain boundaries of diamond nanocrystals. The set of observations can be qualitatively explained consistently through the concept of a superconducting transition with a non-s wave order parameter in the diamond heterostructures.


title: Fingerprints of spin-orbital polarons and of their disorder in the photoemission spectra of doped Mott insulators with orbital degeneracy

abstract: We explore the effects of disordered charged defects on the electronic excitations observed in the photoemission spectra of doped transition metal oxides in the Mott insulating regime by the example of the $R_{1-x}$Ca$_x$VO$_3$ perovskites, where $R=$La,$\dots$,Lu. A fundamental characteristic of these vanadium $d^2$ compounds with partly filled $t_{2g}$ valence orbitals is the persistence of spin and orbital order up to high doping, in contrast to the loss of magnetic order in high-$T_c$ cuprates at low defect concentration. We demonstrate that the disordered electronic structure of doped Mott-Hubbard insulators can be obtained with high precision within the unrestricted Hartree-Fock approximation. In particular: (i) the atomic multiplet excitations in the inverse photoemission spectra and the various defect-related states and satellites are well reproduced, (ii) a robust Mott gap survives up to large doping, and (iii) we show that the defect states inside the Mott gap develop a soft gap at the Fermi energy. The soft defect states gap can be characterized by a shape and a scale parameter extracted from a Weibull statistical sampling of the density of states near the chemical potential. We demonstrate that charge defects trigger small spin-orbital polarons, with their internal kinetic energy responsible for the opening of the soft defect states gap. The small size of spin-orbital polarons is inferred by an analysis of the inverse participation ratio which explains the origin of the robustness of spin and orbital order. Using realistic parameters for La$_{1-x}$Ca$_x$VO$_3$, we show that its soft gap is well reproduced as well as the marginal doping dependence of the position of the chemical potential relative to the center of the lower Hubbard band.


title: Co-saliency Detection for RGBD Images Based on Multi-constraint Feature Matching and Cross Label Propagation

abstract: Co-saliency detection aims at extracting the common salient regions from an image group containing two or more relevant images. It is a newly emerging topic in computer vision community. Different from the most existing co-saliency methods focusing on RGB images, this paper proposes a novel co-saliency detection model for RGBD images, which utilizes the depth information to enhance identification of co-saliency. First, the intra saliency map for each image is generated by the single image saliency model, while the inter saliency map is calculated based on the multi-constraint feature matching, which represents the constraint relationship among multiple images. Then, the optimization scheme, namely Cross Label Propagation (CLP), is used to refine the intra and inter saliency maps in a cross way. Finally, all the original and optimized saliency maps are integrated to generate the final co-saliency result. The proposed method introduces the depth information and multi-constraint feature matching to improve the performance of co-saliency detection. Moreover, the proposed method can effectively exploit any existing single image saliency model to work well in co-saliency scenarios. Experiments on two RGBD co-saliency datasets demonstrate the effectiveness of our proposed model.


title: An overview of $\eta$ and $\eta^{\prime}$ decays at BESIII

abstract: The world's largest sample of 1.31 billion $J/\psi$ events accumulated at the BESIII detector, provides a unique opportunity to investigate $\eta$ and $\eta^\prime$ physics via two body $J/\psi$ radiative or hadronic decays. For many $\eta^\prime$ decay channels the low background data samples are up to three orders of magnitude larger than collected in any previous experiment. Here we review the most significant results on $\eta$ and $\eta^\prime$ obtained at BESIII so far. The analyses range from detailed studies of the common decays dynamics, observations of new radiative and Dalitz decays, and search for rare/forbidden decays with sensitivity up to $\cal{B}\sim 10^{-5}$. Finally, prospects of the forthcoming runs at $J/\psi$ peak for the $\eta$ and $\eta^\prime$ physics are discussed.


title: Saliency Detection for Stereoscopic Images Based on Depth Confidence Analysis and Multiple Cues Fusion

abstract: Stereoscopic perception is an important part of human visual system that allows the brain to perceive depth. However, depth information has not been well explored in existing saliency detection models. In this letter, a novel saliency detection method for stereoscopic images is proposed. Firstly, we propose a measure to evaluate the reliability of depth map, and use it to reduce the influence of poor depth map on saliency detection. Then, the input image is represented as a graph, and the depth information is introduced into graph construction. After that, a new definition of compactness using color and depth cues is put forward to compute the compactness saliency map. In order to compensate the detection errors of compactness saliency when the salient regions have similar appearances with background, foreground saliency map is calculated based on depth-refined foreground seeds selection mechanism and multiple cues contrast. Finally, these two saliency maps are integrated into a final saliency map through weighted-sum method according to their importance. Experiments on two publicly available stereo datasets demonstrate that the proposed method performs better than other 10 state-of-the-art approaches.


title: Self-oscillation excitation under condition of positive dissipation in a state-dependent potential well

abstract: The self-oscillatory dynamics is considered as motion of a particle in a potential field in the presence of dissipation. Described mechanism of self-oscillation excitation is not associated with peculiarities of a dissipation function, but results from properties of a potential, whose shape depends on a system state. Moreover, features of a potential function allow to realize the self-oscillation excitation in a case of the dissipation function being positive at each point of the phase space. The phenomenon is explored both numerically and experimentally on the example of a double-well oscillator with a state-dependent potential and dissipation. After that a simplified single-well model is studied.


title: Exoplanet phase curves at large phase angles. Diagnostics for extended hazy atmospheres

abstract: At optical wavelengths, Titan's brightness for large Sun-Titan-observer phase angles significantly exceeds its dayside brightness. The brightening that occurs near back-illumination is due to moderately large haze particles in the moon's extended atmosphere that forward-scatter the incident sunlight. Motivated by this phenomenon, here we investigate the forward scattering from currently known exoplanets, its diagnostics possibilities, the observational requirements to resolve it, and potential implications. An analytical expression is derived for the amount of starlight forward-scattered by an exponential atmosphere that takes into account the finite angular size of the star. We use this expression to tentatively estimate how prevalent this phenomenon may be. Based on numerical calculations that consider exoplanet visibility, we identify numerous planets with predicted out-of-transit forward scattering signals of up to tens of parts-per-million provided that aerosols of $>$1 micron size form over an extended vertical region near the optical radius level. We propose that the interpretation of available optical phase curves should be revised to constrain the strength of this phenomenon that might provide insight into aerosol scale heights and particle sizes. For the relatively general atmospheres considered here, forward scattering reduces the transmission-only transit depth by typically less than the equivalent to a scale height. For short-period exoplanets the finite angular size of the star severely affects the amount of radiation scattered towards the observer at mid-transit.


title: The Order on the Light Cone and its induced Topology

abstract: In this article we first correct a recent misconception about a topology that was suggested by Zeeman as a possible alternative to his Fine topology. This misconception appeared while trying to establish the causality in the ambient boundary-ambient space cosmological model. We then show that this topology is actually the intersection topology (in the sense of G.M. Reed) between the Euclidean topology on $\mathbb{R}^4$ and the order topology whose order, namely horismos, is defined on the light cone. Last, but not least, we show that the order topology from horismos belongs to the class of Zeeman topologies. These results accelerate the need for a deeper and more systematic study of the global topological properties of spacetime manifolds.


title: Checking the Compatibility of the Cold Kuiper Belt with a Planetary Instability Migration Model

abstract: The origin of the orbital structure of the cold component of the Kuiper belt is still a hot subject of investigation. Several features of the solar system suggest that the giant planets underwent a phase of global dynamical instability, but the actual dynamical evolution of the planets during the instability is still debated. To explain the structure of the cold Kuiper belt, Nesvorny (2015, AJ 150,68) argued for a soft instability, during which Neptune never achieved a very eccentric orbit. Here we investigate the possibility of a more violent instability, from an initially more compact fully resonant configuration of 5 giant planets. We show that the orbital structure of the cold Kuiper belt can be reproduced quite well provided that the cold population formed in situ, with an outer edge between 44 and 45 au and never had a large mass.


title: Regularizing Deep Neural Networks by Noise: Its Interpretation and Optimization

abstract: Overfitting is one of the most critical challenges in deep neural networks, and there are various types of regularization methods to improve generalization performance. Injecting noises to hidden units during training, e.g., dropout, is known as a successful regularizer, but it is still not clear enough why such training techniques work well in practice and how we can maximize their benefit in the presence of two conflicting objectives---optimizing to true data distribution and preventing overfitting by regularization. This paper addresses the above issues by 1) interpreting that the conventional training methods with regularization by noise injection optimize the lower bound of the true objective and 2) proposing a technique to achieve a tighter lower bound using multiple noise samples per training example in a stochastic gradient descent iteration. We demonstrate the effectiveness of our idea in several computer vision applications.


title: Space-time $L^2$ estimates, regularity and almost global existence for elastic waves

abstract: In this paper, we first establish a kind of weighted space-time $L^2$ estimate, which belongs to Keel-Smith-Sogge type estimates, for perturbed linear elastic wave equations. This estimate refines the corresponding one established by the second author [J. Differential Equations 263(2017), 1947--1965] and is proved by combining the methods in the former paper, the first author, Wang and Yokoyama's paper [Adv. Differential Equations 17 (2012), 267--306], and some new ingredients. Then together with some weighted Sobolev inequalities, this estimate is used to show a refined version of almost global existence of classical solutions for nonlinear elastic waves with small initial data. Compared with former almost global existence results for nonlinear elastic waves due to John [Comm. Pure Appl. Math. 41 (1988) 615--666], Klaierman-Sideris [Comm. Pure Appl. Math. 49 (1996) 307--321], the main innovation of our one is that it considerably improves the amount of regularity of initial data, i.e., the Sobolev regularity of initial data is assumed to be the smallest among all the admissible Sobolev spaces of integer order in the standard local existence theory. Finally, in the radially symmetric case, we establish the almost global existence of a low regularity solution for every small initial data in $H^3\times H^2$.


title: Neutrino spin oscillations in curved space-time under the influence of external fields

abstract: We study neutrino spin oscillations in background matter under the influence of strong electromagnetic and gravitational fields. The neutrino spin evolution is treated quasiclassically. We derive the effective Hamiltonian governing spin oscillations of a neutrino moving in the vicinity of a black hole and interacting with a relativistic magnetized accretion disk around a black hole. Applications for the studies of spin oscillations of astrophysical neutrinos are considered.


title: Performance of the upgraded ultracold neutron source at Los Alamos National Laboratory and its implication for a possible neutron electric dipole moment experiment

abstract: The ultracold neutron (UCN) source at Los Alamos National Laboratory (LANL), which uses solid deuterium as the UCN converter and is driven by accelerator spallation neutrons, has been successfully operated for over 10 years, providing UCN to various experiments, as the first production UCN source based on the superthermal process. It has recently undergone a major upgrade. This paper describes the design and performance of the upgraded LANL UCN source. Measurements of the cold neutron spectrum and UCN density are presented and compared to Monte Carlo predictions. The source is shown to perform as modeled. The UCN density measured at the exit of the biological shield was $184(32)$ UCN/cm$^3$, a four-fold increase from the highest previously reported. The polarized UCN density stored in an external chamber was measured to be $39(7)$ UCN/cm$^3$, which is sufficient to perform an experiment to search for the nonzero neutron electric dipole moment with a one-standard-deviation sensitivity of $\sigma(d_n) = 3\times 10^{-27}$ $e\cdot$cm.


title: Inferring Mesoscale Models of Neural Computation

abstract: Recent years have seen dramatic progress in the development of techniques for measuring the activity and connectivity of large populations of neurons in the brain. However, as these techniques grow ever more powerful---allowing us to even contemplate measuring every neuron in entire brain---a new problem arises: how do we make sense of the mountains of data that these techniques produce? Here, we argue that the time is ripe for building an intermediate or "mesoscale" computational theory that can bridge between single-cell (microscale) accounts of neural function and behavioral (macroscale) accounts of animal cognition and environmental complexity. Just as digital accounts of computation in conventional computers abstract away the non-essential dynamics of the analog circuits that implementing gates and registers, so too a computational account of animal cognition can afford to abstract from the non-essential dynamics of neurons. We argue that the geometry of neural circuits is essential in explaining the computational limitations and technological innovations inherent in biological information processing. We propose a blueprint for how to employ tools from modern machine learning to automatically infer a satisfying mesoscale account of neural computation that combines functional and structural data, with an emphasis on learning and exploiting regularity and repeating motifs in neuronal circuits. Rather than suggest a specific theory, we present a new class of scientific instruments that can enable neuroscientists to design, propose, implement and test mesoscale theories of neural computation.


title: A Model Predictive Approach to Preventing Cascading Failures of Power Systems

abstract: Power system blackouts are usually triggered by the initial contingency and then deteriorate as the branch outage spreads quickly. Thus, it is crucial to eliminate the propagation of cascading outages in its infancy. In this paper, a model predictive approach is proposed to protect power grids against cascading blackout by timely shedding load on buses. The cascading dynamics of power grids is described by a cascading outage model of transmission lines coupled with the DC power flow equation. In addition, a nonlinear convex optimization formulation is established to characterize the optimal load shedding for the mitigation of cascading failures. As a result, two protection schemes are designed on the basis of the optimization formulation. One scheme carries out the remedial action once for all, while the other focuses on the consecutive protection measures. Saddle point dynamics is employed to provide a numerical solution to the proposed optimization problem, and its global convergence is guaranteed in theory. Finally, numerical simulations on IEEE 57 Bus Systems have been implemented to validate the proposed approach in terms of preventing the degradation of power grids.


title: Approximate Hotspots of Orthogonal Trajectories

abstract: In this paper we study the problem of finding hotspots, i.e. regions in which a moving entity has spent a significant amount of time, for polygonal trajectories. The fastest optimal algorithm, due to Gudmundsson, van Kreveld, and Staals (2013) finds an axis-parallel square hotspot of fixed side length in $O(n^2)$. Limiting ourselves to the case in which the entity moves in a direction parallel either to the $x$ or the $y$-axis, We present an approximation algorithm with the time complexity $O(n \log^3 n)$ and approximation factor $1/2$.


title: Search for the Standard Model Higgs boson decaying to b-quark pairs with the ATLAS detector at the LHC

abstract: After many years of searches, the Higgs boson was observed by the ATLAS and CMS collaborations in July 2012. Since then, its properties have been measured using primarily the bosonic decay channels: $H \to \gamma\gamma $, $H \to Z Z$, $H \to W^+ W^-$. In order to probe the predictions of the Standard Model and the nature of the Higgs boson, it is also fundamental to measure its couplings to fermions and, in particular, to quarks. In this paper we present the ATLAS results in the search for the Higgs boson decaying to b-quark pairs using two different production channels: associated production with a vector boson and vector boson fusion with an additional hard photon, using 36.1 $\mathrm{fb^{-1}}$ the $pp$ collisions delivered by the LHC at a center of mass energy of 13 $\mathrm{TeV}$.


title: Towards mechanical characterization of soft digital materials for multimaterial 3D-printing

abstract: We study mechanical behavior of soft rubber-like digital materials used in Polyjet multi-material 3D-printing to create deformable composite materials and flexible structures. These soft digital materials are frequently treated as linear elastic materials in the literature. However, our experiments clearly show that these materials exhibit significant non-linearities under large strain regime. Moreover, the materials demonstrate pronounced rate-dependent behavior. In particular, their instantaneous moduli as well as ultimate strain and stress significantly depend on the strain rate. To take into account both hyper- and viscoelasticity phenomena, we employ the Quasi-Linear Viscoelastic (QLV) model with instantaneous Yeoh strain-energy density function. We show that the QLV-Yeoh model accurately describes the mechanical behavior of the majority of the soft digital materials under uniaxial tension.


title: Signatures of extra dimensions in gravitational waves from black hole quasi-normal modes

abstract: In this work we have derived the evolution equation for gravitational perturbation in four dimensional spacetime in presence of a spatial extra dimension. The evolution equation is derived by perturbing the \emph{effective} gravitational field equations on the four dimensional spacetime, which inherits non-trivial higher dimensional effects. Note that this is \emph{different} from the perturbation of the five dimensional gravitational field equations, existing in literature, and possess quantitatively new features. The gravitational perturbation has further been decomposed into a purely four dimensional part and another piece that depends on extra dimensions. The four dimensional gravitational perturbation now admits massive propagating degrees of freedom, owing to the existence of higher dimensions. We have also studied the influence of these massive propagating modes on the quasi-normal mode frequencies, signaling the higher dimensional nature of the spacetime, and have contrasted these massive modes with the massless modes in general relativity. Surprisingly, it turns out that the massive modes experience much smaller damping compared to the massless modes in general relativity and may even \emph{dominate} over and above the general relativity contribution if one observes the ringdown phase of a black hole merger event at sufficiently late times. Furthermore, the whole analytical framework has been supplemented by the fully numerical Cauchy evolution problem as well. In this context we have shown that except for minute details the overall features of the gravitational perturbations are captured in both the Cauchy evolution as well as in the analysis of quasi-normal modes. The implications on observations of black holes with LIGO and proposed space missions like LISA are also discussed.


title: A graphical, scalable and intuitive method for the placement and the connection of biological cells

abstract: We introduce a graphical method originating from the computer graphics domain that is used for the arbitrary and intuitive placement of cells over a two-dimensional manifold. Using a bitmap image as input, where the color indicates the identity of the different structures and the alpha channel indicates the local cell density, this method guarantees a discrete distribution of cell position respecting the local density function. This method scales to any number of cells, allows to specify several different structures at once with arbitrary shapes and provides a scalable and versatile alternative to the more classical assumption of a uniform non-spatial distribution. Furthermore, several connection schemes can be derived from the paired distances between cells using either an automatic mapping or a user-defined local reference frame, providing new computational properties for the underlying model. The method is illustrated on a discrete homogeneous neural field, on the distribution of cones and rods in the retina and on a coronal view of the basal ganglia.


title: Neutrino telescope searches for dark matter in the Sun

abstract: I give a brief review of a few recent developments and future directions in the search for dark matter using high-energy neutrinos from the Sun. This includes the ability to recast neutrino telescope limits on nuclear scattering of dark matter to arbitrary new theories, and new calculations of the solar atmospheric background relevant to such searches. I also touch on applications to global searches for new physics, and prospects for improving searches for asymmetric dark matter in the Sun.


title: Microaneurysm Detection in Fundus Images Using a Two-step Convolutional Neural Networks

abstract: Diabetic Retinopathy (DR) is a prominent cause of blindness in the world. The early treatment of DR can be conducted from detection of microaneurysms (MAs) which appears as reddish spots in retinal images. An automated microaneurysm detection can be a helpful system for ophthalmologists. In this paper, deep learning, in particular convolutional neural network (CNN), is used as a powerful tool to efficiently detect MAs from fundus images. In our method a new technique is used to utilise a two-stage training process which results in an accurate detection, while decreasing computational complexity in comparison with previous works. To validate our proposed method, an experiment is conducted using Keras library to implement our proposed CNN on two standard publicly available datasets. Our results show a promising sensitivity value of about 0.8 at the average number of false positive per image greater than 6 which is a competitive value with the state-of-the-art approaches.


title: Kernel-based Approximation Methods for Generalized Interpolations: A Deterministic or Stochastic Problem?

abstract: In this article, we solve a deterministically generalized interpolation problem by a stochastic approach. We introduce a kernel-based probability measure on a Banach space by a covariance kernel which is defined on the dual space of the Banach space. The kernel-based probability measure provides a numerical tool to construct and analyze the kernel-based estimators conditioned on non-noise data or noisy data including algorithms and error analysis. Same as meshfree methods, we can also obtain the kernel-based approximate solutions of elliptic partial differential equations by the kernel-based probability measure.


title: K-means clustering for efficient and robust registration of multi-view point sets

abstract: Generally, there are three main factors that determine the practical usability of registration, i.e., accuracy, robustness, and efficiency. In real-time applications, efficiency and robustness are more important. To promote these two abilities, we cast the multi-view registration into a clustering task. All the centroids are uniformly sampled from the initially aligned point sets involved in the multi-view registration, which makes it rather efficient and effective for the clustering. Then, each point is assigned to a single cluster and each cluster centroid is updated accordingly. Subsequently, the shape comprised by all cluster centroids is used to sequentially estimate the rigid transformation for each point set. For accuracy and stability, clustering and transformation estimation are alternately and iteratively applied to all point sets. We tested our proposed approach on several benchmark datasets and compared it with state-of-the-art approaches. Experimental results validate its efficiency and robustness for the registration of multi-view point sets.


title: Interstellar Pickup Ion Observations to 38 AU

abstract: We provide the first direct observations of interstellar H+ and He+ pickup ions in the solar wind from 22 AU to 38 AU. We use the Vasyliunas and Siscoe model functional form to quantify the pickup ion distributions, and while the fit parameters generally lie outside their physically expected ranges, this form allows fits that quantify variations in the pickup H+ properties with distance. By ~20 AU, the pickup ions already provide the dominant internal pressure in the solar wind. We determine the radial trends and extrapolate them to the termination shock at ~90 AU, where the pickup H+ to core solar wind density reaches ~0.14. The pickup H+ temperature and thermal pressure increase from 22-38 AU, indicating additional heating of the pickup ions. This produces very large extrapolated ratios of pickup H+ to solar wind temperature and pressure and an extrapolated ratio of the pickup ion pressure to the solar wind dynamic pressure at the termination shock of ~0.16. Such a large ratio has profound implications for moderating the termination shock and the overall outer heliospheric interaction. We also identify suprathermal tails in the H+ spectra and complex features in the He+ spectra, likely indicating variations in the pickup ion history and processing. Finally, we discover enhancements in both H+ and He+ populations just below their cutoff energies, which may be associated with enhanced local pickup. This study serves to document the release and as the citable reference of these pickup ion data for broad community use and analysis.


title: Thermomechanical properties of zirconium tungstate/hydrogenated nitrile butadiene rubber (HNBR) composites for low-temperature applications

abstract: Rubber compounds for pressure sealing application typically have inferior dimensional stability with temperature fluctuations compared with their steel counterparts. This effect may result in seal leakage failures when subjected to decreases in temperature. Composites of hydrogenated nitrile butadiene rubber (HNBR) and zirconium tungstate as a negative thermal expansion filler were prepared in order to control the thermal expansivity of the material. The amount of zirconium tungstate (ZrW2O8) was varied in the range of 0 to about 40 vol%. The coefficient of thermal expansion (CTE), bulk modulus, uniaxial extension and compression set properties were measured. The CTE of the ZrW2O8-filled HNBR decreases with the filler content and it is reduced by a factor of 2 at the highest filler concentration used. The filler effect on CTE is found to be stronger when HNBR is below the glass transition temperature. The experimental thermal expansion data of the composites are compared with the theoretical estimates and predictions given by FEA. The effect of ZrW2O8 on the mechanical characteristics and compression set of these materials is also discussed.


title: Quantum Locker Using a Novel Verification Algorithm and Its Experimental Realization in IBM Quantum Computer

abstract: It is well known that Grover's algorithm asymptotically transforms an equal superposition state into an eigenstate (of a given basis). Here, we demonstrate a verification algorithm based on weak measurement which can achieve the same purpose even if the qubit is not in an equal superposition state. The proposed algorithm highlights the distinguishability between any arbitrary single qubit superposition state and an eigenstate. We apply this algorithm to propose the scheme of a Quantum Locker, a protocol in which any legitimate party can verify his/her authenticity by using a newly developed quantum One-Time Password (OTP) and retrieve the necessary message from the locker. We formally explicate the working of quantum locker in association with the quantum OTP, which theoretically offers a much higher security against any adversary, as compared to any classical security device.


title: Hyperbolic structures on groups

abstract: For every group $G$, we introduce the set of hyperbolic structures on $G$, denoted $\mathcal{H}(G)$, which consists of equivalence classes of (possibly infinite) generating sets of $G$ such that the corresponding Cayley graph is hyperbolic; two generating sets of $G$ are equivalent if the corresponding word metrics on $G$ are bi-Lipschitz equivalent. Alternatively, one can define hyperbolic structures in terms of cobounded $G$-actions on hyperbolic spaces. We are especially interested in the subset $\mathcal{AH}(G)\subseteq \mathcal{H}(G)$ of acylindrically hyperbolic structures on $G$, i.e., hyperbolic structures corresponding to acylindrical actions. Elements of $\mathcal{H}(G)$ can be ordered in a natural way according to the amount of information they provide about the group $G$. The main goal of this paper is to initiate the study of the posets $\mathcal{H}(G)$ and $\mathcal{AH}(G)$ for various groups $G$. We discuss basic properties of these posets such as cardinality and existence of extremal elements, obtain several results about hyperbolic structures induced from hyperbolically embedded subgroups of $G$, and study to what extent a hyperbolic structure is determined by the set of loxodromic elements and their translation lengths.


title: Universal scaling in first-order phase transitions mixed with nucleation and growth

abstract: Matter exhibits phases and their transitions. These transitions are classified as first-order phase transitions (FOPTs) and continuous ones. While the latter has a well-established theory of the renormalization group, the former is only qualitatively accounted for by classical theories of nucleation, since their predictions often disagree with experiments by orders of magnitude. A theory to integrate FOPTs into the framework of the renormalization-group theory has been proposed but seems to contradict with extant wisdom and lacks numerical evidence. Here we show that universal hysteresis scaling as predicted by the renormalization-group theory emerges unambiguously when the theory is combined intimately with the theory of nucleation and growth in the FOPTs of the paradigmatic two-dimensional Ising model driven by a linearly varying externally applied field below its critical point. This not only provides a new method to rectify the nucleation theories, but also unifies the theories for both classes of transitions and FOPTs can be studied using universality and scaling similar to their continuous counterpart.


title: Community Aware Random Walk for Network Embedding

abstract: Social network analysis provides meaningful information about behavior of network members that can be used for diverse applications such as classification, link prediction. However, network analysis is computationally expensive because of feature learning for different applications. In recent years, many researches have focused on feature learning methods in social networks. Network embedding represents the network in a lower dimensional representation space with the same properties which presents a compressed representation of the network. In this paper, we introduce a novel algorithm named "CARE" for network embedding that can be used for different types of networks including weighted, directed and complex. Current methods try to preserve local neighborhood information of nodes, whereas the proposed method utilizes local neighborhood and community information of network nodes to cover both local and global structure of social networks. CARE builds customized paths, which are consisted of local and global structure of network nodes, as a basis for network embedding and uses the Skip-gram model to learn representation vector of nodes. Subsequently, stochastic gradient descent is applied to optimize our objective function and learn the final representation of nodes. Our method can be scalable when new nodes are appended to network without information loss. Parallelize generation of customized random walks is also used for speeding up CARE. We evaluate the performance of CARE on multi label classification and link prediction tasks. Experimental results on various networks indicate that the proposed method outperforms others in both Micro and Macro-f1 measures for different size of training data.


title: Objective acceleration for unconstrained optimization

abstract: Acceleration schemes can dramatically improve existing optimization procedures. In most of the work on these schemes, such as nonlinear Generalized Minimal Residual (N-GMRES), acceleration is based on minimizing the $\ell_2$ norm of some target on subspaces of $\mathbb{R}^n$. There are many numerical examples that show how accelerating general purpose and domain-specific optimizers with N-GMRES results in large improvements. We propose a natural modification to N-GMRES, which significantly improves the performance in a testing environment originally used to advocate N-GMRES. Our proposed approach, which we refer to as O-ACCEL (Objective Acceleration), is novel in that it minimizes an approximation to the \emph{objective function} on subspaces of $\mathbb{R}^n$. We prove that O-ACCEL reduces to the Full Orthogonalization Method for linear systems when the objective is quadratic, which differentiates our proposed approach from existing acceleration methods. Comparisons with L-BFGS and N-CG indicate the competitiveness of O-ACCEL. As it can be combined with domain-specific optimizers, it may also be beneficial in areas where L-BFGS or N-CG are not suitable.


title: Emergence of a Balanced Core through Dynamical Competition in Heterogeneous Neuronal Networks

abstract: The balance between excitation and inhibition is crucial for neuronal computation. It is observed that the balanced state of neuronal networks exists in many experiments, yet its underlying mechanism remains to be fully clarified. Theoretical studies of the balanced state mainly focus on the analysis of the homogeneous Erd$\ddot{\text{o}}$s-R\'enyi network. However, neuronal networks have been found to be inhomogeneous in many cortical areas. In particular, the connectivity of neuronal networks can be of the type of scale-free, small-world, or even with specific motifs. In this work, we examine the questions of whether the balanced state is universal with respect to network topology and what characteristics the balanced state possesses in inhomogeneous networks such as scale-free and small-world networks. We discover that, for a sparsely but strongly connected inhomogeneous network, despite that the whole network receives external inputs, there is a small active subnetwork (active core) inherently embedded within it. The neurons in this active core have relatively high firing rates while the neurons in the rest of the network are quiescent. Surprisingly, the active core possesses a balanced state and this state is independent of the model of single-neuron dynamics. The dynamics of the active core can be well predicted using the Fokker-Planck equation with the mean-field assumption. Our results suggest that, in the presence of inhomogeneous network connectivity, the balanced state may be ubiquitous in the brain, and the network connectivity in the active core is essentially close to the Erd$\ddot{\text{o}}$s-R\'enyi structure. The existence of the small active core embedded in a large network may provide a potential dynamical scenario underlying sparse coding in neuronal networks.


title: Malware Lineage in the Wild

abstract: Malware lineage studies the evolutionary relationships among malware and has important applications for malware analysis. A persistent limitation of prior malware lineage approaches is to consider every input sample a separate malware version. This is problematic since a majority of malware are packed and the packing process produces many polymorphic variants (i.e., executables with different file hash) of the same malware version. Thus, many samples correspond to the same malware version and it is challenging to identify distinct malware versions from polymorphic variants. This problem does not manifest in prior malware lineage approaches because they work on synthetic malware, malware that are not packed, or packed malware for which unpackers are available. In this work, we propose a novel malware lineage approach that works on malware samples collected in the wild. Given a set of malware executables from the same family, for which no source code is available and which may be packed, our approach produces a lineage graph where nodes are versions of the family and edges describe the relationships between versions. To enable our malware lineage approach, we propose the first technique to identify the versions of a malware family and a scalable code indexing technique for determining shared functions between any pair of input samples. We have evaluated the accuracy of our approach on 13 open-source programs and have applied it to produce lineage graphs for 10 popular malware families. Our malware lineage graphs achieve on average a 26 times reduction from number of input samples to number of versions.


title: Comment on the paper by Popova et al. "On a role of quadruple component of magnetic field in defining solar activity in grand cycles"

abstract: The paper by Popova et al. presents an oversimplified mathematical model of solar activity with a claim of predicting/postdicting it for several millennia ahead/backwards. The work contains several flaws devaluating the results: (1) the method is unreliable from the point of view of signal processing (it is impossible to make harmonic predictions for thousands of years based on only 35 years of data) and lacks quality control, (2) the result of post-diction apparently contradicts the observational data. (3) theoretical speculations make little sense, To summarize, a multi-harmonic mathematical model, hardly related to full solar dynamo theory, is presented, which is not applicable to realistic solar conditions because of the significant chaotic/stochastic intrinsic component and strong non-stationarity of solar activity. The obtained result is apparently inconsistent with the data in the past and thus cannot be trusted for the future predictions.


title: Sequential Design and Spatial Modeling for Portfolio Tail Risk Measurement

abstract: We consider calculation of capital requirements when the underlying economic scenarios are determined by simulatable risk factors. In the respective nested simulation framework, the goal is to estimate portfolio tail risk, quantified via VaR or TVaR of a given collection of future economic scenarios representing factor levels at the risk horizon. Traditionally, evaluating portfolio losses of an outer scenario is done by computing a conditional expectation via inner-level Monte Carlo and is computationally expensive. We introduce several inter-related machine learning techniques to speed up this computation, in particular by properly accounting for the simulation noise. Our main workhorse is an advanced Gaussian Process (GP) regression approach which uses nonparametric spatial modeling to efficiently learn the relationship between the stochastic factors defining scenarios and corresponding portfolio value. Leveraging this emulator, we develop sequential algorithms that adaptively allocate inner simulation budgets to target the quantile region. The GP framework also yields better uncertainty quantification for the resulting VaR/TVaR estimators that reduces bias and variance compared to existing methods. We illustrate the proposed strategies with two case-studies in two and six dimensions.


title: An Onsager Singularity Theorem for Leray Solutions of Incompressible Navier-Stokes

abstract: We study in the inviscid limit the global energy dissipation of Leray solutions of incompressible Navier-Stokes on the torus ${\mathbb T}^d$, assuming that the solutions have norms for Besov space $B^{\sigma,\infty}_3({\mathbb T}^d),$ $\sigma\in (0,1],$ that are bounded in the $L^3$-sense in time, uniformly in viscosity. We establish an upper bound on energy dissipation of the form $O(\nu^{(3\sigma-1)/(\sigma+1)}),$ vanishing as $\nu\to0$ if $\sigma>1/3.$ A consequence is that Onsager-type "quasi-singularities" are required in the Leray solutions, even if the total energy dissipation vanishes in the limit $\nu\to 0$, as long as it does so sufficiently slowly. We also give two sufficient conditions which guarantee the existence of limiting weak Euler solutions $u$ which satisfy a local energy balance with possible anomalous dissipation due to inertial-range energy cascade in the Leray solutions. For $\sigma\in (1/3,1)$ the anomalous dissipation vanishes and the weak Euler solutions may be spatially "rough" but conserve energy.


title: Transport in a Single Self-Doped Nanocrystal

abstract: Addressing the optical properties of a single nanoparticle in the infrared is particularly challenging, thus alternative methods for characterizing the conductance spectrum of nanoparticles in this spectral range need to be developed. Here we describe an efficient method of fabricating single nanoparticle tunnel junctions on a chip circuit. We apply this method to narrow band gap nanoparticles of HgSe, which band structure combine the inverted character of the bulk semimetal with quantum confinement and self-doping. Upon tuning the gate bias, measurement reveals the presence of two energy gaps in the spectrum. The wider gap results from the interband gap, while the narrower gap results from intraband transitions. The observation of the latter near zero gate voltage confirms the doped character of the nanoparticle at the single particle level, which is in full agreement with the ensemble optical and transport measurements. Finally we probe the phototransport within a single quantum dot and demonstrate a large photogain mechanism resulting from photogating.


title: Network Model Selection Using Task-Focused Minimum Description Length

abstract: Networks are fundamental models for data used in practically every application domain. In most instances, several implicit or explicit choices about the network definition impact the translation of underlying data to a network representation, and the subsequent question(s) about the underlying system being represented. Users of downstream network data may not even be aware of these choices or their impacts. We propose a task-focused network model selection methodology which addresses several key challenges. Our approach constructs network models from underlying data and uses minimum description length (MDL) criteria for selection. Our methodology measures efficiency, a general and comparable measure of the network's performance of a local (i.e. node-level) predictive task of interest. Selection on efficiency favors parsimonious (e.g. sparse) models to avoid overfitting and can be applied across arbitrary tasks and representations. We show stability, sensitivity, and significance testing in our methodology.


title: Average Spectral Properties of Type Ia Supernova Host Galaxies

abstract: We construct average spectra of host galaxies of slower, faster, bluer, and redder Type Ia Supernovae (SNe Ia) from the SDSS-II supernova survey. The average spectrum of slower declining (broader light-curve width or higher stretch) SN Ia hosts shows stronger emission lines compared to the average spectrum of faster declining (narrower light-curve width or lower stretch) SN Ia hosts. Using pPXF, we find that hosts of slower declining SNe Ia have metallicities that are, on average, 0.24 dex lower than average metallicities of faster declining SN Ia hosts. Similarly, redder SN Ia hosts have slightly higher metallicities than bluer SN Ia hosts. Lick index analysis of metallic lines and Balmer lines show that faster declining SN Ia hosts have relatively higher metal content and have relatively older stellar populations compared with slower declining SN Ia hosts. We calculate average $\rm H_{\alpha}$ Star Formation Rate (SFR), stellar mass, and the specific-SFR (sSFR) of host galaxies in these subgroups of SNe Ia. We find that slower declining SN Ia hosts have significantly higher ($>5\sigma$) sSFR than faster declining SN Ia hosts. A Kolmogorov-Smirnov test shows that these two types of hosts originate from different parent distributions. Our results, when compared with the models of \cite{childress14}, indicate that slower declining SNe Ia, being hosted in actively star-forming galaxies, are young (prompt) SNe Ia, originating from similar progenitor age groups.


title: Near-optimal Sample Complexity Bounds for Robust Learning of Gaussians Mixtures via Compression Schemes

abstract: We prove that $\tilde{\Theta}(k d^2 / \varepsilon^2)$ samples are necessary and sufficient for learning a mixture of $k$ Gaussians in $\mathbb{R}^d$, up to error $\varepsilon$ in total variation distance. This improves both the known upper bounds and lower bounds for this problem. For mixtures of axis-aligned Gaussians, we show that $\tilde{O}(k d / \varepsilon^2)$ samples suffice, matching a known lower bound. Moreover, these results hold in the agnostic-learning/robust-estimation setting as well, where the target distribution is only approximately a mixture of Gaussians. The upper bound is shown using a novel technique for distribution learning based on a notion of `compression.' Any class of distributions that allows such a compression scheme can also be learned with few samples. Moreover, if a class of distributions has such a compression scheme, then so do the classes of products and mixtures of those distributions. The core of our main result is showing that the class of Gaussians in $\mathbb{R}^d$ admits a small-sized compression scheme.


title: Anisotropic flow of identified particles in Pb--Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV

abstract: Anisotropic flow is sensitive to the shear $(\eta/s)$ and bulk ($\zeta/s$) viscosity of the quark-gluon plasma created in heavy-ion collisions, as well as the initial state of such collisions and hadronization mechanisms. In these proceedings, elliptic ($v_2$) and higher harmonic ($v_3, v_4$) flow coefficients of $\pi^{\pm}$, $K^{\pm}$, p$(\overline{\rm{p}})$ and the $\phi$-meson, are presented for Pb-Pb collisions at the highest-ever center-of-mass energy of $\sqrt{s_{\rm NN}}$ = 5.02 TeV. Comparisons to hydrodynamic calculations (IP-Glasma, MUSIC, UrQMD) are shown to constrain the initial conditions and viscosity of the medium.


title: Affine special Kaehler structures in real dimension two

abstract: We review properties of affine special Kaehler structures focusing on singularities of such structures in the simplest case of real dimension two. We describe all possible isolated singularities and compute the monodromy of the flat symplectic connection, which is a part of a special Kaehler structure, near a singularity. Beside numerous local examples, we construct continuous families of special Kaehler structures with isolated singularities on the projective line.


title: On Supporting Digital Journalism: Case Studies in Co-Designing Journalistic Tools

abstract: Since 2013 researchers at University College Dublin in the Insight Centre for Data Analytics have been involved in a significant research programme in digital journalism, specifically targeting tools and social media guidelines to support the work of journalists. Most of this programme was undertaken in collaboration with The Irish Times. This collaboration involved identifying key problems currently faced by digital journalists, developing tools as solutions to these problems, and then iteratively co-designing these tools with feedback from journalists. This paper reports on our experiences and learnings from this research programme, with a view to informing similar efforts in the future.


title: Simultaneous Matrix Diagonalization for Structural Brain Networks Classification

abstract: This paper considers the problem of brain disease classification based on connectome data. A connectome is a network representation of a human brain. The typical connectome classification problem is very challenging because of the small sample size and high dimensionality of the data. We propose to use simultaneous approximate diagonalization of adjacency matrices in order to compute their eigenstructures in more stable way. The obtained approximate eigenvalues are further used as features for classification. The proposed approach is demonstrated to be efficient for detection of Alzheimer's disease, outperforming simple baselines and competing with state-of-the-art approaches to brain disease classification.


title: A non-iterative formula for straightening fillings of Young diagrams

abstract: Young diagrams are fundamental combinatorial objects in representation theory and algebraic geometry. Many constructions that rely on these objects depend on variations of a straightening process that expresses a filling of a Young diagram as a sum of semistandard tableaux subject to certain relations. This paper solves the long standing open problem of giving a non-iterative formula for straightening a filling. We apply our formula to give a complete generalization of a theorem of Gonciulea and Lakshmibai.


title: Relative Perturbation Bounds for the Joint Spectrum of Commuting Tuple of Matrices

abstract: In this paper, we study the relative perturbation bounds for joint eigenvalues of commuting tuples of normal $n \times n$ matrices. Some Hoffman-Wielandt type relative perturbation bounds are proved using the Clifford algebra technique. A result is also extended for diagonalizable matrices which improves a relative perturbation bound for single matrices.


title: Gassmann Theory Applies to Nanoporous Media

abstract: Recent progress in extraction of unconventional hydrocarbon resources has ignited the interest in the studies of nanoporous media. Since many thermodynamic and mechanical properties of nanoscale solids and fluids differ from the analogous bulk materials, it is not obvious whether wave propagation in nanoporous media can be described using the same framework as in macroporous media. Here we test the validity of Gassmann equation using two published sets of ultrasonic measurements for a model nanoporous medium, Vycor glass, saturated with two different fluids, argon and n-hexane. Predictions of the Gassmann theory depend on the bulk and shear moduli of the dry samples, which are known from ultrasonic measurements, and the bulk moduli of the solid and fluid constituents. The solid bulk modulus can be estimated from adsorption-induced deformation or from elastic effective medium theory. The fluid modulus can be calculated according to the Tait-Murnaghan equation at the solvation pressure in the pore. Substitution of these parameters into the Gassmann equation provides predictions consistent with measured data. Our findings set up a theoretical framework for investigation of fluid-saturated nanoporous media using ultrasonic elastic wave propagation.


title: Modular inequalities for the maximal operator in variable Lebesgue spaces

abstract: A now classical result in the theory of variable Lebesgue spaces due to Lerner [A. K. Lerner, On modular inequalities in variable $L^p$ spaces, Archiv der Math. 85 (2005), no. 6, 538-543] is that a modular inequality for the Hardy-Littlewood maximal function in $L^{p(\cdot)}(\mathbb{R}^n)$ holds if and only if the exponent is constant. We generalize this result and give a new and simpler proof. We then find necessary and sufficient conditions for the validity of the weaker modular inequality \[ \int_\Omega Mf(x)^{p(x)}\,dx \ \leq c_1 \int_\Omega |f(x)|^{q(x)}\,dx + c_2, \] where $c_1,\,c_2$ are non-negative constants and $\Omega$ is any measurable subset of $\mathbb{R}^n$. As a corollary we get sufficient conditions for the modular inequality \[ \int_\Omega |Tf(x)|^{p(x)}\,dx \ \leq c_1 \int_\Omega |f(x)|^{q(x)}\,dx + c_2, \] where $T$ is any operator that is bounded on $L^p(\Omega)$, $1<p<\infty$.


title: Two-player incentive compatible outcome functions are affine maximizers

abstract: In mechanism design, for a given type space, there may be incentive compatible outcome functions which are not affine maximizers. Using tools from linear algebra and tropical geometry, we prove that for two-player games on a discrete type space, any given outcome function can be turned into an affine maximizer through a nontrivial perturbation of the type space. Furthermore, our theorems are the strongest possible in this setup.


title: Mental Sampling in Multimodal Representations

abstract: Both resources in the natural environment and concepts in a semantic space are distributed "patchily", with large gaps in between the patches. To describe people's internal and external foraging behavior, various random walk models have been proposed. In particular, internal foraging has been modeled as sampling: in order to gather relevant information for making a decision, people draw samples from a mental representation using random-walk algorithms such as Markov chain Monte Carlo (MCMC). However, two common empirical observations argue against simple sampling algorithms such as MCMC. First, the spatial structure is often best described by a L\'evy flight distribution: the probability of the distance between two successive locations follows a power-law on the distances. Second, the temporal structure of the sampling that humans and other animals produce have long-range, slowly decaying serial correlations characterized as $1/f$-like fluctuations. We propose that mental sampling is not done by simple MCMC, but is instead adapted to multimodal representations and is implemented by Metropolis-coupled Markov chain Monte Carlo (MC$^3$), one of the first algorithms developed for sampling from multimodal distributions. MC$^3$ involves running multiple Markov chains in parallel but with target distributions of different temperatures, and it swaps the states of the chains whenever a better location is found. Heated chains more readily traverse valleys in the probability landscape to propose moves to far-away peaks, while the colder chains make the local steps that explore the current peak or patch. We show that MC$^3$ generates distances between successive samples that follow a L\'evy flight distribution and $1/f$-like serial correlations, providing a single mechanistic account of these two puzzling empirical phenomena.


title: Effect of Pore Geometry on the Compressibility of a Confined Simple Fluid

abstract: Fluids confined in nanopores exhibit properties different from the properties of the same fluids in bulk, among these properties are the isothermal compressibility or elastic modulus. The modulus of a fluid in nanopores can be extracted from ultrasonic experiments or calculated from molecular simulations. Using Monte Carlo simulations in the grand canonical ensemble, we calculated the modulus for liquid argon at its normal boiling point (87.3~K) adsorbed in model silica pores of two different morphologies and various sizes. For spherical pores, for all the pore sizes (diameters) exceeding 2~nm, we obtained a logarithmic dependence of fluid modulus on the vapor pressure. Calculation of the modulus at saturation showed that the modulus of the fluid in spherical pores is a linear function of the reciprocal pore size. The calculation of the modulus of the fluid in cylindrical pores appeared too scattered to make quantitative conclusions. We performed additional simulations at higher temperature (119.6~K), at which Monte Carlo insertions and removals become more efficient. The results of the simulations at higher temperature confirmed both regularities for cylindrical pores and showed quantitative difference between the fluid moduli in pores of different geometries. Both of the observed regularities for the modulus stem from the Tait-Murnaghan equation applied to the confined fluid. Our results, along with the development of the effective medium theories for nanoporous media, set the groundwork for analysis of the experimentally-measured elastic properties of fluid-saturated nanoporous materials.


title: An Adaptive Framework for Missing Depth Inference Using Joint Bilateral Filter

abstract: Depth imaging has largely focused on sensor and intrinsics properties. However, the accuracy of acquire pixel is largely dependent on the capture. We propose a new depth estimation and approximation algorithm which takes an arbitrary 3D point cloud as input, with potentially complex geometric structures, and generates automatically a bounding box which is used to clamp the 3D distribution into a valid range. We then infer the desired compact geometric network from complex 3D geometries by using a series of adaptive joint bilateral filters. Our approach leverages these input depth in the construction of a compact descriptive adaptive filter framework. The built system that allows a user to control the result of capture depth map to fit the target geometry. In addition, it is desirable to visualize structurally problematic areas of the depth data in a dynamic environment. To provide this feature, we investigate a fast update algorithm for the fragility of each pixel's corresponding 3D point using machine learning. We present a new for of feature vector analysis and demonstrate the effectiveness in the dataset. In our experiment, we demonstrate the practicality and benefits of our proposed method by computing accurate solutions captured depth map from different types of sensors and shows better results than existing methods.


title: The intrinsic shape of bulges in the CALIFA survey

abstract: The intrinsic shape of galactic bulges in nearby galaxies provides crucial information to separate bulge types. We intended to derive accurate constraints to the intrinsic shape of bulges to provide new clues on their formation mechanisms and set new limitations for future simulations. We retrieved the intrinsic shape of a sample of CALIFA bulges using a statistical approach. Taking advantage of GalMer numerical simulations of binary mergers we estimated the reliability of the procedure. Analyzing the $i$-band mock images of resulting lenticular remnants, we studied the intrinsic shape of their bulges at different galaxy inclinations. Finally, we introduced a new ($B/A$, $C/A$) diagram to analyze possible correlations between the intrinsic shape and the properties of bulges. We tested the method on simulated lenticular remnants, finding that for galaxies with inclinations $25^{\circ} < \theta < 65^{\circ}$ we can safely derive the intrinsic shape of their bulges. We found that our CALIFA bulges tend to be nearly oblate systems (66%), with a smaller fraction of prolate spheroids (19%) and triaxial ellipsoids (15%). The majority of triaxial bulges are in barred galaxies (75%). Moreover, we found that bulges with low S\'ersic indices or in galaxies with low bulge-to-total luminosity ratios form a heterogeneous class of objects; additionally, also bulges in late-type galaxies or in less massive galaxies have no preference in being oblate, prolate, or triaxial. On the contrary, bulges with high S\'ersic index, in early-type galaxies, or in more massive galaxies are mostly oblate systems. We concluded that various evolutionary pathways may coexist in galaxies, with merging events and dissipative collapse being the main mechanisms driving the formation of the most massive oblate bulges and bar evolution reshaping the less massive triaxial bulges.


title: DPG* Method

abstract: We introduce a cousin of the DPG method - the DPG* method - discuss their relationship and compare the two methods through numerical experiments.


title: Development of the LAGO Project in Chiapas-Mexico

abstract: The Latin American Giant Observatory (LAGO) is an extended astroparticle observatory with the goal of studying Gamma Ray Bursts (among other extreme universe phenomena), space weather and atmospheric radiation at ground level. It consists of a network of several Water Cherenkov Detectors (WCD) located at different sites and different latitudes along the American Continent (from Mexico up to the Antarctic region). Another interest of LAGO is to encourage and support the development of experimental basic research in Latin America, mainly with low cost equipment. In the case of Chiapas, Mexico, the experimental astroparticle physics activity was limited, up to now, to data analysis from other detectors located far away from the region. Thanks to the collaboration within LAGO, the deployment of one WCD is ongoing at the Universidad Aut\'onoma de Chiapas (UNACH). This will allow, for the first time in the region, to train students and researchers in the deployment processes. Till now the setup of the signal-processing electronics has been performed and the characterization of the photomultiplier tube is currently being done. The main, short-term goal is to install one WCD on top of the Tacan\'a volcano in Chiapas in a short term. The status of the work is presented.


title: Popularity of arXiv.org within Computer Science

abstract: It may seem surprising that, out of all areas of science, computer scientists have been slow to post electronic versions of papers on sites like arXiv.org. Instead, computer scientists have tended to place papers on our individual home pages, but this loses the benefits of aggregation, namely notification and browsing. But this is changing. More and more computer scientists are now using the arXiv. At the same time, there is ongoing discussion and controversy about how prepublication affects peer review, especially for double-blind conferences. This discussion is often carried out with precious little evidence of how popular prepublication is. We measure what percentage of papers in computer science are placed on the arXiv, by cross-referencing published papers in DBLP with e-prints on arXiv. We found: * Usage of arXiv.org has risen dramatically among the most selective conferences in computer science. In 2017, fully 23% of papers had e-prints on arXiv, compared to only 1% ten years ago. * Areas of computer science vary widely in e-print prevalence. In theoretical computer science and machine learning, over 60% of published papers are on arXiv, while other areas are essentially zero. In most areas, arXiv usage is rising. * Many researchers use arXiv for posting preprints. Of the 2017 published papers with arXiv e-prints, 56% were preprints that were posted before or during peer review. Our paper describes these results as well as policy implications for researchers and practitioners.


title: MARXS: A modular software to ray-trace X-ray instrumention

abstract: To obtain the best possible scientific result, astronomers must understand the properties of the available instrumentation well. This is important both when designing new instruments and when using existing instruments close to the limits of their specified capabilities or beyond. Ray-tracing is a technique for numerical simulations where the path of many light rays is followed through the system to understand how individual system components influence the observed properties, such as the shape of the point-spread-function (PSF). In instrument design, such simulations can be used to optimize the performance. For observations with existing instruments this helps to discern instrumental artifacts from a true signal. Here, we describe MARXS, a new python package designed to simulate X-ray instruments on satellites and sounding rockets. MARXS uses probability tracking of photons and has polarimetric capabilities.


title: Brownian Motion with Singular Time-Dependent Drift

abstract: In this paper we study weak solutions for the following type of stochastic differential equation \[ dX_{t}=dW_{t}+b(t, X_{t})dt, \quad t\ge s, \quad X_{s}=x, \] where $b: [0,\infty) \times \mathbb{R}^{d} \to \mathbb{R}^{d}$ is a measurable drift, $W=(W_{t})_{t \ge 0}$ is a $d$-dimensional Brownian motion and $(s,x)\in [0,\infty) \times \mathbb{R}^{d}$ is the starting point. A solution $X=(X_t)_{t \ge s}$ for the above SDE is called a Brownian motion with time-dependent drift $b$ starting from $(s,x)$. Under the assumption that $|b|$ belongs to the forward-Kato class $\mathcal{F} \mathcal{K}_{d-1}^{\alpha}$ for some $\alpha \in (0,1/2)$, we prove that the above SDE has a unique weak solution for every starting point $(s,x)\in [0,\infty) \times \mathbb{R}^{d}$.


title: Torsion subgroups of rational elliptic curves over the compositum of all $D_4$ extensions of the rational numbers

abstract: Let $E/\mathbb{Q}$ be an elliptic curve and let $\mathbb{Q}(D_4^\infty)$ be the compositum of all extensions of $\mathbb{Q}$ whose Galois closure has Galois group isomorphic to a quotient of a subdirect product of a finite number of transitive subgroups of $D_4$. In this article we first show that $\mathbb{Q}(D_4^\infty)$ is in fact the compositum of all $D_4$ extensions of $\mathbb{Q}$ and then we prove that the torsion subgroup of $E(\mathbb{Q}(D_4^\infty))$ is finite and determine the 24 possibilities for its structure. We also give a complete classification of the elliptic curves that have each possible torsion structure in terms of their $j$-invariants.


title: Mysl, smysl, svet

abstract: Speech on the occasion of accepting the Dagmar and Vaclav Havel Foundation VIZE 97 Prize for 2017. Delivered at Prague Crossroads, October 5, 2017


title: Searches for diboson resonances at CMS

abstract: A list of searches is presented for massive resonances decaying to WW, WZ, ZZ, WH, ZH and HH boson pairs in proton-proton collision data collected by the CMS experiment at the LHC. The data are taken at centre-of-mass energies of 13TeV, corresponding to respective integrated luminosities of 12.9/fb(ICHEP) and up to 35.9/fb. The results are interpreted in the context of heavy vector triplet and singlet models that mimic properties of composite-Higgs models predicting W' and Z' bosons decaying to WZ, WW, WH, and ZH bosons. A model with a bulk graviton that decays into WW and ZZ is also considered. Upper limits are set on the production cross section of the two models and no evidence is found for a signal.


title: Hybrid DDP in Clutter (CHDDP): Trajectory Optimization for Hybrid Dynamical System in Cluttered Environments

abstract: We present an algorithm for obtaining an optimal control policy for hybrid dynamical systems in cluttered environments. To the best of our knowledge, this is the first attempt to have a locally optimal solution for this specific problem setting. Our approach extends an optimal control algorithm for hybrid dynamical systems in the obstacle-free case to environments with obstacles. Our method does not require any preset mode sequence or heuristics to prune the exponential search of mode sequences. By first solving the relaxed problem of getting an obstacle-free, dynamically feasible trajectory and then solving for both obstacle-avoidance and optimality, we can generate smooth, locally optimal control policies. We demonstrate the performance of our algorithm on a box-pushing example in a number of environments against the baseline of randomly sampling modes and actions with a Kinodynamic RRT.


title: Curved O-operator systems

abstract: We introduce the notion of curved $\mathcal{O}$-operator systems as a generalization of T. Brzezi\'{n}ski's (curved) Rota-Baxter systems, and then investigate their relations with (tri)dendriform systems and pre-Lie algebras. Moreover, we consider the special cases of generalized Rota-Baxter algebras with respect to associative compatible pairs and then double curved Rota-Baxter systems.


title: Learners that Use Little Information

abstract: We study learning algorithms that are restricted to using a small amount of information from their input sample. We introduce a category of learning algorithms we term $d$-bit information learners, which are algorithms whose output conveys at most $d$ bits of information of their input. A central theme in this work is that such algorithms generalize. We focus on the learning capacity of these algorithms, and prove sample complexity bounds with tight dependencies on the confidence and error parameters. We also observe connections with well studied notions such as sample compression schemes, Occam's razor, PAC-Bayes and differential privacy. We discuss an approach that allows us to prove upper bounds on the amount of information that algorithms reveal about their inputs, and also provide a lower bound by showing a simple concept class for which every (possibly randomized) empirical risk minimizer must reveal a lot of information. On the other hand, we show that in the distribution-dependent setting every VC class has empirical risk minimizers that do not reveal a lot of information.


title: Phase Retrieval via Linear Programming: Fundamental Limits and Algorithmic Improvements

abstract: A recently proposed convex formulation of the phase retrieval problem estimates the unknown signal by solving a simple linear program. This new scheme, known as PhaseMax, is computationally efficient compared to standard convex relaxation methods based on lifting techniques. In this paper, we present an exact performance analysis of PhaseMax under Gaussian measurements in the large system limit. In contrast to previously known performance bounds in the literature, our results are asymptotically exact and they also reveal a sharp phase transition phenomenon. Furthermore, the geometrical insights gained from our analysis led us to a novel nonconvex formulation of the phase retrieval problem and an accompanying iterative algorithm based on successive linearization and maximization over a polytope. This new algorithm, which we call PhaseLamp, has provably superior recovery performance over the original PhaseMax method.


title: Moment and tail estimates and Banach space valued Non-Central Limit Theorem (NCLT) for sums of multi-indexed random variables, processes and fields

abstract: We derive in this preprint the moment and exponential tail estimates, sufficient conditions for the Non-Central Limit Theorem (NCLT) in the ordinary one-dimensional space as well as in the space of continuous functions for the properly (natural) normalized multi-indexed sums of function of random variables, processes or fields (r.f.), on the other words V-statistics, parametric, in general case. We construct also some examples in order to show the exactness of obtained estimates. We will use the theory of the so-called degenerate approximation of the functions of several variables as well as the theory of Grand Lebesgue Spaces (GLS) of measurable functions (random variables).


title: On a Minkowski geometric flow in the plane: evolution of curves with lack of scale invariance

abstract: We consider a planar geometric flow in which the normal velocity is a nonlocal variant of the curvature. The flow is not scaling invariant and in fact has different behaviors at different spatial scales, thus producing phenomena that are different with respect to both the classical mean curvature flow and the fractional mean curvature flow. In particular, we give examples of neckpinch singularity formation, and we discuss convexity properties of the evolution. We also take into account traveling waves for this geometric flow, showing that a new family of $C^{1,1}$ and convex traveling sets arises in this setting.


title: Unknotting numbers for prime $\theta$-curves up to seven crossings

abstract: Determining unknotting numbers is a large and widely studied problem. We consider the more general question of the unknotting number of a spatial graph. We show the unknotting number of spatial graphs is subadditive. Let $g$ be an embedding of a planar graph $G$, then we show $u(g) \geq \max\{u(s) |$ $s$ is a non-overlapping set of constituents of $g\}$. Focusing on $\theta$-curves, we determine the exact unknotting numbers of the $\theta$-curves in the Litherland-Moriuchi Table. Additionally, we demonstrate unknotting crossing changes for all of the curves. In doing this we introduce new methods for obstructing unknotting number $1$ in $\theta$-curves.


title: On inequalities for convex functions

abstract: We study some properties convex functions fulfill. Among the conclusions we obtain from such result, we are able to prove some nontrivial inequalities among real numbers, and we give an improvement of the reverse triangle inequality in the particular case where three aligned points are considered in a normed space.


title: An Online Optimization Framework for Distributed Fog Network Formation with Minimal Latency

abstract: Fog computing is emerging as a promising paradigm to perform distributed, low-latency computation by jointly exploiting the radio and computing resources of end-user devices and cloud servers. However, the dynamic and distributed formation of local fog networks is highly challenging due to the unpredictable arrival and departure of neighboring fog nodes. Therefore, a given fog node must properly select a set of neighboring nodes and intelligently offload its computational tasks to this set of neighboring fog nodes and the cloud in order to achieve low-latency transmission and computation. In this paper, the problem of fog network formation and task distribution is jointly investigated while considering a hybrid fog-cloud architecture. The goal is to minimize the maximum computational latency by enabling a given fog node to form a suitable fog network and optimize the task distribution, under uncertainty on the arrival process of neighboring fog nodes. To solve this problem, a novel online optimization framework is proposed in which the neighboring nodes are selected by using a threshold-based online algorithm that uses a target competitive ratio, defined as the ratio between the latency of the online algorithm and the offline optimal latency. The proposed framework repeatedly updates its target competitive ratio and optimizes the distribution of the fog node's computational tasks in order to minimize latency. Simulation results show that the proposed framework can successfully select a set of neighboring nodes while reducing latency by up to 19.25% compared to a baseline approach based on the well-known online secretary framework. The results also show how, using the proposed framework, the computational tasks can be properly offloaded between the fog network and a remote cloud server in different network settings.


title: Double-diffusive erosion of the core of Jupiter

abstract: We present Direct Numerical Simulations of the transport of heat and heavy elements across a double-diffusive interface or a double-diffusive staircase, in conditions that are close to those one may expect to find near the boundary between the heavy-element rich core and the hydrogen-helium envelope of giant planets such as Jupiter. We find that the non-dimensional ratio of the buoyancy flux associated with heavy element transport to the buoyancy flux associated with heat transport lies roughly between 0.5 and 1, which is much larger than previous estimates derived by analogy with geophysical double-diffusive convection. Using these results in combination with a core-erosion model proposed by Guillot et al. (2004), we find that the entire core of Jupiter would be eroded within less than 1Myr assuming that the core-envelope boundary is composed of a single interface. We also propose an alternative model that is more appropriate in the presence of a well-established double-diffusive staircase, and find that in this limit a large fraction of the core could be preserved. These findings are interesting in the context of Juno's recent results, but call for further modeling efforts to better understand the process of core erosion from first principles.


title: Robust Decentralized Learning Using ADMM with Unreliable Agents

abstract: Many machine learning problems can be formulated as consensus optimization problems which can be solved efficiently via a cooperative multi-agent system. However, the agents in the system can be unreliable due to a variety of reasons: noise, faults and attacks. Providing erroneous updates leads the optimization process in a wrong direction, and degrades the performance of distributed machine learning algorithms. This paper considers the problem of decentralized learning using ADMM in the presence of unreliable agents. First, we rigorously analyze the effect of erroneous updates (in ADMM learning iterations) on the convergence behavior of multi-agent system. We show that the algorithm linearly converges to a neighborhood of the optimal solution under certain conditions and characterize the neighborhood size analytically. Next, we provide guidelines for network design to achieve a faster convergence. We also provide conditions on the erroneous updates for exact convergence to the optimal solution. Finally, to mitigate the influence of unreliable agents, we propose \textsf{ROAD}, a robust variant of ADMM, and show its resilience to unreliable agents with an exact convergence to the optimum.


title: Equidistribution of divergent orbits of the diagonal group in the space of lattices

abstract: We consider divergent orbits of the group of diagonal matrices in the space of lattices in Euclidean space. We define two natural numerical invariants of such orbits: The discriminant - an integer - and the type - an integer vector. We then study the question of the limit distributional behaviour of these orbits as the discriminant goes to infinity. Using entropy methods we prove that for divergent orbits of a specific type, virtually any sequence of orbits equidistribute as the discriminant goes to infinity. Using measure rigidity for higher rank diagonal actions we complement this result and show that in dimension 3 or higher only very few of these divergent orbits can spend all of their life-span in a given compact set before they diverge.


title: Eigenvalues of even very nice Toeplitz matrices can be unexpectedly erratic

abstract: It was shown in a series of recent publications that the eigenvalues of $n\times n$ Toeplitz matrices generated by so-called simple-loop symbols admit certain regular asymptotic expansions into negative powers of $n+1$. On the other hand, recently two of the authors considered the pentadiagonal Toeplitz matrices generated by the symbol $g(x)=(2\sin(x/2))^4$, which does not satisfy the simple-loop conditions, and derived asymptotic expansions of a more complicated form. We here use these results to show that the eigenvalues of the pentadiagonal Toeplitz matrices do not admit the expected regular asymptotic expansion. This also delivers a counter-example to a conjecture by Ekstr\"{o}m, Garoni, and Serra-Capizzano and reveals that the simple-loop condition is essential for the existence of the regular asymptotic expansion.


title: Ensemble Kalman Filters (EnKF) for State Estimation and Prediction of Two-time Scale Nonlinear Systems with Application to Gas Turbine Engines

abstract: In this paper, we propose and develop a methodology for nonlinear systems health monitoring by modeling the damage and degradation mechanism dynamics as "slow" states that are augmented with the system "fast" dynamical states. This augmentation results in a two-time scale nonlinear system that is utilized for development of health estimation and prediction modules within a health monitoring framework. Towards this end, a two-time scale filtering approach is developed based on the ensemble Kalman filtering (EnKF) approach by taking advantage of the model reduction concept. The performance of our proposed two-time scale ensemble Kalman filters is shown to be superior and less computationally intensive in terms of the equivalent flop (EF) complexity metric when compared to well-known particle filtering (PF) approaches. Our proposed methodology is then applied to a gas turbine engine that is affected by erosion of the turbine as the degradation phenomenon and damage mechanism. Extensive comparative studies are conducted to validate and demonstrate the advantages and capabilities of our proposed framework and methodology.


title: Simulation for the ATLAS Upgrade Strip Tracker

abstract: ATLAS is making extensive efforts towards preparing a detector upgrade for the high luminosity operations of the LHC (HL-LHC), which will commence operation in about 10 years. The current ATLAS Inner Detector will be replaced by an all-silicon tracker (comprising an inner Pixel tracker and outer Strip tracker). The software currently used for the new silicon tracker is broadly inherited from that used for the LHC Run-1 and Run-2, but many new developments have been made to better fulfill the future detector and operation requirements. One aspect in particular which will be highlighted is the simulation software for the Strip tracker. The available geometry description software (including the detailed description for all the sensitive elements, the services, etc.) did not allow for accurate modelling of the planned detector design. A range of sensors/layouts for the Strip tracker are being considered and must be studied in detailed simulations in order to assess the performance and ascertain that requirements are met. For this, highly flexibility geometry building is required from the simulation software. A new Xml-based detector description framework has been developed to meet the aforementioned challenges. We will present the design of the framework and its validation results.


title: Shared High Value Research Resources: The CamCAN Human Lifespan Neuroimaging Dataset Processed on the Open Science Grid

abstract: The CamCAN Lifespan Neuroimaging Dataset, Cambridge (UK) Centre for Ageing and Neuroscience, was acquired and processed beginning in December, 2016. The referee consensus solver deployed to the Open Science Grid was used for this task. The dataset includes demographic and screening measures, a high-resolution MRI scan of the brain, and whole-head magnetoencephalographic (MEG) recordings during eyes closed rest (560 sec), a simple task (540 sec), and passive listening/viewing (140 sec). The data were collected from 619 neurologically normal individuals, ages 18-87. The processed results from the resting recordings are completed and available online. These constitute 1.7 TBytes of data including the location within the brain (1 mm resolution), time stamp (1 msec resolution), and 80 msec time course for each of 3.7 billion validated neuroelectric events, i.e. mean 6.1 million events for each of the 619 participants. The referee consensus solver provides high yield (mean 11,000 neuroelectric currents/sec; standard deviation (sd): 3500/sec) high confidence (p < 10-12 for each identified current) measures of the neuroelectric currents whose magnetic fields are detected in the MEG recordings. We describe the solver, the implementation of the solver deployed on the Open Science Grid, the workflow management system, the opportunistic use of high performance computing (HPC) resources to add computing capacity to the Open Science Grid reserved for this project, and our initial findings from the recently completed processing of the resting recordings. This required 14 million core hours, i.e. 40 core hours per second of data.


title: On Hashing-Based Approaches to Approximate DNF-Counting

abstract: Propositional model counting is a fundamental problem in artificial intelligence with a wide variety of applications, such as probabilistic inference, decision making under uncertainty, and probabilistic databases. Consequently, the problem is of theoretical as well as practical interest. When the constraints are expressed as DNF formulas, Monte Carlo-based techniques have been shown to provide a fully polynomial randomized approximation scheme (FPRAS). For CNF constraints, hashing-based approximation techniques have been demonstrated to be highly successful. Furthermore, it was shown that hashing-based techniques also yield an FPRAS for DNF counting without usage of Monte Carlo sampling. Our analysis, however, shows that the proposed hashing-based approach to DNF counting provides poor time complexity compared to the Monte Carlo-based DNF counting techniques. Given the success of hashing-based techniques for CNF constraints, it is natural to ask: Can hashing-based techniques provide an efficient FPRAS for DNF counting? In this paper, we provide a positive answer to this question. To this end, we introduce two novel algorithmic techniques: \emph{Symbolic Hashing} and \emph{Stochastic Cell Counting}, along with a new hash family of \emph{Row-Echelon hash functions}. These innovations allow us to design a hashing-based FPRAS for DNF counting of similar complexity (up to polylog factors) as that of prior works. Furthermore, we expect these techniques to have potential applications beyond DNF counting.


title: A Nonparametric Method for Producing Isolines of Bivariate Exceedance Probabilities

abstract: We present a method for drawing isolines indicating regions of equal joint exceedance probability for bivariate data. The method relies on bivariate regular variation, a dependence framework widely used for extremes. This framework enables drawing isolines corresponding to very low exceedance probabilities and these lines may lie beyond the range of the data. The method we utilize for characterizing dependence in the tail is largely nonparametric. Furthermore, we extend this method to the case of asymptotic independence and propose a procedure which smooths the transition from asymptotic independence in the interior to the first-order behavior on the axes. We propose a diagnostic plot for assessing isoline estimate and choice of smoothing, and a bootstrap procedure to visually assess uncertainty.


title: Multi-time correlators in continuous measurement of qubit observables

abstract: We consider multi-time correlators for output signals from linear detectors, continuously measuring several qubit observables at the same time. Using the quantum Bayesian formalism, we show that for unital (symmetric) evolution in the absence of phase backaction, an $N$-time correlator can be expressed as a product of two-time correlators when $N$ is even. For odd $N$, there is a similar factorization, which also includes a single-time average. Theoretical predictions agree well with experimental results for two detectors, which simultaneously measure non-commuting qubit observables.


title: Commuting Graphs of Boundedly Generated Semigroups

abstract: Ara\'ujo, Kinyon and Konieczny (2011) pose several problems concerning the construction of arbitrary commuting graphs of semigroups. We observe that every star-free graph is the commuting graph of some semigroup. Consequently, we suggest modifications for some of the original problems, generalized to the context of families of semigroups with a bounded number of generators, and pose related problems. We construct monomial semigroups with a bounded number of generators, whose commuting graphs have an arbitrary clique number. In contrast to that, we show that the diameter of the commuting graphs of semigroups in a wider class (containing the class of nilpotent semigroups), is bounded by the minimal number of generators plus two. We also address a problem concerning knit degree.


title: Symmetries in the time-averaged dynamics of networks: reducing unnecessary complexity through minimal network models

abstract: Complex networks are the subject of fundamental interest from the scientific community at large. Several metrics have been introduced to characterize the structure of these networks, such as the degree distribution, degree correlation, path length, clustering coefficient, centrality measures etc. Another important feature is the presence of network symmetries. In particular, the effect of these symmetries has been studied in the context of network synchronization, where they have been used to predict the emergence and stability of cluster synchronous states. Here we provide theoretical, numerical, and experimental evidence that network symmetries play a role in a substantially broader class of dynamical models on networks, including epidemics, game theory, communication, and coupled excitable systems. Namely, we see that in all these models, nodes that are related by a symmetry relation show the same time-averaged dynamical properties. This discovery leads us to propose reduction techniques for exact, yet minimal, simulation of complex networks dynamics, which we show are effective in order to optimize the use of computational resources, such as computation time and memory.


title: Modeling and Analysis of SDN Control Applications using Vector Spaces

abstract: Unlike traditional networks which are statically configured, SDN control applications are dynamic and are becoming more heterogeneous and complex. There is a great need for a framework to reason about the behavior of the various SDN applications. To the best of our knowledge, current network modeling frameworks were not designed to incorporate the application logic into their models, and thus can not be used to accurately model the application. In this paper, we suggest the possibility of leveraging the impact which control applications assert on the network information base to reason about the behavior of such applications. Based on that, we propose SDN-VSA, a framework that models SDN control applications as a set of affine transformations in some vector space. Finally, we present an analytical formulation for such framework, and discuss a use-case. For simplicity, we only consider the case of OpenFlow version 1.0.


title: GOE statistics for Anderson models on antitrees and thin boxes in $\mathbb{Z}^3$ with deformed Laplacian

abstract: Sequences of certain finite graphs, antitrees, are constructed along which the Anderson model shows GOE statistics, i.e. a re-scaled eigenvalue process converges to the ${\rm Sine}_1$ process. The Anderson model on the graph is a random matrix being the sum of the adjacency matrix and a random diagonal matrix with independent identically distributed entries along the diagonal. The strength of the randomness stays fixed, there is no re-scaling with matrix size. These considered random matrices giving GOE statistics can also be viewed as random Schr\"odinger operators $\mathcal{P}\Delta+\mathcal{V}$ on thin finite boxes in $\mathbb{Z}^3$ where the Laplacian $\Delta$ is deformed by a projection $\mathcal{P}$ commuting with $\Delta$.


title: Nearly-Kahler 6-Manifolds of Cohomogeneity Two: Local Theory

abstract: We study nearly-Kahler 6-manifolds equipped with a cohomogeneity-two Lie group action for which the principal orbits are coisotropic. If the metric is complete, then we show that this last condition is automatically satisfied, and both the acting Lie group and the principal orbits are finite quotients of $S^3 \times S^1$. We partition the class of such nearly-Ka}hler structures into three types (called I, II, III) and prove a local existence and generality result for each type. Metrics of Types I and II are shown to be incomplete. We also derive a quasilinear elliptic PDE system on the 2-dimensional orbit space which nearly-Kahler structures of Type I must satisfy. Finally, we remark on a relatively simple one-parameter family of Type III structures that turn out to be incomplete metrics that are cohomogeneity-one under the action of a larger group.