title: Catalytic quantum randomness

abstract: Randomness is a defining element of mixing processes in nature and an essential ingredient to many protocols in quantum information. In this work, we investigate how much randomness is required to transform a given quantum state into another one. Specifically, we ask whether there is a gap between the power of a classical source of randomness compared to that of a quantum one. We provide a complete answer to these questions, by identifying provably optimal protocols for both classical and quantum sources of randomness, based on a dephasing construction. We find that in order to implement any noisy transition on a $d$-dimensional quantum system it is necessary and sufficient to have a quantum source of randomness of dimension $\sqrt{d}$ or a classical one of dimension $d$. Interestingly, coherences provided by quantum states in a source of randomness offer a quadratic advantage. The process we construct has the additional features to be robust and catalytic, i.e., the source of randomness can be re-used. Building upon this formal framework, we illustrate that this dephasing construction can serve as a useful primitive in both equilibration and quantum information theory: We discuss applications describing the smallest measurement device, capturing the smallest equilibrating environment allowed by quantum mechanics, or forming the basis for a cryptographic private quantum channel. We complement the exact analysis with a discussion of approximate protocols based on quantum expanders deriving from discrete Weyl systems. This gives rise to equilibrating environments of remarkably small dimension. Our results highlight the curious feature of randomness that residual correlations and dimension can be traded against each other.


title: Reduced order multi switching hybrid synchronization of chaotic systems

abstract: In this article, a new synchronization scheme is presented by combining the concept of reduced-order synchronization with multi-switching synchronization schemes. The presented scheme, reduced-order multi-switching hybrid synchronization, is notable addition to the earlier multi-switching schemes providing enhanced security in applications of secure communication. Based on the Lyapunov stability theory, the active control method is used to design the controllers and derive sufficient condition for achieving reduced-order multi-switching hybrid synchronization between a new hyperchaotic system taken as drive system and Qi chaotic system serving as response system. Numerical simulations are performed in MATLAB using the Runge-Kutta method to verify the effectiveness of the proposed method. The results show the utility and suitability of the active control method for achieving the reduced-order multi-switching hybrid synchronization among dynamical chaotic systems.


title: Estimation in a simple linear regression model with measurement error

abstract: This paper deals with the problem of estimating a slope parameter in a simple linear regression model, where independent variables have functional measurement errors. Measurement errors in independent variables, as is well known, cause biasedness of the ordinary least squares estimator. A general procedure for the bias reduction is presented in a finite sample situation, and some exact bias-reduced estimators are proposed. Also, it is shown that certain truncation procedures improve the mean square errors of the ordinary least squares and the bias-reduced estimators.


title: SSTRED: Data- and metadata-processing pipeline for CHROMIS and CRISP

abstract: Context: Data from ground-based, high-resolution solar telescopes can only be used for science with calibrations and processing, which requires detailed knowledge about the instrumentation. [...] Aims: We aim to provide observers with a user-friendly data pipeline for data from the Swedish 1-meter Solar Telescope (SST) that delivers science-ready data together with the metadata needed for proper interpretation and archiving. Methods: We briefly describe the [instrumentation]. We summarize the processing steps from raw data to science-ready data cubes in FITS files. We report calibrations and compensations for data imperfections in detail. Misalignment of \ion{Ca}{ii} data due to wavelength-dependent dispersion is identified, characterized, and compensated for. We describe intensity calibrations that remove or reduce the effects of filter transmission profiles as well as solar elevation changes. We present REDUX, a new version of the MOMFBD image restoration code, with multiple enhancements and new features. [...] We describe how image restoration is used [...]. The science-ready output is delivered in FITS files, with metadata compliant with the SOLARNET recommendations. Data cube coordinates are specified within the World Coordinate System (WCS). Cavity errors are specified as distortions of the WCS wavelength coordinate with an extension of existing WCS notation. We establish notation for specifying the reference system for Stokes vectors [...]. [CRISPEX] has been extended to accept SSTRED output [...]. Results: SSTRED is a mature data-processing pipeline for imaging instruments, developed and used for the SST/CHROMIS imaging spectrometer and the SST/CRISP spectropolarimeter. SSTRED delivers well-characterized, science-ready, archival-quality FITS files with well-defined metadata. The SSTRED code, as well as REDUX and CRISPEX, is freely available through git repositories.


title: Reconstructing the impact parameter of proton-nucleus and nucleus-nucleus collisions

abstract: In proton-nucleus and nucleus-nucleus collision experiments, one determines the centrality of a collision according to the multiplicity or energy deposited in a detector. This serves as a proxy for the true collision centrality, as defined by the impact parameter. We show that the probability distribution of impact parameter in a given bin of experiment-defined centrality can be reconstructed without assuming any specific model for the collision dynamics, in both proton-nucleus and nucleus-nucleus systems. The reconstruction is reliable up to about 10\% centrality, and is more accurate for nucleus-nucleus collisions. We perform an application of our procedure to experimental data from all the CERN Large Hadron Collider (LHC) collaborations, from which we extract, in Pb+Pb and $p$+Pb collisions, the corresponding distributions of impact parameter.


title: Approximate k-NN Graph Construction: a Generic Online Approach

abstract: Nearest neighbor search and k-nearest neighbor graph construction are two fundamental issues arise from many disciplines such as multimedia information retrieval, data-mining and machine learning. They become more and more imminent given the big data emerge in various fields in recent years. In this paper, a simple but effective solution both for approximate k-nearest neighbor search and approximate k-nearest neighbor graph construction is presented. These two issues are addressed jointly in our solution. On the one hand, the approximate k-nearest neighbor graph construction is treated as a search task. Each sample along with its k-nearest neighbors are joined into the k-nearest neighbor graph by performing the nearest neighbor search sequentially on the graph under construction. On the other hand, the built k-nearest neighbor graph is used to support k-nearest neighbor search. Since the graph is built online, the dynamic update on the graph, which is not possible from most of the existing solutions, is supported. This solution is feasible for various distance measures. Its effectiveness both as k-nearest neighbor construction and k-nearest neighbor search approaches is verified across different types of data in different scales, various dimensions and under different metrics.


title: A reduced finite element formulation for space fractional partial differential equation

abstract: Applying proper orthogonal decomposition to a usual finite element (FE) formulation for space fractional partial differential equation, we get a reduced FE model, which greatly reduces the complexity of computation. Then, the stability analysis and error estimate for the reduced model are presented. Finally, we verify the effectiveness of the algorithm by numerical experiments.


title: Fenchel-Moreau Conjugation Inequalities with Three Couplings and Application to Stochastic Bellman Equation

abstract: Given two couplings between "primal" and "dual" sets, we prove a general implication that relates an inequality involving "primal" sets to a reverse inequality involving the "dual" sets.% More precisely, let be given two "primal" sets $\PRIMAL$, $\PRIMALBIS$and two "dual" sets $\DUAL$, $\DUALBIS$, together with two {coupling} functions \(\PRIMAL \overset{\coupling}{\leftrightarrow} \DUAL \) and \(\PRIMALBIS \overset{\couplingbis}{\leftrightarrow} \DUALBIS \). We define a new coupling \(\SumCoupling{\coupling}{\couplingbis} \) between the "primal" product set~$\PRIMAL \times \PRIMALBIS$ and the "dual" product set $\DUAL \times \DUALBIS$. Then, we consider any bivariate function \(\kernel : \PRIMAL \times \PRIMALBIS \to \barRR \) and univariate functions \(\fonctionprimal : \PRIMAL \to \barRR \) and \(\fonctionprimalbis : \PRIMALBIS \to \barRR \), all defined on the "primal" sets. We prove that \(\fonctionprimal\np{\primal} \geq \inf\_{\primalbis \in \PRIMALBIS} \Bp{\kernel\np{\primal, \primalbis} \UppPlus \fonctionprimalbis\np{\primalbis}} \) \( \Rightarrow \SFM{\fonctionprimal}{\coupling}\np{\dual} \leq \inf\_{\dualbis \in \DUALBIS} \Bp{\SFM{\kernel}{\SumCoupling{\coupling}{\couplingbis}}\np{\dual,\dualbis} \UppPlus \SFM{\fonctionprimalbis}{-\couplingbis}\np{\dualbis}} \), where we stress that the Fenchel-Moreau conjugates \(\SFM{\fonctionprimal}{\coupling} \) and \(\SFM{\fonctionprimalbis}{-\couplingbis}\) are not necessarily taken with the same coupling. We study the equality case, after having established the classical Fenchel inequality but with a general coupling. % We display several applications. We provide a new formula for the Fenchel-Moreau conjugate of a generalized inf-convolution. We obtain formulas with partial Fenchel-Moreau conjugates. Finally, we consider the Bellman equation in stochastic dynamic programming and we provide a "Bellman-like" equation for the Fenchel conjugates of the value functions.


title: Universal phase diagram and scaling functions of imbalanced Fermi gases

abstract: We discuss the phase diagram and the universal scaling functions of attractive Fermi gases at finite imbalance. The existence of a quantum multicritical point for the unitary gas at vanishing chemical potential $\mu$ and effective magnetic field $h$, first discussed by Nikoli\'{c} and Sachdev, gives rise to three different phase diagrams, depending on whether the inverse scattering length $1/a$ is negative, positive or zero. Within a Luttinger-Ward formalism, the phase diagram and pressure of the unitary gas is calculated as a function of the dimensionless scaling variables $T/\mu$ and $h/\mu$. The results indicate that beyond the Clogston-Chandrasekhar limit at $(h/\mu)_c\simeq 1.09$, the unitary gas exhibits an inhomogeneous superfluid phase with FFLO order that can reach critical temperatures near unitarity of $\simeq 0.03\, T_F$.


title: Image Moment Models for Extended Object Tracking

