title: A General Family of Trimmed Estimators for Robust High-dimensional Data Analysis

abstract: We consider the problem of robustifying high-dimensional structured estimation. Robust techniques are key in real-world applications which often involve outliers and data corruption. We focus on trimmed versions of structurally regularized M-estimators in the high-dimensional setting, including the popular Least Trimmed Squares estimator, as well as analogous estimators for generalized linear models and graphical models, using possibly non-convex loss functions. We present a general analysis of their statistical convergence rates and consistency, and then take a closer look at the trimmed versions of the Lasso and Graphical Lasso estimators as special cases. On the optimization side, we show how to extend algorithms for M-estimators to fit trimmed variants and provide guarantees on their numerical convergence. The generality and competitive performance of high-dimensional trimmed estimators are illustrated numerically on both simulated and real-world genomics data.


title: Secure Repairable Fountain Codes

abstract: In this letter, we provide the construction of repairable fountain codes (RFCs) for distributed storage systems that are information-theoretically secure against an eavesdropper that has access to the data stored in a subset of the storage nodes and the data downloaded to repair an additional subset of storage nodes. The security is achieved by adding random symbols to the message, which is then encoded by the concatenation of a Gabidulin code and an RFC. We compare the achievable code rates of the proposed codes with those of secure minimum storage regenerating codes and secure locally repairable codes.


title: Predictive Coarse-Graining

abstract: We propose a data-driven, coarse-graining formulation in the context of equilibrium statistical mechanics. In contrast to existing techniques which are based on a fine-to-coarse map, we adopt the opposite strategy by prescribing a probabilistic coarse-to-fine map. This corresponds to a directed probabilistic model where the coarse variables play the role of latent generators of the fine scale (all-atom) data. From an information-theoretic perspective, the framework proposed provides an improvement upon the relative entropy method and is capable of quantifying the uncertainty due to the information loss that unavoidably takes place during the CG process. Furthermore, it can be readily extended to a fully Bayesian model where various sources of uncertainties are reflected in the posterior of the model parameters. The latter can be used to produce not only point estimates of fine-scale reconstructions or macroscopic observables, but more importantly, predictive posterior distributions on these quantities. Predictive posterior distributions reflect the confidence of the model as a function of the amount of data and the level of coarse-graining. The issues of model complexity and model selection are seamlessly addressed by employing a hierarchical prior that favors the discovery of sparse solutions, revealing the most prominent features in the coarse-grained model. A flexible and parallelizable Monte Carlo - Expectation-Maximization (MC-EM) scheme is proposed for carrying out inference and learning tasks. A comparative assessment of the proposed methodology is presented for a lattice spin system and the SPC/E water model.


title: Parton Distributions at a 100 TeV Hadron Collider

abstract: The determination of the parton distribution functions (PDFs) of the proton will be an essential input for the physics program of a future 100 TeV hadron collider. The unprecedented center-of-mass energy will require knowledge of PDFs in currently unexplored kinematical regions such as the ultra low-x region or the region of multi-TeV momentum transfers. In this contribution we briefly summarise the studies presented in the PDF section of the upcoming report on "Physics at a 100 TeV pp collider: Standard Model processes". First we map the PDF kinematical coverage in the $(x,Q^2)$ plane, quantify PDF uncertainties, and compute ratios of PDF luminosities between 100 TeV and 14 TeV. Then we show how the extreme kinematics of such collider lead to a number of remarkable PDF-related phenomena such as the top quark as a massless parton, an increased role of photon-initiated processes and the possible need of PDFs with high-energy resummation.


title: Symmetric cumulants and event-plane correlations in Pb+Pb collisions

abstract: The ALICE Collaboration has recently measured the correlations between amplitudes of anisotropic flow in different Fourier harmonics, referred to as symmetric cumulants. We derive approximate relations between symmetric cumulants involving $v_4$ and $v_5$ and the event-plane correlations measured by ATLAS. The validity of these relations is tested using event-by-event hydrodynamic calculations. The corresponding results are in better agreement with ALICE data than existing hydrodynamic predictions. We make quantitative predictions for three symmetric cumulants which are not yet measured.


title: Isoperimetric equalities for rosettes

abstract: In this paper we study the isoperimetric-type equalities for rosettes, i.e. regular closed planar curves with non-vanishing curvature. We find the exact relations between the length and the oriented area of rosettes based on the oriented areas of the Wigner caustic, the Constant Width Measure Set and the Spherical Measure Set. We also study and find new results about the geometry of affine equidistants of rosettes and of the union of rosettes.


title: Spaces of directed paths on pre-cubical sets

abstract: The spaces of directed paths on the geometric realizations of pre-cubical sets, called also $\square$--sets, can be interpreted as the spaces of possible executions of Higher Dimensional Automata, which are models for concurrent computations. In this paper we construct, for a sufficiently good pre-cubical set $K$, a CW-complex $W(K)_v^w$ that is homotopy equivalent to the space of directed paths between given vertices $v$, $w$ of $K$. This construction is functorial with respect to $K$, and minimal among all functorial constructions. Furthermore, explicit formulas for incidence numbers of the cells of $W(K)_v^w$ are provided.


title: Enhanced flexural wave sensing by adaptive gradient-index metamaterials

abstract: Increasing sensitivity and signal to noise ratios of conventional wave sensors is an interesting topic in structural health monitoring, medical imaging, aerospace and nuclear instrumentation. Here, we report the concept of a gradient piezoelectric self-sensing system by integrating shunting circuitry into conventional sensors. By tuning circuit elements properly, both the quality and quantity of the flexural wave measurement data can be significantly increased for new adaptive sensing applications. Through analytical, numerical and experimental studies, we demonstrate that the metamaterial-based sensing system (MBSS) with gradient bending stiffness can be designed by connecting gradient negative capacitance circuits to an array of piezoelectric patches (sensors). We further demonstrate that the proposed system can achieve more than two orders of magnitude amplification of flexural wave signals to overcome the detection limit. This research encompasses fundamental advancements in the MBSS with improved performance and functionalities, and will yield significant advances for a range of applications.


title: On singular Fano varieties with a divisor of Picard number one

abstract: In this paper we study the geometry of mildly singular Fano varieties on which there is an effective prime divisor of Picard number one. Afterwards, we address the case of toric varieties. Finally, we treat the lifting of extremal contractions to universal covering spaces in codimension 1.


title: Efficient Multiple Line-Based Intra Prediction for HEVC

abstract: Traditional intra prediction usually utilizes the nearest reference line to generate the predicted block when considering strong spatial correlation. However, this kind of single line-based method does not always work well due to at least two issues. One is the incoherence caused by the signal noise or the texture of other object, where this texture deviates from the inherent texture of the current block. The other reason is that the nearest reference line usually has worse reconstruction quality in block-based video coding. Due to these two issues, this paper proposes an efficient multiple line-based intra prediction scheme to improve coding efficiency. Besides the nearest reference line, further reference lines are also utilized. The further reference lines with relatively higher quality can provide potential better prediction. At the same time, the residue compensation is introduced to calibrate the prediction of boundary regions in a block when we utilize further reference lines. To speed up the encoding process, this paper designs several fast algorithms. Experimental results show that, compared with HM-16.9, the proposed fast search method achieves 2.0% bit saving on average and up to 3.7%, with increasing the encoding time by 112%.


title: Overview of ISM bands and Software-defined Radio Experimentation

abstract: Wireless systems using low-power wireless communication protocol are rapidly gain popularity in the license-free industrial scientific, and medical (ISM) frequency bands. One such emerging trend in ISM frequency bands is home automation. Historically, all the home devices were once unconnected, today are now being connected either by a wired or wireless connection. The low-power wireless communication protocols enable integration of all the digital home devices into a single system and enhance end user product experience. The rapid prototyping of these end user product using programmable radio system such as software-defined radio (SDR) has been the game-changer in the field of home automation. In this article, we present an overview of popular ISM bands in different regions of the world. Furthermore, we analyze the signals from the wireless electric switch using a software-defined radio in 433 MHz ISM band.


title: Commutative and non-commutative bialgebras of quasi-posets and applications to Ehrhart polynomials

abstract: To any poset or quasi-poset is attached a lattice polytope, whose Ehrhart polynomial we study from a Hopf-algebraic point of view. We use for this two interacting bialgebras on quasi-posets. The Ehrhart polynomial defines a Hopf algebra morphism taking its values in Q[X], we deduce from the interacting bialgebras an algebraic proof of the duality principle, a generalization and a new proof of a result on B-series due to Whright and Zhao, using a monoid of characters on quasi-posets, and a generalization of Faulhaber's formula. We also give non-commutative versions of these results: polynomials are replaced by packed words. We obtain in particular a non-commutative duality principle.


title: 3D Stochastic Geometry Model for Large-Scale Molecular Communication Systems

abstract: Information delivery using chemical molecules is an integral part of biology at multiple distance scales and has attracted recent interest in bioengineering and communication. The collective signal strength at the receiver (i.e., the expected number of observed molecules inside the receiver), resulting from a large number of transmitters at random distances (e.g., due to mobility), can have a major impact on the reliability and efficiency of the molecular communication system. Modeling the collective signal from multiple diffusion sources can be computationally and analytically challenging. In this paper, we present the first tractable analytical model for the collective signal strength due to randomly-placed transmitters, whose positions are modelled as a homogeneous Poisson point process in three-dimensional (3D) space. By applying stochastic geometry, we derive analytical expressions for the expected number of observed molecules at a fully absorbing receiver and a passive receiver. Our results reveal that the collective signal strength at both types of receivers increases proportionally with increasing transmitter density. The proposed framework dramatically simplifies the analysis of large-scale molecular systems in both communication and biological applications.


title: Homogenization for A-quasiconvexity with variable coefficients

