title: Hysteretic superconducting heat-flux quantum modulator

abstract: We discuss heat transport in a thermally-biased SQUID in the presence of an external magnetic flux, when a non-negligible inductance of the SQUID ring is taken into account. A properly sweeping driving flux causes the thermal current to modulate and behave hysteretically. The response of this device is analysed as a function of both the hysteresis parameter and degree of asymmetry of the SQUID, highlighting the parameter range over which hysteretic behavior is observable. Markedly, also the temperature of the SQUID shows hysteretic evolution, with sharp transitions characterized by temperature jumps up to, e.g., ~ 38mK for a realistic Al-based setup. In view of these results, the proposed device can effectively find application as a temperature-based superconducting memory element, working even at GHz frequencies by suitably choosing the superconductor on which the device is based.


title: Different origins or different evolutions? Decoding the spectral diversity among C-type asteroids

abstract: Anhydrous pyroxene-rich interplanetary dust particles (IDPs) have been proposed as surface analogs for about two-thirds of all C-complex asteroids. However, this suggestion appears to be inconsistent with the presence of hydrated silicates on the surfaces of some of these asteroids including Ceres. Here we report the presence of enstatite (pyroxene) on the surface of two C-type asteroids (Ceres and Eugenia) based on their spectral properties in the mid-infrared range. The presence of this component is particularly unexpected in the case of Ceres because most thermal evolution models predict a surface consisting of hydrated compounds only. The most plausible scenario is that Ceres' surface has been partially contaminated by exogenous enstatite-rich material, possibly coming from the Beagle asteroid family. This scenario questions a similar origin for Ceres and the remaining C-types, and it possibly supports recent results obtained by the Dawn mission (NASA) that Ceres may have formed in the very outer solar system. Concerning the smaller C-types such as Eugenia, both their derived surface composition (enstatite and amorphous silicates) and low density suggest that these bodies accreted from the same building blocks, namely chondritic porous, pyroxene-rich IDPs and volatiles (mostly water ice), and that a significant volume fraction of these bodies has remained unaffected by hydrothermal activity likely implying a late accretion. In addition, their current heliocentric distance may best explain the presence or absence of water ice at their surfaces. Finally, we raise the possibility that CI chondrites, Tagish Lake-like material, or hydrated IDPs may be representative samples of the cores of these bodies.


title: Ultra-relativistic double explosions

abstract: We consider fluid dynamics of relativistic double explosion - when a point explosion with energy $E_1$ is followed by a second explosion with energy $E_2$ after time $t_d$ (the second explosion could be in a form of a long lasting wind). The primary explosion creates a self-similar relativistic blast wave propagating with \Lf\ $\Gamma_1(t)$. A sufficiently strong second explosion, with total energy $E_2 \geq 10^{-2} E_1$, creates a fast second shock in the external fluid previously shocked by the primary shock. At times longer than the interval between the explosions $t_d$, yet short compared with the time when the second shock catches up the primary shock at $\sim t_d \Gamma_1^2$, the structure of the second shock is approximately self-similar. Self-similar structure of the second shock exist for the case of constant external density (in this case $\Gamma_2 \propto t^{-7/3}$), but not for the wind environment. At early times the \Lf\ of the second shock may exceed that of the primary shock and may boost the synchrotron emission of locally accelerated electrons into the Fermi LAT range.


title: Identifying Nonlinear 1-Step Causal Influences in Presence of Latent Variables

abstract: We propose an approach for learning the causal structure in stochastic dynamical systems with a $1$-step functional dependency in the presence of latent variables. We propose an information-theoretic approach that allows us to recover the causal relations among the observed variables as long as the latent variables evolve without exogenous noise. We further propose an efficient learning method based on linear regression for the special sub-case when the dynamics are restricted to be linear. We validate the performance of our approach via numerical simulations.


title: Intersection Cuts with Infinite Split Rank

abstract: We consider mixed integer linear programs where free integer variables are expressed in terms of nonnegative continuous variables. When this model only has two integer variables, Dey and Louveaux characterized the intersection cuts that have infinite split rank. We show that, for any number of integer variables, the split rank of an intersection cut generated from a bounded convex set $P$ is finite if and only if the integer points on the boundary of $P$ satisfy a certain "2-hyperplane property". The Dey-Louveaux characterization is a consequence of this more general result.


title: Stable Recovery Of Sparse Vectors From Random Sinusoidal Feature Maps

abstract: Random sinusoidal features are a popular approach for speeding up kernel-based inference in large datasets. Prior to the inference stage, the approach suggests performing dimensionality reduction by first multiplying each data vector by a random Gaussian matrix, and then computing an element-wise sinusoid. Theoretical analysis shows that collecting a sufficient number of such features can be reliably used for subsequent inference in kernel classification and regression. In this work, we demonstrate that with a mild increase in the dimension of the embedding, it is also possible to reconstruct the data vector from such random sinusoidal features, provided that the underlying data is sparse enough. In particular, we propose a numerically stable algorithm for reconstructing the data vector given the nonlinear features, and analyze its sample complexity. Our algorithm can be extended to other types of structured inverse problems, such as demixing a pair of sparse (but incoherent) vectors. We support the efficacy of our approach via numerical experiments.


title: Linear correlations of the divisor function

abstract: Motivated by arithmetic applications on the number of points in a bihomogeneous variety and on moments of Dirichlet $L$-functions, we provide analytic continuation for the series $\mathcal A_{\boldsymbol{a}}(s):=\sum_{n_1,\dots,n_k\geq1}\frac{d(n_1)\cdots d(n_k)}{(n_1\cdots n_k)^{s}}$ with the sum restricted to solutions of a non-trivial linear equation $a_1n_1+\cdots+a_kn_k=0$. The series $\mathcal A_{\boldsymbol{a}}(s)$ converges absolutely for $\Re(s)>1-\frac1k$ and we show it can be meromorphically continued to $\Re(s)>1-\frac 2{k+1}$ with poles at $s=1-\frac1{k-j}$ only, for $1\leq j< (k-1)/2$. As an application, we obtain an asymptotic formula with power saving error term for the number of points in the variety $a_1x_1y_1+\cdots+a_kx_ky_k=0$ in $\mathbb P^{k-1}(\mathbb Q)\times \mathbb P^{k-1}(\mathbb Q)$.


title: On attainability of optimal controls in coefficients for system of Hammerstein type with anisotropic p-Laplacia

abstract: In this paper we consider an optimal control problem for the coupled system of a nonlinear monotone Dirichlet problem with anisotropic p-Laplacian and matrix-valued nonsmooth controls in its coefficients and a nonlinear equation of Hammerstein type. Using the direct method in calculus of variations, we prove the existence of an optimal control in considered problem and provide sensitivity analysis for a specific case of considered problem with respect to two-parameter regularization.


title: The Augustin Center and The Sphere Packing Bound For Memoryless Channels

abstract: For any channel with a convex constraint set and finite Augustin capacity, existence of a unique Augustin center and associated Erven-Harremoes bound are established. Augustin-Legendre capacity, center, and radius are introduced and proved to be equal to the corresponding Renyi-Gallager entities. Sphere packing bounds with polynomial prefactors are derived for codes on two families of channels: (possibly non-stationary) memoryless channels with multiple additive cost constraints and stationary memoryless channels with convex constraints on the empirical distribution of the input codewords.


title: Shape stability of optimal control problems in coefficients for coupled system of Hammerstein type

abstract: In this paper we consider an optimal control problem (OCP) for the coupled system of a nonlinear monotone Dirichlet problem with matrix- valued non-smooth controls in coefficients and a nonlinear equation of Ham- merstein type. Since problems of this type have no solutions in general, we make a special assumption on the coefficients of the state equation and in- troduce the class of so-called solenoidal admissible controls. Using the direct method in calculus of variations, we prove the existence of an optimal control. We also study the stability of the optimal control problem with respect to the domain perturbation. In particular, we derive the sufficient conditions of the Mosco-stability for the given class of OCPs.


title: Evaluating QBF Solvers: Quantifier Alternations Matter

abstract: We present an experimental study of the effects of quantifier alternations on the evaluation of quantified Boolean formula (QBF) solvers. The number of quantifier alternations in a QBF in prenex conjunctive normal form (PCNF) is directly related to the theoretical hardness of the respective QBF satisfiability problem in the polynomial hierarchy. We show empirically that the performance of solvers based on different solving paradigms substantially varies depending on the numbers of alternations in PCNFs. In related theoretical work, quantifier alternations have become the focus of understanding the strengths and weaknesses of various QBF proof systems implemented in solvers. Our results motivate the development of methods to evaluate orthogonal solving paradigms by taking quantifier alternations into account. This is necessary to showcase the broad range of existing QBF solving paradigms for practical QBF applications. Moreover, we highlight the potential of combining different approaches and QBF proof systems in solvers.


title: Polaron resonances in two vertically stacked quantum dots

abstract: In this work, we present a theoretical study of polaron states in a double quantum dot system. We present realistic calculations which combine 8 band \kp model, configuration interaction approach and collective modes method. We investigate the dependence of polaron energy branches on axial electric field. We show that coupling between carriers and longitudinal optical phonons via Fr\"{o}hlich interaction leads to qualitative and quantitative reconstruction of the optical spectra. In particular, we study the structure of resonances between the states localized in different dots. We show that $p$-shell states are strongly coupled to the phonon replicas of $s$-shell states, in contrast to the weak direct $s$-$p$ coupling. We discuss also the dependence of the phonon-assisted tunnel coupling strength on the separation between the dots.


title: A congruence sum and rational approximations

abstract: We give a reciprocity formula for a two-variable sum where the variables satisfy a linear congruence condition. We also prove that such sum is a measure of how well a rational is approximable from below and show that the reciprocity formula is a simple consequence of this fact.


title: Evolutionary games on cycles with strong selection

