title: The superpotential method in cosmological inflation

abstract: Scalar field cosmological inflation has a first integral relating the Hubble function and the lagrangian of the scalar field(s), which is known under the names of "Hamilton-Jacobi approach" or "superpotential equation". Here we exploit the simplicity of this superpotential equation and use it as an alternative but equivalent cosmological evolution equation during inflation, replacing the Friedman-Robertson-Walker (FRW) equations. It turns out that all inflationary observables can be calculated directly from its solution (the superpotential). Further, the superpotential equation allows for a simple and direct calculation of the slow-roll expansion to arbitrary order and, in many cases, for an exact determination of the slow-roll attractor. It also allows for a power series expansion in the inflaton field which permits to estimate the radius of convergence of the slow-roll expansion. We consider several examples of single-field inflationary models to demonstrate the simplicity and usefulness of the method.


title: Turbulent dynamo in a weakly ionized medium

abstract: The small-scale turbulent dynamo is an important process contributing to the cosmic magnetization. In partially ionized astrophysical plasmas, the dynamo growth of magnetic energy strongly depends on the coupling state between ions and neutrals and the ion-neutral collisional damping effect. A new damping stage of turbulent dynamo in a weakly ionized medium was theoretically predicted by Xu \& Lazarian (2016). By carrying out a 3D two-fluid dynamo simulation, here we for the first time numerically confirmed the physical conditions and the linear-in-time growth of magnetic field strength of the damping stage of dynamo. The dynamo-amplified magnetic field has a characteristic length as the damping scale, which increases with time and can reach the injection scale of turbulence after around eight largest eddy-turnover times given sufficiently low ionization fraction and weak initial magnetic field. Due to the weak coupling between ions and neutrals, most turbulent energy carried by neutrals cannot be converted to the magnetic energy, resulting in a relatively weak magnetic field at the end of dynamo. This result has important implications for the growth of magnetic fields in the partially ionized interstellar medium and shock acceleration of Galactic cosmic rays.


title: Black Holes That Repel

abstract: The recent observation that black holes in certain Einstein-Maxwell-Dilaton (EMD) theories can violate the entropy super-additivity led to the suggestion that these black holes might repel each other. In this paper, we consider EMD theories with two Maxwell fields $A_i$, with general exponential couplings $\exp(a_i \phi)$ in their kinetic terms. We calculate the gravi-electrostatic force between charged black holes $(m_1,e_1)$ and $(M_2,Q_2)$; the former is sufficiently small and can be treated as a point-like object. We find there is a potential barrier caused by the dilaton coupling at $r_0$ outside the back hole horizon $r_+$, provided that $-a_1 a_2> 2(D-3)/(D-2)$. As the black hole approaches extremality, both $r_+$ and $r_0$ vanish, the barrier becomes infinitesimally thin but infinitely high, and the two black holes repel each other in the whole space. There is no electrostatic force between them; the dilaton is the antigravity agent. Furthermore we find that the exact constraint on $a_1 a_2$ can be derived from the requirements that two-charged extremal black holes have a fusion bomb like mass formula and the violation of entropy super-additivity can occur. The two very different approaches give a consistent picture of the black hole repulsion.


title: Dwarf galaxies might not be the birth sites of supermassive black holes

abstract: Intermediate-mass black holes (BHs) in local dwarf galaxies are considered the relics of the early seed BHs. However, their growth might have been impacted by galaxy mergers and BH feedback so that they cannot be treated as tracers of the early seed BH population.


title: The inner dark matter distribution of the Cosmic Horseshoe (J1148+1930) with gravitational lensing and dynamics

abstract: We present a detailed analysis of the inner mass structure of the Cosmic Horseshoe (J1148+1930) strong gravitational lens system observed with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3). In addition to the spectacular Einstein ring, this systems shows a radial arc. We obtained the redshift of the radial arc counter image $z_\text{s,r} = 1.961 \pm 0.001$ from Gemini observations. To disentangle the dark and luminous matter, we consider three different profiles for the dark matter distribution: a power-law profile, the NFW, and a generalized version of the NFW profile. For the luminous matter distribution, we base it on the observed light distribution that is fitted with three components: a point mass for the central light component resembling an active galactic nucleus, and the remaining two extended light components scaled by a constant M/L. To constrain the model further, we include published velocity dispersion measurements of the lens galaxy and perform a self-consistent lensing and axisymmetric Jeans dynamical modeling. Our model fits well to the observations including the radial arc, independent of the dark matter profile. Depending on the dark matter profile, we get a dark matter fraction between 60 % and 70 %. With our composite mass model we find that the radial arc helps to constrain the inner dark matter distribution of the Cosmic Hoseshoe independently of the dark matter profile.


title: Climate sensitivity to ozone and its relevance on the habitability of Earth-like planets

abstract: Atmospheric ozone plays an important role on the temperature structure of the atmosphere. However, it has not been included in previous studies on the effect of an increasing solar radiation on the Earth's climate. Here we study the climate sensitivity to the presence/absence of ozone with an increasing solar forcing for the first time with a global climate model. We show that the warming effect of ozone increases both the humidity of the lower atmosphere and the surface temperature. Under the same solar irradiance, the mean surface temperature is 7 K higher than in an analog planet without ozone. Therefore, the moist greenhouse threshold, the state at which water vapor becomes abundant in the stratosphere, is reached at a lower solar irradiance (1572 Wm^-2 with respect to 1647 Wm^-2 in the case without ozone). Our results imply that ozone reduces the maximum solar irradiance at which Earth-like planets would remain habitable.


title: Four-loop quark form factor with quartic fundamental colour factor

abstract: We analytically compute the four-loop QCD corrections for the colour structure $(d_F^{abcd})^2$ to the massless non-singlet quark form factor. The computation involves non-trivial non-planar integral families which have master integrals in the top sector. We compute the master integrals by introducing a second mass scale and solving differential equations with respect to the ratio of the two scales. We present details of our calculational procedure. Analytical results for the cusp and collinear anomalous dimensions, and the finite part of the form factor are presented. We also provide analytic results for all master integrals expanded up to weight eight.


title: Quantum acousto-optic control of light-matter interactions in nanophotonic networks

abstract: We analyze the coupling of atoms or atom-like emitters to nanophotonic waveguides in the presence of propagating acoustic waves. Specifically, we show that strong index modulations induced by such waves can drastically modify the effective photonic density of states and thereby influence the strength, the directionality, as well as the overall characteristics of photon emission and absorption processes. These effects enable a complete dynamical control of light-matter interactions in waveguide structures, which even in a two dimensional system can be used to efficiently exchange individual photons along selected directions and with a very high fidelity. Such a quantum acousto-optical control provides a versatile tool for various quantum networking applications ranging from the distribution of entanglement via directional emitter-emitter interactions to the routing of individual photonic quantum states via acoustic conveyor belts.


title: The birth of the Milky Way as uncovered by accurate stellar ages with Gaia

abstract: Knowledge of ages for stars formed over a galaxy's lifetime is fundamental to understand its formation and evolution. However, stellar ages are difficult to obtain since they cannot be measured from observations, being comparison with stellar models (Soderblom 2010) required. Alternatively, age distributions can be derived applying the robust technique of colour-magnitude diagram fitting (Gallart et al. 2005), till now mainly employed to study nearby galaxies. The new distances to individual Milky Way stars from the Gaia mission (Brown et al. 2018) have allowed us to use this technique to derive ages from a thick disk colour-magnitude diagram, and from the enigmatic, two-sequenced colour-magnitude diagram of the kinematically hot local halo (Babusiaux et al. 2018), which blue-sequence has been linked to a major accretion event (Haywood et al. 2018, Helmi et al. 2018). Because accurate ages were lacking, the time of the merger and its role on our Galaxy's early evolution remained unclear. We show that the stars in both halo sequences share identical age distributions, and are older than the bulk of thick disc stars. The sharp halo age cut 10 Gyr ago can be identified with the accretion of Gaia-Enceladus. Along with state-of-the-art cosmological simulations of galaxy formation (Brook et al. 2012), these robust ages allow us to order the early sequence of events that shaped our Galaxy, identifying the red-sequence as the first stars formed within the Milky Way progenitor which, because of their kinematics, can be described as its long sought in-situ halo.


title: Climate Sensitivity to Carbon Dioxide and Moist Greenhouse threshold of Earth-like planets under an increasing solar forcing

abstract: Carbon dioxide is one of the major contributors to the radiative forcing, increasing both the temperature and the humidity of Earth's atmosphere. If the stellar irradiance increases and water becomes abundant in the stratosphere of an Earth-like planet, it will be dissociated and the resultant hydrogen will escape from the atmosphere. This state is called the moist greenhouse threshold (MGT). Using a global climate model (GCM) of intermediate complexity, we explore how to identify this state for different CO2 concentrations and including the radiative effect of atmospheric ozone for the first time. We show that the moist greenhouse threshold correlates with the inflection point in the water vapor mixing ratio in the stratosphere and a peak in the climate sensitivity. For CO2 concentrations between 560 ppm and 200 ppm, the moist greenhouse threshold is reached at a surface temperature of 320 K. Despite the higher simplicity of our model, our results are consistent with similar simulations without ozone by complex GCMs, suggesting that they are robust indicators of the MGT. We discuss the implications for inner edge of the habitable zone as well as the water loss timescales for Earth analog planets.


title: Length scales in the many-body localized phase and their spectral signatures

abstract: We compute and compare the decay lengths of several correlation functions and effective coupling constants in the many-body localized (MBL) phase. To this end, we consider the distribution of the logarithms of these couplings and correlators: in each case the log-coupling follows a normal distribution with mean and variance that grow linearly with separation. Thus, a localization length is asymptotically sharply defined for each of these quantities. These localization lengths differ numerically from one another, but all of them remain short up to the numerically observed MBL transition, indicating stability of the MBL phase against isolated ergodic inclusions. We also show how these broad distributions may be extracted using interferometric probes such as double electron-electron resonance (DEER) and the statistics of local spin precession frequencies.


title: Multiphase gas in the circumgalactic medium: relative role of $t_{\rm cool}/t_{\rm ff}$ and density fluctuations