abstract: A homogenization result for a family of oscillating integral energies is presented, where the fields under consideration are subjected to first order linear differential constraints depending on the space variable x. The work is based on the theory of A-quasiconvexity with variable coefficients and on two-scale convergence techniques, and generalizes the previously obtained results in the case in which the differential constraints are imposed by means of a linear first order differential operator with constant coefficients. The identification of the relaxed energy in the framework of A-quasiconvexity with variable coefficients is also recovered as a corollary of the homogenization result.


title: A Light-powered, Always-On, Smart Camera with Compressed Domain Gesture Detection

abstract: In this paper we propose an energy-efficient camera-based gesture recognition system powered by light energy for "always on" applications. Low energy consumption is achieved by directly extracting gesture features from the compressed measurements, which are the block averages and the linear combinations of the image sensor's pixel values. The gestures are recognized using a nearest-neighbour (NN) classifier followed by Dynamic Time Warping (DTW). The system has been implemented on an Analog Devices Black Fin ULP vision processor and powered by PV cells whose output is regulated by TI's DC-DC buck converter with Maximum Power Point Tracking (MPPT). Measured data reveals that with only 400 compressed measurements (768x compression ratio) per frame, the system is able to recognize key wake-up gestures with greater than 80% accuracy and only 95mJ of energy per frame. Owing to its fully self-powered operation, the proposed system can find wide applications in "always-on" vision systems such as in surveillance, robotics and consumer electronics with touch-less operation.


title: Statistical mechanics of the Huxley-Simmons model

abstract: The chemomechanical model of Huxley and Simmons (HS) [A. F. Huxley and R. M. Simmons, Nature 233, 533 (1971)] provides a paradigmatic description of mechanically induced collective conformational changes relevant in a variety of biological contexts, from muscles power-stroke and hair cell gating to integrin binding and hairpin unzipping. We develop a statistical mechanical perspective on the HS model by exploiting a formal analogy with a paramagnetic Ising model. We first study the equilibrium HS model with a finite number of elements and compute explicitly its mechanical and thermal properties. To model kinetics, we derive a master equation and solve it for several loading protocols. The developed formalism is applicable to a broad range of allosteric systems with mean-field interactions.


title: A second order minimality condition for a free-boundary problem

abstract: The goal of this paper is to derive in the two-dimensional case necessary and sufficient minimality conditions in terms of the second variation for the functional \[ v\mapsto\int_{\Omega}\big(|\nabla v|^{2}+\chi_{\{v>0\}}Q^{2} \big)\,dx, \] introduced in a classical paper of Alt and Caffarelli. For a special choice of $Q$ this includes water waves. The second variation is obtained by computing the second derivative of the functional along suitable variations of the free boundary. It is proved that the strict positivity of the second variation gives a sufficient condition for local minimality. Also, it is shown that smooth critical points are local minimizers in a small tubular neighborhood of the free-boundary.


title: Spin gravitational resonance and graviton detection

abstract: We develop a gravitational analogue of spin magnetic resonance, called spin gravitational resonance, whereby a gravitational wave interacts with a magnetic field to produce a spin transition. In particular, an external magnetic field separates the energy spin states of a spin-1/2 particle, and the presence of the gravitational wave produces a perturbation in the components of the magnetic field orthogonal to the gravitational wave propagation. In this framework we test Dyson's conjecture that individual gravitons cannot be detected. Although we find no fundamental laws preventing single gravitons being detected with spin gravitational resonance, we show that it cannot be used in practice, in support of Dyson's conjecture.


title: Leavitt path algebras are B\'ezout

abstract: Let $E$ be a directed graph, $K$ any field, and let $L_K(E)$ denote the Leavitt path algebra of $E$ with coefficients in $K$. We show that $L_K(E)$ is a B\'{e}zout ring, i.e., that every finitely generated one-sided ideal of $L_K(E)$ is principal.


title: Abundance of HOCO+ and CO2 in the outer layers of the L1544 prestellar core

abstract: The L1544 prestellar core has been observed as part of the ASAI IRAM Large Program at 3 mm. These observations led to the detection of many complex molecules. In this Letter, we report the detection of two lines, at 85.5 GHz (4,0,4-3,0,3) and 106.9 GHz (5,0,5-4,0,4), respectively, of the protonated carbon dioxide ion, HOCO+. We also report the tentative detection of the line at 100.4 GHz (5,0,5-4,0,4) of DOCO+. The non-LTE analysis of the detected lines shows that the HOCO+ emission originates in the external layer where non-thermal desorption of other species has previously been observed. Its abundance is (5 +/- 2) e-11. Modelling of the chemistry involved in the formation and destruction of HOCO+ provides a gaseous CO2 abundance of 2e-7 (with respect to H2) with an upper limit of 2e-6.


title: A resource-frugal probabilistic dictionary and applications in (meta)genomics

abstract: Genomic and metagenomic fields, generating huge sets of short genomic sequences, brought their own share of high performance problems. To extract relevant pieces of information from the huge data sets generated by current sequencing techniques, one must rely on extremely scalable methods and solutions. Indexing billions of objects is a task considered too expensive while being a fundamental need in this field. In this paper we propose a straightforward indexing structure that scales to billions of element and we propose two direct applications in genomics and metagenomics. We show that our proposal solves problem instances for which no other known solution scales-up. We believe that many tools and applications could benefit from either the fundamental data structure we provide or from the applications developed from this structure.


title: Stochastic discs that roll

abstract: We study a model of rolling particles subject to stochastic fluctuations, which may be relevant in systems of nano- or micro-scale particles where rolling is an approximation for strong static friction. We consider the simplest possible non-trivial system: a linear polymer of three of discs constrained to remain in contact, and immersed in an equilibrium heat bath so the internal angle of the polymer changes due to stochastic fluctuations. We compare two cases: one where the discs can slide relative to each other, and the other where they are constrained to roll, like gears. Starting from the Langevin equations with arbitrary linear velocity constraints, we use formal homogenization theory to derive the overdamped equations that describe the process in configuration space only. The resulting dynamics have the formal structure of a Brownian motion on a Riemannian or sub-Riemannian manifold, depending on if the velocity constraints are holonomic or non-holonomic. We use this to compute the trimer's equilibrium distribution both with, and without, the rolling constraints. Surprisingly, the two distributions are different. We suggest two possible interpretations of this result: either (i) dry friction (or other dissipative, nonequilibrium forces) changes basic thermodynamic quantities like the free energy of a system, a statement that could be tested experimentally, or (ii) as a lesson in modeling rolling or friction more generally as a velocity constraint when stochastic fluctuations are present. In the latter case, we speculate there could be a "roughness" entropy whose inclusion as an effective force could compensate the constraint and preserve classical Boltzmann statistics. Regardless of the interpretation, our calculation shows the word "rolling" must be used with care when stochastic fluctuations are present.

categories: math-ph cond-mat.mes-hall cond-mat.stat-mech math.MP


title: Near Optimal Channel Assignment for Interference Mitigation in Wireless Mesh Networks

abstract: In multi-radio multi-channel (MRMC) WMNs, interference alleviation is affected through several network design techniques e.g., channel assignment (CA), link scheduling, routing etc., intelligent CA schemes being the most effective tool for interference mitigation. CA in WMNs is an NP-Hard problem, and makes optimality a desired yet elusive goal in real-time deployments which are characterized by fast transmission and switching times and minimal end-to-end latency. The trade-off between optimal performance and minimal response times is often achieved through CA schemes that employ heuristics to propose efficient solutions. WMN configuration and physical layout are also crucial factors which decide network performance, and it has been demonstrated in numerous research works that rectangular/square grid WMNs outperform random or unplanned WMN deployments in terms of network capacity, latency, and network resilience. In this work, we propose a smart heuristic approach to devise a near-optimal CA algorithm for grid WMNs (NOCAG). We demonstrate the efficacy of NOCAG by evaluating its performance against the minimal-interference CA generated through a rudimentary brute-force technique (BFCA), for the same WMN configuration. We assess its ability to mitigate interference both, theoretically (through interference estimation metrics) and experimentally (by running rigorous simulations in NS-3). We demonstrate that the performance of NOCAG is almost as good as the BFCA, at a minimal computational overhead of O(n) compared to the exponential of BFCA.


title: Targets, local weak $\sigma$-Gibbs measures and a generalized Bowen dimension formula

abstract: For a dynamical system, we study the set of points $\cal W$ whose orbit approximates any chosen point at certain specified rates. Our basic setting is that of left shift acting on topological Markov chains endowed with a local weak Gibbs measure. Our rates of recurrence are so fast that the corresponding set $\cal W$ has measure zero, but we obtain a generalized Bowen formula for Carath\'eodory dimension. For the case of Markov transformations with countable partition and big image (BI) property a Bowen-type formula is obtained for the Hausdorff dimension of those exceptional sets. In particular, we apply our general results to Gauss and Luroth maps, a not Bernoulli modification of Gauss map and some inner functions. Since we only require the existence of weak Gibbs measures we can deal with non H\"older potentials and we can also consider intermittent systems as the Manneville-Pomeau map.


title: Aerial image geolocalization from recognition and matching of roads and intersections

abstract: Aerial image analysis at a semantic level is important in many applications with strong potential impact in industry and consumer use, such as automated mapping, urban planning, real estate and environment monitoring, or disaster relief. The problem is enjoying a great interest in computer vision and remote sensing, due to increased computer power and improvement in automated image understanding algorithms. In this paper we address the task of automatic geolocalization of aerial images from recognition and matching of roads and intersections. Our proposed method is a novel contribution in the literature that could enable many applications of aerial image analysis when GPS data is not available. We offer a complete pipeline for geolocalization, from the detection of roads and intersections, to the identification of the enclosing geographic region by matching detected intersections to previously learned manually labeled ones, followed by accurate geometric alignment between the detected roads and the manually labeled maps. We test on a novel dataset with aerial images of two European cities and use the publicly available OpenStreetMap project for collecting ground truth roads annotations. We show in extensive experiments that our approach produces highly accurate localizations in the challenging case when we train on images from one city and test on the other and the quality of the aerial images is relatively poor. We also show that the the alignment between detected roads and pre-stored manual annotations can be effectively used for improving the quality of the road detection results.