abstract: Evolutionary games on graphs describe how strategic interactions and population structure determine evolutionary success, quantified by the probability that a single mutant takes over a population. Graph structures, compared to the well-mixed case, can act as amplifiers or suppressors of selection by increasing or decreasing the fixation probability of a beneficial mutant. Properties of the associated mean fixation times can be more intricate, especially when selection is strong. The intuition is that fixation of a beneficial mutant happens fast (in a dominance game), that fixation takes very long (in a coexistence game), and that strong selection eliminates demographic noise. Here we show that these intuitions can be misleading in structured populations. We analyze mean fixation times on the cycle graph under strong frequency-dependent selection for two different microscopic evolutionary update rules (death-birth and birth-death). We establish exact analytical results for fixation times under strong selection, and show that there are coexistence games in which fixation occurs in time polynomial in population size. Depending on the underlying game, we observe inherence of demographic noise even under strong selection, if the process is driven by random death before selection for birth of an offspring (death-birth update). In contrast, if selection for an offspring occurs before random removal (birth-death update), strong selection can remove demographic noise almost entirely.


title: Click-based porous cationic polymers for enhanced carbon dioxide capture

abstract: Imidazolium based porous cationic polymers were synthesized using an innovative and facile approach, which takes advantage of the Debus Radziszewski reaction to obtain meso- and microporous polymers following click chemistry principles. In the obtained set of materials, click based porous cationic polymers have the same cationic backbone, whereas they bear the commonly used anions of imidazolium poly(ionic liquid)s. These materials show hierarchical porosity and a good specific surface area. Furthermore, their chemical structure was extensively characterized using ATR FTIR and SS NMR spectroscopies, and HR MS. These polymers show good performance towards carbon dioxide sorption, especially those possessing the acetate anion. This polymer has an uptake of 2 mmol per g of CO2 at 1 bar and 273 K, a value which is among the highest recorded for imidazolium poly(ionic liquid)s. These polymers were also modified in order to introduce N-heterocyclic carbenes along the backbone. Carbon dioxide loading in the carbene-containing polymer is in the same range as that of the non-modified versions, but the nature of the interaction is substantially different. The combined use of in situ FTIR spectroscopy and microcalorimetry evidenced a chemisorption phenomenon that brings about the formation of an imidazolium carboxylate zwitterion.


title: On the regularity of the covariance matrix of a discretized scalar field on the sphere

abstract: We present a comprehensive study of the regularity of the covariance matrix of a discretized field on the sphere. In a particular situation, the rank of the matrix depends on the number of pixels, the number of spherical harmonics, the symmetries of the pixelization scheme and the presence of a mask. Taking into account the above mentioned components, we provide analytical expressions that constrain the rank of the matrix. They are obtained by expanding the determinant of the covariance matrix as a sum of determinants of matrices made up of spherical harmonics. We investigate these constraints for five different pixelizations that have been used in the context of Cosmic Microwave Background (CMB) data analysis: Cube, Icosahedron, Igloo, GLESP and HEALPix, finding that, at least in the considered cases, the HEALPix pixelization tends to provide a covariance matrix with a rank closer to the maximum expected theoretical value than the other pixelizations. The effect of the propagation of numerical errors in the regularity of the covariance matrix is also studied for different computational precisions, as well as the effect of adding a certain level of noise in order to regularize the matrix. In addition, we investigate the application of the previous results to a particular example that requires the inversion of the covariance matrix: the estimation of the CMB temperature power spectrum through the Quadratic Maximum Likelihood algorithm. Finally, some general considerations in order to achieve a regular covariance matrix are also presented.


title: Subcritical convection in an internally heated layer

abstract: Thermal convection in a horizontal layer with uniform internal heating and stress-free constant-temperature boundaries is analyzed numerically. The work is motivated by the questions arising in development of liquid metal batteries, in which convection is induced by the Joule heating of electrolyte. It is demonstrated that three-dimensional convection cells exist at subcritical Rayleigh numbers.


title: Bayesian Inference in the Presence of Intractable Normalizing Functions

abstract: Models with intractable normalizing functions arise frequently in statistics. Common examples of such models include exponential random graph models for social networks and Markov point processes for ecology and disease modeling. Inference for these models is complicated because the normalizing functions of their probability distributions include the parameters of interest. In Bayesian analysis they result in so-called doubly intractable posterior distributions which pose significant computational challenges. Several Monte Carlo methods have emerged in recent years to address Bayesian inference for such models. We provide a framework for understanding the algorithms and elucidate connections among them. Through multiple simulated and real data examples, we compare and contrast the computational and statistical efficiency of these algorithms and discuss their theoretical bases. Our study provides practical recommendations for practitioners along with directions for future research for MCMC methodologists.


title: Time evolution of the Luttinger model with nonuniform temperature profile

abstract: We study the time evolution of a one-dimensional interacting fermion system described by the Luttinger model starting from a nonequilibrium state defined by a smooth temperature profile $T(x)$. As a specific example we consider the case when $T(x)$ is equal to $T_L$ ($T_R$) far to the left (right). Using a series expansion in $\epsilon = 2(T_{R} - T_{L})/(T_{L}+T_{R})$, we compute the energy density, the heat current density, and the fermion two-point correlation function for all times $t \geq 0$. For local (delta-function) interactions, the first two are computed to all orders, giving simple exact expressions involving the Schwarzian derivative of the integral of $T(x)$. For nonlocal interactions, breaking scale invariance, we compute the nonequilibrium steady state (NESS) to all orders and the evolution to first order in $\epsilon$. The heat current in the NESS is universal even when conformal invariance is broken by the interactions, and its dependence on $T_{L,R}$ agrees with numerical results for the $XXZ$ spin chain. Moreover, our analytical formulas predict peaks at short times in the transition region between different temperatures and show dispersion effects that, even if nonuniversal, are qualitatively similar to ones observed in numerical simulations for related models, such as spin chains and interacting lattice fermions.


title: A Counterexample to a Conjecture of Gomory and Johnson

abstract: In Mathematical Programming 2003, Gomory and Johnson conjecture that the facets of the infinite group problem are always generated by piecewise linear functions. In this paper we give an example showing that the Gomory-Johnson conjecture is false.


title: Radio spectra of bright compact sources at z>4.5

abstract: High-redshift quasars are important to study galaxy and active galactic nuclei (AGN) evolution, test cosmological models, and study supermassive black hole growth. Optical searches for high-redshift sources have been very successful, but radio searches are not hampered by dust obscuration and should be more effective at finding sources at even higher redshifts. Identifying high-redshift sources based on radio data is, however, not trivial. Here we report on new multi-frequency Giant Metrewave Radio Telescope (GMRT) observations of eight z>4.5 sources previously studied at high angular resolution with very long baseline interferometry (VLBI). Combining these observations with those from the literature, we construct broad-band radio spectra of all 30 z>4.5 sources that have been observed with VLBI. In the sample we found flat, steep and peaked spectra in approximately equal proportions. Despite several selection effects, we conclude that the z>4.5 VLBI (and likely also non-VLBI) sources have diverse spectra and that only about a quarter of the sources in the sample have flat spectra. Previously, the majority of high-redshift radio sources were identified based on their ultra-steep spectra (USS). Recently a new method has been proposed to identify these objects based on their megahertz-peaked spectra (MPS). Neither method would have identified more than 18% of the high-redshift sources in this sample. More effective methods are necessary to reliably identify complete samples of high-redshift sources based on radio data.


title: $H^\infty$-calculus for semigroup generators on BMO

abstract: We prove that the negative infinitesimal generator $L$ of a semigroup of positive contractions on $L^\infty$ has a bounded $H^\infty(S_\eta^0)$-calculus on the associated Poisson semigroup-BMO space for any angle $\eta>\pi/2$, provided the semigroup satisfies Bakry-Emry's $\Gamma_2 $ criterion. Our arguments only rely on the properties of the underlying semigroup and works well in the noncommutative setting. A key ingredient of our argument is a quasi monotone property for the subordinated semigroup $T_{t,\alpha}=e^{-tL^\alpha},0<\alpha<1$, that is proved in the first half of the article.


title: Revenue Forecasting for Enterprise Products

abstract: For any business, planning is a continuous process, and typically business-owners focus on making both long-term planning aligned with a particular strategy as well as short-term planning that accommodates the dynamic market situations. An ability to perform an accurate financial forecast is crucial for effective planning. In this paper, we focus on providing an intelligent and efficient solution that will help in forecasting revenue using machine learning algorithms. We experiment with three different revenue forecasting models, and here we provide detailed insights into the methodology and their relative performance measured on real finance data. As a real-world application of our models, we partner with Microsoft's Finance organization (department that reports Microsoft's finances) to provide them a guidance on the projected revenue for upcoming quarters.


title: Econophysics Macroeconomic Model

abstract: This paper presents macroeconomic model that is based on parallels between macroeconomic multi-agent systems and multi-particle systems. We use risk ratings of economic agents as their coordinates on economic space. Aggregates of economic or financial variables like Investment, Assets, Demand, Credits and etc. of economic agents near point x define corresponding macroeconomic variables as functions of time t and coordinates x on economic space. Parallels between multi-agent and multi-particle systems on economic space allow describe transition from economic kinetic-like to economic hydrodynamic-like approximation and derive macroeconomic hydrodynamic-like equations on economic space. Economic or financial transactions between economic agents determine evolution of macroeconomic variables This paper describes local macroeconomic approximation that takes into account transactions between economic agents with coordinates near same point x on economic space only and describes interaction between macroeconomic variables by linear differential operators. For simple model of interaction between macroeconomic variables as Demand on Investment and Interest Rate we derive hydrodynamic-like equations in a closed form. For perturbations of these macroeconomic variables we derive macroeconomic wave equations. Macroeconomic waves on economic space can propagate with exponential growth of amplitude and cause irregular time fluctuations of macroeconomic variables or induce economic crises.


title: A New Formulation of the $3D$ Compressible Euler Equations With Dynamic Entropy: Remarkable Null Structures and Regularity Properties