abstract: In this paper, a novel image moments based model for shape estimation and tracking of an object moving with a complex trajectory is presented. The camera is assumed to be stationary looking at a moving object. Point features inside the object are sampled as measurements. An ellipsoidal approximation of the shape is assumed as a primitive shape. The shape of an ellipse is estimated using a combination of image moments. Dynamic model of image moments when the object moves under the constant velocity or coordinated turn motion model is derived as a function for the shape estimation of the object. An Unscented Kalman Filter-Interacting Multiple Model (UKF-IMM) filter algorithm is applied to estimate the shape of the object (approximated as an ellipse) and track its position and velocity. A likelihood function based on average log-likelihood is derived for the IMM filter. Simulation results of the proposed UKF-IMM algorithm with the image moments based models are presented that show the estimations of the shape of the object moving in complex trajectories. Comparison results, using intersection over union (IOU), and position and velocity root mean square errors (RMSE) as metrics, with a benchmark algorithm from literature are presented. Results on real image data captured from the quadcopter are also presented.


title: 3D Fluid Flow Estimation with Integrated Particle Reconstruction

abstract: The standard approach to densely reconstruct the motion in a volume of fluid is to inject high-contrast tracer particles and record their motion with multiple high-speed cameras. Almost all existing work processes the acquired multi-view video in two separate steps, utilizing either a pure Eulerian or pure Lagrangian approach. Eulerian methods perform a voxel-based reconstruction of particles per time step, followed by 3D motion estimation, with some form of dense matching between the precomputed voxel grids from different time steps. In this sequential procedure, the first step cannot use temporal consistency considerations to support the reconstruction, while the second step has no access to the original, high-resolution image data. Alternatively, Lagrangian methods reconstruct an explicit, sparse set of particles and track the individual particles over time. Physical constraints can only be incorporated in a post-processing step when interpolating the particle tracks to a dense motion field. We show, for the first time, how to jointly reconstruct both the individual tracer particles and a dense 3D fluid motion field from the image data, using an integrated energy minimization. Our hybrid Lagrangian/Eulerian model reconstructs individual particles, and at the same time recovers a dense 3D motion field in the entire domain. Making particles explicit greatly reduces the memory consumption and allows one to use the high-res input images for matching. Whereas the dense motion field makes it possible to include physical a-priori constraints and account for the incompressibility and viscosity of the fluid. The method exhibits greatly (~70%) improved results over our recently published baseline with two separate steps for 3D reconstruction and motion estimation. Our results with only two time steps are comparable to those of sota tracking-based methods that require much longer sequences.


title: Optimal Control of AGC Systems Considering Non-Gaussian Wind Power Uncertainty

abstract: Wind power uncertainty poses significant challenges for automatic generation control (AGC) systems. It can enhance control performances to explicitly consider wind power uncertainty distributions within controller design. However, widely accepted wind uncertainties usually follow non-Gaussian distributions which may lead to complicated stochastic AGC modeling and high computational burdens. To overcome the issue, this paper presents a novel Ito-theory-based model for the stochastic control problem (SCP) of AGC systems, which reduces the computational burden of optimization considering non-Gaussian wind power uncertainty to the same scale as that for deterministic control problems. We present an Ito process model to exactly describe non-Gaussian wind power uncertainty, and then propose an SCP based on the concept of stochastic assessment functions (SAFs). Based on a convergent series expansion of the SAF, the SCP is reformulated as a certain deterministic control problem without sacrificing performance under non-Gaussian wind power uncertainty. The reformulated control problem is proven as a convex optimization which can be solved efficiently. A case study demonstrates the efficiency and accuracy of the proposed approach compared with several conventional approaches.


title: An Infinite Family of Maximally Superintegrable Systems in a Magnetic Field with Higher Order Integrals

abstract: We construct an additional independent integral of motion for a class of three dimensional minimally superintegrable systems with constant magnetic field. This class was introduced in [J. Phys. A: Math. Theor. 50 (2017), 245202, 24 pages] and it is known to possess periodic closed orbits. In the present paper we demonstrate that it is maximally superintegrable. Depending on the values of the parameters of the system, the newly found integral can be of arbitrarily high polynomial order in momenta.


title: Zero-temperature limit and statistical quasiparticles in many-body perturbation theory

abstract: The order-by-order renormalization of the self-consistent mean-field potential in many-body perturbation theory for normal Fermi systems is investigated in detail. Building on previous work mainly by Balian and de Dominicis, as a key result we derive a thermodynamic perturbation series that manifests the consistency of the adiabatic zero-temperature formalism with perturbative statistical mechanics---for both isotropic and anisotropic systems---and satisfies at each order and for all temperatures the thermodynamic relations associated with Fermiliquid theory. These properties are proved to all orders.


title: How to Synthesize Exceptional Points with Three Resonators

abstract: In non-Hermitian coulped-resonator networks, the eigenvectors of degenerate eigenmodes may become parallel due to the singularity at so-called Exceptional Points (EP). To exploit the parametric sensitivity at EPs, an important problem is, given an arbitrary set of coupled resonators, how to generate a desired EP by properly coupling them together. This paper provides the solution for the case of three resonators. We show that all physically admissible EPs can be realized with either weakly coupled linear networks or strongly coupled circular networks, and the latter type of EPs has not been reported in the literature. Each admissible EP eigenvalue can be realized by two and only two resonator networks, and the formulas for calculating the required coupling constants are provided. The characteristics of these EPs are illustrated by the change of transmission spectra near them, which verify the enhanced sensitivity induced by the singularity of EPs.


title: Spatial search for a general multi-vertex state on graph by continuous-time quantum walks

abstract: In this work, we consider the spatial search for a general marked state on graphs by continuous time quantum walks. As a simplest case, we compute the amplitude expression of the search for the multi-vertex uniform superposition state on hypercube, and find that the spatial search algorithm is optimal for the two-vertex uniform state. However, on general graphs, a common formula can't be obtained for searching a general non-uniform superposition state. Fortunately, a Laplacian spectrum condition which determines whether the associated graph could be appropriate for performing the optimal spatial search is presented. The condition implies that if the proportion of the maximum and the non-zero minimum Laplacian eigenvalues is less or equal to 1+sqrt(1/2), then the spatial search is optimal for any general state. At last, we apply this condition to three kind graphs, the induced complete graph, the strongly regular graph and the regular complete multi-partite graphs. By the condition, one can conclude that these graphs will become suitable for optimal search with properly setting their graph parameters.


title: Multiresolution on n-dimensional spheres

abstract: In the present paper, multiscale systems of polynomial wavelets on an n-dimensional sphere are constructed. Scaling functions and wavelets are investigated,and their reproducing and localization properties and positive definiteness are examined. Decomposition and reconstruction algorithms for the wavelet transform are presented. Formulae for variances in space and momentum domain, as well as for the uncertainty product, of zonal functions over n-dimensional spheres are derived and applied to the scaling functions.


title: Angular multiselectivity with spherical wavelets

abstract: We construct spherical wavelets based on approximate identities that are directional, i.e. not rotation-invariant, and have an adaptive angular selectivity. The problem of how to find a proper representation of distinct kinds of details of real images, ranging from highly directional to fully isotropic ones, was quite intensively studied for the case of signals over the Euclidean space. However, the present paper is the first attempt to deal with this task in the case of spherical signals. A multiselectivity scheme, similar to that proposed for $\mathbb R^2$-functions, is presented.


title: Uncertainty of Poisson wavelets

abstract: Poisson wavelets are a powerful tool in analysis of spherical signals. In order to have a deeper characterization of them, we compute their uncertainty product, a quantity introduced for the first time by Narcowich and Ward in~\cite{NW96} and used to measure the trade-off between the space and frequency localization of a function. Surprisingly, the uncertainty product of Poisson wavelets tends to the minimal value in some limiting cases. This shows that in the case of spherical functions, not only Gauss kernel has this property.


title: Uncertainty product of spherical wavelets

abstract: In the paper, asymptotic behavior of the uncertainty product for a family of zonal spherical wavelets is computed. The family contains the most popular wavelets, such as Gauss-Weierstrass, Abel-Poisson and Poisson wavelets and Mexican needlets. Boundedness of the uncertainty constant is in general not given, but it is a property of some of the wavelets from this class.


title: Positive definite functions on semilattices

abstract: We introduce a notion of positive definiteness for functions $f\!:P\to\mathbb{R}$ defined on meet semilattices $(P,\preceq,\wedge)$ and prove several properties for these functions. In addition, we utilize the $LDL^{\rm T}$ decomposition of meet matrices in order to explore the properties of multivariate positive definite arithmetic functions $f\!:\mathbb{Z}_+^d\to\mathbb{R}$. Finally, we give a series of examples and counterexamples of positive definite functions.


title: Clustrophile 2: Guided Visual Clustering Analysis

abstract: Data clustering is a common unsupervised learning method frequently used in exploratory data analysis. However, identifying relevant structures in unlabeled, high-dimensional data is nontrivial, requiring iterative experimentation with clustering parameters as well as data features and instances. The number of possible clusterings for a typical dataset is vast, and navigating in this vast space is also challenging. The absence of ground-truth labels makes it impossible to define an optimal solution, thus requiring user judgment to establish what can be considered a satisfiable clustering result. Data scientists need adequate interactive tools to effectively explore and navigate the large clustering space so as to improve the effectiveness of exploratory clustering analysis. We introduce \textit{Clustrophile~2}, a new interactive tool for guided clustering analysis. \textit{Clustrophile~2} guides users in clustering-based exploratory analysis, adapts user feedback to improve user guidance, facilitates the interpretation of clusters, and helps quickly reason about differences between clusterings. To this end, \textit{Clustrophile~2} contributes a novel feature, the Clustering Tour, to help users choose clustering parameters and assess the quality of different clustering results in relation to current analysis goals and user expectations. We evaluate \textit{Clustrophile~2} through a user study with 12 data scientists, who used our tool to explore and interpret sub-cohorts in a dataset of Parkinson's disease patients. Results suggest that \textit{Clustrophile~2} improves the speed and effectiveness of exploratory clustering analysis for both experts and non-experts.