title: The 750 GeV diphoton excess and SUSY

abstract: The LHC experiments ATLAS and CMS have reported an excess in the diphoton spectrum at \sim 750 GeV. At the same time the motivation for Supersymmetry (SUSY) remains unbowed. Consequently, we review briefly the proposals to explain this excess in SUSY, focusing on "pure" (N)MSSM solutions. We then review in more detail a proposal to realize this excess within the NMSSM. In this particular scenario a Higgs boson with mass around 750 GeV decays to two light pseudo-scalar Higgs bosons. Via mixing with the pion these pseudo-scalars decay into a pair of highly collimated photons, which are identified as one photon, thus resulting in the observed signal.


title: Theano-MPI: a Theano-based Distributed Training Framework

abstract: We develop a scalable and extendable training framework that can utilize GPUs across nodes in a cluster and accelerate the training of deep learning models based on data parallelism. Both synchronous and asynchronous training are implemented in our framework, where parameter exchange among GPUs is based on CUDA-aware MPI. In this report, we analyze the convergence and capability of the framework to reduce training time when scaling the synchronous training of AlexNet and GoogLeNet from 2 GPUs to 8 GPUs. In addition, we explore novel ways to reduce the communication overhead caused by exchanging parameters. Finally, we release the framework as open-source for further research on distributed deep learning


title: The BMO-Dirichlet problem for elliptic systems in the upper-half space and quantitative characterizations of VMO