abstract: We derive a new formulation of the $3D$ compressible Euler equations with dynamic entropy exhibiting remarkable null structures and regularity properties. Our results hold for an arbitrary equation of state (which yields the pressure in terms of the density and the entropy) in non-vacuum regions where the speed of sound is positive. Our work is an extension of our prior joint work with J. Luk, in which we derived a similar new formulation in the special case of a barotropic fluid, that is, when the equation of state depends only on the density. The new formulation comprises covariant wave equations for the Cartesian components of the velocity and the logarithmic density coupled to a transport equation for the specific vorticity (defined to be vorticity divided by density), a transport equation for the entropy, and some additional transport-divergence-curl-type equations involving special combinations of the derivatives of the solution variables. The good geometric structures in the equations allow one to use the full power of the vectorfield method in treating the "wave part" of the system. In a forthcoming application, we will use the new formulation to give a sharp, constructive proof of finite-time shock formation, tied to the intersection of acoustic "wave" characteristics, for solutions with nontrivial vorticity and entropy at the singularity. In the present article, we derive the new formulation and overview the central role that it plays in the proof of shock formation.


title: Representation of light pressure resultant force and moment as a tensor series

abstract: In this article we address a problem of determination of light pressure upon space structures of the complex geometric shape. For the surface element, we wrote a condition that this element can interact with light only from the front side, not from the back side. This condition To in the form of series of Chebyshev polynomials of the first kind. Chebyshev expansion lets us move to the series of tensors of increasing rank for the problem of determination of force and moment. We obtained expressions for the fiber method for determination of light pressure on space structures of complex geometry taking into account self--shadowing and reflections within the structure. We also give the expressions for tensor parametrization by the specularity coefficient in case of specular-diffuse reflection. For these expressions we calculated the principal moment and force upon two-sided flat solar sail, spherical and cylindrical bodies, and approximated light pressure upon a perspective space observatory Millimetron. The proposed expressions can be used in the ballistic analysis of solar sails and other space objects, which are significantly affected by the radiation pressure. Also, these results can be used to analyze the dynamics of movement of large-scale space structures around the center of gravity under the light pressure.


title: Unique lifting of integer variables in minimal inequalities

abstract: This paper contributes to the theory of cutting planes for mixed integer linear programs (MILPs). Minimal valid inequalities are well understood for a relaxation of an MILP in tableau form where all the nonbasic variables are continuous; they are derived using the gauge function of maximal lattice-free convex sets. In this paper we study lifting functions for the nonbasic integer variables starting from such minimal valid inequalities. We characterize precisely when the lifted coefficient is equal to the coefficient of the corresponding continuous variable in every minimal lifting. The answer is a nonconvex region that can be obtained as a finite union of convex polyhedra. We then establish a necessary and sufficient condition for the uniqueness of the lifting function.


title: On exponential suppression of the cosmological constant in non-SUSY strings at two loops and beyond

abstract: Two independent criteria are presented that together guarantee exponential suppression of the two-loop cosmological constant in non-supersymmetric heterotic strings. They are derived by performing calculations in both the full string theory and in its underlying effective field theory, and come respectively from contributions that involve only physical untwisted states, and contributions that include orbifold twisted states. The criteria depend purely on the spectrum and charges, so a model that satisfies them will do so with no fine-tuning. An additional consistency condition (emerging from the so-called separating degeneration limit of the two-loop diagram) is that the one-loop cosmological constant must also be suppressed, by Bose-Fermi degeneracy in the massless spectrum. We remark that theories of this kind, that have continued exponential suppression to higher orders, can form the basis for a string implementation of the "naturalness without supersymmetry" idea.


title: L\'{e}vy Processes and Infinitely Divisible Measures in the Dual of a Nuclear Space

abstract: Let $\Phi$ be a nuclear space and let $\Phi'_{\beta}$ denote its strong dual. In this work we establish the one-to-one correspondence between infinitely divisible measures on $\Phi'_{\beta}$ and L\'{e}vy processes taking values in $\Phi'_{\beta}$. Moreover, we prove the L\'{e}vy-It\^{o} decomposition, the L\'{e}vy-Khintchine formula and the existence of c\`{a}dl\`{a}g versions for $\Phi'_{\beta}$-valued L\'{e}vy processes. A characterization for L\'{e}vy measures on $\Phi'_{\beta}$ is also established. Finally, we prove the L\'{e}vy-Khintchine formula for infinitely divisible measures on $\Phi'_{\beta}$.


title: Block-Diagonal Coding for Distributed Computing With Straggling Servers

abstract: We consider the distributed computing problem of multiplying a set of vectors with a matrix. For this scenario, Li et al. recently presented a unified coding framework and showed a fundamental tradeoff between computational delay and communication load. This coding framework is based on maximum distance separable (MDS) codes of code length proportional to the number of rows of the matrix, which can be very large. We propose a block-diagonal coding scheme consisting of partitioning the matrix into submatrices and encoding each submatrix using a shorter MDS code. We show that the assignment of coded matrix rows to servers to minimize the communication load can be formulated as an integer program with a nonlinear cost function, and propose an algorithm to solve it. We further prove that, up to a level of partitioning, the proposed scheme does not incur any loss in terms of computational delay (as defined by Li et al.) and communication load compared to the scheme by Li et al.. We also show numerically that, when the decoding time is also taken into account, the proposed scheme significantly lowers the overall computational delay with respect to the scheme by Li et al.. For heavy partitioning, this is achieved at the expense of a slight increase in communication load.


title: Shatter functions with polynomial growth rates

abstract: We study how a single value of the shatter function of a set system restricts its asymptotic growth. Along the way, we refute a conjecture of Bondy and Hajnal which generalizes Sauer's Lemma.


title: Nanoscale electrical measurements in liquids using AFM - progress and outlook

abstract: Fundamental mechanisms of energy storage, corrosion, sensing, and multiple biological functionalities are directly coupled to electrical processes and ionic dynamics at solid-liquid interfaces. In many cases, these processes are spatially inhomogeneous taking place at grain boundaries, step edges, point defects, ion channels, etc. and possess complex time and voltage dependent dynamics. This necessitates time-resolved and real-space probing of these phenomena. In this review, we discuss the applications of force-sensitive voltage modulated scanning probe microscopy (SPM) for probing electrical phenomena at solid-liquid interfaces. We first describe the working principles behind electrostatic and Kelvin Probe Force microscopies (EFM & KPFM) at the gas-solid interface, review the state of the art in advanced KPFM methods and developments to (i) overcome limitations of classical KPFM, (i) expand the information accessible from KPFM, and (iii) extend KPFM operation to liquid environments. We briefly discuss the theoretical framework of the electrical double layer (EDL) forces and dynamics, the implications and breakdown of classical EDL models for highly charged interfaces, or under high ion concentrations, and briefly describe recent modifications of the classical EDL theory relevant for understanding nanoscale electrical measurements at the solid-liquid interface. We further review the latest achievements in mapping surface charge, dielectric constants, and electrodynamic and electrochemical processes in liquids. Finally, we outline the key challenges and opportunities that exist in the field of nanoscale electrical measurements in liquid as well as provide a roadmap for the future development of liquid KPFM.


title: On the migration-induced formation of the 9:7 mean motion resonance

abstract: We study formation of 9:7 mean motion resonance (MMR) as a result of convergent migration of two planets embedded in a disc. Depending on the migration parameters, initial orbits and planets' masses the system may pass through the resonance or enter it (permanently or temporarily). We illustrate that a stable equilibrium of the averaged system (a periodic orbit of the N-body model) is surrounded by the region of permanent resonance capture, whose size depends on the migration parameters and the planets mass ratio. A system located inside this region tends towards the equilibrium (the capture is permanent), while a system located outside the region evolves away from the equilibrium and leaves the resonance. We verify recent results of Delisle et al. and Xu & Lai where they show that for $m_1 \lesssim m_2$ ($m_1, m_2$ are the inner and outer planets' masses) the equilibrium is unstable when the migration is added, so the capture cannot be permanent. We show that for particular migration parameters the situation may be reversed (the equilibria are unstable for $m_1 \gtrsim m_2$). We illustrate that 9:7 MMR consists of two modes separated with a separatrix. The inner one is centred at the equilibrium and the outer one has no equilibrium in its centre. A system located outside the region of stable capture evolves from the inner into the outer mode. The evolution occurs along families of periodic orbits of the averaged system, that play a crucial role in the dynamics after the resonance capture.


title: Space-Time Graph Modeling of Ride Requests Based on Real-World Data

abstract: This paper focuses on modeling ride requests and their variations over location and time, based on analyzing extensive real-world data from a ride-sharing service. We introduce a graph model that captures the spatial and temporal variability of ride requests and the potentials for ride pooling. We discover these ride request graphs exhibit a well known property called densification power law often found in real graphs modelling human behaviors. We show the pattern of ride requests and the potential of ride pooling for a city can be characterized by the densification factor of the ride request graphs. Previous works have shown that it is possible to automatically generate synthetic versions of these graphs that exhibit a given densification factor. We present an algorithm for automatic generation of synthetic ride request graphs that match quite well the densification factor of ride request graphs from actual ride request data.


title: Primordial perturbations generated by Higgs field and $R^2$ operator

abstract: If the very early Universe is dominated by the non-minimally coupled Higgs field and Starobinsky's curvature-squared term together, the potential diagram would mimic the landscape of a valley, serving as a cosmological attractor. The inflationary dynamics along this valley is studied, model parameters are constrained against observational data, and the isocurvature perturbation is evaluated.


title: Some Problems in Harmonic Analysis

abstract: In May 2016, we organized a conference in harmonic analysis in honor of Professor Michael Christ, on the campus of the University of Wisconsin in Madison. We are happy to present sixteen open problems, almost all of which were contributed by participants of a problem session held in the afternoon of May 19, 2016.


title: Stability of casein micelles cross-linked with genipin: a physicochemical study as a function of pH