abstract: We perform a suite of simulations with realistic gravity and thermal balance in shells to quantify the role of the ratio of cooling time to the free-fall time ($t_{\rm cool}/t_{\rm ff}$) and the amplitude of density perturbations ($\delta \rho/\rho$) in the production of multiphase gas in the circumgalactic medium (CGM). Previous idealized simulations, focussing on small amplitude perturbations in the intracluster medium (ICM), found that cold gas can condense out of the hot ICM in global thermal balance when the background $t_{\rm cool}/t_{\rm ff} \lesssim 10$. Recent observations suggest the presence of cold gas even when the background profiles have somewhat large values of ${t_{\rm cool}/t_{\rm ff}}$. This partly motivates a better understanding of additional factors such as large density perturbations that can enhance the propensity for cooling and condensation even when the background ${t_{\rm cool}/t_{\rm ff}}$ is high. Such large density contrasts can be seeded by galaxy wakes or dense cosmological filaments. From our simulations, we introduce a condensation curve in the $(\delta \rho/\rho)$ - min$(t_{\rm cool}/t_{\rm ff})$ space, that defines the threshold for condensation of multiphase gas in the CGM. We show that this condensation curve corresponds to ${(t_{\rm cool}/t_{\rm ff})}_{\rm blob} \lesssim 10$ applied to the overdense blob instead of the background for which $t_{\rm cool}/t_{\rm ff}$ can be higher. We also study the modification in the condensation curve by varying entropy stratification. Steeper (positive) entropy gradients shift the condensation curve to higher amplitudes of perturbations (i.e., make condensation difficult). A constant entropy core, applicable to the CGM in smaller halos, shows condensation over a larger range of radii as compared to the steeper entropy profiles in the ICM.


title: A Novel Approach to Constrain Rotational Mixing & Convective-Core Overshoot in Stars Using the Initial-Final Mass Relation

abstract: The semi-empirical initial-final mass relation (IFMR) connects spectroscopically analyzed white dwarfs in star clusters to the initial masses of the stars that formed them. Most current stellar evolution models, however, predict that stars will evolve to white dwarfs $\sim$0.1 M$_\odot$ less massive than that found in the IFMR. We first look at how varying theoretical mass-loss rates, third dredge-up efficiencies, and convective-core overshoot may help explain the differences between models and observations. These parameters play an important role at the lowest masses (M$_{\rm initial}$ $<$ 3 M$_\odot$). At higher masses, only convective-core overshoot meaningfully affects white dwarf mass, but alone it likely cannot explain the observed white dwarf masses nor why the IFMR scatter is larger than observational errors predict. These higher masses, however, are also where rotational mixing in main sequence stars begins to create more massive cores, and hence more massive white dwarfs. This rotational mixing also extends a star's lifetime, making faster rotating progenitors appear like less massive stars in their semi-empirical age analysis. Applying the observed range of young B-dwarf rotations to the MIST or SYCLIST rotational models demonstrates a marked improvement in reproducing both the observed IFMR data and its scatter. The incorporation of both rotation and efficient convective-core overshoot significantly improves the match with observations. This work shows that the IFMR provides a valuable observational constraint on how rotation and convective-core overshoot affect the core evolution of a star.


title: Observation of Interfacial Antiferromagnetic Coupling between Magnetic Topological Insulator and Antiferromagnetic Insulator

abstract: Inducing magnetic orders in a topological insulator (TI) to break its time reversal symmetry has been predicted to reveal many exotic topological quantum phenomena. The manipulation of magnetic orders in a TI layer can play a key role in harnessing these quantum phenomena towards technological applications. Here we fabricated a thin magnetic TI film on an antiferromagnetic (AFM) insulator Cr2O3 layer and found that the magnetic moments of the magnetic TI layer and the surface spins of the Cr2O3 layers favor interfacial AFM coupling. Field cooling studies show a crossover from negative to positive exchange bias clarifying the competition between the interfacial AFM coupling energy and the Zeeman energy in the AFM insulator layer. The interfacial exchange coupling also enhances the Curie temperature of the magnetic TI layer. The unique interfacial AFM alignment in magnetic TI on AFM insulator heterostructures opens a new route toward manipulating the interplay between topological states and magnetic orders in spin-engineered heterostructures, facilitating the exploration of proof-of-concept TI-based spintronic and electronic devices with multi-functionality and low power consumption.


title: Fixed Points of Parking Functions

abstract: We define an action of words in $[m]^n$ on $\mathbb{R}^m$ to give a new characterization of rational parking functions -- they are exactly those words whose action has a fixed point. We use this viewpoint to give a simple definition of Gorsky, Mazin, and Vazirani's zeta map on rational parking functions when m and n are coprime, and prove that this zeta map is invertible. A specialization recovers Loehr and Warrington's sweep map on rational Dyck paths.


title: Learning by Fictitious Play in Large Populations

abstract: We consider learning by fictitious play in a large population of agents engaged in single-play, two-person rounds of a symmetric game, and derive a mean-filed type model for the corresponding stochastic process. Using this model, we describe qualitative properties of the learning process and discuss its asymptotic behavior. Of the special interest is the comparative characteristics of the fictitious play learning with and without a memory factor. As a part of the analysis, we show that the model leads to the continuous, best-response dynamics equation of Gilboa and Matsui (1991), when all agents have similar empirical probabilities.


title: Adaptive Learning in Large Populations

abstract: We consider the adaptive learning rule of Harley (1981) for behavior selection in symmetric conflict games in large populations. The rule uses organisms' past, accumulated rewards as the predictor for the future behavior, and can be traced in many life forms from bacteria to humans. We derive a partial differential equation (PDE) that describes the stochastic learning in a population of agents. The equation has simple structure of the `conservation of mass'-type equation in the space of stimuli to engage in a particular type of behavior. We analyze the solutions of the PDE model for typical 2x2 games. It is found that in games with small residual stimuli, adaptive learning rules with faster memory decay have an evolutionary advantage.


title: Combining Dynamical Quantum Typicality and Numerical Linked Cluster Expansions

abstract: We demonstrate that numerical linked cluster expansions (NLCE) yield a powerful approach to calculate time-dependent correlation functions for quantum many-body systems in one dimension. As a paradigmatic example, we study the dynamics of the spin current in the spin-1/2 XXZ chain for different values of anisotropy, as well as in the presence of an integrability-breaking next-nearest neighbor interaction. For short to intermediate time scales, we unveil that NLCE yields a convergence towards the thermodynamic limit already for small cluster sizes, which is much faster than in direct calculations of the autocorrelation function for systems with open or periodic boundary conditions. Most importantly, we show that the range of accessible cluster sizes in NLCE can be extended by evaluating the contributions of larger clusters by means of a pure-state approach based on the concept of dynamical quantum typicality (DQT). Even for moderate computational effort, this combination of DQT and NLCE provides a competitive alternative to existing state-of-the-art techniques, which may be applied in higher dimensions as well.


title: The space density of z>4 blazars

abstract: High redshift blazars are an important class of Active Galactic Nuclei (AGN) that can provide an independent estimate of the supermassive black-hole mass function in high redshift radio-loud AGN without the bias due to absorption along the line-of-sight. Using the Cosmic Lens All Sky Survey (CLASS) we built a complete radio flux-limited sample of high redshift (z>4) blazars suitable for statistical studies. By combining dedicated optical observations and the SDSS spectroscopic database, we obtained a sample of 26 blazar candidates with a spectroscopic redshift above 4. On the basis of their radio spectrum we distinguish between blazars and QSO with a Gigahertz Peaked Spectrum (GPS) like spectrum. Out of the 18 confirmed blazars 14 constitute a completely identified, flux-limited sample down to a magnitude of 21 (AB). Using this complete sample we derive a space density of blazars with 4<z<5.5 of rho=0.13 (+0.05,-0.03) Gpc^-3. This is the first actual estimate of the blazar space density in this range of redshift. This value is in good agreement with the extrapolation of the luminosity function and cosmological evolution based on a sample of flat-spectrum radio quasars selected at lower redshifts and it is consistent with a cosmological evolution peaking at z$\sim$2 similar to radio-quiet QSO. We do not confirm, instead, the presence of a peak at z~4 in the space density evolution, recently suggested using an X-ray selected sample of blazars. It is possible that this extreme peak of the evolution is present only among the most luminous blazars.


title: Myocardial Infarction Quantification From Late Gadolinium Enhancement MRI Using Top-hat Transforms and Neural Networks

abstract: Significance: Late gadolinium enhanced magnetic resonance imaging (LGE-MRI) is the gold standard technique for myocardial viability assessment. Although the technique accurately reflects the damaged tissue, there is no clinical standard for quantifying myocardial infarction (MI), demanding most algorithms to be expert dependent. Objectives and Methods: In this work a new automatic method for MI quantification from LGE-MRI is proposed. Our novel segmentation approach is devised for accurately detecting not only hyper-enhanced lesions, but also microvascular-obstructed areas. Moreover, it includes a myocardial disease detection step which extends the algorithm for working under healthy scans. The method is based on a cascade approach where firstly, diseased slices are identified by a convolutional neural network (CNN). Secondly, by means of morphological operations a fast coarse scar segmentation is obtained. Thirdly, the segmentation is refined by a boundary-voxel reclassification strategy using an ensemble of CNNs. For its validation, reproducibility and further comparison against other methods, we tested the method on a big multi-field expert annotated LGE-MRI database including healthy and diseased cases. Results and Conclusion: In an exhaustive comparison against nine reference algorithms, the proposal achieved state-of-the-art segmentation performances and showed to be the only method agreeing in volumetric scar quantification with the expert delineations. Moreover, the method was able to reproduce the intra- and inter-observer variability ranges. It is concluded that the method could suitably be transferred to clinical scenarios.


title: Generating functions for finite sums involving higher powers of binomial coefficients: Analysis of hypergeometric functions including new families of polynomials and numbers