title: Growth in linear algebraic groups and permutation groups: towards a unified perspective

abstract: By now, we have a product theorem in every finite simple group $G$ of Lie type, with the strength of the bound depending only in the rank of $G$. Such theorems have numerous consequences: bounds on the diameters of Cayley graphs, spectral gaps, and so forth. For the alternating group Alt_n, we have a quasipolylogarithmic diameter bound (Helfgott-Seress 2014), but it does not rest on a product theorem. We shall revisit the proof of the bound for Alt_n, bringing it closer to the proof for linear algebraic groups, and making some common themes clearer. As a result, we will show how to prove a product theorem for Alt_n -- not of full strength, as that would be impossible, but strong enough to imply the diameter bound.


title: Boiling crisis dynamics: low gravity experiments at high pressure

abstract: To understand the boiling crisis mechanism, one can take advantage of the slowing down of boiling at high pressures, in the close vicinity of the liquid-vapor critical point of the given fluid. To preserve conventional bubble geometry, such experiments need to be carried out in low gravity. We report here two kinds of saturated boiling experiments. First we discuss the spatial experiments with SF$_6$ at 46$^\circ$ C. Next we address two ground-based experiments under magnetic gravity compensation with H$_2$ at 33 K. We compare both kinds of experiments and show their complementarity. The dry spots under vapor bubbles are visualized by using transparent heaters made with metal oxide films. We evidence two regimes of the dry spots growth: the regime of circular dry spots and the regime of chain coalescence of dry spots that immediately precedes the heater dryout. A recent H$_2$ experiment is shown to bridge the gap between the near-critical and low pressure boiling experiments.


title: Gromov Product Decomposition of 7-point Metric Spaces

abstract: Let $X$ be a finite metric space with elements $P_i$, $i=1,\dots,n$ and with distance functions $d_{ij}$. The Gromov product of the triangle with vertices $P_i$, $P_j$ and $P_k$ at the vertex $P_i$ is defined by $\Delta_{ijk}=\frac{1}{2}(d_{ij}+d_{ik}-d_{jk})$. A metric space is called $\Delta$-generic, if the set of Gromov products at each $P_i$ has a unique smallest element $\Delta_{ijk}$. For a $\Delta$-generic metric space, the map $P_i\to (P_jP_k)$, where $(P_jP_k)$ is the edge joining $P_j$ to $P_k$ is a well defined map called the "Gromov product structure" [Bilge, Celik and Kocak, "An equivalence class decomposition of finite metric spaces", Discrete Mathmetics, Vol 340, (2017) 1928-1932]. For n=5, the 3 $\Delta$-equivalence classes coincide with the classification of $5$-point metrics. For $n=6$, there are 26 $\Delta$-equivalence classes refined by 339 metric classes. In this paper, we present the first systematic treatment of 7-point spaces and we obtain the $\Delta$-equivalence decomposition of 7-point metric spaces that consist of 431 equivalence classes.


title: Vision as an Interlingua: Learning Multilingual Semantic Embeddings of Untranscribed Speech

abstract: In this paper, we explore the learning of neural network embeddings for natural images and speech waveforms describing the content of those images. These embeddings are learned directly from the waveforms without the use of linguistic transcriptions or conventional speech recognition technology. While prior work has investigated this setting in the monolingual case using English speech data, this work represents the first effort to apply these techniques to languages beyond English. Using spoken captions collected in English and Hindi, we show that the same model architecture can be successfully applied to both languages. Further, we demonstrate that training a multilingual model simultaneously on both languages offers improved performance over the monolingual models. Finally, we show that these models are capable of performing semantic cross-lingual speech-to-speech retrieval.


title: Radiation mediated shocks in gamma-ray bursts: Subshock photon production

abstract: Internal shocks provide a plausible heating mechanism in the jets of gamma-ray bursts (GRBs). Shocks occurring below the jet photosphere are mediated by radiation. It was previously found that radiation mediated shocks (RMSs) inside GRB jets are inefficient photon producers, and the photons that mediate the RMS must originate from an earlier stage of the explosion. We show that this conclusion is valid only for non-magnetized jets. RMSs that propagate in moderately magnetized plasma develop a collisionless subshock which locally heats the plasma to a relativistic temperature, and the hot electrons emit copious synchrotron photons inside the RMS. We find that this mechanism is generally effective for mildly relativistic shocks and may be the main source of photons observed in GRBs. We derive a simple analytical formula for the generated photon number per proton, $Z$, which gives $Z=10^5$-$10^6$, consistent with observations. The photons are initially injected with low energies, well below the observed GRB peak. Their number is controlled by two main factors: (1) the abundance of electron-positron pairs created in the shock, which is self-consistently calculated, and (2) the upper limit on the brightness temperature of soft radiation set by induced downscattering. The injected soft photons that survive induced downscattering gain energy in the RMS via bulk Comptonization and shape its nonthermal spectrum.


title: Set Similarity Search for Skewed Data

abstract: Set similarity join, as well as the corresponding indexing problem set similarity search, are fundamental primitives for managing noisy or uncertain data. For example, these primitives can be used in data cleaning to identify different representations of the same object. In many cases one can represent an object as a sparse 0-1 vector, or equivalently as the set of nonzero entries in such a vector. A set similarity join can then be used to identify those pairs that have an exceptionally large dot product (or intersection, when viewed as sets). We choose to focus on identifying vectors with large Pearson correlation, but results extend to other similarity measures. In particular, we consider the indexing problem of identifying correlated vectors in a set S of vectors sampled from {0,1}^d. Given a query vector y and a parameter alpha in (0,1), we need to search for an alpha-correlated vector x in a data structure representing the vectors of S. This kind of similarity search has been intensely studied in worst-case (non-random data) settings. Existing theoretically well-founded methods for set similarity search are often inferior to heuristics that take advantage of skew in the data distribution, i.e., widely differing frequencies of 1s across the d dimensions. The main contribution of this paper is to analyze the set similarity problem under a random data model that reflects the kind of skewed data distributions seen in practice, allowing theoretical results much stronger than what is possible in worst-case settings. Our indexing data structure is a recursive, data-dependent partitioning of vectors inspired by recent advances in set similarity search. Previous data-dependent methods do not seem to allow us to exploit skew in item frequencies, so we believe that our work sheds further light on the power of data dependence.


title: Geometry and the Imagination in Minneapolis