abstract: Chemical or enzymatic cross-linking of casein micelles (CMs) increases their stability against dissociating agents. In this paper, a comparative study of stability between native CMs and CMs cross-linked with genipin (CMs-GP) as a function of pH is described. Stability to temperature and ethanol were investigated in the pH range 2.0-7.0. The size and the charge ($\zeta$-potential) of the particles were determined by dynamic light scattering. Native CMs precipitated below pH 5.5, CMs-GP precipitated from pH 3.5 to 4.5, whereas no precipitation was observed at pH 2.0-3.0 or pH 4.5-7.0. The isoelectric point of CMs-GP was determined to be pH 3.7. Highest stability against heat and ethanol was observed for CMs-GP at pH 2, where visible coagulation was determined only after 800 s at 140 $^\circ$C or 87.5% (v/v) of ethanol. These results confirmed the hypothesis that cross-linking by GP increased the stability of CMs.


title: On the complexity of generalized chromatic polynomials

abstract: J. Makowsky and B. Zilber (2004) showed that many variations of graph colorings, called CP-colorings in the sequel, give rise to graph polynomials. This is true in particular for harmonious colorings, convex colorings, mcc_t-colorings, and rainbow colorings, and many more. N. Linial (1986) showed that the chromatic polynomial $\chi(G;X)$ is #P-hard to evaluate for all but three values X=0,1,2, where evaluation is in P. This dichotomy includes evaluation at real or complex values, and has the further property that the set of points for which evaluation is in P is finite. We investigate how the complexity of evaluating univariate graph polynomials that arise from CP-colorings varies for different evaluation points. We show that for some CP-colorings (harmonious, convex) the complexity of evaluation follows a similar pattern to the chromatic polynomial. However, in other cases (proper edge colorings, mcc_t-colorings, H-free colorings) we could only obtain a dichotomy for evaluations at non-negative integer points. We also discuss some CP-colorings where we only have very partial results.


title: On one fractal property of the Minkowski function

abstract: It is shown in the present article that the Minkowski function does not preserve the Hausdorff-Besicovitch dimension. Results of this article were presented by the author of this article on the Second Interuniversity Scientific Conference on Mathematics and Physics for Young Scientists in April, 2011 (Links to the conference paper (in Ukrainian): www.researchgate.net/publication/301637057, www.researchgate.net/publication/311665828. A link to the presentation (in Ukrainian): www.researchgate.net/publication/303752636).


title: Perceptually Optimized Image Rendering

abstract: We develop a framework for rendering photographic images, taking into account display limitations, so as to optimize perceptual similarity between the rendered image and the original scene. We formulate this as a constrained optimization problem, in which we minimize a measure of perceptual dissimilarity, the Normalized Laplacian Pyramid Distance (NLPD), which mimics the early stage transformations of the human visual system. When rendering images acquired with higher dynamic range than that of the display, we find that the optimized solution boosts the contrast of low-contrast features without introducing significant artifacts, yielding results of comparable visual quality to current state-of-the art methods with no manual intervention or parameter settings. We also examine a variety of other display constraints, including limitations on minimum luminance (black point), mean luminance (as a proxy for energy consumption), and quantized luminance levels (halftoning). Finally, we show that the method may be used to enhance details and contrast of images degraded by optical scattering (e.g. fog).


title: Outliers and related problems

abstract: We define outliers as a set of observations which contradicts the proposed mathematical (statistical) model and we discuss the frequently observed types of the outliers. Further we explore what changes in the model have to be made in order to avoid the occurance of the outliers. We observe that some variants of the outliers lead to classical results in probability, such as the law of large numbers and the concept of heavy tailed distributions. Key words: outlier; the law of large numbers; heavy tailed distributions; model rejection.


title: Residual and Plain Convolutional Neural Networks for 3D Brain MRI Classification

abstract: In the recent years there have been a number of studies that applied deep learning algorithms to neuroimaging data. Pipelines used in those studies mostly require multiple processing steps for feature extraction, although modern advancements in deep learning for image classification can provide a powerful framework for automatic feature generation and more straightforward analysis. In this paper, we show how similar performance can be achieved skipping these feature extraction steps with the residual and plain 3D convolutional neural network architectures. We demonstrate the performance of the proposed approach for classification of Alzheimer's disease versus mild cognitive impairment and normal controls on the Alzheimer's Disease National Initiative (ADNI) dataset of 3D structural MRI brain scans.


title: Polynomial-time Algorithms for Computing Distances of Fuzzy Transition Systems

abstract: Behaviour distances to measure the resemblance of two states in a (nondeterministic) fuzzy transition system have been proposed recently in the literature. Such a distance, defined as a pseudo-ultrametric over the state space of the model, provides a quantitative analogue of bisimilarity. In this paper, we focus on the problem of computing these distances. We first extend the definition of the pseudo-ultrametric by introducing discount such that the discounting factor being equal to 1 captures the original definition. We then provide polynomial-time algorithms to calculate the behavioural distances, in both the non-discounted and the discounted setting. The algorithm is strongly polynomial in the former case. Furthermore, we give a polynomial-time algorithm to compute bisimulation over fuzzy transition systems which captures the distance being equal to 0.


title: Semi-Compact Skyrmion-like Structures

abstract: We study three distinct types of planar, spherically symmetric and localized structures, one of them having non-topological behavior and the two others being of topological nature. The non-topological structures have energy density localized in a compact region in the plane, but are unstable against spherically symmetric fluctuations. The topological structures are stable and behave as vortices and skyrmions at larger distances, but they engender interesting compact behavior as one approaches their inner cores. They are semi-compact skyrmion-like spin textures generated from models that allow to control the internal behavior of such topological structures.


title: Synergistic effects in threshold models on networks

abstract: Network structure can have significant effects on the propagation of diseases, memes, and information on social networks. Such effects depend on the specific type of dynamical process that affects the nodes and edges of a network, and it is important to develop tractable models of spreading processes on networks to explore how network structure affects dynamics. In this paper, we incorporate the idea of \emph{synergy} into a two-state ("active" or "passive") threshold model of social influence on networks. Our model's update rule is deterministic, and the influence of each meme-carrying (i.e., active) neighbor can --- depending on a parameter --- either be enhanced or inhibited by an amount that depends on the number of active neighbors of a node. Such a synergistic system models social behavior in which the willingness to adopt either accelerates or saturates depending on the number of neighbors who have adopted that behavior. We illustrate that the synergy parameter in our model has a crucial effect on system dynamics, as it determines whether degree-$k$ nodes are possible or impossible to activate. We simulate synergistic meme spreading on both random-graph models and networks constructed from empirical data. Using a local-tree approximation, we examine the spreading of synergistic memes and find good agreement on all but one of the networks on which we simulate spreading. We find for any network and for a broad family of synergistic models that one can predict which synergy-parameter values allow degree-$k$ nodes to be activated.


title: Geometric Structure of Pseudo-plane Quadratic Flows

abstract: Quadratic flows have the unique property of uniform strain and are commonly used in turbulence modeling and hydrodynamic analysis. While previous application focused on two-dimensional homogeneous fluid, this study examines the geometric structure of three-dimensional quadratic flows in stratified fluid by solving a steady-state pseudo-plane flow model. The complete set of exact solutions reveals that steady quadratic flows have invariant conic type in non-rotating frame and non-rotatory vertical structure in rotating frame. Three baroclinic solutions with vertically non-aligned structure disprove an earlier conjecture. The rich topology of quadratic flows stands in contrast to the depleted geometry of high-degree polynomial flows. A paradox in the steady solutions of shallow-water reduced-gravity models is also explained.


title: Weight recursions for any rotation symmetric Boolean functions

abstract: Let $f_n(x_1, x_2, \ldots, x_n)$ denote the algebraic normal form (polynomial form) of a rotation symmetric Boolean function of degree $d$ in $n \geq d$ variables and let $wt(f_n)$ denote the Hamming weight of this function. Let $(1, a_2, \ldots, a_d)_n$ denote the function $f_n$ of degree $d$ in $n$ variables generated by the monomial $x_1x_{a_2} \cdots x_{a_d}.$ Such a function $f_n$ is called {\em monomial rotation symmetric} (MRS). It was proved in a $2012$ paper that for any MRS $f_n$ with $d=3,$ the sequence of weights $\{w_k = wt(f_k):~k = 3, 4, \ldots\}$ satisfies a homogeneous linear recursion with integer coefficients. In this paper it is proved that such recursions exist for any rotation symmetric function $f_n;$ such a function is generated by some sum of $t$ monomials of various degrees. The last section of the paper gives a Mathematica program which explicitly computes the homogeneous linear recursion for the weights, given any rotation symmetric $f_n.$ The reader who is only interested in finding some recursions can use the program and not be concerned with the details of the rather complicated proofs in this paper.


title: Constant Size Molecular Descriptors For Use With Machine Learning

abstract: A set of molecular descriptors whose length is independent of molecular size is developed for machine learning models that target thermodynamic and electronic properties of molecules. These features are evaluated by monitoring performance of kernel ridge regression models on well-studied data sets of small organic molecules. The features include connectivity counts, which require only the bonding pattern of the molecule, and encoded distances, which summarize distances between both bonded and non-bonded atoms and so require the full molecular geometry. In addition to having constant size, these features summarize information regarding the local environment of atoms and bonds, such that models can take advantage of similarities resulting from the presence of similar chemical fragments across molecules. Combining these two types of features leads to models whose performance is comparable to or better than the current state of the art. The features introduced here have the advantage of leading to models that may be trained on smaller molecules and then used successfully on larger molecules.


title: Electrically driving nuclear spin qubits with microwave photonic bangap resonators

abstract: The electronic and nuclear spin degrees of freedom for donor impurities in semiconductors form ultra coherent two-level systems that are useful for quantum information applications. Spins naturally have magnetic dipoles, so alternating current (AC) magnetic fields are frequently used to drive spin transitions and perform quantum gates. These fields can be difficult to spatially confine to single donor qubits so alternative methods of control such as AC electric field driven spin resonance are desirable. However, donor spin qubits do not have electric dipole moments so that they can not normally be driven by electric fields. In this work we challenge that notion by demonstrating a new, all-electric-field method for controlling neutral $^{31}$P and $^{75}$As donor nuclear spins in silicon through modulation of their donor-bound electrons. This method has major advantages over magnetic field control since electric fields are easy to confine at the nanoscale. This leads to lower power requirements, higher qubit densities, and faster gate times. We also show that this form of control allows for driving nuclear spin qubits at either their resonance frequency or the first subharmonic of that frequency, thus reducing device bandwidth requirements. Interestingly, as we relax the bandwidth requirements, we demonstrate that the computational Hilbert space is expanded to include double quantum transitions, making it feasible to use all four nuclear spin states to implement nuclear-spin-based qudits in Si:As. Based on these results, one can envision novel high-density, low-power quantum computing architectures using nuclear spins in silicon.