abstract: The origin of this study is based on not only explicit formulas of finite sums involving higher powers of binomial coefficients, but also explicit evaluations of generating functions for this sums. It should be emphasized that this study contains both new results and literature surveys about some of the related results that have existed so far. With the aid of hypergeometric function, generating functions for a new family of the combinatorial numbers, related to finite sums of powers of binomial coefficients, are constructed. By using these generating functions, a number of new identities have been obtained and at the same time previously well-known formulas and identities have been generalized. Moreover, on this occasion, we identify new families of polynomials including some families of well-known numbers such as Bernoulli numbers, Euler numbers, Stirling numbers, Franel numbers, Catalan numbers, Changhee numbers, Daehee numbers and the others, and also for the polynomials such as the Legendre polynomials, Michael Vowe polynomial, the Mirimanoff polynomial, Golombek type polynomials, and the others. We also give both Riemann and $p$-adic integral representations of these polynomials. Finally, we give combinatorial interpretations of these new families of numbers, polynomials and finite sums of the powers of binomial coefficients. We also give open questions for ordinary generating functions for these numbers.


title: Hybrid Codes

abstract: A hybrid code can simultaneously encode classical and quantum information into quantum digits such that the information is protected against errors when transmitted through a quantum channel. It is shown that a hybrid code has the remarkable feature that it can detect more errors than a comparable quantum code that is able to encode the classical and quantum information. Weight enumerators are introduced for hybrid codes that allow to characterize the minimum distance of hybrid codes. Surprisingly, the weight enumerators for hybrid codes do not obey the usual MacWilliams identity.


title: Binary Message Passing Decoding of Product Codes Based on Generalized Minimum Distance Decoding

abstract: We propose a binary message passing decoding algorithm for product codes based on generalized minimum distance decoding (GMDD) of the component codes, where the last stage of the GMDD makes a decision based on the Hamming distance metric. The proposed algorithm closes half of the gap between conventional iterative bounded distance decoding (iBDD) and turbo product decoding based on the Chase--Pyndiah algorithm, at the expense of some increase in complexity. Furthermore, the proposed algorithm entails only a limited increase in data flow compared to iBDD.


title: Revealing interpretable object representations from human behavior

abstract: To study how mental object representations are related to behavior, we estimated sparse, non-negative representations of objects using human behavioral judgments on images representative of 1,854 object categories. These representations predicted a latent similarity structure between objects, which captured most of the explainable variance in human behavioral judgments. Individual dimensions in the low-dimensional embedding were found to be highly reproducible and interpretable as conveying degrees of taxonomic membership, functionality, and perceptual attributes. We further demonstrated the predictive power of the embeddings for explaining other forms of human behavior, including categorization, typicality judgments, and feature ratings, suggesting that the dimensions reflect human conceptual representations of objects beyond the specific task.


title: Inversion of walkaway VSP data in the presence of lateral velocity heterogeneity

abstract: Multi-azimuth walkaway vertical seismic profiling (VSP) is an established technique for the estimation of in situ slowness surfaces and inferring anisotropy parameters. Normally, this the technique requires the assumption of lateral homogeneity, which makes the horizontal slowness components at depths of downhole receivers equal to those measured at the surface. Any violations of this assumption, such as lateral heterogeneity or a nonzero dip of intermediate interfaces, lead to distortions in reconstructed slowness surfaces and, consequently, to errors in estimated anisotropic parameters. Here, we relax the assumption of lateral homogeneity and discuss how to correct VSP data for weak, lateral heterogeneity (LH). We describe a procedure of downward continuation of recorded travel times that accounts for the presence of both vertical inhomogeneity and weak lateral heterogeneity, which produces correct slowness surfaces at depths of downhole receivers. Once the slowness surfaces are found and the desired type of anisotropic model to be inverted is selected, the corresponding anisotropic parameters, providing the best fit to the estimated slowness can be obtained. We invert the slowness surfaces of $P$-waves for parameters of the simplest anisotropic model describing dipping fractures -- transversely isotropic medium with a tilted symmetry axis. Five parameters of this model -- the $P$-wave velocity $V_0$ in the direction of the symmetry axis, Thomsen's anisotropic coefficients $\epsilon$ and $\delta$, the tilt $\nu$, and the azimuth $\beta$ of the symmetry axis can be estimated in a stable manner when maximum source offset is greater than half of the receiver depth.


title: Revealing the Physico-Chemical Basis of Organic Solid-Solid Wetting Deposition: Casimir-Like Forces, Hydrophobic Collapse, and the Role of the Zeta Potential

abstract: Supramolecular self-assembly at the solid-solid interface enables the deposition and monolayer formation of insoluble organic semiconductors under ambient conditions. The underlying process, termed as the Organic Solid-Solid Wetting Deposition (OSWD), generates two-dimensional adsorbates directly from dispersed three-dimensional organic crystals. This straightforward process has important implications in various fields of research and technology, such as in the domains of low-dimensional crystal engineering, the chemical doping and band-gap engineering of graphene, and in the area of field-effect transistor fabrication. However, till date, lack of an in-depth understanding of the physico-chemical basis of the OSWD prevented the identification of important parameters, essential to achieve a better control of the growth of monolayers and supramolecular assemblies with defined structures, sizes, and coverage areas. Here we propose a detailed model for the OSWD, derived from experimental and theoretical results that have been acquired by using the organic semiconductor quinacridone as an example system. The model reveals the vital role of the zeta potential and includes Casimir-like fluctuation-induced forces and the effect of dewetting in hydrophobic nano-confinements. Based on our results, the OSWD of insoluble organic molecules can hence be applied to environmental friendly and low-cost dispersing agents, such as water. In addition, the model substantially enhances the ability to control the OSWD in terms of adsorbate structure and substrate coverage.


title: Making AI meaningful again

abstract: Artificial intelligence (AI) research enjoyed an initial period of enthusiasm in the 1970s and 80s. But this enthusiasm was tempered by a long interlude of frustration when genuinely useful AI applications failed to be forthcoming. Today, we are experiencing once again a period of enthusiasm, fired above all by the successes of the technology of deep neural networks or deep machine learning. In this paper we draw attention to what we take to be serious problems underlying current views of artificial intelligence encouraged by these successes, especially in the domain of language processing. We then show an alternative approach to language-centric AI, in which we identify a role for philosophy.


title: A population dynamics model of cell-cell adhesion incorporating population pressure and density saturation

abstract: We discuss several continuum cell-cell adhesion models based on the underlying microscopic assumptions. We propose an improvement on these models leading to sharp fronts and intermingling invasion fronts between different cell type populations. The model is based on basic principles of localized repulsion and nonlocal attraction due to adhesion forces at the microscopic level. The new model is able to capture both qualitatively and quantitatively experiments by Katsunuma et al. (2016) [J. Cell Biol. 212(5), pp. 561--575]. We also review some of the applications of these models in other areas of tissue growth in developmental biology. We finally explore the resulting qualitative behavior due to cell-cell repulsion.


title: TraceCaps: A Capsule-based Neural Network for Semantic Segmentation

abstract: In this paper, we propose a capsule-based neural network model to solve the semantic segmentation problem. By taking advantage of the extractable part-whole dependencies available in capsule layers, we derive the probabilities of the class labels for individual capsules through a recursive, layer-by-layer procedure. We model this procedure as a traceback pipeline and take it as a central piece to build an end-to-end segmentation network. Under the proposed framework, image-level class labels and object boundaries are jointly sought in an explicit manner, which poses a significant advantage over the state-of-the-art fully convolutional network (FCN) solutions. With the capability to extracted part-whole information, our traceback pipeline can potentially be utilized as the building blocks to design interpretable neural networks. Experiments conducted on modified MNIST and neuroimages demonstrate that our model considerably enhance the segmentation performance compared to the leading FCN variants.


title: Tensor-based subspace learning for tracking salt-dome boundaries

abstract: The exploration of petroleum reservoirs has a close relationship with the identification of salt domes. To efficiently interpret salt-dome structures, in this paper, we propose a method that tracks salt-dome boundaries through seismic volumes using a tensor-based subspace learning algorithm. We build texture tensors by classifying image patches acquired along the boundary regions of seismic sections and contrast maps. With features extracted from the subspaces of texture tensors, we can identify tracked points in neighboring sections and label salt-dome boundaries by optimally connecting these points. Experimental results show that the proposed method outperforms the state-of-the-art salt-dome detection method by employing texture information and tensor-based analysis.


title: Degrees of compression and inertia for free-abelian times free groups

abstract: We introduce the concepts of degree of inertia, $\text{di}_G(H)$, and degree of compression, $\text{dc}_G(H)$, of a finitely generated subgroup $H$ of a given group $G$. For the case of direct products of free-abelian and free groups, we compute the degree of compression and give an upper bound for the degree of inertia. Imposing some technical assumptions to the supremum involved in the definition of degree of inertia, we introduce the notion called restricted degree of inertia, $\text{di}'_G(H)$, and, again for the case $\mathbb{Z}^m \times F_n$, we provide an explicit formula relating it to the restricted degree of inertia of its projection to the free part, $\text{di}'_{F_n}(H\pi)$.


title: A Minimal Incentive-based Demand Response Program With Self Reported Baseline Mechanism

abstract: In this paper, we propose a novel incentive based Demand Response (DR) program with a self reported baseline mechanism. The System Operator (SO) managing the DR program recruits consumers or aggregators of DR resources. The recruited consumers are required to only report their baseline, which is the minimal information necessary for any DR program. During a DR event, a set of consumers, from this pool of recruited consumers, are randomly selected. The consumers are selected such that the required load reduction is delivered. The selected consumers, who reduce their load, are rewarded for their services and other recruited consumers, who deviate from their reported baseline, are penalized. The randomization in selection and penalty ensure that the baseline inflation is controlled. We also justify that the selection probability can be simultaneously used to control SO's cost. This allows the SO to design the mechanism such that its cost is almost optimal when there are no recruitment costs or at least significantly reduced otherwise. Finally, we also show that the proposed method of self-reported baseline outperforms other baseline estimation methods commonly used in practice.


title: On the translation invariant operators in $\ell^p(\mathbb{Z}^d)$