abstract: This document consists of the collection of handouts for a two-week summer workshop entitled 'Geometry and the Imagination', led by John Conway, Peter Doyle, Jane Gilman and Bill Thurston at the Geometry Center in Minneapolis, June 17-28, 1991. The workshop was based on a course `Geometry and the Imagination' which we had taught twice before at Princeton.


title: Automated Discovery of Internet Censorship by Web Crawling

abstract: Censorship of the Internet is widespread around the world. As access to the web becomes increasingly ubiquitous, filtering of this resource becomes more pervasive. Transparency about specific content that citizens are denied access to is atypical. To counter this, numerous techniques for maintaining URL filter lists have been proposed by various individuals and organisations that aim to empirical data on censorship for benefit of the public and wider censorship research community. We present a new approach for discovering filtered domains in different countries. This method is fully automated and requires no human interaction. The system uses web crawling techniques to traverse between filtered sites and implements a robust method for determining if a domain is filtered. We demonstrate the effectiveness of the approach by running experiments to search for filtered content in four different censorship regimes. Our results show that we perform better than the current state of the art and have built domain filter lists an order of magnitude larger than the most widely available public lists as of Jan 2018. Further, we build a dataset mapping the interlinking nature of blocked content between domains and exhibit the tightly networked nature of censored web resources.


title: Convex fair partitions into an arbitrary number of pieces

abstract: We prove that any convex body in the plane can be partitioned into $m$ convex parts of equal areas and perimeters for any integer $m\ge 2$; this result was previously known for prime powers $m=p^k$. We also discuss possible higher-dimensional generalizations and difficulties of extending our technique to equalizing more than one non-additive function.


title: C^+ distribution around S1 in rho Ophiuchi

abstract: We analyze a [C II] 158 micron map obtained with the L2 GREAT receiver on SOFIA of the emission/reflection nebula illuminated by the early B star S1 in the rho-OphA cloud core. This data set has been complemented with maps of CO(3-2), 13CO(3-2) and C18O(3-2), observed as a part of the JCMT Gould Belt Survey, with archival HCO^+(4-3) JCMT data, as well as with [O I] 63 and 145 micron imaging with Herschel/PACS. The [C II] emission is completely dominated by the strong PDR emission from the nebula surrounding S1 expanding into the dense Oph A molecular cloud west and south of S1. The [C II] emission is significantly blue shifted relative to the CO spectra and also relative to the systemic velocity, particularly in the northwestern part of the nebula. The [C II] lines are broader towards the center of the S1 nebula and narrower towards the PDR shell. The [C II] lines are strongly self-absorbed over an extended region in the S1 PDR. Based on the strength of the [13C II] F = 2-1 hyperfine component, [C II] is significantly optically thick over most of the nebula. CO and 13CO(3-2) spectra are strongly self-absorbed, while C18O(3-2) is single peaked and centered in the middle of the self-absorption. We have used a simple two-layer LTE model to characterize the background and foreground cloud contributing to the [C II] emission. From this analysis we estimate the extinction due to the foreground cloud to be ~9.9 mag, which is slightly less than the reddening estimated towards S1. Since some of the hot gas in the PDR is not traced by low J CO emission, this result appears quite plausible. Using a plane parallel PDR model with the observed [OI(145)]/[C II] brightness ratio and an estimated FUV intensity of 3100-5000 G0 suggests that the density of the [C II] emitting gas is ~3-4x10^3 cm^-3.


title: Metric cones, N-body collisions, and Marchal's lemma

abstract: Marchal's lemma is the basic tool for eliminating collisions when using the direct method of the calculus of variations to establish existence of "designer" solutions to the classical N-body problem. Our goal here is to understand why Marchal's lemma holds, by taking a metric geometry perspective and employing the Jacobi-Maupertuis [JM] metric reformulation of mechanics. Using analysis inspired by the conical metric nature of the standard Kepler problem at zero-energy, we are able to manufacture potentials, or "counterexamples", for which Marchal's lemma fails. These counterexamples overlap significantly with results obtained by Barutello et al. A novel feature in our proof for the counterexample is the use of piecewise constant potentials, and the resulting piecewise constant metrics.


title: Superconductivity from Meissner Effect and Zero Resistivity in a Phenyl Molecule

abstract: Recently, phenyl molecules have been reported to exhibit Meissner effect mainly from magnetization measurements. Realizing zero-resistivity state in these materials seems a challenge due to many practical difficulties but is required to characterize the existence of superconductivity. By choosing potassium-doped tris(2-methylphenyl)bismuthine as an example, we perform temperature-dependent magnetic susceptibility and resistivity measurements at different magnetic fields and pressures. The solid evidence for supporting superconductivity is achieved from the obtained Meissner effect and zero resistivity with the critical temperature ($T_c$) of 3.6 K at atmosphere pressure. Upon compression, we observe the gradual evolution of superconductivity from its initial phase with a parabolic behavior of $T_{c}$ to the second one with almost constant value of $T_{c}$ of 7 K. The 7 K phase seems a common feature for these newly discovered phenyl-based superconductors.


title: Time-Resolved Mechanical Spectroscopy of Soft Materials via Optimally Windowed Chirps

abstract: The ability to measure the bulk dynamic behavior of soft materials with combined time- and frequency-resolution is instrumental for improving our fundamental understanding of connections between the microstructural dynamics and the macroscopic mechanical response. Current state-of-the-art techniques are often limited by a compromise between resolution in the time and frequency domain, mainly due to the use of elementary input signals that have not been designed for fast time-evolving systems such as materials undergoing gelation, curing or self-healing. In this work, we develop an optimized and robust excitation signal for time-resolved mechanical spectroscopy through the introduction of joint frequency- and amplitude-modulated exponential chirps. Inspired by the biosonar signals of bats and dolphins, we optimize the signal profile to maximize the signal-to-noise ratio while minimizing spectral leakage with a carefully-designed modulation of the envelope of the chirp. A combined experimental and numerical investigation reveals that there exists an optimal range of window profiles that minimizes the error with respect to standard single frequency sweep methods. The minimum error is set by the noise floor of the instrument, suggesting that the accuracy of an optimally windowed chirp signal is directly comparable to that achievable with a standard frequency sweep, while the acquisition time can be reduced by up to two orders of magnitude, for comparable spectral content. Finally, we demonstrate the ability of this optimized signal to provide time- and frequency-resolved rheometric data by studying the fast gelation process of an acid-induced protein gel. The use of optimally windowed chirps enables a robust rheological characterization of a wide range of soft materials undergoing rapid mutation and has the potential to become an invaluable tool for researchers across different disciplines.


title: On marginal and conditional parameters in logistic regression models

abstract: A fundamental research question is how much a variation in a covariate influences a binary response variable in a logistic regression model, both directly or through mediators. We derive the exact formula linking the parameters of marginal and conditional regression models with binary mediators when no conditional independence assumptions can be made. The formula has the appealing property of being the sum of terms that vanish whenever parameters of the conditional models vanish, thereby recovering well-known results as particular cases. It also permits to quantify the distortion induced by omission of some relevant covariates, opening the way to sensitivity analysis. Also in this case, as the parameters of the conditional models are multiplied by terms that are always positive or bounded, the formula may be used to construct reasonable bounds on the parameters of interest. We assume that, conditionally on a set of covariates, the data-generating process can be represented by a Directed Acyclic Graph. We also show how the results here presented lead to the extension of path analysis to a system of binary random variables.


title: Charge Regulating Macro-ions in Salt Solutions: Screening Properties and Electrostatic Interactions

abstract: We revisit the charge-regulation mechanism of macro-ions and apply it to mobile macro-ions in a bathing salt solution. In particular, we examine the effects of correlation between various adsorption/desorption sites and analyze the collective behavior in terms of the solution effective screening properties. We show that such a behavior can be quantified in terms of the charge {\em asymmetry} of the macro-ions, defined by their preference for a non-zero effective charge, and their {\em donor/acceptor} propensity for exchanging salt ions with the bathing solution. Asymmetric macro-ions tend to increase the screening, while symmetric macro-ions can in some cases decrease it. Macro-ions that are classified as donors display a rather regular behavior, while those that behave as acceptors exhibit an anomalous non-monotonic Debye length. The screening properties, in their turn, engender important modifications to the disjoining pressure between two charged surfaces. Our findings are in particular relevant for solutions of proteins, whose exposed amino acids can undergo charge dissociation/association processes to/from the bathing solution, and can be considered as a solution of charged regulated macro-ions, as analyzed here.

categories: cond-mat.soft cond-mat.mtrl-sci cond-mat.stat-mech physics.chem-ph


title: Disk stars in the Milky Way detected beyond 25 kpc from its center

abstract: CONTEXT. The maximum size of the Galactic stellar disk is not yet known. Some studies have suggested an abrupt drop-off of the stellar density of the disk at Galactocentric distances $R\gtrsim 15$ kpc, which means that in practice no disk stars or only very few of them should be found beyond this limit. However, stars in the Milky Way plane are detected at larger distances. In addition to the halo component, star counts have placed the end of the disk beyond 20 kpc, although this has not been spectroscopically confirmed so far. AIMS. Here, we aim to spectroscopically confirm the presence of the disk stars up to much larger distances. METHODS. With data from the LAMOST and SDSS-APOGEE spectroscopic surveys, we statistically derived the maximum distance at which the metallicity distribution of stars in the Galactic plane is distinct from that of the halo populations. RESULTS. Our analysis reveals the presence of disk stars at R>26 kpc (99.7% C.L.) and even at R>31 kpc (95.4% C.L.).


title: Efficient Anomaly Detection via Matrix Sketching

abstract: We consider the problem of finding anomalies in high-dimensional data using popular PCA based anomaly scores. The naive algorithms for computing these scores explicitly compute the PCA of the covariance matrix which uses space quadratic in the dimensionality of the data. We give the first streaming algorithms that use space that is linear or sublinear in the dimension. We prove general results showing that \emph{any} sketch of a matrix that satisfies a certain operator norm guarantee can be used to approximate these scores. We instantiate these results with powerful matrix sketching techniques such as Frequent Directions and random projections to derive efficient and practical algorithms for these problems, which we validate over real-world data sets. Our main technical contribution is to prove matrix perturbation inequalities for operators arising in the computation of these measures.


title: Constraining the time variation of Newton's constant $G$ with gravitational-wave standard sirens and supernovae

abstract: The intrinsic peak luminosity of Type Ia supernovae (SNIa) depends on the value of Newton's gravitational constant $G$, through the Chandrasekhar mass $M_{\rm Ch}\propto G^{-3/2}$. If the luminosity distance can be independently determined, the SNIa can be treated as a tracker to constrain the possible time variation of $G$ in different redshift ranges. The gravitational-wave (GW) standard sirens, caused by the coalescence of binary neutron stars, provide a model-independent way to measure the distance of GW events, which can be used to determine the luminosity distances of SNIa by interpolation, provided the GW and SNIa samples have similar redshift ranges. We demonstrate that combining the GW observations of third-generation detectors with SNIa data provides a powerful and model-independent way to measure $G$ in a wide redshift range, which can constrain the ratio $G/G_0$, where $G$ and $G_0$ are respectively the values in the redshift ranges $z>0.1$ and $z<0.1$, at the level of $1.5\%$.


title: Cosmological structure formation with negative mass

abstract: We construct a family of models with negative gravitational mass in the context of Newtonian gravity. We focus in particular on a model that reproduces the features of the so-called Dirac-Milne universe, a matter-antimatter symmetric universe that was recently proposed as an alternative cosmological scenario [A. Benoit-Levy and G. Chardin, A&A 537, A78 (2012)]. We perform one-dimensional N-body simulations of these negative-mass models for an expanding universe and study the associated formation of gravitational structures. The similarities and differences with the standard cosmological model are highlighted and discussed.


title: Robust fusion algorithms for unsupervised change detection between multi-band optical images - A comprehensive case study

abstract: Unsupervised change detection techniques are generally constrained to two multi-band optical images acquired at different times through sensors sharing the same spatial and spectral resolution. This scenario is suitable for a straight comparison of homologous pixels such as pixel-wise differencing. However, in some specific cases such as emergency situations, the only available images may be those acquired through different kinds of sensors with different resolutions. Recently some change detection techniques dealing with images with different spatial and spectral resolutions, have been proposed. Nevertheless, they are focused on a specific scenario where one image has a high spatial and low spectral resolution while the other has a low spatial and high spectral resolution. This paper addresses the problem of detecting changes between any two multi-band optical images disregarding their spatial and spectral resolution disparities. We propose a method that effectively uses the available information by modeling the two observed images as spatially and spectrally degraded versions of two (unobserved) latent images characterized by the same high spatial and high spectral resolutions. Covering the same scene, the latent images are expected to be globally similar except for possible changes in spatially sparse locations. Thus, the change detection task is envisioned through a robust fusion task which enforces the differences between the estimated latent images to be spatially sparse. We show that this robust fusion can be formulated as an inverse problem which is iteratively solved using an alternate minimization strategy. The proposed framework is implemented for an exhaustive list of applicative scenarios and applied to real multi-band optical images. A comparison with state-of-the-art change detection methods evidences the accuracy of the proposed robust fusion-based strategy.


title: K-cut on paths and some trees

abstract: We define the (random) $k$-cut number of a rooted graph to model the difficulty of the destruction of a resilient network. The process is as the cut model of Meir and Moon except now a node must be cut $k$ times before it is destroyed. The first order terms of the expectation and variance of $\mathcal{X}_{n}$, the $k$-cut number of a path of length $n$, are proved. We also show that $\mathcal{X}_{n}$, after rescaling, converges in distribution to a limit $\mathcal{B}_{k}$, which has a complicated representation. The paper then briefly discusses the $k$-cut number of some trees and general graphs. We conclude by some analytic results which may be of interest.


title: Glassy dynamics due to a trajectory phase transition in dissipative Rydberg gases

abstract: The physics of highly excited Rydberg atoms is governed by blockade or exclusion interactions that hinder the excitation of atoms in the proximity of a previously excited one. This leads to cooperative effects and a relaxation dynamics displaying space-time heterogeneity similar to what is observed in the relaxation of glass-forming systems. Here we establish theoretically the existence of a glassy dynamical regime in an open Rydberg gas, associated with phase coexistence at a first-order transition in dynamical large deviation functions. This transition occurs between an active phase of low density in which dynamical processes take place on short timescales, and an inactive phase in which excited atoms are dense and the dynamics is highly arrested. We perform a numerically exact study and develop a mean-field approach that allows to understand the mechanics of this phase transition. We show that radiative decay --- which becomes experimentally relevant for long times --- moves the system away from dynamical phase coexistence. Nevertheless, the dynamical phase transition persists and causes strong fluctuations in the observed dynamics.


title: The correlation constant of a field

abstract: We study the correlation of edges, vectors or elements to be in a randomly chosen spanning tree or a basis, respectively. Here we follow the guideline of Huh and Wang and introduce as a measure an invariant that is called the correlation constant of a graph, vector configuration, matroid or field. It follows from one of their results that these correlation constants are numbers between $0$ and $2$. Here, we show that the correlation constant of every field is at least $\frac{8}{7}$. In our proof we explicitly construct vector configurations and matroids with positively correlated elements.


title: Snapping of Bistable, Prestressed Cylindrical Shells

abstract: Bistable shells can reversibly change between two stable configurations with very little energetic input. Understanding what governs the shape and snap-through criteria of these structures is crucial for designing devices that utilize instability for functionality. Bistable cylindrical shells fabricated by stretching and bonding multiple layers of elastic plates will contain residual stress that will impact the shell's shape and the magnitude of stimulus necessary to induce snapping. Using the framework of non-Euclidean shell theory, we first predict the mean curvature of a nearly cylindrical shell formed by arbitrarily prestretching one layer of a bilayer plate with respect to another. Then, beginning with a residually stressed cylinder, we determine the amount of the stimuli needed to trigger the snapping between two configurations through a combination of numerical simulations and theory. We demonstrate the role of prestress on the snap-through criteria, and highlight the important role that the Gaussian curvature in the boundary layer of the shell plays in dictating shell stability.


title: Constraints on Dark Matter with a moderately large and velocity-dependent DM-nucleon cross-section

abstract: We derive constraints on a possible velocity-dependent DM-nucleon scattering cross section, for Dark Matter in the 10 MeV -- 100 GeV mass range, using the XQC, DAMIC, and CRESST 2017 Surface Run experiments. We report the limits on cross sections of the form $\sigma=\sigma_0\,v^n$, for a range of velocity dependencies with $n\in\{-4,-2,-1,0,1,2\}$. We point out the need to measure the efficiency with which nuclear recoil energy in the sub-keV range thermalizes, rather than being stored as Frenkel pairs in the semi-conductor lattice. The possibility of a significant inefficiency leaves open a considerable `hole' in the limits for mass in the $\sim$ 0.2 -- 2 GeV range, which XQC and CRESST can potentially fill when the thermalization efficiency is measured. We call attention to the asymmetry between a conventional lower limit cross section and the `upper-reach cross section' imposed by attenuation in an overburden -- an upper boundary being extremely sharp but quite insensitive to the statistics of the experiment. Considering the recent interest to use dark matter-baryon interaction with velocity dependence $n=-4$ to explain the EDGES 21 cm anomaly, we also derive the limits on milli-charged DM that scatters off protons and electrons under a Coulomb-like interaction. We find that much but not all of the region of interest for the EDGES anomaly can be excluded.