abstract: We prove that for any homogeneous, second order, constant complex coefficient elliptic system $L$, the Dirichlet problem in $\mathbb{R}^{n}_{+}$ with boundary data in BMO is well-posed in the class of functions $u$ with $d\mu_u(x',t):=|\nabla u(x',t)|^2\,t\,dx'dt$ being a Carleson measure. We establish a Fatou type theorem guaranteeing the existence of the pointwise nontangential boundary trace for smooth null-solutions $u$ of such systems satisfying the said Carleson measure condition. These imply that BMO can be characterized as the collection of nontangential pointwise traces of smooth null-solutions $u$ to the elliptic system $L$ with the property that $\mu_u$ is a Carleson measure. We establish a regularity result for the BMO-Dirichlet problem in the upper-half space: the nontangential pointwise trace of any given smooth null-solutions of $L$ satisfying the above Carleson measure condition belongs to Sarason's space VMO if and only if $\mu_u$ satsifies a vanishing Carleson measure condition. Moreover, we obtain the well-posedness of the Dirichlet problems when the boundary data are prescribed in Morrey-Campanato and solutions are required to satisfy a vanishing Carleson measure condition of fractional order. As a consequence, we characterize the space VMO as the closure in BMO of classes of smooth functions contained in BMO within which uniform continuity may be suitably quantified (such as the class of smooth functions satisfying a H\"older or Lipschitz condition), improving Sarason's classical result describing VMO as the closure in BMO of the space of uniformly continuous functions with bounded mean oscillations. Finally, we show that any Calder\'on-Zygmund operator $T$ satisfying $T(1)=0$ extends as a bounded linear operator in $\mathrm{VMO}$, and characterize the membership to $\mathrm{VMO}$ via the action of various classes of singular integral operators.


title: Quantum Communication Between Remote Mechanical Resonators

abstract: Mechanical resonators represent one of the most promising candidates to mediate the interaction between different quantum technologies, bridging the gap between efficient quantum computation and long-distance quantum communication. In this letter, we introduce a novel interferometric scheme where the interaction of a mechanical resonator with input/output quantum pulses is controlled by an independent classical drive. We design protocols for state teleportation and direct quantum state transfer, between distant mechanical resonators. The proposed device, feasible with state-of-the-art technology, can serve as building block for the implementation of long-distance quantum networks of mechanical resonators.


title: Discussion of "Causal inference using invariant prediction: identification and confidence intervals" by Peters, B\"uhlmann and Meinshausen

abstract: Contribution to the discussion of the paper "Causal inference using invariant prediction: identification and confidence intervals" by Peters, B\"uhlmann and Meinshausen, to appear in the Journal of the Royal Statistical Society, Series B.


title: Planets, debris and their host metallicity correlations

abstract: Recent observations of debris discs, believed to be made up of remnant planetesimals, brought a number of surprises. Debris disc presence does not correlate with the host star's metallicity, and may anti-correlate with the presence of gas giant planets. These observations contradict both assumptions and predictions of the highly successful Core Accretion model of planet formation. Here we explore predictions of the alternative Tidal Downsizing (TD) scenario of planet formation. In TD, small planets and planetesimal debris is made only when gas fragments, predecessors of giant planets, are tidally disrupted. We show that these disruptions are rare in discs around high metallicity stars but release more debris per disruption than their low [M/H] analogs. This predicts no simple relation between debris disc presence and host star's [M/H], as observed. A detected gas giant planet implies in TD that its predecessor fragment was not disputed, potentially explaining why DDs are less likely to be found around stars with gas giants. Less massive planets should correlate with DD presence, and sub-Saturn planets ($M_{p} \sim 50 \M_{\oplus}$) should correlate with DD presence stronger than sub-Neptunes ($M_{p} \leq 15 \M_{\oplus}$). These predicted planet-DD correlations will be diluted and weakened in observations by planetary systems' long term evolution and multi-fragment effects neglected here. Finally, although presently difficult to observe, DDs around M dwarf stars should be more prevalent than around Solar type stars.


title: Sharp Degree Bounds for Sum-of-Squares Certificates on Projective Curves

abstract: Given a real projective curve with homogeneous coordinate ring R and a nonnegative homogeneous element f in R, we bound the degree of a nonzero homogeneous sum-of-squares g in R such that the product fg is again a sum of squares. Better yet, our degree bounds only depend on geometric invariants of the curve and we show that there exist smooth curves and nonnegative elements for which our bounds are sharp. We deduce the existence of a multiplier g from a new Bertini Theorem in convex algebraic geometry and prove sharpness by deforming rational Harnack curves on toric surfaces. Our techniques also yield similar bounds for multipliers on surfaces of minimal degree, generalizing Hilbert's work on ternary forms.


title: Forming first-ranked early-type galaxies through hierarchical dissipationless merging

abstract: We have developed a computationally competitive N-body model of a previrialized aggregation of galaxies in a flat LambdaCDM universe to assess the role of the multiple mergers that take place during the formation stage of such systems in the configuration of the remnants assembled at their centres. An analysis of a suite of 48 simulations of low-mass forming groups (of about 1E13 solar masses) demonstrates that the gravitational dynamics involved in their hierarchical collapse is capable of creating realistic first-ranked galaxies without the aid of dissipative processes. Our simulations indicate that the brightest group galaxies (BGGs) constitute a distinct population from other group members, sketching a scenario in which the assembly path of these objects is dictated largely by the formation of their host system. We detect significant differences in the distribution of Sersic indices and total magnitudes, as well as a luminosity gap between BGGs and the next brightest galaxy that is positively correlated with the total luminosity of the parent group. Such gaps arise from both the grow of BGGs at the expense of lesser companions and the decrease in the relevance of second-ranked objects in equal measure. This results in a dearth of intermediate-mass galaxies which explains the characteristic central dip detected in their luminosity functions in dynamically young galaxy aggregations. The fact that the basic global properties of our BGGs define a thin mass fundamental plane strikingly similar to that followed giant early-type galaxies in the local universe reinforces confidence in the results obtained.


title: A note on implementation of boundary variation diminishing algorithm to high-order local polynomial-based schemes

abstract: A novel approach for selecting appropriate reconstructions is implemented to the hyperbolic conservation laws in the high-order local polynomial-based framework, e.g., the discontinuous Galerkin (DG) and flux reconstruction (FR) schemes. The high-order polynomial approximation generally fails to correctly capture a strong discontinuity inside a cell due to the Runge phenomenon, which is replaced by more stable approximation on the basis of a troubled-cell indicator such as that used in the total variation bounded (TVB) limiter. This paper examines the applicability of a new algorithm, so-called boundary variation diminishing (BVD) reconstruction, to the weighted essentially non-oscillatory (WENO) methodology in the FR framework including the nodal type DG method. The BVD reconstruction adaptively chooses a proper approximation for the solution function so as to minimize the jump between values at the left- and right-side of cell boundaries without any ad hoc constant such as the TVB parameter. Several numerical tests are conducted for a linear advection equation as well as the selection of an appropriate reconstruction, where the results of the BVD algorithm are comparable to those using the conventional TVB limiter. Note that the present work is limited to third or lower order polynomials, leaving the implementations for higher-order schemes a future work.


title: Polarization of the Sunyaev-Zel'dovich effect: relativistic imprint of thermal and non-thermal plasma

abstract: [Abridged] Inverse Compton scattering of CMB fluctuations off cosmic electron plasma generates a polarization of the associated Sunyaev-Zel'dovich (SZ) effect. This signal has been studied so far mostly in the non-relativistic regime and for a thermal electron population and, as such, has limited astrophysical applications. Partial attempts to extend this calculation for a thermal electron plasma in the relativistic regime have been done but cannot be applied to a general relativistic electron distribution. Here we derive a general form of the SZ effect polarization valid in the full relativistic approach for both thermal and non-thermal electron plasmas, as well as for a generic combination of various electron population co-spatially distributed in the environments of galaxy clusters or radiogalaxy lobes. We derive the spectral shape of the Stokes parameters induced by the IC scattering of every CMB multipole, focusing on the CMB quadrupole and octupole that provide the largest detectable signals in galaxy clusters. We found that the CMB quadrupole induced Stoke parameter Q is always positive with a maximum amplitude at 216 GHz which increases slightly with increasing cluster temperature. The CMB octupole induced Q spectrum shows, instead, a cross-over frequency which depends on the cluster electron temperature, or on the minimum momentum p_1 as well as on the power-law spectral index of a non-thermal electron population. We discuss some possibilities to disentangle the quadrupole-induced Q spectrum from the octupole-induced one which allow to measure these quantities through the SZ effect polarization. We finally apply our model to the realistic case of the Bullet cluster and derive the visibility windows of the total, quandrupole-induced and octupole-induced Stoke parameter Q in the frequency ranges accessible to SKA, ALMA, MILLIMETRON and CORE++ experiments.


title: Ultra-forward particle production from CGC+Lund fragmentation

abstract: We present an analysis of data on single inclusive pion production measured by the LHCf collaboration in high-energy proton-proton and proton-nucleus at ultra-forward rapidities, $8.8 \leq\! y\leq \!10.8$. We also analyse forward RHIC data for calibration purposes. Our analysis relies on the use of a Monte Carlo event generator that combines a perturbative description of the elementary scattering process at partonic level based on the hybrid formalism of the Color Glass Condensate with an implementation of hadronization in the framework of the Lund string fragmentation model. This procedure allows us to reach values of the momenta of the produced particles as low as detected experimentally $p_t\sim0.1$ GeV. We achieve a good description of single inclusive spectra of charged particles and neutral pions at RHIC and the LHC respectively, and nuclear modification factors for proton-lead collisions at the LHC. Our results add evidence to the idea that particle production in the domain of very small Bjorken-$x$ is dominated by the saturation effects encoded in the unintegrated gluon distribution of the target. Being forward particle production of key importance in the development of air showers, we stress that this approach allows for a theoretically controlled extrapolation of our results to the scale of ultra-high energy cosmic rays, thus serving as starting point for future works on this topic.


title: Gauge Invariance and the Quantum Metric Tensor

abstract: The quantum metric tensor was introduced for defining the distance in the parameter space of a system. However, it is also useful for other purposes, like predicting quantum phase transitions. Due to the physical information this tensor provides, its gauge independence sounds reasonable. More over, its original construction was made by looking for this gauge independence. The aim of this paper, however, is to prove that the quantum metric tensor does depend on the gauge. In addition, a real gauge invariant quantum metric tensor is introduced. In this paper, the gauge dependence is explicitly shown by computing the quantum metric tensor of the Landau problem in different gauges. Then, a real gauge independent metric tensor is proposed and computed for the same Landau problem. Since the gauge dependence has not been observed before, the results of this paper might lead to a new study of topics that were believed to be completely understood.

categories: quant-ph cond-mat.stat-mech math-ph math.MP


title: Distributed Gauss-Newton Method for AC State Estimation: A Belief Propagation Approach

abstract: In this paper, we propose a solution to an AC state estimation problem in electric power systems using a fully distributed Gauss-Newton method. The proposed method is placed within the context of factor graphs and belief propagation algorithms and closed-form expressions for belief propagation messages exchanged along the factor graph are derived. The obtained algorithm provides the same solution as the conventional weighted least-squares state estimation. Using a simple example, we provide a step-by-step presentation of the proposed algorithm. Finally, we discuss the convergence behaviour using the IEEE 14 bus test case.


title: Compactifications of 5d SCFTs with a twist

abstract: We study the compactification of 5d SCFTs to 4d on a circle with a twist in a discrete global symmetry element of these SCFTs. We present evidence that this leads to various 4d N=2 isolated SCFTs. These include many known theories as well as seemingly new ones. The known theories include the recently discovered rank 1 SU(4) SCFT and its mass deformations. One application of the new SCFTs is in the dual descriptions of the 4d gauge theory SU(N)+1S+(N-2)F. Also interesting is the appearance of a theory with rank 1 and $F_4$ global symmetry.


title: Primordial inhomogeneities from massive defects during inflation

abstract: We consider the imprints of local massive defects, such as a black hole or a massive monopole, during inflation. The massive defect breaks the background homogeneity. We consider the limit that the physical Schwarzschild radius of the defect is much smaller than the inflationary Hubble radius so a perturbative analysis is allowed. The inhomogeneities induced in scalar and gravitational wave power spectrum are calculated. We obtain the amplitudes of dipole, quadrupole and octupole anisotropies in curvature perturbation power spectrum and identify the relative configuration of the defect to CMB sphere in which large observable dipole asymmetry can be generated. We observe a curious reflection symmetry in which the configuration where the defect is inside the CMB comoving sphere has the same inhomogeneous variance as its mirror configuration where the defect is outside the CMB sphere.


title: Random walks on hyperplane arrangements and stopping times

abstract: Consider a real hyperplane arrangement and let $\mathcal{C}$ denote the occurring chambers. Bidigare, Hanlon and Rockmore introduced a Markov chain on $\mathcal{C}$ which is a generalization of some card shuffling models used in computer science, biology and card games. This paper introduces strong stationary arguments for this Markov chain, which provide explicit bounds for the separation distance.


title: Optimal Pricing under Mixed Logit Choice

abstract: In this paper we consider the problem of pricing multiple differentiated products. This is challenging as a price change in one product, not only changes the demand of that particular product, but also the demand for the other products. To address this problem, customer choice models have recently been introduced as these are capable of describing customer choice behavior across differentiated products. In the present paper the objective is to obtain the revenue-maximizing prices when the customer's decision making process is modelled according to a particular customer choice model, namely the mixed logit model. The main advantage of using the mixed logit model, also known as the random coefficients logit model, for this purpose is its flexibility. In the single-product case we establish log-concavity of the optimization problem under certain regularity conditions. In addition, in the multi-product case, we present the results of our extensive numerical experiments. These suggest that the mixed logit model, by taking unobserved customer heterogeneity and flexible substitution patterns into account, can significantly improve the attainable revenue.


title: Observational selection biases in time-delay strong lensing and their impact on cosmography

abstract: Inferring cosmological parameters from time-delay strong lenses requires a significant investment of telescope time; it is therefore tempting to focus on the systems with the brightest sources, the highest image multiplicities and the widest image separations. We investigate if this selection bias can influence the properties of the lenses studied and the cosmological parameters that are inferred. Using a population of lenses with ellipsoidal powerlaw density profiles, we build a sample of double and quadruple image systems. Assuming reasonable thresholds on image separation and flux, based on current lens monitoring campaigns, we find that the typical density profile slopes of monitorable lenses are significantly shallower than the input ensemble. From a sample of quadruple image lenses we find that this selection function can introduce a 3.5% bias on the inferred time-delay distances if the ensemble of deflector properties is used as a prior for a cosmographical analysis. This bias remains at the 2.4% level when high resolution imaging of the quasar host is used to precisely infer the density profiles of individual lenses. We also investigate if the lines-of-sight for monitorable strong lenses are biased. After adding external convergence, $\kappa$, and shear to our lens population we find that the expectation value for $\kappa$ is increased by 0.004 and 0.009 for doubles and quads respectively. $\kappa$ is degenerate with the value of $H_0$ inferred from time delays; fortunately the shift in $\kappa$ only induces a 0.9 (0.4) percent bias on $H_0$ for quads (doubles). We therefore conclude that whilst the properties of typical quasar lenses and their lines-of-sight do deviate from the global population, the total magnitude of this effect is likely a subdominant effect for current analyses, but has the potential to be a major systematic for samples of $\sim$25 or more lenses.


title: Unveiling the formation route of the largest galaxies in the universe

abstract: Observational evidence indicates that the role of gas is secondary to that of gravity in the formation of the most luminous spheroids inhabiting the centres of galaxy associations, as originally conjectured in the late 80's/early 90's. However, attempts to explain the origin of the Fundamental Plane (FP) of massive early-type galaxies (ETGs) -- a tilted version of the scaling relation connecting the size, velocity dispersion and mass of virialized homologous systems -- based on sequences of pairwise mergers, have systematically concluded that dissipation cannot be ignored. We use controlled simulations of the previrialization stage of galaxy groups to show that multiple collisionless merging is capable of creating realistic first-ranked galaxies. Our mock remnants define a thin FP that perfectly fits data from all kinds of giant ETGs in the local volume, showing the existence of a unified relationship for these systems. High-ranked galaxies occupy in the FP different areas than standard objects, a segregation which is viewed essentially as zero-point offsets in the 2D correlations arising from standard projections of this plane. Our findings make a strong case for considering hierarchical dissipationless merging a viable route for the formation of the largest galaxies in the universe.


title: Fundamental Groups of Simplicial Complexes

abstract: We define two different simplicial complexes, the common divisor simplicial complex and the prime divisor simplicial complex, from a set of integers, and explore their similarities. We will define a map between the two simplicial complexes, and use this map to show that for any set of integers, the fundamental groups of the resulting simplicial complexes are isomorphic.


title: Cohomological Hasse principle for schemes over valuation rings of higher dimensional local fields

abstract: K. Kato's conjecture about the cohomological Hasse principle for regular connected schemes $\mathfrak X$ which are flat and proper over the complete discrete valuation rings $\mathcal O_N$ of higher local fields $F_N$ is proven. This generalizes the work of M. Kerz, S. Saito and U. Jannsen for finite fields to the case of all higher local fields. For that purpose a $p$-alteration theorem for the local uniformization of schemes over valuation rings of arbitrary finite rank and a corresponding Bertini theorem is developed extending the results of O. Gabber, J. deJong, L. Illusie, M. Temkin, S. Saito, U. Jannsen to the non-noetherian world. As an application it is shown that certain motivic cohomology groups of varieties over higher local fields are finite. This is one of the rare cases where such a result could be shown for schemes without finite or separably closed residue fields. Furthermore, it will be derived that the kernels of the reciprocity map $\rho^X : \text{SK}_N(X) \to \pi_1^\text{ab}(X)$ and norm map $N_{X|F}: \text{SK}_N(X) \to K_N^M(F_N)$ modulo maximal $p'$-divisible subgroups are finite for regular $X$ which are proper over a higher local field $F_N$ with final residue characteristic $p$. This generalizes results of S. Bloch, K. Kato, U. Jannsen, S. Saito from varieties over finite and local fields to varieties over higher local fields, both of arbitrary dimensions.


title: Modeling non-stationary, non-axisymmetric heat patterns in DIII-D tokamak

abstract: Non-axisymmetric stationary magnetic perturbations lead to the formation of homoclinic tangles near the divertor magnetic saddle in tokamak discharges. These tangles intersect the divertor plates in static helical structures that delimit the regions reached by open magnetic field lines reaching the plasma column and leading the charged particles to the strike surfaces by parallel transport. In this article we introduce a non-axisymmetric rotating magnetic perturbation to model the time development of the three-dimensional magnetic field of a single-null DIII-D tokamak discharge developing a rotating tearing mode. The stable and unstable manifolds of the asymmetric magnetic saddle are calculated through an adaptive method providing the manifold cuts at a given poloidal plane and the strike surfaces. For the modeled shot, the experimental heat pattern and its time development are well described by the rotating unstable manifold, indicating the emergence of homoclinic lobes in a rotating frame due to the plasma instabilities. In the model it is assumed that the magnetic field is created by a stationary axisymmetric plasma current and a set of rotating internal helical filamentary currents. The currents in the filaments are adjusted to match the waveforms of the magnetic probes at the mid-plane and the rotating magnetic field is introduced as a perturbation to the axisymmetric field obtained from a Grad-Shafranov equilibrium reconstruction code.


title: Distributed Sequence Memory of Multidimensional Inputs in Recurrent Networks

abstract: Recurrent neural networks (RNNs) have drawn interest from machine learning researchers because of their effectiveness at preserving past inputs for time-varying data processing tasks. To understand the success and limitations of RNNs, it is critical that we advance our analysis of their fundamental memory properties. We focus on echo state networks (ESNs), which are RNNs with simple memoryless nodes and random connectivity. In most existing analyses, the short-term memory (STM) capacity results conclude that the ESN network size must scale linearly with the input size for unstructured inputs. The main contribution of this paper is to provide general results characterizing the STM capacity for linear ESNs with multidimensional input streams when the inputs have common low-dimensional structure: sparsity in a basis or significant statistical dependence between inputs. In both cases, we show that the number of nodes in the network must scale linearly with the information rate and poly-logarithmically with the ambient input dimension. The analysis relies on advanced applications of random matrix theory and results in explicit non-asymptotic bounds on the recovery error. Taken together, this analysis provides a significant step forward in our understanding of the STM properties in RNNs.


title: Enabling science with Gaia observations of naked-eye stars

abstract: ESA's Gaia space astrometry mission is performing an all-sky survey of stellar objects. At the beginning of the nominal mission in July 2014, an operation scheme was adopted that enabled Gaia to routinely acquire observations of all stars brighter than the original limit of G~6, i.e. the naked-eye stars. Here, we describe the current status and extent of those observations and their on-ground processing. We present an overview of the data products generated for G<6 stars and the potential scientific applications. Finally, we discuss how the Gaia survey could be enhanced by further exploiting the techniques we developed.


title: Existence of Ground State Eigenvalues for the Spin-Boson Model with Critical Infrared Divergence and Multiscale Analysis

abstract: A two-level atom coupled to the radiation field is studied. First principles in physics suggest that the coupling function, representing the interaction between the atom and the radiation field, behaves like $\vert k \vert^{- 1/2}$, as the photon momentum k tends to zero. Previous results on non-existence of ground state eigenvalues suggest that in the most general case binding does not occur in the spin-boson model, i.e., the minimal energy of the atom-photon system is not an eigenvalue of the energy operator. Hasler and Herbst have shown [12], however, that under the additional hypothesis that the coupling function be off-diagonal -which is customary to assume-binding does indeed occur. In this paper an alternative proof of binding in case of off-diagonal coupling is given, i.e., it is proven that, if the coupling function is off-diagonal, the ground state energy of the spin-boson model is an eigenvalue of the Hamiltonian. We develop a multiscale method that can be applied in the situation we study, identifying a new key symmetry operator which we use to demonstrate that the most singular terms appearing in the multiscale analysis vanish.


title: Enhancement of fusion at near and sub-barrier energies for neutron-rich light nuclei

abstract: Measurement of the fusion cross-section for neutron-rich light nuclei is crucial in ascertaining if fusion of these nuclei occurs in the outer crust of a neutron star. We have therefore measured the fusion excitation function at near-barrier energies for the 19O + 12C system and compared the experimental results with the fusion excitation function of 18O + 12C and 16O + 12C. In the experiment a beam of 19O, produced via the 18O(d,p) reaction, was incident on a 12C target at energies near the Coulomb barrier. Evaporation residues produced in fusion of 18,19O ions with 12C target nuclei were detected with good geometric efficiency and identified by measuring their energy and time-of-flight. A significant enhancement in the fusion probability of 19O ions with a 12C target as compared to 18O ions is observed. The significantly larger cross-sections observed at near barrier energies are not predicted by a static model of fusion for 19O + 12C indicating that dynamics play an important role in the fusion of neutron-rich light nuclei.


title: Benign-Malignant Lung Nodule Classification with Geometric and Appearance Histogram Features

abstract: Lung cancer accounts for the highest number of cancer deaths globally. Early diagnosis of lung nodules is very important to reduce the mortality rate of patients by improving the diagnosis and treatment of lung cancer. This work proposes an automated system to classify lung nodules as malignant and benign in CT images. It presents extensive experimental results using a combination of geometric and histogram lung nodule image features and different linear and non-linear discriminant classifiers. The proposed approach is experimentally validated on the LIDC-IDRI public lung cancer screening thoracic computed tomography (CT) dataset containing nodule level diagnostic data. The obtained results are very encouraging correctly classifying 82% of malignant and 93% of benign nodules on unseen test data at best.


title: 'Spectro-temporal' variabilities and possible physical mechanism for Jet ejections

abstract: In this paper, we attempt to find out the `spectro-temporal' characteristics during the jet ejection, of several outbursting Galactic black hole sources based on RXTE-PCA/HEXTE data in the energy band of 2 - 100 keV. We present results of detailed analysis of these sources during the rising phase of their outburst, whenever simultaneous or near-simultaneous X-ray and Radio observations are `available'. We find that before the peak radio flare (transient jet) a few of the sources (in addition to those reported earlier) exhibit `local' softening within the soft intermediate state itself. Except the duration, all the properties of the `local' softening (QPO not observed, reduction in total rms, soft spectra) are observed to be similar to the canonical soft state. We find similar `local' softening for the recent outburst of V404 Cyg also based on SWIFT observations. Fast changes in the `spectro-temporal' properties during the `local' softening implies that it may not be occurring due to change in Keplerian accretion rate. We discuss these results in the framework of the magnetized two component advective flow model.


title: Automated Parton-Shower Variations in Pythia 8

abstract: In the era of precision physics measurements at the LHC, efficient and exhaustive estimations of theoretical uncertainties play an increasingly crucial role. In the context of Monte Carlo (MC) event generators, the estimation of such uncertainties traditionally requires independent MC runs for each variation, for a linear increase in total run time. In this work, we report on an automated evaluation of the dominant (renormalization-scale and non-singular) perturbative uncertainties in the PYTHIA 8 event generator, with only a modest computational overhead. Each generated event is accompanied by a vector of alternative weights (one for each uncertainty variation), with each set separately preserving the total cross section. Explicit scale-compensating terms can be included, reflecting known coefficients of higher-order splitting terms and reducing the effect of the variations. The formalism also allows for the enhancement of rare partonic splittings, such as $g \to b \bar{b}$ and $q\to q \gamma$, to obtain weighted samples enriched in these splittings while preserving the correct physical Sudakov factors.


title: User Performance in Small Cells Networks with Inter-Cell Mobility

abstract: We analyze the impact of intra-cell mobility on user performance in dense networks such as that enabled by LTE-A and 5G. To this end, we consider a homogeneous network of small cells and first show how to reduce the evaluation of user performance to the case of a single representative cell. We then propose simple analytical models that capture mobility through the distribution of the residual sojourn time of mobile users in the cell. An approximate model, based on Quasi-Stationary (QS) assumptions, is developed in order to speed up computation in the Markovian framework. We use these models to derive the average throughput of both mobile and static users, along with the probability of handover for mobile users. Numerical evaluation and simulation results are provided to assess the accuracy of the proposed models. We show, in particular, that both classes of users benefit from a throughput gain induced by the "opportunistic" displacement of mobile users among cells.


title: Can an interdisciplinary field contribute to one of the parent disciplines from which it emerged?

abstract: In the light of contemporary discussions of inter and transdisciplinarity, this paper approaches econophysics and sociophysics to seek a response to the question -- whether these interdisciplinary fields could contribute to physics and economics. Drawing upon the literature on history and philosophy of science, the paper argues that the two way traffic between physics and economics has a long history and this is likely to continue in the future.


title: On the 4D generalized Proca action for an Abelian vector field

abstract: We summarize previous results on the most general Proca theory in 4 dimensions containing only first-order derivatives in the vector field (second-order at most in the associated St\"uckelberg scalar) and having only three propagating degrees of freedom with dynamics controlled by second-order equations of motion. Discussing the Hessian condition used in previous works, we conjecture that, as in the scalar galileon case, the most complete action contains only a finite number of terms with second-order derivatives of the St\"uckelberg field describing the longitudinal mode, which is in agreement with the results of JCAP 1405, 015 (2014) and Phys. Lett. B 757, 405 (2016) and complements those of JCAP 1602, 004 (2016). We also correct and complete the parity violating sector, obtaining an extra term on top of the arbitrary function of the field $A_\mu$, the Faraday tensor $F_{\mu \nu}$ and its Hodge dual $\tilde{F}_{\mu \nu}$.


title: Two-gap superconductivity in LaNiGa$_2$ with non-unitary triplet pairing and even parity gap symmetry

abstract: The nature of the pairing states of superconducting LaNiC$_2$ and LaNiGa$_2$ has to date remained a puzzling question. Broken time reversal symmetry has been observed in both compounds and a group theoretical analysis implies a non-unitary triplet pairing state. However all the allowed non-unitary triplet states have nodal gap functions but most thermodynamic and NMR measurements indicate fully gapped superconductivity in LaNiC$_2$. Here we probe the gap symmetry of LaNiGa$_2$ by measuring the London penetration depth, specific heat and upper critical field. These measurements demonstrate two-gap nodeless superconductivity in LaNiGa$_2$, suggesting that this is a common feature of both compounds. These results allow us to propose a novel triplet superconducting state, where the pairing occurs between electrons of the same spin, but on different orbitals. In this case the superconducting wavefunction has a triplet spin component but isotropic even parity gap symmetry, yet the overall wavefunction remains antisymmetric under particle exchange. This model leads to a nodeless two-gap superconducting state which breaks time reversal symmetry, and therefore accounts well for the seemingly contradictory experimental results.


title: MobileAppScrutinator: A Simple yet Efficient Dynamic Analysis Approach for Detecting Privacy Leaks across Mobile OSs

abstract: Smartphones, the devices we carry everywhere with us, are being heavily tracked and have undoubtedly become a major threat to our privacy. As "tracking the trackers" has become a necessity, various static and dynamic analysis tools have been developed in the past. However, today, we still lack suitable tools to detect, measure and compare the ongoing tracking across mobile OSs. To this end, we propose MobileAppScrutinator, based on a simple yet efficient dynamic analysis approach, that works on both Android and iOS (the two most popular OSs today). To demonstrate the current trend in tracking, we select 140 most representative Apps available on both Android and iOS AppStores and test them with MobileAppScrutinator. In fact, choosing the same set of apps on both Android and iOS also enables us to compare the ongoing tracking on these two OSs. Finally, we also discuss the effectiveness of privacy safeguards available on Android and iOS. We show that neither Android nor iOS privacy safeguards in their present state are completely satisfying.


title: On $M$ -- terms approximations Besov classes in Lorentz spaces

abstract: In this paper we consider Lorentz space with a mixed norm of periodic functions of many variables. We obtain the exact estimation of the best M-term approximations of Nikol'ski's, Besov's classes in the Lorentz space with the mixed norm.


title: Pairwise Decomposition of Image Sequences for Active Multi-View Recognition

abstract: A multi-view image sequence provides a much richer capacity for object recognition than from a single image. However, most existing solutions to multi-view recognition typically adopt hand-crafted, model-based geometric methods, which do not readily embrace recent trends in deep learning. We propose to bring Convolutional Neural Networks to generic multi-view recognition, by decomposing an image sequence into a set of image pairs, classifying each pair independently, and then learning an object classifier by weighting the contribution of each pair. This allows for recognition over arbitrary camera trajectories, without requiring explicit training over the potentially infinite number of camera paths and lengths. Building these pairwise relationships then naturally extends to the next-best-view problem in an active recognition framework. To achieve this, we train a second Convolutional Neural Network to map directly from an observed image to next viewpoint. Finally, we incorporate this into a trajectory optimisation task, whereby the best recognition confidence is sought for a given trajectory length. We present state-of-the-art results in both guided and unguided multi-view recognition on the ModelNet dataset, and show how our method can be used with depth images, greyscale images, or both.


title: Probing interstellar turbulence in cirrus with deep optical imaging: no sign of energy dissipation at 0.01 pc scale

abstract: Diffuse Galactic light has been observed in the optical since the 1930s. We propose that, when observed in the optical with deep imaging surveys, it can be used as a tracer of the turbulent cascade in the diffuse interstellar medium (ISM), down to scales of about 1 arcsec. Here we present a power spectrum analysis of the dust column density of a diffuse cirrus at high Galactic latitude (l ~ 198 deg, b ~ 32 deg) as derived from the combination of a MegaCam g-band image, obtained as part of the MATLAS large programme at the CFHT, with Planck radiance and Wise 12 micron data. The combination of these three datasets have allowed us to compute the density power spectrum of the HI over scales of more than three orders of magnitude. We found that the density field is well described by a single power law over scales ranging from 0.01 to 50 pc. The exponent of the power spectrum, gamma=-2.9 +- 0.1, is compatible with what is expected for thermally bi-stable and turbulent HI. We did not find any steepening of the power spectrum at small scales indicating that the typical scale at which turbulent energy is dissipated in this medium is smaller than 0.01pc. The ambipolar diffusion scenario that is usually proposed as the main dissipative agent, is consistent with our data only if the density of the cloud observed is higher than the typical values assumed for the cold neutral medium gas. We discuss the new avenue offered by deep optical imaging surveys for the study of the low density ISM structure and turbulence.


title: No bad local minima: Data independent training error guarantees for multilayer neural networks

abstract: We use smoothed analysis techniques to provide guarantees on the training loss of Multilayer Neural Networks (MNNs) at differentiable local minima. Specifically, we examine MNNs with piecewise linear activation functions, quadratic loss and a single output, under mild over-parametrization. We prove that for a MNN with one hidden layer, the training error is zero at every differentiable local minimum, for almost every dataset and dropout-like noise realization. We then extend these results to the case of more than one hidden layer. Our theoretical guarantees assume essentially nothing on the training data, and are verified numerically. These results suggest why the highly non-convex loss of such MNNs can be easily optimized using local updates (e.g., stochastic gradient descent), as observed empirically.


title: Naturally light neutrinos in $Diracon$ model

abstract: We propose a simple model for Dirac neutrinos where the smallness of neutrino mass follows from a parameter $\kappa$ whose absence enhances the symmetry of the theory. The symmetry breaking is performed by a two-doublet Higgs sector supplemented by an extra gauge singlet scalar, realizing an accidental global $U(1)$ symmetry. Its spontaneous breaking at the few TeV scale leads to a physical Nambu-Goldstone boson - the $Diracon$, denoted $\mathcal{D}$ - which is restricted by astrophysics and induces invisible Higgs decays $h\to \mathcal{D} \mathcal {D}$. The scheme provides a rich, yet very simple reference scenario for symmetry breaking studies at colliders such as the LHC.


title: Asymmetric Quantum Dialogue in Noisy Environment

abstract: A notion of asymmetric quantum dialogue (AQD) is introduced. Conventional protocols of quantum dialogue are essentially symmetric as both the users (Alice and Bob) can encode the same amount of classical information. In contrast, the scheme for AQD introduced here provides different amount of communication powers to Alice and Bob. The proposed scheme, offers an architecture, where the entangled state and the encoding scheme to be shared between Alice and Bob depends on the amount of classical information they want to exchange with each other. The general structure for the AQD scheme has been obtained using a group theoretic structure of the operators introduced in (Shukla et al., Phys. Lett. A, 377 (2013) 518). The effect of different types of noises (e.g., amplitude damping and phase damping noise) on the proposed scheme is investigated, and it is shown that the proposed AQD is robust and uses optimized amount of quantum resources.


title: Duration problem: basic concept and some extensions

abstract: We consider a sequence of independent random variables with the known distribution observed sequentially. The observation $n$ is assumed to be a value of one order statistics such as s:n-th, where 1 is less than s is less than n. It the instances following the $n$th observation it may remain of the s:m or it will be the value of the order statistics r:m (of m> n observations). Changing the rank of the observation, along with expanding a set of observations there is a random phenomenon that is difficult to predict. From practical reasons it is of great interest. Among others, we pose the question of the moment in which the observation appears and whose rank will not change significantly until the end of sampling of a certain size. We also attempt to answer which observation should be kept to have the "good quality observation" as long as possible. This last question was analysed by Ferguson, Hardwick and Tamaki (1991) in the abstract form which they called the problem of duration. This article gives a systematical presentation of the known duration models and a new modifications. We collect results from different papers on the duration of the extremal observation in the no-information (denoted as rank based) case and the full-information case. In the case of non-extremal observation duration models the most appealing are various settings related to the two extremal order statistics. In the no-information case it will be the maximizing duration of owning the relatively best or the second best object. The idea was formulated and the problem was solved by Szajowski and Tamaki (2006). The full-information duration problem with special requirement was presented by Kurushima and Ano (2010).


title: Electrostatic Control over Temperature-Dependent Tunneling across a Single Molecule Junction

abstract: Understanding how the mechanism of charge transport through molecular tunnel junctions depends on temperature is crucial to control electronic function in molecular electronic devices. With just a few systems investigated as a function of bias and temperature so far, thermal effects in molecular tunnel junctions remain poorly understood. Here we report a detailed charge transport study of an individual redox-active ferrocene-based molecule over a wide range of temperatures and applied potentials. The results show the temperature dependence of the current to vary strongly as a function of the gate voltage. Specifically, the current across the molecule exponentially increases in the Coulomb blockade regime and decreases at the charge degeneracy points, while remaining temperature-independent at resonance. Our observations can be well accounted for by a formal single-level tunneling model where the temperature dependence relies on the thermal broadening of the Fermi distributions of the electrons in the leads.


title: Hitting minors, subdivisions, and immersions in tournaments

abstract: The Erd\H{o}s-P\'osa property relates parameters of covering and packing of combinatorial structures and has been mostly studied in the setting of undirected graphs. In this note, we use results of Chudnovsky, Fradkin, Kim, and Seymour to show that, for every directed graph $H$ (resp. strongly-connected directed graph $H$), the class of directed graphs that contain $H$ as a strong minor (resp. butterfly minor, topological minor) has the vertex-Erd\H{o}s-P\'osa property in the class of tournaments. We also prove that if $H$ is a strongly-connected directed graph, the class of directed graphs containing $H$ as an immersion has the edge-Erd\H{o}s-P\'osa property in the class of tournaments.


title: Probabilistic Inference Modulo Theories

abstract: We present SGDPLL(T), an algorithm that solves (among many other problems) probabilistic inference modulo theories, that is, inference problems over probabilistic models defined via a logic theory provided as a parameter (currently, propositional, equalities on discrete sorts, and inequalities, more specifically difference arithmetic, on bounded integers). While many solutions to probabilistic inference over logic representations have been proposed, SGDPLL(T) is simultaneously (1) lifted, (2) exact and (3) modulo theories, that is, parameterized by a background logic theory. This offers a foundation for extending it to rich logic languages such as data structures and relational data. By lifted, we mean algorithms with constant complexity in the domain size (the number of values that variables can take). We also detail a solver for summations with difference arithmetic and show experimental results from a scenario in which SGDPLL(T) is much faster than a state-of-the-art probabilistic solver.


title: Inferring biological networks by sparse identification of nonlinear dynamics

abstract: Inferring the structure and dynamics of network models is critical to understanding the functionality and control of complex systems, such as metabolic and regulatory biological networks. The increasing quality and quantity of experimental data enable statistical approaches based on information theory for model selection and goodness-of-fit metrics. We propose an alternative method to infer networked nonlinear dynamical systems by using sparsity-promoting $\ell_1$ optimization to select a subset of nonlinear interactions representing dynamics on a fully connected network. Our method generalizes the sparse identification of nonlinear dynamics (SINDy) algorithm to dynamical systems with rational function nonlinearities, such as biological networks. We show that dynamical systems with rational nonlinearities may be cast in an implicit form, where the equations may be identified in the null-space of a library of mixed nonlinearities including the state and derivative terms; this approach applies more generally to implicit dynamical systems beyond those containing rational nonlinearities. This method, implicit-SINDy, succeeds in inferring three canonical biological models: Michaelis-Menten enzyme kinetics, the regulatory network for competence in bacteria, and the metabolic network for yeast glycolysis.


title: Semiparametric Estimation of Dynamic Discrete-Choice Models

abstract: We consider the estimation of dynamic discrete choice models in a semiparametric setting, in which the per-period utility functions are known up to a finite number of parameters, but the distribution of utility shocks is left unspecified. This semiparametric setup differs from most of the existing identification and estimation literature for dynamic discrete choice mod- els. To show identification we derive and exploit a new Bellman-like recursive representation for the unknown quantile function of the utility shocks. Our estimators are straightforward to compute; some are simple and require no iteration, and resemble classic estimators from the literature on semiparametric regression and average derivative estimation. Monte Carlo simulations demonstrate that our estimator performs well in small samples. To highlight features of this estimator, we estimate a structural model of dynamic labor supply for New York City taxicab drivers.


title: Provable Efficient Online Matrix Completion via Non-convex Stochastic Gradient Descent

abstract: Matrix completion, where we wish to recover a low rank matrix by observing a few entries from it, is a widely studied problem in both theory and practice with wide applications. Most of the provable algorithms so far on this problem have been restricted to the offline setting where they provide an estimate of the unknown matrix using all observations simultaneously. However, in many applications, the online version, where we observe one entry at a time and dynamically update our estimate, is more appealing. While existing algorithms are efficient for the offline setting, they could be highly inefficient for the online setting. In this paper, we propose the first provable, efficient online algorithm for matrix completion. Our algorithm starts from an initial estimate of the matrix and then performs non-convex stochastic gradient descent (SGD). After every observation, it performs a fast update involving only one row of two tall matrices, giving near linear total runtime. Our algorithm can be naturally used in the offline setting as well, where it gives competitive sample complexity and runtime to state of the art algorithms. Our proofs introduce a general framework to show that SGD updates tend to stay away from saddle surfaces and could be of broader interests for other non-convex problems to prove tight rates.


title: Multi-threaded Geant4 on the Xeon-Phi with Complex High-Energy Physics Geometry

abstract: To study the performance of multi-threaded Geant4 for high-energy physics experiments, an application has been developed which generalizes and extends previous work. A highly-complex detector geometry is used for benchmarking on an Intel Xeon Phi coprocessor. In addition, an implementation of parallel I/O based on Intel SCIF and ROOT technologies is incorporated and studied.


title: On constructing weight structures and extending them to idempotent extensions

abstract: We describe a new method for constructing a weight structure $w$ on a triangulated category $C$. For a given $C$ and $w$ it allow us to give a fairly comprehensive (and new) description of those triangulated categories consisting of retracts of objects of $C$ (i.e., of subcategories of the Karoubi envelope of $C$ that contain $C$; we call them idempotent extensions of $C$) such that $w$ extends to them. In particular, any bounded above or below $w$ extends to any idempotent extension of $C$. We also discuss the applications of our results to certain triangulated categories of ("relative") motives.


title: Measuring cellular traction forces on non-planar substrates

abstract: Animal cells use traction forces to sense the mechanics and geometry of their environment. Measuring these traction forces requires a workflow combining cell experiments, image processing and force reconstruction based on elasticity theory. Such procedures have been established before mainly for planar substrates, in which case one can use the Green's function formalism. Here we introduce a worksflow to measure traction forces of cardiac myofibroblasts on non-planar elastic substrates. Soft elastic substrates with a wave-like topology were micromolded from polydimethylsiloxane (PDMS) and fluorescent marker beads were distributed homogeneously in the substrate. Using feature vector based tracking of these marker beads, we first constructed a hexahedral mesh for the substrate. We then solved the direct elastic boundary volume problem on this mesh using the finite element method (FEM). Using data simulations, we show that the traction forces can be reconstructed from the substrate deformations by solving the corresponding inverse problem with a L1-norm for the residue and a L2-norm for 0th order Tikhonov regularization. Applying this procedure to the experimental data, we find that cardiac myofibroblast cells tend to align both their shapes and their forces with the long axis of the deformable wavy substrate.


title: Kronecker Determinantal Point Processes

abstract: Determinantal Point Processes (DPPs) are probabilistic models over all subsets a ground set of $N$ items. They have recently gained prominence in several applications that rely on "diverse" subsets. However, their applicability to large problems is still limited due to the $\mathcal O(N^3)$ complexity of core tasks such as sampling and learning. We enable efficient sampling and learning for DPPs by introducing KronDPP, a DPP model whose kernel matrix decomposes as a tensor product of multiple smaller kernel matrices. This decomposition immediately enables fast exact sampling. But contrary to what one may expect, leveraging the Kronecker product structure for speeding up DPP learning turns out to be more difficult. We overcome this challenge, and derive batch and stochastic optimization algorithms for efficiently learning the parameters of a KronDPP.


title: Generalization Properties and Implicit Regularization for Multiple Passes SGM

abstract: We study the generalization properties of stochastic gradient methods for learning with convex loss functions and linearly parameterized functions. We show that, in the absence of penalizations or constraints, the stability and approximation properties of the algorithm can be controlled by tuning either the step-size or the number of passes over the data. In this view, these parameters can be seen to control a form of implicit regularization. Numerical results complement the theoretical findings.


title: Warburg's impedance revisited

abstract: The derivation of Warburg's impedance presented in several books and scientific papers is reconsidered. It has been obtained by assuming that the total electric current across the sample is just due to the diffusion, and that the external potential applied to the electrode is responsible for an increase of the bulk density of charge described by Nernst's model. We show that these assumptions are not correct, and hence the proposed derivations questionable. A correct determination of the electrochemical impedance of a cell of an insulating material where are injected external charges of a given sign, when the diffusion and the displacement currents are taken into account, does not predict, in the high frequency region, for the real and imaginary parts of the impedance, the trends predicted by Warburg's impedance in the Nernstian approximation. The presented model can be generalized to the case of asymmetric cell, assuming boundary conditions physically sound.


title: Secondary radiation measurements for particle therapy applications: prompt photons produced by $^{4}$He, $^{12}$C and $^{16}$O ion beams in a PMMA target

abstract: Charged particle beams are used in Particle Therapy (PT) to treat oncological patients due to their selective dose deposition in tissues and to their high biological effect in killing cancer cells with respect to photons and electrons used in conventional radiotherapy. Nowadays, protons and carbon ions are used in PT clinical routine but, recently, the interest on the potential application of helium and oxygen beams is growing due to their reduced multiple scattering inside the body and increased linear energy transfer, relative biological effectiveness and oxygen enhancement ratio. The precision of PT demands for online dose monitoring techniques, crucial to improve the quality assurance of treatments. The beam range confined in the irradiated target can be monitored thanks to the neutral or charged secondary radiation emitted by the interactions of hadron beams with matter. Prompt photons are produced by nuclear de-excitation processes and, at present, different dose monitoring and beam range verification techniques based on the prompt {\gamma} detection have been proposed. It is hence of importance to perform the {\gamma} yield measurement in therapeutical-like conditions. In this paper we report the yields of prompt photons produced by the interaction of helium, carbon and oxygen ion beams with a PMMA target. The measurements were performed at the Heidelberg Ion-beam Therapy center (HIT) with beams of different energies. A LYSO scintillator has been used as photon detector. The obtained {\gamma} yields for $^{12}$C ion beams are compared with results from literature, while no other results from $^{4}$He and $^{16}$O beams have been published yet. A discussion on the expected resolution of a slit camera detector is presented, demonstrating the feasibility of a prompt-{\gamma} based monitoring technique for PT treatments using helium, carbon and oxygen ion beams.


title: Newtonian Cosmology with a Quantum Bounce

abstract: It has been known for some time that the cosmological Friedmann equation deduced from General Relativity can be also obtained within the Newtonian framework under certain assumptions. We use this result together with quantum corrections to the Newtonian potentials to derive a set a of quantum corrected Friedmann equations. We examine the behavior of the solutions of these modified cosmological equations paying special attention to the sign of the quantum corrections. We find different quantum effects crucially depending on this sign. One such a solution displays a qualitative resemblance to other quantum models like Loop Quantum Gravity or non-commutative geometry.


title: Determination of the Total Solar Modulation Factors in the Heliosphere For Cosmic Ray Protons and Electrons by Comparing Interstellar Spectra Deduced from Voyager Measurements and PAMELA Spectra of These Particles at the Earth

abstract: We have determined the interstellar spectra of cosmic ray protons and electrons from a few MeV to ~10 GeV. These interstellar spectra are based on Voyager data and a normalization of specific galactic propagation model calculations of both protons and electrons to PAMELA data at the Earth at 10 GeV, where the solar modulation is small. These resulting interstellar spectra are then compared with spectra of protons and electrons measured at lower energies at the Earth by PAMELA in 2009. The total amount of modulation at lower rigidities (energies) is found to be nearly the same at the same rigidity for both protons and electrons and ranges in magnitude from a factor ~400 at 0.1 GV for electrons, to a factor ~15 at 0.44 GV (100 MeV for protons), to a factor ~3.3 at 1 GV for both components. The magnitude of this total modulation of both components are the same to within + 10% from ~0.3 to ~3 GV in rigidity. The observed total modulation for protons can be matched quite closely using a simple spherically symmetric modulation picture involving a force field model for the modulation and a constant energy loss at all energies. The electrons require a set of parameters to calculate more detailed features of the modulation using the diffusion coefficient and its variation with rigidity and radius at rigidities less than few GV.


title: Active Control of Probability Amplitudes in a Mesoscale System via Feedback-Induced Suppression of Dissipation and Noise

abstract: We introduce a classical potentiostatic feedback mechanism that attenuates the dissipation in a quantum system arising from coupling to the surrounding thermodynamic bath, preserving the inter-state interference in an electronic excitation transfer (EET) process. A three-terminal potentiostat device applies a low-noise voltage bias to the terminals of the EET system and reduces the physical coupling between the quantum system and its environment. We introduce a classical equivalent circuit to model the environment-coupled excitation transfer in an elementary two-state system. This model provides qualitative insight into how classical feedback action affects the transition probabilities between the states and selectively reduces the dissipative coupling for one of the vibronic energy levels of the transfer system. Furthermore, we show that negative feedback results in persistent spectral coherence between the energy level of the decoupled state and the vibronic levels of the complementary state, making the decoupled vibronic channel a probe for characterizing the vibronic structure of the complementary channel of the EET system.


title: On Modeling Coverage and Rate of Random Cellular Networks under Generic Channel Fading

abstract: In this paper we provide an analytic framework for computing the expected downlink coverage probability, and the associated data rate of cellular networks, where base stations are distributed in a random manner. The provided expressions are in computable integral forms that accommodate generic channel fading conditions. We develop these expressions by modelling the cellular interference using stochastic geometry analysis, then we employ them for comparing the coverage resulting from various channel fading conditions namely Rayleigh and Rician fading, in addition to the fading-less channel. Furthermore, we expand the work to accommodate the effects of random frequency reuse on the cellular coverage and rate. Monte-Carlo simulations are conducted to validate the theoretical analysis, where the results show a very close match.


title: Rank parity for congruent supersingular elliptic curves

abstract: A recent paper of Shekhar compares the ranks of elliptic curves $E_1$ and $E_2$ for which there is an isomorphism $E_1[p] \simeq E_2[p]$ as $\mathrm{Gal}(\bar{\mathbf{Q}}/\mathbf{Q})$-modules, where $p$ is a prime of good ordinary reduction for both curves. In this paper we prove an analogous result in the case where $p$ is a prime of good supersingular reduction.


title: The number of cycles in random permutations without long cycles is asymptotically Gaussian

abstract: For uniform random permutations conditioned to have no long cycles, we prove that the total number of cycles satisfies a central limit theorem. Under additional assumptions on the asymptotic behavior of the set of allowed cycle lengths, we derive asymptotic expansions for the corresponding expected value and variance.


title: Quantum Router with Network Coding

abstract: Many protocols of quantum information processing, like quantum key distribution or measurement-based quantum computation, "consume" entangled quantum states during their execution. When participants are located at distant sites, these resource states need to be distributed. Due to transmission losses quantum repeater become necessary for large distances (e.g. $\gtrsim$ 300 km). Here we generalize the concept of the graph state repeater to $D$-dimensional graph states and to repeaters that can perform basic measurement-based quantum computations, which we call quantum routers. This processing of data at intermediate network nodes is called quantum network coding. We describe how a scheme to distribute general two-colorable graph states via quantum routers with network coding can be constructed from classical linear network codes. The robustness of the distribution of graph states against outages of network nodes is analysed by establishing a link to stabilizer error correction codes. Furthermore we show, that for any stabilizer error correction code there exists a corresponding quantum network code with similar error correcting capabilities.


title: Comparative study of gyrokinetic, hybrid-kinetic and fully kinetic wave physics for space plasmas

abstract: A set of numerical solvers for the linear dispersion relations of the gyrokinetic, the hybrid-kinetic, and the fully kinetic model is employed to study the physics of the kinetic Alfv\'en wave and the fast magnetosonic mode in these models. In particular, we focus on parameters that are relevant for solar wind oriented applications (using a homogeneous, isotropic background), which are characterized by wave propagation angles averaging close to 90{\deg}. It is found that the gyrokinetic model, while lacking high-frequency solutions and cyclotron effects, faithfully reproduces the fully kinetic Alfv\'en wave physics close to, and sometimes significantly beyond, the boundaries of its range of validity. The hybrid-kinetic model, on the other hand, is much more complete in terms of high-frequency waves, but owing to its simple electron model it is found to severely underpredict wave damping rates even on ion spatial scales across a large range of parameters, despite containing full kinetic ion physics.


title: Heat-bath random walks with Markov bases

abstract: Graphs on lattice points are studied whose edges come from a finite set of allowed moves of arbitrary length. We show that the diameter of these graphs on fibers of a fixed integer matrix can be bounded from above by a constant. We then study the mixing behaviour of heat-bath random walks on these graphs. We also state explicit conditions on the set of moves so that the heat-bath random walk, a generalization of the Glauber dynamics, is an expander in fixed dimension.


title: Notes on frequencies and timescales in nonequilibrium Green's functions

abstract: We discuss the ringdown behavior of the nonequilibrium Green's function in a strongly coupled theory with the holographic dual with a focus on quasinormal-mode equilibration. We study the time resolved spectral function for a probe scalar in Vaidya-AdS spacetime in detail as a complement to the preceding work arXiv:1603.06935 using further numerical results in very nonadiabatic temperature changes. It is shown that the relaxation of the nonequilibrium spectral function obtained through the Wigner transform is governed by the lowest quasinormal mode frequency. The timescale of the background temperature change is also observed in the frequency analysis. We then consider a toy model motivated by the quasinormal mode behavior and discuss these main features in numerical results are simply realized.


title: Computational Design Of Surfaces, Nanostructures and Optoelectronic Materials

abstract: Properties of engineering materials are generally influenced by defects such as point defects (vacancies, interstitials, substitutional defects), line defects (dislocations), planar defects (grain boundaries, free surfaces/nanostructures, interfaces, stacking faults) and volume defects (voids). Classical physics based molecular dynamics and quantum physics based density functional theory can be useful in designing materials with controlled defect properties. In this thesis, empirical potential based molecular dynamics was used to study the surface modification of polymers due to energetic polyatomic ion, thermodynamics and mechanics of metal-ceramic interfaces and nanostructures, while density functional theory was used to screen substituents in optoelectronic materials.


title: Impact of the neutron matter equation of state on neutron skin and neutron drip lines in chiral effective field theory

abstract: We present predictions of the binding energy per nucleon and the neutron skin thickness in highly neutron-rich isotopes of Oxygen, Magnesium, and Aluminum. The calculations are carried out at and below the neutron drip line as predicted by our model. The nuclear properties are obtained via an energy functional whose input is the equation of state of isospin-asymmetric infinite matter. The latter is based on a microscopic derivation of the energy per particle in neutron matter applying chiral few-nucleon forces together with a phenomenological model for the equation of state of symmetric nuclear matter. We highlight the impact of the neutron matter equation of state at different orders of chiral effective field theory on neutron skins and the binding energy per particle and quantify the uncertainty carried by our predictions.


title: Estimation of Passenger Route Choice Pattern Using Smart Card Data for Complex Metro Systems

abstract: Nowadays, metro systems play an important role in meeting the urban transportation demand in large cities. The understanding of passenger route choice is critical for public transit management. The wide deployment of Automated Fare Collection(AFC) systems opens up a new opportunity. However, only each trip's tap-in and tap-out timestamp and stations can be directly obtained from AFC system records; the train and route chosen by a passenger are unknown, which are necessary to solve our problem. While existing methods work well in some specific situations, they don't work for complicated situations. In this paper, we propose a solution that needs no additional equipment or human involvement than the AFC systems. We develop a probabilistic model that can estimate from empirical analysis how the passenger flows are dispatched to different routes and trains. We validate our approach using a large scale data set collected from the Shenzhen metro system. The measured results provide us with useful inputs when building the passenger path choice model.


title: Absolute semi-deviation risk measure for ordering problem with transportation cost in Supply Chain

abstract: We present a decomposition method for stochastic programs with 0-1 variables in the second-stage with absolute semi-deviation (ASD) risk measure. Traditional stochastic programming models are risk-neutral where expected costs are considered for the second-stage. A common approach to address risk is to include a dispersion statistic in addition with expected costs and weighted appropriately. Due to the lack of block angular structure, stochastic programs with ASD risk-measure possess computational challenges. The proposed decomposition algorithm uses another risk-measure `expected excess', and provides tighter bounds for ASD stochastic models. We perform computational study on a supply chain replenishment problem and standard knapsack instances. The computational results using supply chain instances demonstrate the usefulness of ASD risk-measure in decision making under uncertainty, and knapsack instances indicate that the proposed methodology outperforms a direct solver.