abstract: In this paper we study boundedness of translation invariant operators in the discrete space $\ell^p(\mathbb{Z}^d)$. In this context a Mikhlin type multiplier theorem is given, yielding boundedness for certain known operators . We also give $\ell^p-\ell^q$ boundedness of a discrete wave equation.


title: Convection in a coupled fluid-porous media system

abstract: We perform linear and nonlinear stability analysis for thermal convection in a fluid overlying a saturated porous medium. We use a coupled system, with the Navier-Stokes equations and Darcy's equation governing the free-flow and the porous regions respectively. Incorporating a dynamic pressure term in the Lions interface condition (which specifies the normal force balance across the fluid-medium interface) permits an energy bound on the typically uncooperative nonlinear advection term, enabling new nonlinear stability results. Within certain regimes, the nonlinear stability thresholds agree closely with the linear ones, and we quantify the differences that exist. We then compare stability thresholds produced by several common variants of the tangential interface conditions, using both numerics and asymptotics in the small Darcy number limit. Finally, we investigate the transition between full convection and fluid-dominated convection using both numerics and a heuristic theory. This heuristic theory is based on comparing the ratio of the Rayleigh number in each domain to its corresponding critical value, and it is shown to agree reasonably well with the numerics regarding how the transition depends on the depth ratio, the Darcy number, and the thermal-diffusivity ratio.


title: Three Other Models of Computer System Performance

abstract: This note argues for more use of simple models beyond Amdahl's Law: Bottleneck Analysis, Little's Law, and a M/M/1 Queue.


title: A note on isomorphisms between submonoids of $\mathbb N^k$ and numerical semigroups

abstract: Given a submonoid $H$ of $\mathbb N^k$, we give some characterizations of the minimum $r\in \mathbb N^+$ such that $H$ is isomorphic to a submonoid of $\mathbb N^r$. In the context of submonoids of $\mathbb N$, we prove that if two numerical semigroups are isomorphic submonoids of $\mathbb N$, then they are equal and the identity map is the unique isomorphism between them.


title: Beyond the EM Algorithm: Constrained Optimization Methods for Latent Class Model

abstract: Latent class model (LCM), which is a finite mixture of different categorical distributions, is one of the most widely used models in statistics and machine learning fields. Because of its non-continuous nature and the flexibility in shape, researchers in practice areas such as marketing and social sciences also frequently use LCM to gain insights from their data. One likelihood-based method, the Expectation-Maximization (EM) algorithm, is often used to obtain the model estimators. However, the EM algorithm is well-known for its notoriously slow convergence. In this research, we explore alternative likelihood-based methods that can potential remedy the slow convergence of the EM algorithm. More specifically, we regard likelihood-based approach as a constrained nonlinear optimization problem, and apply quasi-Newton type methods to solve them. We examine two different constrained optimization methods to maximize the log likelihood function. We present simulation study results to show that the proposed methods not only converge in less iterations than the EM algorithm but also produce more accurate model estimators.


title: Report of the Third Global Experimentation for Future Internet (GEFI 2018) Workshop

abstract: The third Global Experimentation for Future Internet (GEFI 2018) workshop was held October 25-26, 2018 in Tokyo, Japan, hosted by the University of Tokyo. A total of forty-four participants attended, representing Belgium, Brazil, China, Denmark, France, Ireland, Japan, the Republic of Korea, and the United States. The workshop employed a mixed format of presentations and open group discussions to advance multi-national coordination and interoperation of research infrastructure for advanced networking and computer science research. Major topic areas included: softwareization and virtualization of radios and networks; testbed support for networking experiments; EdgeNet; a federated testbed of elastic optical networks; and reproducibility in experimentation. Workshop goals included both the formulation of specific new research collaborations and strategies for coordination and interoperation of research testbeds. Workshop outcomes include a variety of new and ongoing collaborative efforts, ranging from an agreement to pursue the development of optical "white boxes" in support of elastic optical testbeds to the identification of strategies for effective use of open-source software and hardware platforms in future research infrastructure.


title: Conditions de stabilit\'e et g\'eom\'etrie birationnelle [d'apr\`es Bridgeland, Bayer-Macr\`i, ...]

abstract: We survey some aspects of stability conditions both in general and on the derived category of coherent sheaves on a surface, with applications to the birational geometry of certain holomorphic symplectic varieties.


title: Multiple propagation paths enhance locating the source of diffusion in complex networks

abstract: We investigate the problem of locating the source of diffusion in complex networks without complete knowledge of nodes' states. Some currently known methods assume the information travels via a single, shortest path, which by assumption is the fastest way. We show that such a method leads to the overestimation of propagation time for synthetic and real networks, where multiple shortest paths as well as longer paths between vertices exist. We propose a new method of source estimation based on maximum likelihood principle, that takes into account existence multiple shortest paths. It shows up to 1.6 times higher accuracy in synthetic and real networks.


title: Inference of Demographic Attributes based on Mobile Phone Usage Patterns and Social Network Topology

abstract: Mobile phone usage provides a wealth of information, which can be used to better understand the demographic structure of a population. In this paper, we focus on the population of Mexican mobile phone users. We first present an observational study of mobile phone usage according to gender and age groups. We are able to detect significant differences in phone usage among different subgroups of the population. We then study the performance of different machine learning (ML) methods to predict demographic features (namely, age and gender) of unlabeled users by leveraging individual calling patterns, as well as the structure of the communication graph. We show how a specific implementation of a diffusion model, harnessing the graph structure, has significantly better performance over other node-based standard ML methods. We provide details of the methodology together with an analysis of the robustness of our results to changes in the model parameters. Furthermore, by carefully examining the topological relations of the training nodes (seed nodes) to the rest of the nodes in the network, we find topological metrics which have a direct influence on the performance of the algorithm.


title: Polynomial-time Capacity Calculation and Scheduling for Half-Duplex 1-2-1 Networks

abstract: This paper studies the 1-2-1 half-duplex network model, where two half-duplex nodes can communicate only if they point `beams' at each other; otherwise, no signal can be exchanged or interference can be generated. The main result of this paper is the design of two polynomial-time algorithms that: (i) compute the approximate capacity of the 1-2-1 half-duplex network and, (ii) find the network schedule optimal for the approximate capacity. The paper starts by expressing the approximate capacity as a linear program with an exponential number of constraints. A core technical component consists of building a polynomial-time separation oracle for this linear program, by using algorithmic tools such as perfect matching polytopes and Gomory-Hu trees.


title: Signal Space Cognitive Cooperation

abstract: In this work, a new transmission scheme, signal space cognitive cooperation, is introduced by applying the idea of signal space diversity in an underlay spectrum sharing decode-and-forward multi-relay cooperative network. In the proposed structure, the secondary source signal is rotated by a certain angle and then the source and the secondary best relay transmit the in-phase and the quadrature components of the rotated signal. As a consequence, two source signals, rather than one, are transmitted in two-time slots which improves data rates considerably, compared to the conventional cognitive cooperative schemes. In this work, proactive relaying mode is used in which the best relay is selected based on the max-min selection criterion before executing the transmission. Considering both statistical and the instantaneous channel state information of the feedback channel between the primary receiver and the secondary network, two power allocation methods are adopted at the source. For both methods, closed-form expressions of error probability are derived. Moreover, asymptotic analysis is performed and diversity gain is obtained to provide further insights about the system performance. Finally, analytical expressions are verified by Monte-Carlo simulations.


title: An optimization framework for simulation and kinematic control of Constrained Collaborative Mobile Agents (CCMA) system

abstract: We present a concept of constrained collaborative mobile agents (CCMA) system, which consists of multiple wheeled mobile agents constrained by a passive kinematic chain. This mobile robotic system is modular in nature, the passive kinematic chain can be easily replaced with different designs and morphologies for different functions and task adaptability. Depending solely on the actuation of the mobile agents, this mobile robotic system can manipulate or position an end-effector. However, the complexity of the system due to presence of several mobile agents, passivity of the kinematic chain and the nature of the constrained collaborative manipulation requires development of an optimization framework. We therefore present an optimization framework for forward simulation and kinematic control of this system. With this optimization framework, the number of deployed mobile agents, actuation schemes, the design and morphology of the passive kinematic chain can be easily changed, which reinforces the modularity and collaborative aspects of the mobile robotic system. We present results, in simulation, for spatial 4-DOF to 6-DOF CCMA system examples. Finally, we present experimental quantitative results for two different fabricated 4-DOF prototypes, which demonstrate different actuation schemes, control and collaborative manipulation of an end-effector.


title: The Mahalanobis kernel for heritability estimation in genome-wide association studies: fixed-effects and random-effects methods

abstract: Linear mixed models (LMMs) are widely used for heritability estimation in genome-wide association studies (GWAS). In standard approaches to heritability estimation with LMMs, a genetic relationship matrix (GRM) must be specified. In GWAS, the GRM is frequently a correlation matrix estimated from the study population's genotypes, which corresponds to a normalized Euclidean distance kernel. In this paper, we show that reliance on the Euclidean distance kernel contributes to several unresolved modeling inconsistencies in heritability estimation for GWAS. These inconsistencies can cause biased heritability estimates in the presence of linkage disequilibrium (LD), depending on the distribution of causal variants. We show that these biases can be resolved (at least at the modeling level) if one adopts a Mahalanobis distance-based GRM for LMM analysis. Additionally, we propose a new definition of partitioned heritability -- the heritability attributable to a subset of genes or single nucleotide polymorphisms (SNPs) -- using the Mahalanobis GRM, and show that it inherits many of the nice consistency properties identified in our original analysis. Partitioned heritability is a relatively new area for GWAS analysis, where inconsistency issues related to LD have previously been known to be especially pernicious.


title: SalSi: A new seismic attribute for salt dome detection

abstract: In this paper, we propose a saliency-based attribute, SalSi, to detect salt dome bodies within seismic volumes. SalSi is based on the saliency theory and modeling of the human vision system (HVS). In this work, we aim to highlight the parts of the seismic volume that receive highest attention from the human interpreter, and based on the salient features of a seismic image, we detect the salt domes. Experimental results show the effectiveness of SalSi on the real seismic dataset acquired from the North Sea, F3 block. Subjectively, we have used the ground truth and the output of different salt dome delineation algorithms to validate the results of SalSi. For the objective evaluation of results, we have used the receiver operating characteristics (ROC) curves and area under the curves (AUC) to demonstrate SalSi is a promising and an effective attribute for seismic interpretation.