title: Moments of zeta and correlations of divisor-sums: V

abstract: In this series of papers we examine the calculation of the $2k$th moment and shifted moments of the Riemann zeta-function on the critical line using long Dirichlet polynomials and divisor correlations. The present paper completes the general study of what we call Type II sums which utilize a circle method framework and a convolution of shifted convolution sums to obtain all of the lower order terms in the asymptotic formula for the mean square along $[T,2T]$ of a Dirichlet polynomial of arbitrary length with divisor functions as coefficients.


title: Convex Parameterizations and Fidelity Bounds for Nonlinear Identification and Reduced-Order Modelling

abstract: Model instability and poor prediction of long-term behavior are common problems when modeling dynamical systems using nonlinear "black-box" techniques. Direct optimization of the long-term predictions, often called simulation error minimization, leads to optimization problems that are generally non-convex in the model parameters and suffer from multiple local minima. In this work we present methods which address these problems through convex optimization, based on Lagrangian relaxation, dissipation inequalities, contraction theory, and semidefinite programming. We demonstrate the proposed methods with a model order reduction task for electronic circuit design and the identification of a pneumatic actuator from experiment.


title: Candidate X-ray-Emitting OB Stars in the MYStIX Massive Star-Forming Regions

abstract: Massive, O and early B-type (OB) stars remain incompletely catalogued in the nearby Galaxy due to high extinction, bright visible and infrared nebular emission in H II regions, and high field star contamination. These difficulties are alleviated by restricting the search to stars with X-ray emission. Using the X-ray point sources from the Massive Young star-forming complex Study in Infrared and X-rays (MYStIX) survey of OB-dominated regions, we identify 98 MYStIX candidate OB (MOBc) stars by fitting their 1-8 micron spectral energy distributions (SEDs) with reddened stellar atmosphere models. We identify 27 additional MOBc stars based on JHK photometry of X-ray stars lacking SED fitting. These candidate OB stars indicate that the current census of stars earlier than B1, taken across the 18 MYStIX regions studied, is less than 50% complete. We also fit the SEDs of 239 previously-published OB stars to measure interstellar extinction and bolometric luminosities, revealing six candidate massive binary systems and five candidate O-type (super)giants. As expected, candidate OB stars have systematically higher extinction than previously-published OB stars. Notable results for individual regions include: identification of the OB population of a recently discovered massive cluster in NGC 6357; an older OB association in the M17 complex; and new massive luminous O stars near the Trifid Nebula. In several relatively poorly-studied regions (RCW 38, NGC 6334, NGC 6357, Trifid, and NGC 3576), the OB populations may increase by factors of >2.


title: Elastic Strain Fields In Lateral Double Ge/Si Quantum Dots

abstract: Simulations of the elastic strain fields for double Ge/Si quantum dots located on the same wetting layer are carried out. The cluster approximation is used for the atomistic model based on the Keating potential. The spatial distributions of the strain energy density and electron potential energy are calculated using clusters containing atoms of 150 coordination spheres. It is shown that the used cluster boundary conditions are close to the periodic boundary conditions.


title: Patchwork Kriging for Large-scale Gaussian Process Regression

abstract: This paper presents a new approach for Gaussian process (GP) regression for large datasets. The approach involves partitioning the regression input domain into multiple local regions with a different local GP model fitted in each region. Unlike existing local partitioned GP approaches, we introduce a technique for patching together the local GP models nearly seamlessly to ensure that the local GP models for two neighboring regions produce nearly the same response prediction and prediction error variance on the boundary between the two regions. This largely mitigates the well-known discontinuity problem that degrades the boundary accuracy of existing local partitioned GP methods. Our main innovation is to represent the continuity conditions as additional pseudo-observations that the differences between neighboring GP responses are identically zero at an appropriately chosen set of boundary input locations. To predict the response at any input location, we simply augment the actual response observations with the pseudo-observations and apply standard GP prediction methods to the augmented data. In contrast to heuristic continuity adjustments, this has an advantage of working within a formal GP framework, so that the GP-based predictive uncertainty quantification remains valid. Our approach also inherits a sparse block-like structure for the sample covariance matrix, which results in computationally efficient closed-form expressions for the predictive mean and variance. In addition, we provide a new spatial partitioning scheme based on a recursive space partitioning along local principal component directions, which makes the proposed approach applicable for regression domains having more than two dimensions. Using three spatial datasets and three higher dimensional datasets, we investigate the numerical performance of the approach and compare it to several state-of-the-art approaches.


title: A multiphase Cahn-Hilliard-Darcy model for tumour growth with necrosis

abstract: We derive a Cahn-Hilliard-Darcy model to describe multiphase tumour growth taking interactions with multiple chemical species into account as well as the simultaneous occurrence of proliferating, quiescent and necrotic regions. Via a coupling of the Cahn-Hilliard-Darcy equations to a system of reaction-diffusion equations a multitude of phenomena such as nutrient diffusion and consumption, angiogenesis, hypoxia, blood vessel growth, and inhibition by toxic agents, which are released for example by the necrotic cells, can be included. A new feature of the modelling approach is that a volume-averaged velocity is used, which dramatically simplifies the resulting equations. With the help of formally matched asymptotic analysis we develop new sharp interface models. Finite element numerical computations are performed and in particular the effects of necrosis on tumour growth is investigated numerically.


title: Global $\Lambda$ hyperon polarization in nuclear collisions: evidence for the most vortical fluid

abstract: The extreme temperatures and energy densities generated by ultra-relativistic collisions between heavy nuclei produce a state of matter with surprising fluid properties. Non-central collisions have angular momentum on the order of 1000$\hbar$, and the resulting fluid may have a strong vortical structure that must be understood to properly describe the fluid. It is also of particular interest because the restoration of fundamental symmetries of quantum chromodynamics is expected to produce novel physical effects in the presence of strong vorticity. However, no experimental indications of fluid vorticity in heavy ion collisions have so far been found. Here we present the first measurement of an alignment between the angular momentum of a non-central collision and the spin of emitted particles, revealing that the fluid produced in heavy ion collisions is by far the most vortical system ever observed. We find that $\Lambda$ and $\overline{\Lambda}$ hyperons show a positive polarization of the order of a few percent, consistent with some hydrodynamic predictions. A previous measurement that reported a null result at higher collision energies is seen to be consistent with the trend of our new observations, though with larger statistical uncertainties. These data provide the first experimental access to the vortical structure of the "perfect fluid" created in a heavy ion collision. They should prove valuable in the development of hydrodynamic models that quantitatively connect observations to the theory of the Strong Force. Our results extend the recent discovery of hydrodynamic spin alignment to the subatomic realm.


title: A highly efficient 3D level-set grain growth algorithm tailored for ccNUMA architecture

abstract: A highly efficient simulation model for 2D and 3D grain growth and recrystallization was developed based on the level-set method. The model introduces modern computational concepts to achieve excellent performance on parallel computer architectures. Strong scalability was measured on ccNUMA architectures. To achieve this, the proposed approach considers the application of local level-set functions at the grain level. Ideal and non-ideal grain growth was simulated in 3D with the objective to study the evolution of statistical representative volume elements in polycrystals. In addition, microstructure evolution in an anisotropic magnetic material affected by an external magnetic field was simulated.


title: DSSD : Deconvolutional Single Shot Detector

abstract: The main contribution of this paper is an approach for introducing additional context into state-of-the-art general object detection. To achieve this we first combine a state-of-the-art classifier (Residual-101[14]) with a fast detection framework (SSD[18]). We then augment SSD+Residual-101 with deconvolution layers to introduce additional large-scale context in object detection and improve accuracy, especially for small objects, calling our resulting system DSSD for deconvolutional single shot detector. While these two contributions are easily described at a high-level, a naive implementation does not succeed. Instead we show that carefully adding additional stages of learned transformations, specifically a module for feed-forward connections in deconvolution and a new output module, enables this new approach and forms a potential way forward for further detection research. Results are shown on both PASCAL VOC and COCO detection. Our DSSD with $513 \times 513$ input achieves 81.5% mAP on VOC2007 test, 80.0% mAP on VOC2012 test, and 33.2% mAP on COCO, outperforming a state-of-the-art method R-FCN[3] on each dataset.


title: Yangian and SUSY symmetry of High Spin Parton Splitting Amplitudes in Generalised Yang-Mills Theory

abstract: We have calculated the high spin parton splitting amplitudes postulating the Yangian symmetry of the scattering amplitudes for tensorgluons. The resulting splitting amplitudes coincide with the earlier calculations, which were based on the BCFW recursion relations. The resulting formula unifies all known splitting probabilities found earlier in gauge field theories. It describes splitting probabilities for integer and half-integer spin particles. We also checked that the splitting probabilities fulfil generalised Kounnas-Ross N=1 supersymmetry relations hinting to the fact that the underlying theory can be formulated in an explicit supersymmetric manner.


title: Mean Field Stochastic Games with Binary Action Spaces and Monotone Costs

abstract: This paper considers mean field games in a multi-agent Markov decision process (MDP) framework. Each player has a continuum state and binary action. By active control, a player can bring its state to a resetting point. All players are coupled through their cost functions. The structural property of the individual strategies is characterized in terms of threshold policies when the mean field game admits a solution. We further introduce a stationary equation system of the mean field game and analyze uniqueness of its solution under positive externalities.


title: Duality between k-essence and Rastall gravity