title: Standard map-like models for single and multiple walkers in an annular cavity

abstract: Recent experiments on walking droplets in an annular cavity showed the existence of complex dynamics including chaotically changing velocity. This article presents models, influenced by the kicked rotator/standard map, for both single and multiple droplets. The models are shown to achieve both qualitative and quantitative agreement with the experiments, and makes predictions about heretofore unobserved behavior. Using dynamical systems techniques and bifurcation theory, the single droplet model is analyzed to prove dynamics suggested by the numerical simulations.


title: Investigating hot-Jupiter inflated radii with hierarchical Bayesian modelling

abstract: As of today, hundreds of hot Jupiters have been found, yet the inflated radii of a large fraction of them remain unexplained. It is still unclear whether a single inflation mechanism is enough to explain the entire distribution of radii, or whether a combination of them is needed. We seek to understand the relationship of hot Jupiter radii with stellar irradiation and mass. We also aim to find the intrinsic physical scatter in their radii caused by unobservable parameters. By constructing a hierarchical Bayesian model, we infer the probabilistic relation between planet radius, mass and incident flux for a sample of 286 gas giants. We separately incorporate the observational uncertainties of the data and the intrinsic physical scatter in the population. This allows us to treat the intrinsic physical scatter in radii (due to latent parameters such as the heavy element fraction) as a parameter to be inferred. We find that the planetary mass plays a key role in the inflation extent, with planets in the range between 0.37 - 0.98 MJ showing the most inflated radii. At higher masses, the radius response to incident flux begins to decrease. Below a threshold of 0.37 +- 0.03 MJ we find that giant exoplanets as a population are unable to maintain inflated radii greater than 1.4 RJ but instead exhibit smaller sizes as the incident flux is increased beyond 10^6 Wm-2. We also find that below 1 MJ, there is a cutoff point at high incident flux beyond which we find no more inflated planets, and that this cutoff point decreases as the mass decreases. At incident fluxes higher than 1.6*10^6 Wm -2 and in a mass range 0.37 - 0.98 MJ, we find no evidence for a population of non-inflated hot Jupiters. Our study demonstrates the importance of population-level analysis to grasp the underlying properties of exoplanets.


title: Applying a Particle-only Model to the HL Tau Disk

abstract: Observations have revealed rich structures in protoplanetary disks, offering clues about their embedded planets. Due to the complexities introduced by the abundance of gas in these disks, modeling their structure in detail is computationally intensive, requiring complex hydrodynamic codes and substantial computing power. It would be advantageous if computationally simpler models could provide some preliminary information on these disks. Here we apply a particle-only model (that we developed for gas-poor debris disks) to the gas-rich disk, HL Tauri, to address the question of whether such simple models can inform the study of these systems. Assuming three potentially embedded planets, we match HL Tau's radial profile fairly well and derive best-fit planetary masses and orbital radii (0.40, 0.02, 0.21 Jupiter masses for the planets orbiting a 0.55 solar-mass star at 11.22, 29.67, 64.23 AU). Our derived parameters are comparable to those estimated by others, except for the mass of the second planet. Our simulations also reproduce some narrower gaps seen in the ALMA image away from the orbits of the planets. The nature of these gaps is debated but, based on our simulations, we argue they could result from planet-disk interactions via mean-motion resonances, and need not contain planets. Our results suggest that a simple particle-only model can be used as a first step to understanding dynamical structures in gas disks, particularly those formed by planets, and determine some parameters of their hidden planets, serving as useful initial inputs to hydrodynamic models which are needed to investigate disk and planet properties more thoroughly.


title: A plug-in approach to maximising precision at the top and recall at the top

abstract: For information retrieval and binary classification, we show that precision at the top (or precision at k) and recall at the top (or recall at k) are maximised by thresholding the posterior probability of the positive class. This finding is a consequence of a result on constrained minimisation of the cost-sensitive expected classification error which generalises an earlier related result from the literature.


title: Unconventional field induced phases in a quantum magnet formed by free radical tetramers

abstract: We report experimental and theoretical studies on the magnetic and thermodynamic properties of NIT-2Py, a free radical-based organic magnet. From magnetization and specific heat measurements we establish the temperature versus magnetic field phase diagram which includes two Bose-Einstein condensates (BEC) and an infrequent half magnetization plateau. Calculations based on density functional theory demonstrates that magnetically this system can be mapped to a quasi-two-dimensional structure of weakly coupled tetramers. Density matrix renormalization group calculations show the unusual characteristics of the BECs where the spins forming the low-field condensate are different than those participating in the high-field one.


title: User Scheduling for Millimeter Wave Hybrid Beamforming Systems with Low-Resolution ADCs

abstract: We investigate uplink user scheduling for millimeter wave (mmWave) hybrid analog/digital beamforming systems with low-resolution analog-to-digital converters (ADCs). Deriving new scheduling criteria for the mmWave systems, we show that the channel structure in the beamspace, in addition to the channel magnitude and orthogonality, plays a key role in maximizing the achievable rates of scheduled users due to quantization error. The criteria show that to maximize the achievable rate for a given channel gain, the channels of the scheduled users need to have (1) as many propagation paths as possible with unique angle-of-arrivals (AoAs) and (2) even power distribution in the beamspace. Leveraging the derived criteria, we propose an efficient scheduling algorithm for mmWave zero-forcing receivers with low-resolution ADCs. We further propose a chordal distance-based scheduling algorithm that exploits only the AoA knowledge and analyze the performance by deriving ergodic rates in closed form. Based on the derived rates, we show that the beamspace channel leakage resulting from phase offsets between AoAs and quantized angles of analog combiners can lead to sum rate gain by reducing quantization error compared to the channel without leakage. Simulation results validate the sum rate performance of the proposed algorithms and derived ergodic rates expressions.


title: Binge Watching: Scaling Affordance Learning from Sitcoms