title: Private Information Retrieval from Locally Repairable Databases with Colluding Servers

abstract: We consider information-theoretical private information retrieval (PIR) from a coded database with colluding servers. We target, for the first time, locally repairable storage codes (LRCs). We consider any number of local groups $ g $, locality $ r $, local distance $ \delta $ and dimension $ k $. Our main contribution is a PIR scheme for maximally recoverable (MR) LRCs based on linearized Reed--Solomon codes, which achieve the smallest field sizes among MR-LRCs for many parameter regimes. In our scheme, nodes are identified with codeword symbols and servers are identified with local groups of nodes. Only locally non-redundant information is downloaded from each server, that is, only $ r $ nodes (out of $ r+\delta-1 $) are downloaded per server. The PIR scheme achieves the (download) rate $ R = (N - k - rt + 1)/N $, where $ N = gr $ is the length of the MDS code obtained after removing the local parities, and for any $ t $ colluding servers such that $ k + rt \leq N $. For an unbounded number of stored files, the obtained rate is strictly larger than those of known PIR schemes that work for any MDS code. Finally, the obtained PIR scheme can also be adapted when communication between the user and each server is performed via linear network coding, achieving the same rate as previous PIR schemes for this scenario but with polynomial finite field sizes, instead of exponential. Our rates are equal to those of PIR schemes for Reed--Solomon codes, but Reed--Solomon codes are incompatible with the MR-LRC property or linear network coding, thus our PIR scheme is less restrictive in its applications.


title: Quantum Brownian motion in a quasiperiodic potential

abstract: We consider a quantum particle subject to Ohmic dissipation, moving in a bichromatic quasiperiodic potential. In a periodic potential the particle undergoes a zero-temperature localization-delocalization transition as dissipation strength is decreased. We show that the delocalized phase is absent in the quasiperiodic case, even when the deviation from periodicity is infinitesimal. Using the renormalization group, we determine how the effective localization length depends on the dissipation. We show that {a similar problem can emerge in} the strong-coupling limit of a mobile impurity moving in a periodic lattice and immersed in a one-dimensional quantum gas.


title: A note on the eightfold way

abstract: Assuming the existence of a Mahlo cardinal, we construct a model in which there exists an $\omega_2$-Aronszajn tree, the $\omega_1$-approachability property fails, and every stationary subset of $\omega_2 \cap \mathrm{cof}(\omega)$ reflects.


title: The Cholewinski-Fock space in the Slice Hyperholomorphic Setting

abstract: Inspired from the Cholewinski approach see [5], we investigate a family of Fock spaces in the quaternionic slice hyperholomorphic setting as well as some associated quaternionic linear operators. In a particular case, we reobtain the slice hyperholomorphic Fock space introduced and studied in [2].


title: DASPS: A Database for Anxious States based on a Psychological Stimulation

abstract: Anxiety affects human capabilities and behavior as much as it affects productivity and quality of life. It can be considered as the main cause of depression and suicide. Anxious states are easily detectable by humans due to their acquired cognition, humans interpret the interlocutor's tone of speech, gesture, facial expressions and recognize their mental state. There is a need for non-invasive reliable techniques that performs the complex task of anxiety detection. In this paper, we present DASPS database containing recorded Electroencephalogram (EEG) signals of 23 participants during anxiety elicitation by means of face-to-face psychological stimuli. EEG signals were captured with Emotiv Epoc headset as it's a wireless wearable low-cost equipment. In our study, we investigate the impact of different parameters, notably: trial duration, feature type, feature combination and anxiety levels number. Our findings showed that anxiety is well elicited in 1 second. For instance, stacked sparse autoencoder with different type of features achieves 83.50% and 74.60% for 2 and 4 anxiety levels detection, respectively. The presented results prove the benefits of the use of a low-cost EEG headset instead of medical non-wireless devices and create a starting point for new researches in the field of anxiety detection.


title: Brain Electrical Stimulation for Animal Navigation

abstract: The brain stimulation and its widespread use is one of the most important subjects in studies of neurophysiology. In brain electrical stimulation methods, following the surgery and electrode implantation, electrodes send electrical impulses to the specific targets in the brain. The use of this stimulation method is provided therapeutic benefits for treatment chronic pain, essential tremor, Parkinsons disease, major depression, and neurological movement disorder syndrome (dystonia). One area in which advancements have been recently made is in controlling the movement and navigation of animals in a specific pathway. It is important to identify brain targets in order to stimulate appropriate brain regions for all the applications listed above. An animal navigation system based on brain electrical stimulation is used to develop new behavioral models for the aim of creating a platform for interacting with the animal nervous system in the spatial learning task. In the context of animal navigation the electrical stimulation has been used either as creating virtual sensation for movement guidance or virtual reward for movement motivation. In this paper, different approaches and techniques of brain electrical stimulation for this application has been reviewed. Keywords: Rat Robot, Brain Computer Interface, Electrical Stimulation, Cyborg Intelligence, Brain to Brain Interface


title: Maximally entangling states and dynamics in one dimensional nearest neighbor Floquet systems

abstract: We describe conditions for generating entanglement between two regions at the optimal rate in a class of one-dimensional quantum circuits with Floquet dynamics. The optimal value follows from subadditivity and Araki-Lieb inequalities. A quantum circuit composed of parallel SWAP gates that act periodically on entangled pairs is a simple system that saturates the bound. We show that any other system that entangles at this maximal rate must act as a generalized SWAP gate dynamics on the relevant states of the Hilbert space. We further discuss some characterizations of states according to entropy generation. States with multipartite entanglement generically fail to entangle efficiently as time evolves. This suggests that chaos, which tend to produce such entanglement patterns, is expected to work against the process of spreading information efficiently. It also provides a simple intuition for why the entangling tsunami velocity must be slower than the Lieb-Robinson velocity.


title: Algorithmic Bayesian Group Gibbs Selection

abstract: Bayesian model selection, with precedents in George and McCulloch (1993) and Abramovich et al. (1998), support credibility measures that relate model uncertainty, but computation can be costly when sparse priors are approximate. We design an exact selection engine suitable for Gauss noise, t-distributed noise, and logistic learning, benefiting from data-structures derived from coordinate descent lasso. Gibbs sampler chains are stored in a compressed binary format compatible with Equi-Energy (Kou et al., 2006) tempering. We achieve a grouped-effects selection model, similar to the setting for group lasso, to determine co-entry of coefficients into the model. We derive a functional integrand for group inclusion, and introduce a new MCMC switching step to avoid numerical integration. Theorems show this step has exponential convergence to target distribution. We demonstrate a role for group selection to inform on genetic decomposition in a diallel experiment, and identify potential quantitative trait loci in p > 40K Heterogenous Stock haplotype/phenotype studies.


title: Anapole Dark Matter Interactions as Soft Hidden Photons

abstract: We propose a model where the anapole appears as a hidden photon that is coupled to visible matter through a kinetic mixing. For low momentum $|{\bf p}| \ll M$ where $M$ is the cutoff the model (soft hidden photons limit) is reduced to the Ho-Scherrer description. We show that the hidden gauge boson is stable and therefore the hidden photons, indeed, are candidates for dark matter. Our approach shows that anapole and kinetic mixing terms are equivalent descriptions seen from different scales of energy.


title: An interval-valued GARCH model for range-measured return processes

abstract: Range-measured return contains more information than the traditional scalar-valued return. In this paper, we propose to model the [low, high] price range as a random interval and suggest an interval-valued GARCH (Int-GARCH) model for the corresponding range-measured return process. Under the general framework of random sets, the model properties are investigated. Parameters are estimated by the maximum likelihood method, and the asymptotic properties are established. Empirical application to stocks and financial indices data sets suggests that our Int-GARCH model overall outperforms the traditional GARCH for both in-sample estimation and out-of-sample prediction of volatility.


title: IC 361:- Near infrared and UBVRI photometric Analysis

abstract: We present here the detailed optical and infra-red photometric analysis of the open star cluster IC 361. On studying the radial density profile, radial extent of the cluster is found to be 8.0 +/- 0.5 arcmin. The basic physical parameters of the cluster such as E(B-V) = 0.56 +/- 0.10 mag, E(V-K) = 1.72+/-0.12 mag, log(Age)=9.10+/-0.05, and (m-M)0 = 12.54 +/-.05 mag are obtained using the color-color and color-magnitude diagrams. IC 361 is found to be located at a distance of 3.22 +/- 0.07 kpc. Using the archival proper motion catalogues, we estimate mean proper motions of IC 361 as 4.97+/-0.17 mas yr-1 and -5.80+/-0.18 mas yr-1 in the direction of RA and DEC, respectively. We derive the luminosity and mass functions for the cluster main sequence stars. The mass function slope is found to be -1.06+/-0.09 which is too low compare than Salpeter value.


title: A Bayesian Cognition Approach to Improve Data Visualization

abstract: People naturally bring their prior beliefs to bear on how they interpret the new information, yet few formal models exist for accounting for the influence of users' prior beliefs in interactions with data presentations like visualizations. We demonstrate a Bayesian cognitive model for understanding how people interpret visualizations in light of prior beliefs and show how this model provides a guide for improving visualization evaluation. In a first study, we show how applying a Bayesian cognition model to a simple visualization scenario indicates that people's judgments are consistent with a hypothesis that they are doing approximate Bayesian inference. In a second study, we evaluate how sensitive our observations of Bayesian behavior are to different techniques for eliciting people subjective distributions, and to different datasets. We find that people don't behave consistently with Bayesian predictions for large sample size datasets, and this difference cannot be explained by elicitation technique. In a final study, we show how normative Bayesian inference can be used as an evaluation framework for visualizations, including of uncertainty.


title: Spectral Approach to Verifying Non-linear Arithmetic Circuits