abstract: The k-essence theory with a power-law function of $(\partial\phi)^2$ and Rastall's non-conservative theory of gravity with a scalar field are shown to have the same solutions for the metric under the assumption that both the metric and the scalar fields depend on a single coordinate. This equivalence (called k-R duality) holds for static configurations with various symmetries (spherical, plane, cylindrical, etc.) and all homogeneous cosmologies. In the presence of matter, Rastall's theory requires additional assumptions on how the stress-energy tensor non-conservation is distributed between different contributions. Two versions of such non-conservation are considered in the case of isotropic spatially flat cosmological models with a perfect fluid: one (R1) in which there is no coupling between the scalar field and the fluid, and another (R2) in which the fluid separately obeys the usual conservation law. In version R1 it is shown that k-R duality holds not only for the cosmological models themselves but also for their adiabatic perturbations. In version R2, among other results, a particular model is singled out that reproduces the same cosmological expansion history as the standard $\Lambda$CDM model but predicts different behaviors of small fluctuations in the k-essence and Rastall frameworks.


title: The $L^2$-cutoffs for reversible Markov chains

abstract: In this article, we considers reversible Markov chains of which $L^2$-distances can be expressed in terms of Laplace transforms. The cutoff of Laplace transforms was first discussed by Chen and Saloff-Coste in [8], while we provide here a completely different pathway to analyze the $L^2$-distance. Consequently, we obtain several considerably simplified criteria and this allows us to proceed advanced theoretical studies, including the comparison of cutoffs between discrete time lazy chains and continuous time chains. For an illustration, we consider product chains, a rather complicated model which could be involved to analyze using the method in [8], and derive the equivalence of their $L^2$-cutoffs.


title: Repair Duality with Locally Repairable and Locally Regenerating Codes

abstract: We construct an explicit family of locally repairable and locally regenerating codes whose existence was proven in a recent work by Kamath et al. about codes with local regeneration but no explicit construction was given. This explicit family of codes is based on HashTag codes. HashTag codes are recently defined vector codes with different vector length $\alpha$ (also called a sub-packetization level) that achieve the optimal repair bandwidth of MSR codes or near-optimal repair bandwidth depending on the sub-packetization level. We applied the technique of parity-splitting code construction. We show that the lower bound on the size of the finite field for the presented explicit code constructions can be lower than the one given in the work of Kamath et al. Finally, we discuss the importance of having two ways for node repair with locally regenerating HashTag codes: repair only with local parity nodes or repair with both local and global parity nodes. To the best of the authors' knowledge, this is the first work where this duality in repair process is discussed. We give a practical example and experimental results in Hadoop where we show the benefits of having this repair duality.


title: Cutoffs for product chains

abstract: In this article, we consider products of ergodic Markov chains and discuss their cutoffs in the total variation. Through a new inequality relating the total variation and the Hellinger distance, we may identify the total variation cutoffs with cutoffs in the Hellinger distance. This provides a new scheme to study the total variation mixing of Markov chains, in particular, product chains. In the theoretical framework, a series of criteria are introduced to examine cutoffs and a comparison of mixing between the product chain and its coordinate chains is made in detail. For illustration, we consider products of two-state chains, cycles and other typical examples.


title: Diffusion of Oxygen Isotopes in Thermally Evolving Planetesimals and Size Ranges of Presolar Silicate Grains

abstract: Presolar grains are small particles found in meteorites through their isotopic compositions which are considerably different from those of materials in the Solar System. If some isotopes in presolar grains diffused out beyond their grain sizes when they were embedded in parent bodies of meteorites, their isotopic compositions could be washed out, and hence the grains cannot be identified as presolar grains any more. We explore this possibility for the first time by self-consistently simulating the thermal evolution of planetesimals and the diffusion length of $^{18}$O in presolar silicate grains. Our results show that presolar silicate grains smaller than $\sim$ 0.03 $\mu m$ cannot keep their original isotopic compositions even if the host planetesimals experienced maximum temperature as low as 600 $^{\circ}$C. Since this temperature corresponds to the one experienced by petrologic type 3 chondrites, the isotopic diffusion can constrain the size of presolar silicate grains discovered in such chondrites to be larger than $\sim$ 0.03 $\mu m$. We also find that the diffusion lengths of $^{18}$O reach $\sim$ 0.3-2 $\mu m$ in planetesimals that were heated up to 700-800 $^{\circ}$C. This indicates that, if the original size of presolar grains spans a range from $\sim$ 0.001 $\mu m$ to $\sim$ 0.3 $\mu m$ like that in the interstellar medium, the isotopic records of the presolar grains may be almost completely lost in such highly thermalized parent bodies. We propose that isotopic diffusion could be a key process to control the size distribution and abundance of presolar grains in some types of chondrites.


title: Three-dimensionality of the bulk electronic structure in WTe2

abstract: We use temperature- and field-dependent resistivity measurements [Shubnikov--de Haas (SdH) quantum oscillations] and ultrahigh resolution, tunable, vacuum ultraviolet (VUV) laser-based angle-resolved photoemission spectroscopy (ARPES) to study the three-dimensionality (3D) of the bulk electronic structure in WTe2, a type-II Weyl semimetal. The bulk Fermi surface (FS) consists of two pairs of electron pockets and two pairs of hole pockets along the X-Gamma-X direction as detected by using an incident photon energy of 6.7 eV, which is consistent with the previously reported data. However, if using an incident photon energy of 6.36 eV, another pair of tiny electron pockets is detected on both sides of the Gamma point, which is in agreement with the small quantum oscillation frequency peak observed in the magnetoresistance. Therefore, the bulk, 3D FS consists of three pairs of electron pockets and two pairs of hole pockets in total. With the ability of fine tuning the incident photon energy, we demonstrate the strong three-dimensionality of the bulk electronic structure in WTe2. The combination of resistivity and ARPES measurements reveal the complete, and consistent, picture of the bulk electronic structure of this material.


title: Linear Instability of the Reissner-Nordstr\"om Cauchy Horizon

abstract: This work studies solutions of the scalar wave equation \[\Box_g\phi=0\] on a fixed subextremal Reissner-Nordstr\"{o}m spacetime with non-vanishing charge $q$ and mass $M$. In a recent paper, Luk and Oh established that generic smooth and compactly supported initial data on a Cauchy hypersurface lead to solutions which are singular in the $W^{1,2}_{loc}$ sense near the Cauchy horizon in the black hole interior, and it follows easily that they are also singular in the $W^{1,p}_{loc}$ sense for $p>2$. On the other hand, the work of Franzen shows that such solutions are non-singular near the Cauchy horizon in the $W^{1,1}_{loc}$ sense. Motivated by these results, we investigate the strength of the singularity at the Cauchy horizon. We identify a sufficient condition on the black hole interior (which holds for the full subextremal parameter range $0<|q|<M$) ensuring $W^{1,p}_{loc}$ blow up near the Cauchy horizon of solutions arising from generic smooth and compactly supported data for every $1<p<2$. We moreover prove that provided the spacetime parameters satisfy $\frac{2\sqrt e}{e+1}<\frac{|q|}{M}<1$, we in fact have $W^{1,p}_{loc}$ blow up near the Cauchy horizon for such solutions for every $1<p<2$. This shows that the singularity is even stronger than was implied by the work of Luk and Oh for this restricted parameter range.


title: Algorithmic complexity of proper labeling problems

abstract: A proper labeling of a graph is an assignment of integers to some elements of a graph, which may be the vertices, the edges, or both of them, such that we obtain a proper vertex coloring via the labeling subject to some conditions. The problem of proper labeling offers many variants and received a great interest during recent years. We consider the algorithmic complexity of some variants of the proper labeling problems, we present some polynomial time algorithms and $ \mathbf{NP} $-completeness results for them.


title: Arbitrary order 2D virtual elements for polygonal meshes: Part I, elastic problem

abstract: The present work deals with the formulation of a Virtual Element Method (VEM) for two dimensional structural problems. The contribution is split in two parts: in part I, the elastic problem is discussed, while in part II [3] the method is extended to material nonlinearity, considering different inelastic responses of the material. In particular, in part I a standardized procedure for the construction of all the terms required for the implementation of the method in a code is explained. The procedure is initially illustrated for the simplest case of quadrilateral virtual elements with linear approximation of displacement variables on the boundary of the element. Then, the case of polygonal elements with quadratic and, even, higher order interpolation is considered. The construction of the method is detailed, deriving the approximation of the consistent term, the required stabilization term and the loading term for all the considered virtual elements. A wide numerical investigation is performed to assess the performances of the developed virtual elements, considering different number of edges describing the elements and different order of approximations of the unknown field. Numerical results are also compared with the one recovered using the classical finite element method.


title: Quantum phase transitions in the driven dissipative Jaynes-Cummings oscillator: from the dispersive regime to resonance

abstract: We follow the passage from complex amplitude bistability to phase bistability in the driven dissipative Jaynes-Cummings oscillator. Quasidistribution functions in the steady state are employed, for varying qubit-cavity detuning and drive parameters, in order to track a first-order dissipative quantum phase transition up to the critical point marking a second-order transition and spontaneous symmetry breaking. We demonstrate the photon blockade breakdown in the dispersive regime, and find that the coexistence of cavity states in the regime of quantum bistability is accompanied by pronounced qubit-cavity entanglement. Focusing on the role of quantum fluctuations in the response of both coupled quantum degrees of freedom (cavity and qubit), we move from a region of minimal entanglement in the dispersive regime, where we derive analytical perturbative results, to the threshold behaviour of spontaneous dressed-state polarization at resonance.


title: Additive cyclic codes over finite commutative chain rings

abstract: Additive cyclic codes over Galois rings were investigated in previous works. In this paper, we investigate the same problem but over a more general ring family, finite commutative chain rings. When we focus on non-Galois finite commutative chain rings, we observe two different kinds of additivity. One of them is a natural generalization of the previous studies, whereas the other one has some unusual properties especially while constructing dual codes. We interpret the reasons of such properties and illustrate our results giving concrete examples.


title: Decentralized Coded Caching in Wireless Networks: Trade-off between Storage and Latency