abstract: In recent years, there has been a renewed interest in jointly modeling perception and action. At the core of this investigation is the idea of modeling affordances(Affordances are opportunities of interaction in the scene. In other words, it represents what actions can the object be used for). However, when it comes to predicting affordances, even the state of the art approaches still do not use any ConvNets. Why is that? Unlike semantic or 3D tasks, there still does not exist any large-scale dataset for affordances. In this paper, we tackle the challenge of creating one of the biggest dataset for learning affordances. We use seven sitcoms to extract a diverse set of scenes and how actors interact with different objects in the scenes. Our dataset consists of more than 10K scenes and 28K ways humans can interact with these 10K images. We also propose a two-step approach to predict affordances in a new scene. In the first step, given a location in the scene we classify which of the 30 pose classes is the likely affordance pose. Given the pose class and the scene, we then use a Variational Autoencoder (VAE) to extract the scale and deformation of the pose. The VAE allows us to sample the distribution of possible poses at test time. Finally, we show the importance of large-scale data in learning a generalizable and robust model of affordances.


title: Routing Game on Parallel Networks: the Convergence of Atomic to Nonatomic

abstract: We consider an instance of a nonatomic routing game. We assume that the network is parallel, that is, constituted of only two nodes, an origin and a destination. We consider infinitesimal players that have a symmetric network cost, but are heterogeneous through their set of feasible strategies and their individual utilities. We show that if an atomic routing game instance is correctly defined to approximate the nonatomic instance, then an atomic Nash Equilibrium will approximate the nonatomic Wardrop Equilibrium. We give explicit bounds on the distance between the equilibria according to the parameters of the atomic instance. This approximation gives a method to compute the Wardrop equilibrium at an arbitrary precision.


title: Attribute-Centered Loss for Soft-Biometrics Guided Face Sketch-Photo Recognition

abstract: Face sketches are able to capture the spatial topology of a face while lacking some facial attributes such as race, skin, or hair color. Existing sketch-photo recognition approaches have mostly ignored the importance of facial attributes. In this paper, we propose a new loss function, called attribute-centered loss, to train a Deep Coupled Convolutional Neural Network (DCCNN) for the facial attribute guided sketch to photo matching. Specifically, an attribute-centered loss is proposed which learns several distinct centers, in a shared embedding space, for photos and sketches with different combinations of attributes. The DCCNN simultaneously is trained to map photos and pairs of testified attributes and corresponding forensic sketches around their associated centers, while preserving the spatial topology information. Importantly, the centers learn to keep a relative distance from each other, related to their number of contradictory attributes. Extensive experiments are performed on composite (E-PRIP) and semi-forensic (IIIT-D Semi-forensic) databases. The proposed method significantly outperforms the state-of-the-art.


title: Multi-Higgs models with CP symmetries of increasingly high order

abstract: When building CP-symmetric models beyond the Standard Model, one can impose CP-symmetry of higher order. This means that one needs to apply the CP-transformation more than two times to get the identity transformation, but still the model is perfectly CP-conserving. A multi-Higgs-doublet model based on CP-symmetry of order 4, dubbed CP4, was recently proposed and its phenomenology is being explored. Here, we show that the construction does not stop at CP4. We build examples of renormalizable multi-Higgs-doublet potentials which are symmetric under CP8 or CP16, without leading to any accidental symmetry. If the vacuum conserves CP-symmetry of order 2k, then the neutral scalars become CP-eigenstates, which are characterized not by CP-parities but by CP-charges defined modulo 2k. One or more lightest states can be the dark matter candidates, which are protected against decay not by the internal symmetry but by the exotic CP. We briefly discuss their mass spectra and interaction patterns for CP8 and CP16.


title: A ZX-Calculus with Triangles for Toffoli-Hadamard, Clifford+T, and Beyond

abstract: We consider a ZX-calculus augmented with triangle nodes which is well-suited to reason on the so-called Toffoli-Hadamard fragment of quantum mechanics. We precisely show the form of the matrices it represents, and we provide an axiomatisation which makes the language complete for the Toffoli-Hadamard quantum mechanics. We extend the language with arbitrary angles and show that any true equation involving linear diagrams which constant angles are multiple of Pi are derivable. We show that a single axiom is then necessary and sufficient to make the language equivalent to the ZX-calculus which is known to be complete for Clifford+T quantum mechanics. As a by-product, it leads to a new and simple complete axiomatisation for Clifford+T quantum mechanics.


title: How costly punishment, diversity, and density of connectivity influence cooperation in a biological network

abstract: It has been an old unsolved puzzle to evolutionary theorists on which mechanisms would increase large-scale cooperation in human societies. Thus, how such mechanisms operate in a biological network is still not very understood. This study addresses these questions with empirical evidence from agent-based models designed to understand these network interactions. Repeated Prisoner's Dilemma games were designed to study how costly punishment, diversity, and density of connectivity interact to influence cooperation in a biological network. There are 1000 rounds in each game made up of 18 players engaged in pairwise relationship with their neighbors. This study shows three important interactions. (1) Introducing diversity to costly punishment favors both cooperation and defection, but not vice versa. Introducing costly punishment to diversity disfavors defection but favors cooperation. (2) Costly Punishment, alone, disfavors defection but decreases average payoff. Decreasing the density of connectivity, Dc, when there is no costly punishment applied, increases average payoff. (3) A synergy of diversity and decreasing density of connectivity favors cooperation in a biological network. Furthermore, this study also suggests a likelihood from empirical findings that spatial structures may not be favoring cooperation, as is the widely-accepted notion, but rather disfavoring defection in the global scale.


title: Localization transitions and mobility edges in coupled Aubry-Andr\'e chains

abstract: We study the localization transitions for coupled one-dimensional lattices with quasiperiodic potential. Besides the localized and extended phases there is an intermediate mixed phase which can be easily explained decoupling the system so as to deal with effective uncoupled Aubry-Andr\'e chains with different transition points. We clarify, therefore, the origin of such an intermediate phase finding the conditions for getting a uniquely defined mobility edge for such coupled systems. Finally we consider many coupled chains with an energy shift which compose an extension of the Aubry-Andr\'e model in two dimensions. We study the localization behavior in this case comparing the results with those obtained for a truly aperiodic two-dimensional (2D) Aubry-Andr\'e model, with quasiperiodic potentials in any directions, and for the 2D Anderson model.


title: Testing the Cosmic Anisotropy with Supernovae Data: Hemisphere Comparison and Dipole Fitting

abstract: The cosmological principle is one of the cornerstones in modern cosmology. It assumes that the universe is homogeneous and isotropic on cosmic scales. Both the homogeneity and the isotropy of the universe should be tested carefully. In the present work, we are interested in probing the possible preferred direction in the distribution of type Ia supernovae (SNIa). To our best knowledge, two main methods have been used in almost all of the relevant works in the literature, namely the hemisphere comparison (HC) method and the dipole fitting (DF) method. However, the results from these two methods are not always approximately coincident with each other. In this work, we test the cosmic anisotropy by using these two methods with the Joint Light-Curve Analysis (JLA) and simulated SNIa datasets. In many cases, both methods work well, and their results are consistent with each other. However, in the cases with two (or even more) preferred directions, the DF method fails while the HC method still works well. This might shed new light on our understanding of these two methods.


title: Application of instruments of nuclear physics to the calculation of theoretical dose distributions in various organs of the human body for beams used in hadrontherapy

abstract: The area of interests of nuclear physics are studies of reactions, wherein atomic nuclei of projectile collide with target nuclei. An amount of energy lost by projectile nucleus during its passing through the target is a major issue - it is important to know how charged particles interact with matter. It is possible to afford this knowledge by using theoretical programs, which calculate energy loss applying Bethe-Bloch equation. Hadrontherapy, which is a field of still growing interest, is based on interactions of charged particles with matter. Therefore exists a need of creating a simple model, which could be used to the calculation of dose distributions in biological matter. Two programs (SRIM, Xeloss), used to the calculation of energy loss by nuclear physicist, have been adapted to determine dose distributions in analogs of human tissues. Results of the calculations with those programs for beams used in hadrontherapy (e.g. $^1$H, $^{12}$C) will be compared with experimental data available in references.


title: Quantifying precision loss in local quantum thermometry via diagonal discord

abstract: When quantum information is spread over a system through nonclassical correlation, it makes retrieving information by local measurements difficult---making global measurement necessary for optimal parameter estimation. In this paper, we consider temperature estimation of a system in a Gibbs state and quantify the separation between the estimation performance of the global optimal measurement scheme and a greedy local measurement scheme by diagonal quantum discord. In a greedy local scheme, instead of global measurements, one performs sequential local measurement on subsystems, which is potentially enhanced by feed-forward communication. We show that, for finite-dimensional systems, diagonal discord quantifies the difference in the quantum Fisher information quantifying the precision limits for temperature estimation of these two schemes, and we analytically obtain the relation in the high-temperature limit. We further verify this result by employing the examples of spins with Heisenberg's interaction.


title: Mean Field Games in the weak noise limit : A WKB approach to the Fokker-Planck equation

abstract: Motivated by the study of a Mean Field Game toy model called the "seminar problem", we consider the Fokker-Planck equation in the small noise regime for a specific drift field. This gives us the opportunity to discuss the application to diffusion problem of the WKB approach "a la Maslov", making it possible to solve directly the time dependant problem in an especially transparent way.


title: Polarizable Molecular Simulations Reveal How Silicon-containing Functional Groups Govern the Desalination Mechanism in Nanoporous Graphene

abstract: We report a molecular dynamics (MD) simulation study of reverse osmosis desalination using nanoporous monolayer graphene passivated by SiH$_2$ and Si(OH)$_2$ functional groups. A highly accurate and detailed polarizable molecular mechanics force field model was developed for simulating graphene nanopores of various sizes and geometries. The simulated water fluxes and ion rejection percentages are explained using detailed atomistic mechanisms derived from analysis of the simulation trajectories. Our main findings are: (1) The Si(OH)$_2$ pores possess superior ion rejection rates due to selective electrostatic repulsion of Cl$^-$ ions, but Na$^+$ ions are attracted to the pore and block water transfer. (2) By contrast, the SiH$_2$ pores operate via a steric mechanism that excludes ions based on the size and flexibility of their hydration layers. (3) In absence of ions, water flux is directly proportional to the solvent accessible area within the pore; however, simulated fluxes are lower than those inferred from recent experimental work. We also provide some hypotheses that could resolve the differences between simulation and experiment.