abstract: This paper presents a fast and effective computer algebraic method for analyzing and verifying non-linear integer arithmetic circuits using a novel algebraic spectral model. It introduces a concept of algebraic spectrum, a numerical form of polynomial expression; it uses the distribution of coefficients of the monomials to determine the type of arithmetic function under verification. In contrast to previous works, the proof of functional correctness is achieved by computing an algebraic spectrum combined with a local rewriting of word-level polynomials. The speedup is achieved by propagating coefficients through the circuit using And-Inverter Graph (AIG) datastructure. The effectiveness of the method is demonstrated with experiments including standard and Booth multipliers, and other synthesized non-linear arithmetic circuits up to 1024 bits containing over 12 million gates.


title: Factorization of composed polynomials and applications

abstract: Let $\mathbb{F}_q$ be the finite field with $q$ elements, where $q$ is a prime power and $n$ be a positive integer. In this paper, we explore the factorization of $f(x^{n})$ over $\mathbb{F}_q$, where $f(x)$ is an irreducible polynomial over $\mathbb F_q$. Our main results provide generalizations of recent works on the factorization of binomials $x^n-1$. As an application, we provide an explicit formula for the number of irreducible factors of $f(x^n)$ under some generic conditions on $f$ and $n$.


title: Two-dimensional nuclear magnetic resonance spectroscopy with a microfluidic diamond quantum sensor

abstract: Quantum sensors based on nitrogen-vacancy centers in diamond have emerged as a promising detection modality for nuclear magnetic resonance (NMR) spectroscopy owing to their micron-scale detection volume and non-inductive based detection. A remaining challenge is to realize sufficiently high spectral resolution and concentration sensitivity for multidimensional NMR analysis of picoliter sample volumes. Here, we address this challenge by spatially separating the polarization and detection phases of the experiment in a microfluidic platform. We realize a spectral resolution of 0.65 +/- 0.05 Hz, an order-of-magnitude improvement over previous diamond NMR studies. We use the platform to perform two-dimensional correlation spectroscopy of liquid analytes within an effective ~20 picoliter detection volume. The use of diamond quantum sensors as in-line microfluidic NMR detectors is a significant step towards applications in mass-limited chemical analysis and single cell biology.

categories: physics.ins-det cond-mat.mes-hall physics.app-ph physics.chem-ph


title: Black hole scalarisation from the breakdown of scale-invariance

abstract: Electro-vacuum black holes are scale-invariant; their energy-momentum tensor is traceless. Quantum corrections of various sorts, however, can often produce a trace anomaly and a breakdown of scale-invariance. The (quantum-corrected) black hole solutions of the corresponding gravitational effective field theory (EFT) have a non-vanishing Ricci scalar. Then, the presence of a scalar field with the standard non-minimal coupling $\xi \phi^2 R$ naturally triggers a spontaneous scalarisation of the corresponding black holes. This scalarisation phenomenon occurs for an (infinite) discrete set of $\xi$. We illustrate the occurrence of this phenomenon for two examples of static, spherically symmetric, asymptotically flat black hole solution of EFTs. In one example the trace anomaly comes from the matter sector -- a novel, closed form, generalisation of the Reissner-Nordstr\"om solution with an $F^4$ correction -- whereas in the other example it comes from the geometry sector -- a noncommutative geometry generalization of the Schwarzschild black hole. For comparison, we also consider the scalarisation of a black hole surrounded by (non-conformally invariant) classical matter (Einstein-Maxwell-dilaton black holes). We find that the scalarised solutions are, generically, entropically favoured.


title: Evolution of instabilities in filament buckling processes

abstract: In this work we study the dynamical buckling process of a thin filament immersed in a high viscous medium. We perform an experimental study to track the shape evolution of the filament during a constant velocity compression. Numerical simulations reproduce the dynamical features observed for the experimental data and allow quantifying the filament's load. We observe that both the filament's load and the wavenumber evolve in a step-wise manner. To achieve a physical insight of the process we apply a theoretical model to describe the buckling of a filament in a viscous medium. We solve a hydrodynamic equation in terms of normal modes for clamped-clamped boundary conditions and applied constant load. We find a good agreement with experimental data and simulations, suggesting that the proposed mechanistic model captures the essential features underlying the dynamical buckling process.


title: Simulating Non-hydrostatic atmospheres on Planets (SNAP): formulation, validation and application to the Jovian atmosphere

abstract: A new non-hydrostatic and cloud-resolving atmospheric model is developed for studying moist convection and cloud formation in planetary atmospheres. It is built on top of the Athena++ framework, utilizing its static/adaptive mesh-refinement, parallelization, curvilinear geometry, and dynamic task scheduling. We extend the original hydrodynamic solver to vapors, clouds, and precipitation. Microphysics is formulated generically so that it can be applied to both Earth and Jovian planets. We implemented the Low Mach number Approximate Riemann Solver (LMARS) for simulating low speed atmospheric flows in addition to the usual Roe and HLLC Riemann solvers. Coupled with a fifth-order Weighted Essentially Nonoscillatory (WENO) subgrid-reconstruction method, the sharpness of critical fields such as clouds is well-preserved, and no extra hyperviscosity or spatial filter is needed to stabilize the model. Unlike many atmospheric models, total energy is used as the prognostic variable of the thermodynamic equation. One significant advantage of using total energy as a prognostic variable is that the entropy production due to irreversible mixing process can be properly captured. The model is designed to provide a unified framework for exploring planetary atmospheres across various conditions, both terrestrial and Jovian. First, a series of standard numerical tests for Earth's atmosphere is carried out to demonstrate the performance and robustness of the new model. Second, simulation of an idealized Jovian atmosphere in radiative-convective equilibrium shows that 1) the temperature gradient is superadiabatic near the water condensation level because of the changing of the mean molecular weight, and 2) the mean profile of ammonia gas shows a depletion in the subcloud layer down to nearly 10 bars. Relevance to the recent Juno observations is discussed.


title: The Bishop--Phelps--Bollob\'as property for Lipschitz maps

abstract: In this paper, we introduce and study a Lipschitz version of the Bishop-Phelps-Bollob\'as property (Lip-BPB property). This property deals with the possibility of making a uniformly simultaneous approximation of a Lipschitz map $F$ and a pair of points at which $F$ almost attains its norm by a Lipschitz map $G$ and a pair of points such that $G$ strongly attains its norm at the new pair of points. We first show that if $M$ is a finite pointed metric space and $Y$ is a finite-dimensional Banach space, then the pair $(M,Y)$ has the Lip-BPB property, and that both finiteness assumptions are needed. Next, we show that if $M$ is a uniformly Gromov concave pointed metric space (i.e.\ the molecules of $M$ form a set of uniformly strongly exposed points), then $(M,Y)$ has the Lip-BPB property for every Banach space $Y$. We further prove that this is the case for finite concave metric spaces, ultrametric spaces, and H\"older metric spaces. The extension of the Lip-BPB property from $(M,\mathbb R)$ to some Banach spaces $Y$ and some results for compact Lipschitz maps are also discussed.


title: Decision-Making Under Uncertainty in Research Synthesis: Designing for the Garden of Forking Paths

abstract: To make evidence-based recommendations to decision-makers, researchers conducting systematic reviews and meta-analyses must navigate a garden of forking paths: a series of analytical decision-points, each of which has the potential to influence findings. To identify challenges and opportunities related to designing systems to help researchers manage uncertainty around which of multiple analyses is best, we interviewed 11 professional researchers who conduct research synthesis to inform decision-making within three organizations. We conducted a qualitative analysis identifying 480 analytical decisions made by researchers throughout the scientific process. We present descriptions of current practices in applied research synthesis and corresponding design challenges: making it more feasible for researchers to try and compare analyses, shifting researchers' attention from rationales for decisions to impacts on results, and supporting communication techniques that acknowledge decision-makers' aversions to uncertainty. We identify opportunities to design systems which help researchers explore, reason about, and communicate uncertainty in decision-making about possible analyses in research synthesis.


title: Asymptotic estimates of solutions to time-fractional diffusion equations with space-dependent variable order

abstract: We examine the short and long-time behaviors of time-fractional diffusion equations with variable space-dependent order. More precisely, we describe the time-evolution of the solution to these equations as the time parameter goes either to zero or to infinity.


title: Disentangling the spatial substructure of Cygnus OB2 from Gaia DR2

abstract: For the first time, we have explored the spatial substructure of the Cygnus OB2 association using parallaxes from the recent second Gaia data release. We find significant line-of-sight substructure within the association, which we quantify using a parameterised model that reproduces the observed parallax distribution. This inference approach is necessary due to the non-linearity of the parallax-distance transformation and the asymmetry of the resulting probability distribution. By using a Markov Chain Monte Carlo ensemble sampler and an unbinned maximum likelihood test we identify two different stellar groups superposed on the association. We find the main Cygnus OB2 group at ~1760 pc, further away than recent estimates have envisaged, and a foreground group at ~1350 pc. We also calculate individual membership probabilities and identify outliers as possible non-members of the association.


title: TOROS Optical follow-up of the Advanced LIGO-VIRGO O2 second observational campaign

abstract: We present the results of the optical follow-up, conducted by the TOROS collaboration, of gravitational wave events detected during the Advanced LIGO-Virgo second observing run (Nov 2016 -- Aug 2017). Given the limited field of view ($\sim100\arcmin$) of our observational instrumentation we targeted galaxies within the area of high localization probability that were observable from our sites. We analyzed the observations using difference imaging, followed by a Random Forest algorithm to discriminate between real and bogus transients. For all three events that we respond to, except GW170817, we did not find any bona fide optical transient that was plausibly linked with the observed gravitational wave event. Our observations were conducted using telescopes at Estaci\'{o}n Astrof\'{\i}sica de Bosque Alegre, Cerro Tololo Inter-American Observatory, and the Dr. Cristina V. Torres Memorial Astronomical Observatory. Our results are consistent with the LIGO-Virgo detections of a binary black hole merger (GW170104) for which no electromagnetic counterparts were expected, as well as a binary neutron star merger (GW170817) for which an optical transient was found as expected.


title: Thermal Conductance across beta-Ga2O3-diamond Van der Waals Heterogeneous Interfaces