abstract: This paper studies the decentralized coded caching for a Fog Radio Access Network (F-RAN), whereby two edge-nodes (ENs) connected to a cloud server via fronthaul links with limited capacity are serving the requests of $K_r$ users. We consider all ENs and users are equipped with caches. A decentralized content placement is proposed to independently store contents at each network node during the off-peak hours. After that, we design a coded delivery scheme in order to deliver the user demands during the peak-hours under the objective of minimizing the normalized delivery time (NDT), which refers to the worst case delivery latency. An information-theoretic lower bound on the minimum NDT is derived for arbitrary number of ENs and users. We evaluate numerically the performance of the decentralized scheme. Additionally, we prove the approximate optimality of the decentralized scheme for a special case when the caches are only available at the ENs.


title: Frustrated Spin One On A Diamond Lattice

abstract: We report the discovery of a spin one diamond lattice in NiRh2O4. This spinel undergoes a cubic to tetragonal phase transition at T = 440 K that leaves all nearest neighbor interactions equivalent. In the tetragonal phase, magnetization measurements show a Ni2+ effective moment of p(eff) = 3.3(1) and dominant antiferromagnetic interactions with {\Theta}(CW) = -11.3(7) K. No phase transition to a long-range magnetically ordered state is observed by specific heat measurements down to T = 0.1 K. Inelastic neutron scattering measurements on sub-stoichiometric NiRh2O4 reveal possible valence-bond behavior and show no visible signs of magnetic ordering. NiRh2O4 provides a platform on which to explore the previously unknown and potentially rich physics of spin one interacting on the diamond lattice, including the realization of theoretically predicted quantum spin liquid and topological paramagnet states.


title: Dynamic Mortality Risk Predictions in Pediatric Critical Care Using Recurrent Neural Networks

abstract: Viewing the trajectory of a patient as a dynamical system, a recurrent neural network was developed to learn the course of patient encounters in the Pediatric Intensive Care Unit (PICU) of a major tertiary care center. Data extracted from Electronic Medical Records (EMR) of about 12000 patients who were admitted to the PICU over a period of more than 10 years were leveraged. The RNN model ingests a sequence of measurements which include physiologic observations, laboratory results, administered drugs and interventions, and generates temporally dynamic predictions for in-ICU mortality at user-specified times. The RNN's ICU mortality predictions offer significant improvements over those from two clinically-used scores and static machine learning algorithms.


title: Arbitrary order 2D virtual elements for polygonal meshes: Part II, inelastic problem

abstract: The present paper is the second part of a twofold work, whose first part is reported in [3], concerning a newly developed Virtual Element Method (VEM) for 2D continuum problems. The first part of the work proposed a study for linear elastic problem. The aim of this part is to explore the features of the VEM formulation when material nonlinearity is considered, showing that the accuracy and easiness of implementation discovered in the analysis inherent to the first part of the work are still retained. Three different nonlinear constitutive laws are considered in the VEM formulation. In particular, the generalized viscoplastic model, the classical Mises plasticity with isotropic/kinematic hardening and a shape memory alloy (SMA) constitutive law are implemented. The versatility with respect to all the considered nonlinear material constitutive laws is demonstrated through several numerical examples, also remarking that the proposed 2D VEM formulation can be straightforwardly implemented as in a standard nonlinear structural finite element method (FEM) framework.


title: Twisted monodromy homomorphisms and Massey products

abstract: Let $\phi: M\to M$ be a diffeomorphism of a $C^\infty$ compact connected manifold, and $X$ its mapping torus. There is a natural fibration $p:X\to S^1$, denote by $\xi\in H^1(X, \mathbb{Z})$ the corresponding cohomology class. Let $\rho:\pi_1(X)\to GL(n,\mathbb{C})$ be a representation, denote by $H^*(X,\rho)$ the corresponding twisted cohomology of $X$. Denote by $\rho_0$ the restriction of $\rho$ to $\pi_1(M)$, and by $\rho^*_0$ the antirepresentation conjugate to $\rho_0$. We construct from these data an automorphism of the group $H_*(M,\rho^*_0)$, that we call the twisted monodromy homomorphism $\phi_*$. The aim of the present work is to establish a relation between Massey products in $H^*(X,\rho)$ and Jordan blocks of $\phi_*$. Given a non-zero complex number $\lambda$ define a representation $\rho_\lambda:\pi_1(X)\to GL(n,\mathbb{C})$ as follows: $\rho_\lambda(g)=\lambda^{\xi(g)}\cdot\rho(g)$. Denote by $J_k(\phi_*, \lambda)$ the maximal size of a Jordan block of eigenvalue $\lambda$ of the automorphism $\phi_*$ in the homology of degree $k$. The main result of the paper says that $J_k(\phi_*, \lambda)$ is equal to the maximal length of a non-zero Massey product of the form $\langle \xi, \ldots , \xi, x\rangle$ where $x\in H^k(X,\rho)$ (here the length means the number of entries of $\xi$). In particular, $\phi_*$ is diagonalizable, if a suitable formality condition holds for the manifold $X$. This is the case if $X$ a compact K\"ahler manifold and $\rho$ is a semisimple representation. The proof of the main theorem is based on the fact that the above Massey products can be identified with differentials in a Massey spectral sequence, which in turn can be explicitly computed in terms of the Jordan normal form of $\phi_*$.


title: Variable-length codes for channels with memory and feedback: error-exponent upper bounds

abstract: The reliability function of memoryless channels with noiseless feedback and variable-length coding has been found to be a linear function of the average rate in the classic work of Burnashev. In this work we consider unifilar channels with noiseless feedback and study upper bounds for the channel reliability function with variable length codes. In unifilar channels the channel state is known to the transmitter but is unknown to the receiver. We generalize Burnashev's analysis and derive a similar expression which is linear in average rate and depends on the channel capacity, as well as an additional parameter which relates to a sequential binary hypothesis testing problem over this channel. This parameter is evaluated by setting up an appropriate Markov decision process (MDP). Furthermore, an upper bound for this parameter is derived using a simplified MDP. Numerical evaluation of the parameter for several binary input/state/output unifilar channels hints at the optimal transmission strategies. Such strategies are studied in a companion paper to provide lower (achievable) bounds on the channel reliability function.


title: Characterization of the Split Closure via Geometric Lifting

abstract: We analyze split cuts from the perspective of cut generating functions via geometric lifting. We show that $\alpha$-cuts, a natural higher-dimensional generalization of the $k$-cuts of Cornu\'{e}jols et al., gives all the split cuts for the mixed-integer corner relaxation. As an immediate consequence we obtain that the $k$-cuts are equivalent to split cuts for the 1-row mixed-integer relaxation. Further, we show that split cuts for finite-dimensional corner relaxations are restrictions of split cuts for the infinite-dimensional relaxation. In a final application of this equivalence, we exhibit a family of pure-integer programs whose split closures have arbitrarily bad integrality gap. This complements the mixed-integer example provided by Basu et al [On the relative strength of split, triangle and quadrilateral cuts, Math. Program. 126(2):281--314, 2011].


title: Growth Models for Tree Stems and Vines

abstract: The paper introduces a PDE model for the growth of a tree stem or a vine. The equations describe the elongation due to cell growth, and the response to gravity and to external obstacles. An additional term accounts for the tendency of a vine to curl around branches of other plants. When obstacles are present, the model takes the form of a differential inclusion with unilateral constraints. At each time t, a cone of admissible reactions is determined by the minimization of an elastic deformation energy. The main theorem shows that local solutions exist and can be prolonged globally in time, except when a specific "breakdown configuration" is reached. Approximate solutions are constructed by an operator-splitting technique. Some numerical simulations are provided at the end of the paper.


title: Variable-length codes for channels with memory and feedback: error-exponent lower bounds

abstract: The reliability function of memoryless channels with noiseless feedback and variable-length coding has been found to be a linear function of the average rate in the classic work of Burnashev. In this work we consider unifilar channels with noiseless feedback and study specific transmission schemes, the performance of which provides lower bounds for the channel reliability function. In unifilar channels the channel state evolves in a deterministic fashion based on the previous state, input, and output, and is known to the transmitter but is unknown to the receiver. We consider a two-stage transmission scheme. In the first stage, both transmitter and receiver summarize their common information in an M-dimensional vector with elements in the state space of the unifilar channel and an M-dimensional probability mass function, with M being the number of messages. The second stage, which is entered when one of the messages is sufficiently reliable, is resolving a binary hypothesis testing problem. The analysis assumes the presence of some common randomness shared by the transmitter and receiver, and is based on the study of the log-likelihood ratio of the transmitted message posterior belief, and in particular on the study of its multi-step drift. Simulation results confirm that the bounds are tight compared to the upper bounds derived in a companion paper.


title: Medium effects on pion production in heavy ion collisions

abstract: Within the framework of the relativistic Vlasov-Uehling-Uhlenbeck transport model based on the relativistic nonlinear NL$\rho$ interaction, we study pion in-medium effects on the $\pi^-/\pi^+$ ratio in Au+Au collisions at the energy of $E/A=400~\mathrm{MeV}$. These effects include the isospin-dependent pion s-wave and p-wave potentials, which are taken from calculations based on the chiral perturbation theory and the $\Delta$-hole model, respectively. We find that the $\pi^-/\pi^+$ ratio in this collision is suppressed by the pion s-wave potential but enhanced by the p-wave potential, with a net effect of a significantly suppressed $\pi^-/\pi^+$ ratio. Including also the in-medium threshold effects on $\Delta$ resonance production and decay and using a nuclear symmetry energy with a slope parameter $L=59~\mathrm{MeV}$ by reducing the coupling of isovector-vector $\rho$ meson to nucleon, our result is in good agreement with measured $\pi^-/\pi^+$ ratio from the FOPI Collaboration. We further investigate the pion in-medium effects on the ratio of charged pions as a function of their kinetic energies.


title: Dynamics of Fragmented Condensates and Macroscopic Entanglement