title: Discrete resonant Rossby/drift wave triads: an explicit parameterisation and a fast direct numerical search algorithm

abstract: We report results on the explicit parameterisation of discrete Rossby-wave resonant triads of the Charney-Hasegawa-Mima equation in the small-scale limit (i.e. large Rossby deformation radius), following up from our previous solution in terms of elliptic curves (Bustamante and Hayat, 2013). We find an explicit parameterisation of the discrete resonant wavevectors in terms of two rational variables. We show that these new variables are restricted to a bounded region and find this region explicitly. We argue that this can be used to reduce the complexity of a direct numerical search for discrete triad resonances. Also, we introduce a new direct numerical method to search for discrete resonances. This numerical method has complexity ${\mathcal{O}}(N^3)$, where $N$ is the largest wavenumber in the search. We apply this new method to find all discrete irreducible resonant triads in the wavevector box of size $5000$, in a calculation that took about $10.5$ days on a $16$-core machine. Finally, based on our method of mapping to elliptic curves, we discuss some dynamical implications regarding the spread of quadratic invariants across scales via resonant triad interactions, in the form of sharp bounds on the size of the interacting wavevectors.


title: BANYAN. XII. New Members of Nearby Young Associations from Gaia-Tycho Data

abstract: We present a search for stellar members of young associations within 150 pc of the Sun based on Gaia-Tycho and an updated version of the BANYAN $\Sigma$ software to determine Bayesian membership probabilities that includes Gaia-2MASS color-magnitude diagrams. We identify 32 new F0-M3-type bona fide members of the 10-200 Myr-old Sco-Cen, Carina, Tucana-Horologium, Columba and Octans associations and the AB Doradus, $\beta$ Pictoris and Carina-Near moving groups. These new bona fide members have measurements of their full kinematics and literature data consistent with a young age. We also confirm the membership of 66 previously known candidate members using their Gaia-Tycho trigonometric distances or new literature radial velocities, and identify 219 additional new candidate members, most of which do not yet have a radial velocity measurement. This work is the first step towards a completeness-corrected survey of young association members based on Gaia-DR2 in the near future.


title: Baryon masses and \sigma-terms in SU(3) BChPT x 1/Nc

abstract: Baryon masses and nucleon \sigma-terms are studied with the effective theory that combines the chiral and 1/Nc expansions for three flavors. In particular the connection between the deviation of the Gell-Mann-Okubo relation and the \sigma-term corresponding to the scalar density associated with the hypercharge is emphasized. The latter is at lowest order related to a mass combination whose low value has given rise to a \sigma- term puzzle. It is shown that while the nucleon \sigma-terms have a well behaved low energy expansion, that mass combination is affected by large higher order corrections non-analytic in quark masses. Adding to the analysis lattice QCD baryon masses, it is found that {\sigma}{\pi}N=69(10) MeV and {\sigma}s has natural magnitude within its relative large uncertainty.


title: Exact non-Markovian dynamics of Gaussian quantum channels: Finite-time and asymptotic regimes

abstract: We investigate the Markovian and non-Markovian dynamics of Gaussian quantum channels, exploiting a recently introduced necessary and sufficient criterion and the ensuing measure of non-Markovianity based on the violation of the divisibility property of the dynamical map. We compare the paradigmatic instances of Quantum Brownian motion (QBM) and Pure Damping (PD) channels, and for the former we find that the exact dynamical evolution is always non-Markovian in the finite-time as well as in the asymptotic regimes, for any nonvanishing value of the non-Markovianity parameter. If one resorts to the rotating wave approximated (RWA) form of the QBM, that neglects the anomalous diffusion contribution to the system dynamics, we show that such approximation fails to detect the non-Markovian nature of the dynamics. Finally, for the exact dynamics of the QBM in the asymptotic regime, we show that the quantifiers of non-Markovianity based on the distinguishability between quantum states fail to detect the non-Markovian nature of the dynamics.

categories: quant-ph cond-mat.other math-ph math.MP physics.optics


title: Atomic Physics Constraints on the X Boson

abstract: Recently, a peak in the light fermion pair spectrum at invariant q^2 approximately equal to 16.7 Me^2. has been observed in the bombardment of Li-7 by protons. This peak has been interpreted in terms of a protophobic interaction of fermions with a gauge boson (X boson) of invariant mass of approximately 16.7 MeV which couples mainly to neutrons. High-precision atomic physics experiments aimed at observing the protophobic interaction need to separate the X boson effect from the nuclear-size effect, which is a problem because of the short range of the interaction 11.8 fm, which is commensurate with a "nuclear halo". Here, we analyze the X boson in terms of its consequences for both electronic atoms as well as muonic hydrogen and deuterium. We find that the most promising atomic systems where the X boson has an appreciable effect, distinguishable from a finite-nuclear-size effect, are muonic atoms of low and intermediate nuclear charge numbers.


title: The Connection between Galaxies and their Dark Matter Halos

abstract: In our modern understanding of galaxy formation, every galaxy forms within a dark matter halo. The formation and growth of galaxies over time is connected to the growth of the halos in which they form. The advent of large galaxy surveys as well as high-resolution cosmological simulations has provided a new window into the statistical relationship between galaxies and halos and its evolution. Here we define this galaxy-halo connection as the multi-variate distribution of galaxy and halo properties that can be derived from observations and simulations. This connection provides a key test of physical galaxy formation models; it also plays an essential role in constraints of cosmological models using galaxy surveys and in elucidating the properties of dark matter using galaxies. We review techniques for inferring the galaxy-halo connection and the insights that have arisen from these approaches. Some things we have learned are that galaxy formation efficiency is a strong function of halo mass; at its peak in halos around a pivot halo mass of 10^12 Msun, less than 20% of the available baryons have turned into stars by the present day; the intrinsic scatter in galaxy stellar mass is small, less than 0.2 dex at a given halo mass above this pivot mass; below this pivot mass galaxy stellar mass is a strong function of halo mass; the majority of stars over cosmic time were formed in a narrow region around this pivot mass. We also highlight key open questions about how galaxies and halos are connected, including understanding the correlations with secondary properties and the connection of these properties to galaxy clustering.


title: Stochastic Comparisons of Two-Units Repairable Systems

abstract: We consider two models of two-units repairable systems: cold standby system and warm standby system. We suppose that the lifetimes and repair times of the units are all independent exponentially distributed random variables. Using stochastic orders we compare the lifetimes of systems under different assumptions on the parameters of exponential distributions. We also consider a cold standby system where the lifetimes and repair times of its units are not necessarily exponentially distributed.


title: Hidden entropy production and work fluctuations in an ideal active gas

abstract: Collections of self-propelled particles that move persistently by continuously consuming free energy are a paradigmatic example of active matter. In these systems, unlike Brownian "hot colloids", the breakdown of detailed balance yields a continuous production of entropy at steady state, even for an ideal active gas. We quantify the irreversibility for a non-interacting active particle in two dimensions by treating both conjugated and time-reversed dynamics. By starting with underdamped dynamics, we identify a hidden rate of entropy production required to maintain persistence and prevent the rapidly relaxing momenta from thermalizing, even in the limit of very large friction. Additionally, comparing two popular models of self-propulsion with identical dissipation on average, we find that the fluctuations and large deviations in work done are markedly different, providing thermodynamic insight into the varying extents to which macroscopically similar active matter systems may depart from equilibrium.


title: Incorporating Surface Convection into a 3D Babcock-Leighton Solar Dynamo Model

abstract: The observed convective flows on the photosphere (e.g., supergranulation, granulation) play a key role in the Babcock-Leighton (BL) process to generate large-scale polar fields from sunspots fields. In most surface flux transport (SFT) and BL dynamo models, the dispersal and migration of surface fields is modeled as an effective turbulent diffusion. Recent SFT models have incorporated explicit, realistic convective flows in order to improve the fidelity of convective transport but, to our knowledge, this has not yet been implemented in previous BL models. Since most Flux-Transport (FT)/BL models are axisymmetric, they do not have the capacity to include such flows. We present the first kinematic 3D FT/BL model to explicitly incorporate realistic convective flows based on solar observations. Though we describe a means to generalize these flows to 3D, we find that the kinematic small-scale dynamo action they produce disrupts the operation of the cyclic dynamo. Cyclic solution is found by limiting the convective flow to act only on the vertical radial component of the magnetic field. The results obtained are generally in good agreement with the observed surface flux evolution and with non-convective models that have a turbulent diffusivity on the order of $3 \times 10^{12}$ cm$^2$ s$^{-1}$ (300 km$^2$ s$^{-1}$). However, we find that the use of a turbulent diffusivity underestimates the dynamo efficiency, producing weaker mean fields and shorter cycle than in the convective models. Also, the convective models exhibit mixed polarity bands in the polar regions that have no counterpart in solar observations. Also, the explicitly computed turbulent electromotive force (emf) bears little resemblance to a diffusive flux. We also find that the poleward migration speed of poloidal flux is determined mainly by the meridional flow and the vertical diffusion.


title: Total gravitational mass of the Galactic Double Neutron Star systems: evidence for a bimodal distribution