abstract: Because of its ultrawide bandgap, high breakdown electric field, and large area affordable substrates grown from the melt, beta Ga2O3 has attracted great attention recently for potential applications of power electronics. However, its thermal conductivity is significantly lower than those of other wide bandgap semiconductors, such as AlN, SiC, GaN, and diamond. To ensure reliable operation with minimal selfheating at high power, proper thermal management is even more essential for Ga2O3 devices. Similarly to the past approaches aiming to alleviate selfheating in GaN HEMTs, a possible solution has been to integrate thin Ga2O3 membranes with diamond to fabricate Ga2O3 on diamond lateral MESFET or MOSFET devices by taking advantage of the ultra high thermal conductivity of diamond. Even though the TBC between wide bandgap semiconductor devices such as GaN HEMTs and a diamond substrate is of primary importance for heat dissipation in these devices, fundamental understanding of the Ga2O3 diamond thermal interface is still missing. In this work, we study the thermal transport across the interfaces of Ga2O3 exfoliated onto a single crystal diamond. The Van der Waals bonded Ga2O3 diamond TBC is measured to be 17 MWm2K1, which is comparable to the TBC of several physical vapor deposited metals on diamond. A Landauer approach is used to help understand phonon transport across perfect Ga2O3 diamond interface, which in turn sheds light on the possible TBC one could achieve with an optimized interface. A reduced thermal conductivity of the Ga2O3 nanomembrane is also observed due to additional phonon membrane boundary scattering. The impact of the Ga2O3substrate TBC and substrate thermal conductivity on the thermal performance of a power device are modeled and discussed.


title: CLOUDY view of the warm corona

abstract: Bright active galaxies show a range of properties but many of these properties are correlated which has led to the concept of the Quasar Main Sequence. We test whether our current understanding of the quasar structure allows to reproduce the pattern observed in the optical plane formed by the kinematic line width of H$\beta$ and the relative importance of the Fe II optical emission. We performed simulations of the H$\beta$ and Fe II production using the code CLOUDY and well justified assumptions about the broad band spectra, distance of the emission line region, and the cloud properties. We show that the presence of the warm corona is an important element of the broad band spectrum which decreases the dependence of the relative Fe II emissivity on the Eddington ratio, and allows to reproduce the rare cases of the particularly strong Fe II emitters. Results are sensitive to the adopted cloud distance, and strong Fe II emission can be obtain either by adopting strongly super-solar metallicity, or much shorter distance than traditionally obtained from reverberation mapping. We modeled in a similar way the UV plane defined by the Mg II line and Fe II UV pseudo-continuum, but here our approach is less successful, in general overproducing the Fe II strength. We found that the Fe II optical and UV emissivity depend in a different way on the turbulent velocity and metallicity, and the best extension of the model in order to cover both planes is to allow very large turbulent velocities in the Broad Line Region clouds.


title: The effects of competition between random sequential nucleation of point-sized seeds and island growth by adsorption of finite-sized grains

abstract: We study random sequential adsorption of particles from pool onto a one dimensional substrate following ballistic deposition rules, with separate nucleation and growth processes occurring simultaneously. Nucleation describes the formation of point-sized seeds, and after a seed is sown, it acts as an attractor and grows in size by the addition of grains of a fixed-sized. At each time step either an already-nucleated seed can increase in size, or a new seed may be nucleated. We incorporate a parameter $m$, to describe the relative rates of growth to nucleation. We solve the model analytically to obtain gap size distribution function and a general expression for the jamming coverage as a function of $m$. We show that the jamming coverage $\theta(m)$ reaches its maximum value $\theta(m)=1$ in the limit $m\rightarrow \infty$ following a power-law $\theta(\infty) - \theta(m) \sim m^{-1/2}$. We also perform extensive Monte Carlo simulation and find excellent agreement between analytic and numerical results.


title: On global stability of optimal rearrangement maps

abstract: We study the nonlocal vectorial transport equation $\partial_ty+ (\mathbb{P} y \cdot \nabla) y=0$ on bounded domains of $\mathbb{R}^d$ where $\mathbb{P}$ denotes the Leray projector. This equation was introduced to obtain the unique optimal rearrangement of a given map $y_0$ as the infinite time limit of the solution with initial data $y_0$ (\cite{AHT, Macthesis, Brenier09}). We rigorously justify this expectation by proving that for initial maps $y_0$ sufficiently close to maps with strictly convex potential, the solutions $y$ are global in time and converge exponentially fast to the optimal rearrangement of $y_0$ as time tends to infinity.


title: The regular black hole in four dimensional Born-Infeld gravity

abstract: In the context of Born-Infeld gravity theories we report the existence of a regular black hole interior representing a spherically symmetric vacuum solution of the theory. It reduces to the Schwarzschild interior metric in the weak field region. In particular, there is a new length scale which is related to the Born-Infeld parameter. This endows the spacetime with an inner (i.e. well inside the event horizon) asymptotic region which is unattainable for observers. The central curvature singularity is replaced by an infinitely long cosmic string with constant curvature invariants related to the Born-Infeld constant. The presence of this limiting curvature spacetime renders the black hole timelike and null geodesically complete, free from the classical Schwarzschild singularity. The transition between the usual black hole interior and this maximum curvature space is achieved without introducing any kind of matter content nor topological changes.


title: Stone-Wales graphene: A Two Dimensional Carbon Semi-Metal with Magic Stability

abstract: A two-dimensional carbon allotrope, Stone-Wales graphene, is identified in stochastic group and graph constrained searches and systematically investigated by first-principles calculations. Stone-Wales graphene consists of well-arranged Stone-Wales defects, and it can be constructed through a 90$^\circ$ bond-rotation in a $\sqrt{8}$$\times$$\sqrt{8}$ super-cell of graphene. Its calculated energy relative to graphene, +149 meV/atom, makes it more stable than the most competitive previously suggested graphene allotropes. We find that Stone-Wales graphene based on a $\sqrt{8}$ super-cell is more stable than those based on $\sqrt{9} \times \sqrt{9}$, $\sqrt{12} \times \sqrt{12}$ and $\sqrt{13} \times \sqrt{13}$ super-cells, and is a "magic size" that can be further understood through a simple "energy splitting and inversion" model. The calculated vibrational properties and molecular dynamics of SW-graphene confirm that it is dynamically stable. The electronic structure shows SW-graphene is a semimetal with distorted, strongly anisotropic Dirac cones.


title: Elpasolite Planetary Ice and Composition Spectrometer (EPICS): A Low-Resource Combined Gamma-Ray and Neutron Spectrometer for Planetary Science

abstract: Neutron and gamma-ray spectroscopy (NGRS) is a well established technique for studying the geochemical composition and volatile abundance relevant to planetary structure and evolution of planetary bodies. Previous NGRS instruments have used separate gamma-ray and neutron spectrometers. The Elpasolite Planetary Ice and Composition Spectrometer (EPICS) instrument is an innovative and fully integrated NGRS with low resource requirements. EPICS utilizes elpasolite scintillator read out by silicon photomultipliers to combine the gamma-ray and neutron spectrometer into a single instrument, leading to a significant reduction in instrument size, weight, and power. An overview and motivation for the EPICS instrument, current status of the EPICS development, and a discussion of the expected sensitivity and performance are presented.


title: Composite Shape Modeling via Latent Space Factorization

abstract: We present a novel neural network architecture, termed Decomposer-Composer, for semantic structure-aware 3D shape modeling. Our method utilizes an auto-encoder-based pipeline, and produces a novel factorized shape embedding space, where the semantic structure of the shape collection translates into a data-dependent sub-space factorization, and where shape composition and decomposition become simple linear operations on the embedding coordinates. We further propose to model shape assembly using an explicit learned part deformation module, which utilizes a 3D spatial transformer network to perform an in-network volumetric grid deformation, and which allows us to train the whole system end-to-end. The resulting network allows us to perform part-level shape manipulation, unattainable by existing approaches. Our extensive ablation study, comparison to baseline methods and qualitative analysis demonstrate the improved performance of the proposed method.


title: $L^2$-type contraction for shocks of scalar viscous conservation laws with strictly convex flux

abstract: We study the $L^2$-type contraction property of large perturbations around shock waves of scalar viscous conservation laws with strictly convex fluxes in one space dimension. The contraction holds up to a shift, and it is measured by a weighted related entropy, for which we choose an appropriate entropy associated with the strictly convex flux. In particular, we handle shocks with small amplitude. This result improves the recent article [18] of the author and Vasseur on $L^2$-contraction property of shocks to scalar viscous conservation laws with a special flux, that is almost the Burgers flux.


title: Normalized Object Coordinate Space for Category-Level 6D Object Pose and Size Estimation

abstract: The goal of this paper is to estimate the 6D pose and dimensions of unseen object instances in an RGB-D image. Contrary to "instance-level" 6D pose estimation tasks, our problem assumes that no exact object CAD models are available during either training or testing time. To handle different and unseen object instances in a given category, we introduce a Normalized Object Coordinate Space (NOCS)---a shared canonical representation for all possible object instances within a category. Our region-based neural network is then trained to directly infer the correspondence from observed pixels to this shared object representation (NOCS) along with other object information such as class label and instance mask. These predictions can be combined with the depth map to jointly estimate the metric 6D pose and dimensions of multiple objects in a cluttered scene. To train our network, we present a new context-aware technique to generate large amounts of fully annotated mixed reality data. To further improve our model and evaluate its performance on real data, we also provide a fully annotated real-world dataset with large environment and instance variation. Extensive experiments demonstrate that the proposed method is able to robustly estimate the pose and size of unseen object instances in real environments while also achieving state-of-the-art performance on standard 6D pose estimation benchmarks.


title: An Algebraic Sparsified Nested Dissection Algorithm Using Low-Rank Approximations