abstract: The relative phase of the order parameters in the collision of two condensates can influence the outcome of their collision in the case of weak coupling. With increasing interaction strength however, the initially independent phases of the two order parameters in the colliding partners quickly become phase locked, as the strong coupling favors an overall phase rigidity of the entire condensate, and upon their separation the emerging superfluid fragments become entangled.


title: Transference of Fermi Contour Anisotropy to Composite Fermions

abstract: There has been a surge of recent interest in the role of anisotropy in interaction-induced phenomena in two-dimensional (2D) charged carrier systems. A fundamental question is how an anisotropy in the energy-band structure of the carriers at zero magnetic field affects the properties of the interacting particles at high fields, in particular of the composite fermions (CFs) and the fractional quantum Hall states (FQHSs). We demonstrate here tunable anisotropy for holes and hole-flux CFs confined to GaAs quantum wells, via applying \textit{in situ} in-plane strain and measuring their Fermi wavevector anisotropy through commensurability oscillations. For strains on the order of $10^{-4}$ we observe significant deformations of the shapes of the Fermi contours for both holes and CFs. The measured Fermi contour anisotropy for CFs at high magnetic field ($\alpha_\mathrm{CF}$) is less than the anisotropy of their low-field hole (fermion) counterparts ($\alpha_\mathrm{F}$), and closely follows the relation: $\alpha_\mathrm{CF} = \sqrt{\alpha_\mathrm{F}}$. The energy gap measured for the $\nu = 2/3$ FQHS, on the other hand, is nearly unaffected by the Fermi contour anisotropy up to $\alpha_\mathrm{F} \sim 3.3$, the highest anisotropy achieved in our experiments.


title: Is there any polynomial upper bound for the universal labeling of graphs?

abstract: A {\it universal labeling} of a graph $G$ is a labeling of the edge set in $G$ such that in every orientation $\ell$ of $G$ for every two adjacent vertices $v$ and $u$, the sum of incoming edges of $v$ and $u$ in the oriented graph are different from each other. The {\it universal labeling number} of a graph $G$ is the minimum number $k$ such that $G$ has {\it universal labeling} from $\{1,2,\ldots, k\}$ denoted it by $\overrightarrow{\chi_{u}}(G) $. We have $2\Delta(G)-2 \leq \overrightarrow{\chi_{u}} (G)\leq 2^{\Delta(G)}$, where $\Delta(G)$ denotes the maximum degree of $G$. In this work, we offer a provocative question that is:" Is there any polynomial function $f$ such that for every graph $G$, $\overrightarrow{\chi_{u}} (G)\leq f(\Delta(G))$?". Towards this question, we introduce some lower and upper bounds on their parameter of interest. Also, we prove that for every tree $T$, $\overrightarrow{\chi_{u}}(T)=\mathcal{O}(\Delta^3) $. Next, we show that for a given 3-regular graph $G$, the universal labeling number of $G$ is 4 if and only if $G$ belongs to Class 1. Therefore, for a given 3-regular graph $G$, it is an $ \mathbf{NP} $-complete to determine whether the universal labeling number of $G$ is 4. Finally, using probabilistic methods, we almost confirm a weaker version of the problem.


title: Nested Markov Properties for Acyclic Directed Mixed Graphs

abstract: Conditional independence models associated with directed acyclic graphs (DAGs) may be characterized in at least three different ways: via a factorization, the global Markov property (given by the d-separation criterion), and the local Markov property. Marginals of DAG models also imply equality constraints that are not conditional independences; the well-known ``Verma constraint'' is an example. Constraints of this type are used for testing edges, and in a computationally efficient marginalization scheme via variable elimination. We show that equality constraints like the ``Verma constraint'' can be viewed as conditional independences in kernel objects obtained from joint distributions via a fixing operation that generalizes conditioning and marginalization. We use these constraints to define, via ordered local and global Markov properties, and a factorization, a graphical model associated with acyclic directed mixed graphs (ADMGs). We prove that marginal distributions of DAG models lie in this model, and that a set of these constraints given by Tian provides an alternative definition of the model. Finally, we show that the fixing operation used to define the model leads to a particularly simple characterization of identifiable causal effects in hidden variable causal DAG models.


title: On the Average Locality of Locally Repairable Codes

abstract: A linear block code with dimension $k$, length $n$, and minimum distance $d$ is called a locally repairable code (LRC) with locality $r$ if it can retrieve any coded symbol by at most $r$ other coded symbols. LRCs have been recently proposed and used in practice in distributed storage systems (DSSs) such as Windows Azure storage and Facebook HDFS-RAID. Theoretical bounds on the maximum locality of LRCs ($r$) have been established. The \textit{average} locality of an LRC ($\overline{r}$) directly affects the costly repair bandwidth, disk I/O, and number of nodes involved in the repair process of a missing data block. There is a gap in the literature studying $\overline{r}$. In this paper, we establish a lower bound on $\overline{r}$ of arbitrary $(n,k,d)$ LRCs. Furthermore, we obtain a tight lower bound on $\overline{r}$ for a practical case where the code rate $(R=\frac{k}{n})$ is greater than $(1-\frac{1}{\sqrt{n}})^2$. Finally, we design three classes of LRCs that achieve the obtained bounds on $\overline{r}$. Comparing with the existing LRCs, our proposed codes improve the average locality without sacrificing such crucial parameters as the code rate or minimum distance.


title: Quantum Information on Spectral Sets

abstract: For convex optimization problems Bregman divergences appear as regret functions. Such regret functions can be defined on any convex set but if a sufficiency condition is added the regret function must be proportional to information divergence and the convex set must be spectral. Spectral set are sets where different orthogonal decompositions of a state into pure states have unique mixing coefficients. Only on such spectral sets it is possible to define well behaved information theoretic quantities like entropy and divergence. It is only possible to perform measurements in a reversible way if the state space is spectral. The most important spectral sets can be represented as positive elements of Jordan algebras with trace 1. This means that Jordan algebras provide a natural framework for studying quantum information. We compare information theory on Hilbert spaces with information theory in more general Jordan algebras, and conclude that much of the formalism is unchanged but also identify some important differences.


title: Original proofs of Stirling's series for log(n!)

abstract: Transcription into modern notations of the derivation by Stirling and De Moivre of an asymptotic series for $\log(n!)$, usually called Stirling's series. The previous discovery by Wallis of an infinite product for $\pi$, and later results on the divergence of the series are also presented. We conclude that James Stirling has priority over Abraham de Moivre for Stirling's formula and Stirling's series.


title: On the almost Gorenstein property of determinantal rings

abstract: In this paper we investigate the question of when the determinantal ring $R$ over a field $k$ is an almost Gorenstein local/graded ring in the sense of Goto, Takahashi, and the author. As a consequence of the main result, we see that if $R$ is a non-Gorenstein almost Gorenstein local/graded ring, then the ring $R$ has a minimal multiplicity.


title: Maximal immediate extensions of valued differential fields

abstract: We show that every valued differential field has an immediate strict extension that is spherically complete. We also discuss the issue of uniqueness up to isomorphism of such an extension.


title: A geometric approach to cut-generating functions

abstract: The cutting-plane approach to integer programming was initiated more that 40 years ago: Gomory introduced the corner polyhedron as a relaxation of a mixed integer set in tableau form and Balas introduced intersection cuts for the corner polyhedron. This line of research was left dormant for several decades until relatively recently, when a paper of Andersen, Louveaux, Weismantel and Wolsey generated renewed interest in the corner polyhedron and intersection cuts. Recent developments rely on tools drawn from convex analysis, geometry and number theory, and constitute an elegant bridge between these areas and integer programming. We survey these results and highlight recent breakthroughs in this area.


title: Emergence of Bloch oscillations in one-dimensional systems

abstract: Electrons in periodic potentials exhibit oscillatory motion in presence of an electric field. Such oscillations are known as Bloch oscillations. In this article we theoretically investigate the emergence of Bloch oscillations for systems where the electric field is confined to a finite region, like in typical electronic devices. We use a one-dimensional tight-binding model within the single-band approximation to numerically study the dynamics of electrons after a sudden switching-on of the electric field. We find a transition from a regime with direct current to Bloch oscillations when increasing the system size or decreasing the field strength. We propose a pump-probe scheme to observe the oscillations by measuring the accumulated charge as a function of the pulse-length.


title: Grain boundary diffusion in severely deformed Al-based alloy

abstract: Grain boundary diffusion in severely deformed Al-based AA5024 alloy is investigated. Different states are prepared by combination of equal channel angular processing and heat treatments, with the radioisotope $^{57}$Co being employed as a sensitive probe of a given grain boundary state. Its diffusion rates near room temperature (320~K) are utilized to quantify the effects of severe plastic deformation and a presumed formation of a previously reported deformation-modified state of grain boundaries, solute segregation at the interfaces, increased dislocation content after deformation and of the precipitation behavior on the transport phenomena along grain boundaries. The dominant effect of nano-sized Al$_3$Sc-based precipitates is evaluated using density functional theory and the Eshelby model for the determination of elastic stresses around the precipitates.


title: The Rational Sectional Category of Certain Universal Fibrations

abstract: We prove that the sectional category of the universal fibration with fibre X, for X any space that satisfies a well-known conjecture of Halperin, equals one after rationalization.


title: Achromatic Metalens in the Visible and Metalens with Reverse Chromatic Dispersion

abstract: In this letter, we experimentally report an achromatic metalens (AML) operating over a continuous bandwidth in the visible. This is accomplished via dispersion engineering of dielectric phase shifters: titanium dioxide nanopillars tiled on a dielectric spacer layer above a metallic mirror. The AML works in reflection mode with a focal length independent of wavelength from wavelength of 490 nm to 550 nm. We also design a metalens with reverse chromatic dispersion, where the focal length increases as the wavelength increases, contrary to conventional diffractive lenses. The ability to engineer the chromatic dispersion of metalenses at will enables a wide variety of applications that were not previously possible. In particular, for the AML design, we envision applications such as imaging under LED illumination, fluorescence, and photoluminescence spectroscopy.