abstract: So far, in total 15 double neutron star systems (DNSs) with a reliable measurement of the total gravitational mass ($M_{\rm T}$) have been detected in the Galaxy. In this work we study the distribution of $M_{\rm T}$. The data prefer the double Gaussian distribution over a single Gaussian distribution and the low and high mass populations center at $M_{\rm T}\sim 2.58M_\odot$ and $\sim 2.72M_\odot$, respectively. The progenitor stars of GW170817 have a $M_{\rm T}=2.74^{+0.04}_{-0.01}M_\odot$, falling into the high mass population. With a local neutron star merger rate of $\sim 10^{3}~{\rm Gpc^{-3}~yr^{-1}}$, supposing the $M_{\rm T}$ of those merging neutron stars also follow the double Gaussian distribution, the upcoming runs of the advanced LIGO/Virgo will soon detect some events with a $M_{\rm T}\lesssim 2.6M_\odot$ that can effectively probe the equation of state of the neutron stars and the distribution function is expected to be reliably reconstructed in the next decade.


title: Dynamic finite-size scaling after a quench at quantum transitions

abstract: We present a general dynamic finite-size scaling theory for the quantum dynamics after an abrupt quench, at both continuous and first-order quantum transitions. For continuous transitions, the scaling laws are naturally ruled by the critical exponents and the renormalization-group dimension of the perturbation at the transition. In the case of first-order transitions, it is possible to recover a universal scaling behavior, which is controlled by the size behavior of the energy gap between the lowest energy levels. We discuss these findings in the framework of the paradigmatic quantum Ising ring, and support the dynamic scaling laws by numerical evidence.


title: The volumetric rate of calcium-rich transients in the local universe

abstract: We present a measurement of the volumetric rate of `calcium-rich' optical transients in the local universe, using a sample of three events from the Palomar Transient Factory (PTF). This measurement builds on a detailed study of the PTF transient detection efficiencies, and uses a Monte Carlo simulation of the PTF survey. We measure the volumetric rate of calcium-rich transients to be higher than previous estimates: $1.21^{+1.13}_{-0.39}\times10^{-5}$ events yr$^{-1}$ Mpc$^{-3}$. This is equivalent to 33-94% of the local volumetric type Ia supernova rate. This calcium-rich transient rate is sufficient to reproduce the observed calcium abundances in galaxy clusters, assuming an asymptotic calcium yield per calcium-rich event of ~0.05$\mathrm{M}_{\odot}$. We also study the PTF detection efficiency of these transients as a function of position within their candidate host galaxies. We confirm as a real physical effect previous results that suggest calcium-rich transients prefer large physical offsets from their host galaxies.


title: Chromo-Natural inflation in Supergravity

abstract: We present supergravity realizations of chromo-natural inflationary models. We show that by using superpotentials with "imaginary" holomorphic functions of the inflaton one can obtain effective theories of inflation that are also consistent truncations of the original supergravity models.


title: Central limit theorems via Stein's method for randomized experiments under interference

abstract: We study conditions under which treatment effect estimators constructed under the no-interference assumption in randomized experiments are asymptotically normal in the presence of interference. We prove that the standard Horvitz-Thompson estimator is asymptotically normal under a restricted interference condition characterized by limiting the degree of the dependency graph. The amount of interference is allowed to grow with the population size. We then provide a central limit theorem for the difference-in-means estimator that can handle interference that exists between all pairs of units, provided most of the interference is captured by a restricted-degree dependency graph. The asymptotic variance admits a decomposition into two terms: (a) the variance that is expected under no-interference and (b) the additional variance contributed by interference. We propose a conservative variance estimator based on this variance decomposition. The results arise as an application of Stein's method. For practitioners, our results show that standard estimators continue to exhibit normality in large sample sizes and that inference can be made robust to mild forms of interference.


title: Optimal Approximation by $sk$-Splines on the Torus

abstract: Fixed a continuous kernel K on the $d$-dimensional torus, we consider a generalization of the univariate $sk$-spline to the torus, associated with the kernel K. It is proved an estimate which provides the rate of convergence of a given function by its interpolating $sk$-splines, in the norm of $L^q$ for functions of the type $f=K*\varphi$ where $\varphi \in L^p$ and $1\leq p \leq 2 \leq q \leq \infty,\ 1/p - 1/q \geq 1/2$. The rate of convergence is obtained for functions f in Sobolev classes and this rate gives optimal error estimate of the same order as best trigonometric approximation, in a special case.


title: Signature of Fermi arc surface states in Andreev reflection at the WTe$_2$ Weyl semimetal surface

abstract: We experimentally investigate charge transport through the interface between a niobium superconductor and a three-dimensional WTe$_2$ Weyl semimetal. In addition to classical Andreev reflection, we observe sharp non-periodic subgap resistance resonances. From an analysis of their positions, magnetic field and temperature dependencies, we can interpret them as an analog of Tomasch oscillations for transport along the topological surface state across the region of proximity-induced superconductivity at the Nb-WTe$_2$ interface. Observation of distinct geometrical resonances implies a specific transmission direction for carriers, which is a hallmark of the Fermi arc surface states.


title: Optimal Transport over Deterministic Discrete-time Nonlinear Systems using Stochastic Feedback Laws

abstract: This paper considers the relaxed version of the transport problem for general nonlinear control systems, where the objective is to design time-varying feedback laws that transport a given initial probability measure to a target probability measure under the action of the closed-loop system. To make the problem analytically tractable, we consider control laws that are stochastic, i.e., the control laws are maps from the state space of the control system to the space of probability measures on the set of admissible control inputs. Under some controllability assumptions on the control system as defined on the state space, we show that the transport problem, considered as a controllability problem for the lifted control system on the space of probability measures, is well-posed for a large class of initial and target measures. We use this to prove the well-posedness of a fixed-endpoint optimal control problem defined on the space of probability measures, where along with the terminal constraints, the goal is to optimize an objective functional along the trajectory of the control system. This optimization problem can be posed as an infinite-dimensional linear programming problem. This formulation facilitates numerical solutions of the transport problem for low-dimensional control systems, as we show in two numerical examples.


title: Tensor Mixed Effects Model with Applications in Nanomanufacturing Inspection

abstract: Raman mapping technique has been used to perform in-line quality inspections of nanomanufacturing processes. In such an application, massive high-dimensional Raman mapping data with mixed effects is generated. In general, fixed effects and random effects in the multi-array Raman data are associated with different quality characteristics such as fabrication consistency, uniformity, defects, et al. The existing tensor decomposition methods cannot separate mixed effects, and existing mixed effects model can only handle matrix data but not high-dimensional multi-array data. In this paper, we propose a tensor mixed effects (TME) model to analyze massive high-dimensional Raman mapping data with complex structure. The proposed TME model can (i) separate fixed effects and random effects in a tensor domain; (ii) explore the correlations along different dimensions; and (iii) realize efficient parameter estimation by a proposed iterative double Flip-Flop algorithm. We also investigate the properties of the TME model, existence and identifiability of parameter estimation. The numerical analysis demonstrates the efficiency and accuracy of the parameter estimation in the TME model. Convergence and asymptotic properties are discussed in the simulation and surrogate data analysis. The case study shows an application of the TME model in quantifying the influence of alignment on carbon nanotubes buckypaper. Moreover, the TME model can be applied to provide potential solutions for a family of tensor data analytics problems with mixed effects.


title: Varying Random Coefficient Models

abstract: This paper provides a new methodology to analyze unobserved heterogeneity when observed characteristics are modeled nonlinearly. The proposed model builds on varying random coefficients (VRC) that are determined by nonlinear functions of observed regressors and additively separable unobservables. This paper proposes a novel estimator of the VRC density based on weighted sieve minimum distance. The main example of sieve bases are Hermite functions which yield a numerically stable estimation procedure. This paper shows inference results that go beyond what has been shown in ordinary RC models. We provide in each case rates of convergence and also establish pointwise limit theory of linear functionals, where a prominent example is the density of potential outcomes. In addition, a multiplier bootstrap procedure is proposed to construct uniform confidence bands. A Monte Carlo study examines finite sample properties of the estimator and shows that it performs well even when the regressors associated to RC are far from being heavy tailed. Finally, the methodology is applied to analyze heterogeneity in income elasticity of demand for housing.


title: Analysis of the angular dependence of time delay in gravitational lensing

abstract: We consider an alternative formula for time delay in gravitational lensing. Imposing a smoothness condition on the gravitationally deformed paths followed by the photons from the source to the observer, we show that our formula displays the same degrees of freedom as the standard one. In addition to this, it is shown that the standard expression for time delay is recovered when small angles are involved. These two features strongly support the claim that the formula for time delay studied in this paper is the generalization to the arbitrary angles of the standard one, which is valid at small angles. This could therefore result in a useful tool in Astrophysics and Cosmology which may be applied to investigate the discrepancy between the various estimates of the Hubble constant. As an aside, two interesting consequences of our proposal for time delay are discussed: the existence of a constraint on the gravitational potential generated by the lens and a formula for the mass of the lens in the case of central potential.


title: Beating the integrality ratio for s-t-tours in graphs

abstract: Among various variants of the traveling salesman problem, the s-t-path graph TSP has the special feature that we know the exact integrality ratio, 3/2, and an approximation algorithm matching this ratio. In this paper, we go below this threshold: we devise a polynomial-time algorithm for the s-t-path graph TSP with approximation ratio 1.497. Our algorithm can be viewed as a refinement of the 3/2-approximation algorithm by Seb\H{o} and Vygen [2014], but we introduce several completely new techniques. These include a new type of ear-decomposition, an enhanced ear induction that reveals a novel connection to matroid union, a stronger lower bound, and a reduction of general instances to instances in which s and t have small distance (which works for general metrics).


title: From quantum unstable systems to the decaying dark energy: Cosmological implications

abstract: We consider a cosmology with decaying metastable dark energy and assume that a decay process of this metastable dark energy is a quantum decay process. Such an assumption implies among others that the evolution of the Universe is irreversible and violates the time reversal symmetry. We show that if to replace the cosmological time $t$ appearing in the equation describing the evolution of the Universe by the Hubble cosmological scale time, then we obtain time dependent $\Lambda (t)$ in the form of the series of even powers of the Hubble parameter $H$: $\Lambda (t) = \Lambda (H)$. Out special attention is focused on radioactive like exponential form of the decay process of the dark energy and on the consequences of this type decay.