abstract: We propose a new algorithm for the fast solution of large, sparse, symmetric positive-definite linear systems, spaND -- sparsified Nested Dissection. It is based on nested dissection, sparsification and low-rank compression. After eliminating all interiors at a given level of the elimination tree, the algorithm sparsifies all separators corresponding to the interiors. This operation reduces the size of the separators by eliminating some degrees of freedom but without introducing any fill-in. This is done at the expense of a small and controllable approximation error. The result is an approximate factorization that can be used as an efficient preconditioner. We then perform several numerical experiments to evaluate this algorithm. We demonstrate that a version using orthogonal factorization and block-diagonal scaling takes fewer CG iterations to converge than previous similar algorithms on various kinds of problems. Furthermore, this algorithm is provably guaranteed to never break down and the matrix stays symmetric positive-definite throughout the process. We evaluate the algorithm on some large problems and show it exhibits near-linear scaling. The factorization time is roughly O(N) and the number of iterations grows slowly with N.


title: A sequential update algorithm for computing the stationary distribution vector in upper block-Hessenberg Markov chains

abstract: This paper proposes a new algorithm for computing the stationary distribution vector in continuous-time upper block-Hessenberg Markov chains. To this end, we consider the last-block-column-linearly-augmented (LBCL-augmented) truncation of the (infinitesimal) generator of the upper block-Hessenberg Markov chain. The LBCL-augmented truncation is a linearly-augmented truncation such that the augmentation distribution has its probability mass only on the last block column. We first derive an upper bound for the total variation distance between the respective stationary distribution vectors of the original generator and its LBCL-augmented truncation. Based on the upper bound, we then establish a series of linear fractional programming (LFP) problems to obtain augmentation distribution vectors such that the bound converges to zero. Using the optimal solutions of the LFP problems, we construct a matrix-infinite-product (MIP) form of the original (i.e., not approximate) stationary distribution vector and develop a sequential update algorithm for computing the MIP form. Finally, we demonstrate the applicability of our algorithm to BMAP/M/$\infty$ queues and M/M/$s$ retrial queues.


title: Existence of a unique solution and invariant measures for the stochastic Landau--Lifshitz--Bloch equation

abstract: The Landau--Lifshitz--Bloch equation perturbed by a space-dependent noise was proposed in Garanin 1991 as a model for evolution of spins in ferromagnatic materials at the full range of temperatures, including the temperatures higher than the Curie temperature. In the case of a ferromagnet filling a bounded domain $D\subset \mathbb R^d$, $d=1,2,3$, we show the existence of strong (in the sense of PDEs) martingale solutions. Furthermore, in cases $d=1,2$ we prove uniqueness of pathwise solutions and the existence of invariant measures.


title: Complex oscillatory patterns near singular Hopf bifurcation in a two time-scale ecosystem

abstract: We consider an ecological model consisting of two species of predators competing for their common prey with explicit interference competition. With a proper rescaling, the model is portrayed as a singularly perturbed system with one-fast (prey dynamics) and two-slow variables (dynamics of the predators). The model exhibits variety of rich and interesting dynamics, including, but not limited to mixed mode oscillations (MMOs), featuring concatenation of small and large amplitude oscillations, relaxation oscillations and bistability between a semi-trivial equilibrium state and a coexistence oscillatory state. Existence of co-dimenison two bifurcations such as fold-Hopf and generalized Hopf bifurcations make the system further intriguing. More interestingly, in a neighborhood of {\emph{singular Hopf}} bifurcation, long lasting transient dynamics in form of chaotic MMOs or relaxation oscillations are observed as the system approaches the periodic attractor born out of supercritical Hopf bifurcation or a semi-trivial equilibrium state respectively. The transient dynamics could persist for hundreds or thousands of generations before the ecosystem experiences a regime shift. The time series of population cycles with different types of irregular oscillations arising in this model stem from a biological realistic feature, namely, by the variation in the intraspecific competition amongst the predators. To explain these oscillations, we use bifurcation analysis and methods from {\emph{geometric singular perturbation theory}}.


title: A witness function based construction of discriminative models using Hermite polynomials

abstract: In machine learning, we are given a dataset of the form $\{(\mathbf{x}_j,y_j)\}_{j=1}^M$, drawn as i.i.d. samples from an unknown probability distribution $\mu$; the marginal distribution for the $\mathbf{x}_j$'s being $\mu^*$. We propose that rather than using a positive kernel such as the Gaussian for estimation of these measures, using a non-positive kernel that preserves a large number of moments of these measures yields an optimal approximation. We use multi-variate Hermite polynomials for this purpose, and prove optimal and local approximation results in a supremum norm in a probabilistic sense. Together with a permutation test developed with the same kernel, we prove that the kernel estimator serves as a `witness function' in classification problems. Thus, if the value of this estimator at a point $\mathbf{x}$ exceeds a certain threshold, then the point is reliably in a certain class. This approach can be used to modify pretrained algorithms, such as neural networks or nonlinear dimension reduction techniques, to identify in-class vs out-of-class regions for the purposes of generative models, classification uncertainty, or finding robust centroids. This fact is demonstrated in a number of real world data sets including MNIST, CIFAR10, Science News documents, and LaLonde data sets.


title: The square root rule for adaptive importance sampling

abstract: In adaptive importance sampling, and other contexts, we have $K>1$ unbiased and uncorrelated estimates $\hat\mu_k$ of a common quantity $\mu$. The optimal unbiased linear combination weights them inversely to their variances but those weights are unknown and hard to estimate. A simple deterministic square root rule based on a working model that $\mathrm{Var}(\hat\mu_k)\propto k^{-1/2}$ gives an unbisaed estimate of $\mu$ that is nearly optimal under a wide range of alternative variance patterns. We show that if $\mathrm{Var}(\hat\mu_k)\propto k^{-y}$ for an unknown rate parameter $y\in [0,1]$ then the square root rule yields the optimal variance rate with a constant that is too large by at most $9/8$ for any $0\le y\le 1$ and any number $K$ of estimates. Numerical work shows that rule is similarly robust to some other patterns with mildly decreasing variance as $k$ increases.


title: Sharing Energy Storage Between Transmission and Distribution

abstract: This paper addresses the problem of how best to coordinate, or `stack,' energy storage services in systems that lack centralized markets. Specifically, its focus is on how to coordinate transmission-level congestion relief with local, distribution-level objectives. We describe and demonstrate a unified communication and optimization framework for performing this coordination. The congestion relief problem formulation employs a weighted $\ell_{1}$-norm objective. This approach determines a set of corrective actions, i.e., energy storage injections and conventional generation adjustments, that minimize the required deviations from a planned schedule. To exercise this coordination framework, we present two case studies. The first is based on a 3-bus test system, and the second on a realistic representation of the Pacific Northwest region of the United States. The results indicate that the scheduling methodology provides congestion relief, cost savings, and improved renewable energy integration. The large-scale case study informed the design of a live demonstration carried out in partnership with the University of Washington, Doosan GridTech, Snohomish County PUD, and the Bonneville Power Administration. The goal of the demonstration was to test the feasibility of the scheduling framework in a production environment with real-world energy storage assets. The demonstration results were consistent with computational simulations.


title: A Truthful FPTAS Mechanism for Emergency Demand Response in Colocation Data Centers

abstract: Demand response (DR) is not only a crucial solution to the demand side management but also a vital means of electricity market in maintaining power grid reliability, sustainability and stability. DR can enable consumers (e.g. data centers) to reduce their electricity consumption when the supply of electricity is a shortage. The consumers will be rewarded in the case of DR if they reduce or shift some of their energy usage during peak hours. Aiming at solving the efficiency of DR, in this paper, we present MEDR, a mechanism on emergency DR in colocation data center. First, we formalize the MEDR problem and propose a dynamic programming to solve the optimization version of the problem. We then design a deterministic mechanism as a solution to solve the MEDR problem. We show that our proposed mechanism is truthful. Next, we prove that our mechanism is an FPTAS, i.e., it can be approximated within $1 + \epsilon$ for any given $\epsilon > 0$, while the running time of our mechanism is polynomial in $n$ and $1/\epsilon$, where $n$ is the number of tenants in the datacenter. Furthermore, we also give an auction system covering the efficient FPTAS algorithm as bidding decision program for DR in colocation datacenter. Finally, we choose a practical smart grid dataset to build a large number of datasets for simulation in performance evaluation. By evaluating metrics of the approximation ratio of our mechanism, the non-negative utility of tenants and social cost of colocation datacenter, the results demonstrate the effectiveness of our work.


title: Active Willis metamaterials for ultra-compact non-reciprocal linear acoustic devices

abstract: Willis materials are complex media characterized by four macroscopic material parameters, the conventional mass density, and bulk modulus and two additional Willis coupling terms, which have been shown to enable unsurpassed control over the propagation of mechanical waves. However, virtually all previous studies on Willis materials involved passive structures which have been shown to have limitations in terms of achievable Willis coupling terms. In this article, we show experimentally that linear active Willis metamaterials breaking these constraints enable highly non-reciprocal sound transport in very subwavelength structures, a feature unachievable through other methods. Furthermore, we present an experimental procedure to extract the effective material parameters expressed in terms of acoustic polarizabilities for media in which the Willis coupling terms are allowed to vary independently. The approach presented here will enable a new generation of Willis materials for enhanced sound control and improved acoustic imaging and signal processing.


title: The FHD/$\boldsymbol{\varepsilon}$ppsilon Epoch of Reionization Power Spectrum Pipeline

abstract: Epoch of Reionization data analysis requires unprecedented levels of accuracy in radio interferometer pipelines. We have developed an imaging power spectrum analysis to meet these requirements and generate robust 21 cm EoR measurements. In this work, we build a signal path framework to mathematically describe each step in the analysis, from data reduction in the FHD package to power spectrum generation in the $\varepsilon$ppsilon package. In particular, we focus on the distinguishing characteristics of FHD/$\varepsilon$ppsilon: highly accurate spectral calibration, extensive data verification products, and end-to-end error propagation. We present our key data analysis products in detail to facilitate understanding of the prominent systematics in image-based power spectrum analyses. As a verification to our analysis, we also highlight a full-pipeline analysis simulation to demonstrate signal preservation and lack of signal loss. This careful treatment ensures that the FHD/$\varepsilon$ppsilon power spectrum pipeline can reduce radio interferometric data to produce credible 21 cm EoR measurements.