title: Regularity estimates for fully nonlinear elliptic PDEs with general Hamiltonian terms and unbounded ingredients

abstract: We develop an optimal regularity theory for $L^p$-viscosity solutions of fully nonlinear uniformly elliptic equations in nondivergence form whose gradient growth is described through a Hamiltonian function with measurable and possibly unbounded coefficients. Our approach treats both superlinear and sublinear gradient regimes in a unified way. We show $C^{0,\alpha}$, $C^{0,\textrm{Log-Lip}}$, $C^{1,\alpha}$, $ C^{1,\textrm{Log-Lip}}$ and $C^{2,\alpha}$ regularity estimates, by displaying the growth allowed to the Hamiltonian in order to deal with an unbounded nonlinear gradient coefficient, whose integrability in turn gets worse as we approach the quadratic regime. Moreover, we find proper compatibility conditions for which our regularity results depend intrinsically on the integrability of the underlying source term. As a byproduct of our methods, we prove a priori BMO estimates; sharp regularity to associated recession and flat profiles under relaxed convexity assumptions; improved regularity for a class of singular PDEs; and a Perron type result under unbounded ingredients.


title: Deciding whether two quadratic surfaces actually intersect

abstract: We answer an open question proposed by P\'{o}lik and Terlaky in 2007 that: {\it how we can decide whether two quadratic surfaces intersect without actually computing the intersections?}


title: Direct Measurement of the Cosmic-Ray Carbon and Oxygen Spectra from 10 GeV$/n$ to 2.2 TeV$/n$ with the Calorimetric Electron Telescope on the International Space Station

abstract: In this paper, we present the measurement of the energy spectra of carbon and oxygen in cosmic rays based on observations with the Calorimetric Electron Telescope (CALET) on the International Space Station from October 2015 to October 2019. Analysis, including the detailed assessment of systematic uncertainties, and results are reported. The energy spectra are measured in kinetic energy per nucleon from 10 GeV$/n$ to 2.2 TeV$/n$ with an all-calorimetric instrument with a total thickness corresponding to 1.3 nuclear interaction length. The observed carbon and oxygen fluxes show a spectral index change of $\sim$0.15 around 200 GeV$/n$ established with a significance $>3\sigma$. They have the same energy dependence with a constant C/O flux ratio $0.911\pm 0.006$ above 25 GeV$/n$. The spectral hardening is consistent with that measured by AMS-02, but the absolute normalization of the flux is about 27% lower, though in agreement with observations from previous experiments including the PAMELA spectrometer and the calorimetric balloon-borne experiment CREAM.


title: Local Dvoretzky-Kiefer-Wolfowitz confidence bands

abstract: In this paper, we revisit the concentration inequalities for the supremum of the cumulative distribution function (CDF) of a real-valued continuous distribution as established by Dvoretzky, Kiefer, Wolfowitz and revisited later by Massart in two seminal papers. We focus on the concentration of the \emph{local} supremum over a sub-interval, rather than on the full domain. That is, denoting $U$ the CDF of the uniform distribution over $[0,1]$ and $U_n$ its empirical version built from $n$ samples, we study $P(\sup_{u\in[\underline{u},\overline{u}]}U_n(u)-U(u)>\epsilon)$ for different values of $\underline{u},\overline{u}\in[0,1]$. Such local controls naturally appear for instance when studying estimation error of spectral risk-measures (such as the conditional value at risk), where $[\underline{u},\overline{u}]$ is typically $[0,\alpha]$ or $[1-\alpha,1]$ for a risk level $\alpha$, after reshaping the CDF $F$ of the considered distribution into $U$ by the general inverse transform $F^{-1}$. Extending a proof technique from Smirnov, we provide exact expressions of the local quantities $P(\sup_{u\in[\underline{u},\overline{u}]}U_n(u)-U(u)>\epsilon)$ and $P(\sup_{u\in [\underline{u},\overline{u}]}U(u)-U_n(u)>\epsilon)$ for each $n,\epsilon,\underline{u},\overline{u}$. Interestingly these quantities, seen as a function of $\epsilon$, can be easily inverted numerically into functions of the probability level $\delta$. Although not explicit, they can be computed and tabulated. We plot such expressions and compare them to the classical bound $\sqrt{\frac{\ln(1/\delta)}{2n}}$ provided by Massart inequality. Last, we extend the local concentration results holding individually for each $n$ to time-uniform concentration inequalities holding simultaneously for all $n$, revisiting a reflection inequality by James, which is of independent interest for the study of sequential decision making strategies.


title: Moments and saturation properties of eigenstates

abstract: Eigenvalues are defined for any element of an algebra of observables and do not require a representation in terms of wave functions or density matrices. A systematic algebraic derivation based on moments is presented here for the harmonic oscillator, together with a perturbative treatment of anharmonic systems. In this process, a collection of inequalities is uncovered which amount to uncertainty relations for higher-order moments saturated by the harmonic-oscillator excited states. Similar saturation properties hold for anharmonic systems order by order in perturbation theory. The new method, based on recurrence relations for moments of a state combined with positivity conditions, is therefore able to show new physical features.


title: Resonant Shattering Flares as Multimessenger Probes of the Nuclear Symmetry Energy

abstract: The behaviour of the nuclear symmetry energy near saturation density is important for our understanding of dense nuclear matter. This density dependence can be parameterised by the nuclear symmetry energy and its derivatives evaluated at nuclear saturation density. In this work we show that the core-crust interface mode of a neutron star is sensitive to these parameters, through the (density-weighted) shear-speed within the crust, which is in turn dependent on the symmetry energy profile of dense matter. We calculate the frequency at which the neutron star quadrupole ($\ell = 2$) crust-core interface mode must be driven by the tidal field of its binary partner to trigger a Resonant Shattering Flare (RSF). We demonstrate that coincident multimessenger timing of an RSF and gravitational wave chirp from a neutron star merger would enable us to place constraints on the symmetry energy parameters that are competitive with those from current nuclear experiments.


title: Minimal generating sets for matrix monoids

abstract: In this paper, we determine minimal generating sets for several well-known monoids of matrices over semirings. In particular, we find minimal generating sets for the monoids consisting of: all $n\times n$ boolean matrices when $n\leq 8$; the $n\times n$ boolean matrices containing the identity matrix (the reflexive boolean matrices) when $n\leq 7$; the $n\times n$ boolean matrices containing a permutation (the Hall matrices) when $n \leq 8$; the upper, and lower, triangular boolean matrices of every dimension; the $2 \times 2$ matrices over the semiring $\mathbb{N} \cup \{-\infty\}$ with addition $\oplus$ defined by $x\oplus y = \max(x, y)$ and multiplication $\otimes$ given by $x\otimes y = x + y$ (the max-plus semiring); the $2\times 2$ matrices over any quotient of the max-plus semiring by the congruence generated by $t = t + 1$ where $t\in \mathbb{N}$; the $2\times 2$ matrices over the min-plus semiring and its finite quotients by the congruences generated by $t = t + 1$ for all $t\in \mathbb{N}$; and the $n \times n$ matrices over $\mathbb{Z} / n\mathbb{Z}$ relative to their group of units.


title: Sources of $H_0$-tensions in dark energy scenarios

abstract: By focusing on the simple $w\neq-1$ extension to $\Lambda$CDM, we assess which epoch(s) possibly source the $H_0$-tension. We consider Cosmic Microwave Background (CMB) data in three possible ways: $i)$ complete CMB data; $ii)$ excluding the $l<30$ temperature and polarization likelihoods; $iii)$ imposing early universe priors, that disentangle early and late time physics. Through a joint analysis with low-redshift supernovae type-Ia and gravitationally lensed time delay datasets, {and neglecting galaxy clustering Baryonic Acoustic Oscillation (BAO) data}, we find that the inclusion of early universe CMB priors is consistent with the local estimate of $H_0$ while excluding the low-$l$+lowE likelihoods mildly relaxes the tension. This is in contrast to joint analyses with the complete CMB data. Our simple implementation of contrasting the effect of different CMB priors on the $H_0$ estimate shows that the early universe information from the CMB data when decoupled from late-times physics could be in agreement with a higher value of $H_0$. {We also find no evidence for the early dark energy model using only the early universe physics within the CMB data. Finally using the BAO data in different redshift ranges to perform inverse distance ladder analysis, we find that the early universe modifications, while being perfectly capable of alleviating the $H_0$-tension when including the BAO galaxy clustering data, would be at odds with the Ly-$\alpha$ BAO data due to the difference in $r_{\rm d}\, vs.\, H_0$ correlation between the two BAO datasets.} We therefore infer and speculate that source for the $H_0$-tension between CMB and local estimates could possibly originate in the modeling of late-time physics within the CMB analysis. This in turn recasts the $H_0$-tension as an effect of late-time physics in CMB, instead of the current early-time CMB vs. local late-time physics perspective.


title: Exploring instantons with spin-lattice systems

abstract: Instanton processes are present in a variety of quantum field theories relevant to high energy as well as condensed matter physics. While they have led to important theoretical insights and physical applications, their underlying features often remain elusive due to the complicated computational treatment. Here, we address this problem by studying topological as well as non-topological instantons using Monte Carlo methods on lattices of interacting spins. As a proof of principle, we systematically construct instanton solutions in $O(3)$ non-linear sigma models with a Dzyaloshinskii-Moriya interaction in $(1 + 1)$ and $(1 + 2)$ dimensions, thereby resembling an example of a chiral magnet. We demonstrate that, due to their close correspondence, Monte Carlo techniques in spin-lattice systems are well suited to describe topologically non-trivial field configurations in these theories. In particular, by means of simulated annealing, we demonstrate how to obtain domain walls, merons and critical instanton solutions.


title: Utilizing a Fully Optical and Reconfigurable PUF as a Quantum Authentication Mechanism

abstract: In this work the novel usage of a physically unclonable function composed of a network of Mach-Zehnder interferometers for authentication tasks is described. The physically unclonable function hardware is completely reconfigurable, allowing for a large number of seemingly independent devices to be utilized, thus imitating a large array of single-response physically unclonable functions. It is proposed that any reconfigurable array of Mach-Zehnder interferometers can be used as an authentication mechanism, not only for physical objects, but for information transmitted both classically and quantumly. The proposed use-case for a fully-optical physically unclonable function, designed with reconfigurable hardware, is to authenticate messages between a trusted and possibly untrusted party; verifying that the messages received are generated by the holder of the authentic device.


title: Solving a new type of quadratic optimization problem having a joint numerical range constraint

abstract: We propose a new formulation of quadratic optimization problems. The objective function $F(f(x),g(x))$ is given as composition of a quadratic function $F(z)$ with two $n$-variate quadratic functions $z_1=f(x)$ and $z_2=g(x).$ In addition, it incorporates with a set of linear inequality constraints in $z=(z_1,z_2)^T,$ while having an implicit constraint that $z$ belongs to the joint numerical range of $(f,g).$ The formulation is very general in the sense that it covers quadratic programming with a single quadratic constraint of all types, including the inequality-type, the equality-type, and the interval-type. Even more, the composition of "quadratic with quadratics" as well as the joint numerical range constraint all together allow us to formulate existing unsolved (or not solved efficiently) problems into the new model. In this paper, we solve the quadratic hypersurfaces intersection problem (QSIC) proposed by P$\acute{{\rm o}}$lik and Terlaky; and the problem (AQP) to minimize the absolute value of a quadratic function over a quadratic constraint proposed by Ye and Zhang. We show that, when $F(z)$ and the joint numerical range constraint are both convex, the optimal value of the convex optimization problem can be obtained by solving an SDP followed from a new development of the $\mathcal{S}$-procedure. The optimal solution can be approximated by conducting a bisection method on $[0,2\pi].$ On the other hand, if the joint numerical range of $f(x)$ and $g(x)$ is non-convex, the respective quadratic matrices of $f(x)$ and $g(x)$ must be linearly dependent. The linear dependence property enables us to solve (QSIC) and (AQP) accordingly by elementary analysis.


title: Deep learning and high harmonic generation

abstract: Using machine learning, we explore the utility of various deep neural networks (NN) when applied to high harmonic generation (HHG) scenarios. First, we train the NNs to predict the time-dependent dipole and spectra of HHG emission from reduced-dimensionality models of di- and triatomic systems based of on sets of randomly generated parameters (laser pulse intensity, internuclear distance, and molecular orientation). These networks, once trained, are useful tools to rapidly generate the HHG spectra of our systems. Similarly, we have trained the NNs to solve the inverse problem - to determine the molecular parameters based on HHG spectra or dipole acceleration data. These types of networks could then be used as spectroscopic tools to invert HHG spectra in order to recover the underlying physical parameters of a system. Next, we demonstrate that transfer learning can be applied to our networks to expand the range of applicability of the networks with only a small number of new test cases added to our training sets. Finally, we demonstrate NNs that can be used to classify molecules by type: di- or triatomic, symmetric or asymmetric, wherein we can even rely on fairly simple fully connected neural networks. With outlooks toward training with experimental data, these NN topologies offer a novel set of spectroscopic tools that could be incorporated into HHG experiments.


title: Development of Techniques Enabling Suborbital Small Object Capture Around Low Gravity Asteroids

abstract: The exploration of small body asteroids provides direct access to the primitive building blocks of our solar system. Sample and return missions enable chemical and radioisotope studies which not only provide evidence for the formation of the solar system, but also a basic understanding of where resources might be found for future exploitation. The touch-down and sample techniques established by Hayabusa-2 and OSIRIS-REx accomplish this mission by physically touching down on the asteroid and collecting samples into a basket extended via a probe from the bottom of the spacecraft. This technique has been demonstrated to work, but contains a high cost in both mission operations as well as the size and complexity of the collection mechanism itself. This paper explores an alternative sample and return technique by exploiting the recent discovery of regolith particle ejections from Bennu. Particles ejected from the surface of Bennu are typically 1 cm in size and spend several hours in flight, suggesting the possibility that nanospacecraft deployed from the mother-spacecraft could chase down, collect and return with the sample with minimal sensor and delta-V capability. Key aspects of this mission are developed to reduce risk, and an overall mission concept is developed to establish plausibility.


title: On monophonic position sets in graphs

abstract: The general position problem in graph theory asks for the largest set $S$ of vertices of a graph $G$ such that no shortest path of $G$ contains more than two vertices of $S$. In this paper we consider a variant of the general position problem called the \emph{monophonic position problem}, obtained by replacing `shortest path' by `induced path'. We prove some basic properties and bounds for the monophonic position number of a graph and determine the monophonic position number of some graph families, including unicyclic graphs, complements of bipartite graphs and split graphs. We show that the monophonic position number of triangle-free graphs is bounded above by the independence number. We present realisation results for the general position number, monophonic position number and monophonic hull number. Finally we discuss the complexity of the monophonic position problem.


title: Neural network fast-classifies biological images using features selected after their random-forests-importance to power smart microscopy

abstract: Artificial intelligence is nowadays used for cell detection and classification in optical microscopy, during post-acquisition analysis. The microscopes are now fully automated and next expected to be smart, to make acquisition decisions based on the images. It calls for analysing them on the fly. Biology further imposes training on a reduced dataset due to cost and time to prepare the samples and have the datasets annotated by experts. We propose here a real-time image processing, compliant with these specifications by balancing accurate detection and execution performance. We characterised the images using a generic, high-dimensional feature extractor. We then classified the images using machine learning for the sake of understanding the contribution of each feature in decision and execution time. We found that the non-linear-classifier random forests outperformed Fisher's linear discriminant. More importantly, the most discriminant and time-consuming features could be excluded without any significant loss in accuracy, offering a substantial gain in execution time. It suggests a feature-group redundancy likely related to the biology of the observed cells. We offer a method to select fast and discriminant features. In our assay, a 79.6 $\pm$ 2.4 % accurate classification of a cell took 68.7 $\pm$ 3.5 ms (mean $\pm$ SD, 5-fold cross-validation nested in 10 bootstrap repeats), corresponding to 14 cells per second, dispatched into 8 phases of the cell cycle using 12 feature-groups and operating a consumer market ARM-based embedded system. Interestingly, a simple neural network offered similar performances paving the way to faster training and classification, using parallel execution on a general-purpose graphic processing unit. Finally, this strategy is also usable for deep neural networks paving the way to optimising these algorithms for smart microscopy.


title: 2-adic Valuations of Quadratic Sequences

abstract: Properties of 2-adic valuation sequences for general quadratic polynomials with integer coefficients are determined directly from the coefficients. These properties include boundedness or unboundedness, periodicity, and valuations at terminating nodes. We completely describe the periodic sequences in the bounded case. Throughout, we frame results in terms of trees and sequences.


title: Learning from History for Byzantine Robust Optimization

abstract: Byzantine robustness has received significant attention recently given its importance for distributed and federated learning. In spite of this, we identify severe flaws in existing algorithms even when the data across the participants is identically distributed. First, we show realistic examples where current state of the art robust aggregation rules fail to converge even in the absence of any Byzantine attackers. Secondly, we prove that even if the aggregation rules may succeed in limiting the influence of the attackers in a single round, the attackers can couple their attacks across time eventually leading to divergence. To address these issues, we present two surprisingly simple strategies: a new robust iterative clipping procedure, and incorporating worker momentum to overcome time-coupled attacks. This is the first provably robust method for the standard stochastic optimization setting. Our code is open sourced at https://github.com/epfml/byzantine-robust-optimizer.


title: Atomic-scale imaging of laser-driven electron dynamics in solids using subcycle-resolved x-ray-optical wave mixing

abstract: We investigate laser-driven electron dynamics in solids on the atomic scale and in real space within Floquet formalism, and develop a method based on subcycle-resolved x-ray-optical wave mixing to reconstruct those dynamics. We analyze how time-reversal and inversion symmetries influence properties of optically-induced charge distributions and microscopic electron currents. Several examples for the $\mu$th-order microscopic optical response of band-gap crystals are shown and compared for cases when there is either a considerable or a vanishing $\mu$th-order macroscopic response. We then analyze the consequence of crystal symmetries on subcycle-resolved x-ray-optical wave mixing, a process in which an x-ray pulse of a duration shorter than the optical cycle of the driving pulse interacts with a laser-dressed crystal. Based on this analysis, we develop a method to reconstruct amplitudes and phases of Fourier components of optically-induced charge distributions from momentum and delay dependence of x-ray-optical wave-mixing spectra. Subcycle-resolved x-ray-optical wave mixing also reveals phases of temporal oscillations of microscopic optical response and some properties of microscopic laser-driven electron currents.


title: Solving Black-Box Optimization Challenge via Learning Search Space Partition for Local Bayesian Optimization

abstract: Black-box optimization is one of the vital tasks in machine learning, since it approximates real-world conditions, in that we do not always know all the properties of a given system, up to knowing almost nothing but the results. This paper describes our approach to solving the black-box optimization challenge at NeurIPS 2020 through learning search space partition for local Bayesian optimization. We describe the task of the challenge as well as our algorithm for low budget optimization that we named \texttt{SPBOpt}. We optimize the hyper-parameters of our algorithm for the competition finals using multi-task Bayesian optimization on results from the first two evaluation settings. Our approach has ranked third in the competition finals.


title: Motional heating of spatially extended ion crystals

abstract: We study heating of motional modes of a single ion and of extended ion crystals trapped in a linear radio frequency (rf) Paul trap with a precision of $\Delta \dot{\bar{n}} \approx 0.2 $ phonons s$^{-1}$. Single-ion axial and radial heating rates are consistent and electric field noise has been stable over the course of four years. At a secular frequency of $\omega_\mathrm{sec}=2\pi\times620$ kHz, we measure $\dot{\bar{n}} = 0.56(6)$ phonons s$^{-1}$ per ion for the center-of-mass (com) mode of linear chains of up to eleven ions and observe no significant heating of the out-of-phase (oop) modes. By displacing the ions away from the nodal line, inducing excess micromotion, rf noise heats the com mode quadratically as a function of radial displacement $r$ by $\dot{\bar{n}}(r)/ r^2 = 0.89(4)$ phonons s$^{-1}$ $\mu$m$^{-2}$ per ion, while the oop modes are protected from rf-noise induced heating in linear chains. By changing the quality factor of the resonant rf circuit from $Q=542$ to $Q=204$, we observe an increase of rf noise by a factor of up to 3. We show that the rf-noise induced heating of motional modes of extended crystals also depends on the symmetry of the crystal and of the mode itself. As an example, we consider several 2D and 3D crystal configurations. Heating rates of up to 500 phonons s$^{-1}$ are observed for individual modes, giving rise to a total kinetic energy increase and thus a fractional time dilation shift of up to $-0.3\times 10^{-18}$ s$^{-1}$ of the total system. In addition, we detail on how the excitation probability of the individual ions is reduced and decoherence is increased due to the Debye-Waller effect.


title: The half-space Airy stat process

abstract: We study the multipoint distribution of stationary half-space last passage percolation with exponentially weighted times. We derive both finite-size and asymptotic results for this distribution. In the latter case we observe a new one-parameter process we call half-space Airy stat. It is a one-parameter generalization of the Airy stat process of Baik-Ferrari-P\'ech\'e, which is recovered far away from the diagonal. All these results extend the one-point results previously proven by the authors.


title: Cellular and Developmental Basis of Avian Structural Coloration

abstract: Vivid structural colors in birds are a conspicuous and vital part of their phenotype. They are produced by a rich diversity of integumentary photonic nanostructures in skin and feathers. Unlike pigmentary coloration, whose molecular genetic basis is being elucidated, little is known regarding the pathways underpinning organismal structural coloration. Here, we review available data on the development of avian structural colors. In particular, feather photonic nanostructures are understood to be intracellularly self-assembled by physicochemical forces typically seen in soft colloidal systems. We identify promising avenues for future research that can address current knowledge gaps, which is also highly relevant for the sustainable engineering of advanced bioinspired and biomimetic materials.


title: Semirigid GCD domains II

abstract: Let $D$ be an integral domain with quotient field $K,$ throughout$.$ Call two elements $x,y\in D\backslash \{0\}$ $v$-coprime if $xD\cap yD=xyD.$ Call a nonzero non unit $r$ of an integral domain $D$ rigid if for all $x,y|r$ we have $x|y$ or $y|x.$ Also call $D$ semirigid if every nonzero non unit of $D$ is expressible as a finite product of rigid elements. We show that a semirigid domain $D$ is a GCD domain if and only if $D$ satisfies $\ast :$ product of every pair of non-$v$-coprime rigid elements is again rigid. Next call $a\in D$ a valuation element if $aV\cap D=aD$ for some valuation ring $% V $ with $D\subseteq V\subseteq K$ and call $D$ a VFD if every nonzero non unit of $D$ is a finite product of valuation elements. It turns out that a valuation element is what we call a packed element: a rigid element $r$ all of whose powers are rigid and $\sqrt{rD}$ is a prime ideal. Calling $D$ a semi packed domain (SPD) if every nonzero non unit of $D$ is a finite product of packed elements, we study SPDs and explore situations in which an SPD is a semirigid GCD domain.


title: Spontaneous pattern growth on chocolate surface: simulation and experiments

abstract: The natural variation of temperature at ambient conditions produces spontaneous patterns on the surface of chocolate, which result from fat bloom. These patterns are peculiar because of their shape and cannot be obtained by controlled temperature conditions. The formation of these spontaneous grains on the surface of chocolate is studied on experimental and theoretical grounds.Three different kinds of experiments were conducted: observation of formed patterns in time, atomic force microscopy of the initial events on the grain formation and rheology of the melted chocolate. The patterns observed in our experiments follow the trends described by the Avrami model, which considers that is possible to define a characteristic time scale that governs the growth of grains starting from germ nuclei. Through computer simulations, in the NVT ensemble using a coarse-grained model of triacylglycerides molecules, we studied the process of nucleation that starts the pattern growth and that is consistent with the Avrami model.


title: Graceful Labellings of Various Cyclic Snakes

abstract: In this paper, we present a new sufficiency condition to obtain a graceful labelling for every $kC_{4n}$ snake and use this condition to label every such snake for $n=1,2,\ldots,6$. Then, we extend this result to cyclic snakes where the cycles lengths vary. Also, we obtain new results on the (near) graceful labelling of cyclic snakes based on cycles of lengths $n=6, 10, 14$, completely solving the case $n=6$.


title: Exploring and Interrogating Astrophysical Data in Virtual Reality

abstract: Scientists across all disciplines increasingly rely on machine learning algorithms to analyse and sort datasets of ever increasing volume and complexity. Although trends and outliers are easily extracted, careful and close inspection will still be necessary to explore and disentangle detailed behavior, as well as identify systematics and false positives. We must therefore incorporate new technologies to facilitate scientific analysis and exploration. Astrophysical data is inherently multi-parameter, with the spatial-kinematic dimensions at the core of observations and simulations. The arrival of mainstream virtual-reality (VR) headsets and increased GPU power, as well as the availability of versatile development tools for video games, has enabled scientists to deploy such technology to effectively interrogate and interact with complex data. In this paper we present development and results from custom-built interactive VR tools, called the iDaVIE suite, that are informed and driven by research on galaxy evolution, cosmic large-scale structure, galaxy-galaxy interactions, and gas/kinematics of nearby galaxies in survey and targeted observations. In the new era of Big Data ushered in by major facilities such as the SKA and LSST that render past analysis and refinement methods highly constrained, we believe that a paradigm shift to new software, technology and methods that exploit the power of visual perception, will play an increasingly important role in bridging the gap between statistical metrics and new discovery. We have released a beta version of the iDaVIE software system that is free and open to the community.


title: Application of computer simulation results and machine learning in analysis of microwave radiothermometry data

abstract: This work was done with the aim of developing the fundamental breast cancer early differential diagnosis foundations based on modeling the space-time temperature distribution using the microwave radiothermometry method and obtained data intelligent analysis. The article deals with the machine learning application in the microwave radiothermometry data analysis. The problems associated with the construction mammary glands temperature fields computer models for patients with various diagnostics classes, are also discussed. With the help of a computer experiment, based on the machine learning algorithms set (logistic regression, naive Bayesian classifier, support vector machine, decision tree, gradient boosting, K-nearest neighbors, etc.) usage, the mammary glands temperature fields computer models set adequacy.


title: Existence and uniqueness for a viscoelastic Kelvin-Voigt model with nonconvex stored energy

abstract: We consider nonlinear viscoelastic materials of Kelvin-Voigt type with stored energies satisfying an Andrews-Ball condition, allowing for non convexity in a compact set. Existence of weak solutions with deformation gradients in $H^1$ is established for energies of any superquadratic growth. In two space dimensions, weak solutions notably turn out to be unique in this class. Conservation of energy for weak solutions in two and three dimensions, as well as global regularity for smooth initial data in two dimensions are established under additional mild restrictions on the growth of the stored energy.


title: The Simons Observatory: Overview of data acquisition, control, monitoring, and computer infrastructure

abstract: The Simons Observatory (SO) is an upcoming polarized cosmic microwave background (CMB) survey experiment with three small-aperture telescopes and one large-aperture telescope that will observe from the Atacama Desert in Chile. In total, SO will field over 60,000 transition-edge sensor (TES) bolometers in six spectral bands centered between 27 and 280 GHz to achieve the sensitivity necessary to measure or constrain numerous cosmological parameters, including the tensor-to-scalar ratio, effective number of relativistic species, and sum of the neutrino masses. The SO scientific goals require coordination and control of the hardware distributed among the four telescopes on site. To meet this need, we have designed and built an open-sourced platform for distributed system management, called the Observatory Control System (ocs). This control system interfaces with all subsystems including the telescope control units, the microwave multiplexing readout electronics, and the cryogenic thermometry. We have also developed a system for live monitoring of housekeeping data and alerting, both of which are critical for remote observation. We take advantage of existing open source projects, such as crossbar for RPC and PubSub, twisted for asynchronous events, grafana for online remote monitoring, and docker for containerization. We provide an overview of the SO software and computer infrastructure, including the integration of SO-developed code with open source resources and lessons learned while testing at SO labs developing hardware systems as we prepare for deployment.


title: A general approach to state-dependent optical tweezer traps for polar molecules

abstract: State-dependent optical tweezers can be used to trap a pair of molecules with a separation much smaller than the wavelength of the trapping light, greatly enhancing the dipole-dipole interaction between them. Here we describe a general approach to producing these state-dependent potentials using the tensor part of the ac Stark shift and show how it can be used to carry out two-qubit gates between pairs of molecules. The method is applicable to broad classes of molecules including bialkali molecules produced by atom association and those amenable to direct laser cooling.


title: Self-evaporation dynamics of quantum droplets in a 41K-87Rb mixture

abstract: We theoretically investigate the self-evaporation dynamics of quantum droplets in a 41K-87Rb mixture, in free-space. The dynamical formation of the droplet and the effects related to the presence of three-body losses are analyzed by means of numerical simulations. We identify a regime of parameters allowing for the observation of the droplet self-evaporation in a feasible experimental setup.


title: Small Business Classification By Name: Addressing Gender and Geographic Origin Biases

abstract: Small business classification is a difficult and important task within many applications, including customer segmentation. Training on small business names introduces gender and geographic origin biases. A model for predicting one of 66 business types based only upon the business name was developed in this work (top-1 f1-score = 60.2%). Two approaches to removing the bias from this model are explored: replacing given names with a placeholder token, and augmenting the training data with gender-swapped examples. The results for these approaches is reported, and the bias in the model was reduced by hiding given names from the model. However, bias reduction was accomplished at the expense of classification performance (top-1 f1-score = 56.6%). Augmentation of the training data with gender-swapping samples proved less effective at bias reduction than the name hiding approach on the evaluated dataset.


title: Junction conditions in infinite derivative gravity

abstract: The junction conditions for the infinite derivative gravity theory ${R{+}RF(\Box)R}$ are derived under the assumption that the conditions can be imposed by avoiding the `ill-defined expressions' in the theory of distributions term by term in infinite summations. We find that the junction conditions of such non-local theories are much more restrictive than in local theories, since the conditions comprise an infinite number of equations for the Ricci scalar. These conditions can constrain the geometry far beyond the matching hypersurface. Furthermore, we derive the junction field equations which are satisfied by the energy-momentum on the hypersurface. It turns out that the theory still allows some matter content on the hypersurface (without external flux and external tension), but with a traceless energy-momentum tensor. We also discuss the proper matching condition where no matter is concentrated on the hypersurface. Finally, we explore the possible applications and consequences of our results to the braneworld scenarios and star models. Particularly, we find that the internal tension is given purely by the trace of the energy-momentum tensor of the matter confined to the brane. Consequences of the junction conditions are illustrated on two simple examples of static and collapsing stars. It is demonstrated that even without solving the field equations the geometry on one side of the hypersurface can be determined to a great extent by the geometry on the other side if the Ricci scalar is analytic. We further show that some usual star models in the general relativity are no longer solutions of the infinite derivative gravity.


title: Scaling Laws for Gaussian Random Many-Access Channels

abstract: This paper considers a Gaussian multiple-access channel with random user activity where the total number of users $\ell_n$ and the average number of active users $k_n$ may grow with the blocklength $n$. For this channel, it studies the maximum number of bits that can be transmitted reliably per unit-energy as a function of $\ell_n$ and $k_n$. When all users are active with probability one, i.e., $\ell_n = k_n$, it is demonstrated that if $k_n$ is of an order strictly below $n/\log n$, then each user can achieve the single-user capacity per unit-energy $(\log e)/N_0$ (where $N_0/ 2$ is the noise power) by using an orthogonal-access scheme. In contrast, if $k_n$ is of an order strictly above $n/\log n$, then the capacity per unit-energy is zero. Consequently, there is a sharp transition between orders of growth where interference-free communication is feasible and orders of growth where reliable communication at a positive rate per unit-energy is infeasible. It is further demonstrated that orthogonal-access schemes in combination with orthogonal codebooks, which achieve the capacity per unit-energy when the number of users is bounded, can be strictly suboptimal. When the user activity is random, i.e., when $\ell_n$ and $k_n$ are different, it is demonstrated that if $k_n\log \ell_n$ is sublinear in $n$, then each user can achieve the single-user capacity per unit-energy $(\log e)/N_0$. Conversely, if $k_n\log \ell_n$ is superlinear in $n$, then the capacity per unit-energy is zero. Consequently, there is again a sharp transition between orders of growth where interference-free communication is feasible and orders of growth where reliable communication at a positive rate is infeasible that depends on the asymptotic behaviours of both $\ell_n$ and $k_n$. It is further demonstrated that orthogonal-access schemes, which are optimal when $\ell_n = k_n$, can be strictly suboptimal.


title: Universal Approximation in Dropout Neural Networks

abstract: We prove two universal approximation theorems for a range of dropout neural networks. These are feed-forward neural networks in which each edge is given a random $\{0,1\}$-valued filter, that have two modes of operation: in the first each edge output is multiplied by its random filter, resulting in a random output, while in the second each edge output is multiplied by the expectation of its filter, leading to a deterministic output. It is common to use the random mode during training and the deterministic mode during testing and prediction. Both theorems are of the following form: Given a function to approximate and a threshold $\varepsilon>0$, there exists a dropout network that is $\varepsilon$-close in probability and in $L^q$. The first theorem applies to dropout networks in the random mode. It assumes little on the activation function, applies to a wide class of networks, and can even be applied to approximation schemes other than neural networks. The core is an algebraic property that shows that deterministic networks can be exactly matched in expectation by random networks. The second theorem makes stronger assumptions and gives a stronger result. Given a function to approximate, it provides existence of a network that approximates in both modes simultaneously. Proof components are a recursive replacement of edges by independent copies, and a special first-layer replacement that couples the resulting larger network to the input. The functions to be approximated are assumed to be elements of general normed spaces, and the approximations are measured in the corresponding norms. The networks are constructed explicitly. Because of the different methods of proof, the two results give independent insight into the approximation properties of random dropout networks. With this, we establish that dropout neural networks broadly satisfy a universal-approximation property.


title: Probabilistic Aspects of Voting, Intransitivity and Manipulation

abstract: These lecture notes are based on lectures given in 2019 Saint-Flour Probability School.


title: Stable maps to Looijenga pairs: orbifold examples

abstract: In arXiv:2011.08830 we established a series of correspondences relating five enumerative theories of log Calabi-Yau surfaces, i.e. pairs $(Y,D)$ with $Y$ a smooth projective complex surface and $D=D_1+\dots +D_l$ an anticanonical divisor on $Y$ with each $D_i$ smooth and nef. In this paper we explore the generalisation to $Y$ being a smooth Deligne-Mumford stack with projective coarse moduli space of dimension 2, and $D_i$ nef $\mathbb{Q}$-Cartier divisors. We consider in particular three infinite families of orbifold log Calabi-Yau surfaces, and for each of them we provide closed form solutions of the maximal contact log Gromov-Witten theory of the pair $(Y,D)$, the local Gromov-Witten theory of the total space of $\bigoplus_i \mathcal{O}_Y(-D_i)$, and the open Gromov-Witten theory of toric orbi-branes in a Calabi-Yau 3-orbifold associated to $(Y,D)$. We also consider new examples of BPS integral structures underlying these invariants, and relate them to the Donaldson-Thomas theory of a symmetric quiver specified by $(Y,D)$, and to a class of open/closed BPS invariants.


title: A two-component Comptonisation model for the type-B QPO in MAXI J1348-630

abstract: Spectral-timing analysis of the fast variability observed in X-rays is a powerful tool to study the physical and geometrical properties of the accretion/ejection flows in black-hole binaries. The origin of type-B quasi-periodic oscillations (QPO), predominantly observed in black-hole candidates in the soft-intermediate state, has been linked to emission arising from the relativistic jet. In this state, the X-ray spectrum is characterised by a soft-thermal blackbody-like emission due to the accretion disc, an iron emission line (in the 6-7 keV range), and a power-law like hard component due to Inverse-Compton scattering of the soft-photon source by hot electrons in a corona or the relativistic jet itself. The spectral-timing properties of MAXI J1348-630 have been recently studied using observations obtained with the NICER observatory. The data show a strong type-B QPO at ~4.5 Hz with increasing fractional rms amplitude with energy and positive lags with respect to a reference band at 2-2.5 keV. We use a variable-Comptonisation model that assumes a sinusoidal coherent oscillation of the Comptonised X-ray flux and the physical parameters of the corona at the QPO frequency, to fit simultaneously the energy-dependent fractional rms amplitude and phase lags of this QPO. We show that two physically-connected Comptonisation regions can successfully explain the radiative properties of the QPO in the full 0.8-10 keV energy range.


title: Assessing Pattern Recognition Performance of Neuronal Cultures through Accurate Simulation

abstract: Previous work has shown that it is possible to train neuronal cultures on Multi-Electrode Arrays (MEAs), to recognize very simple patterns. However, this work was mainly focused to demonstrate that it is possible to induce plasticity in cultures, rather than performing a rigorous assessment of their pattern recognition performance. In this paper, we address this gap by developing a methodology that allows us to assess the performance of neuronal cultures on a learning task. Specifically, we propose a digital model of the real cultured neuronal networks; we identify biologically plausible simulation parameters that allow us to reliably reproduce the behavior of real cultures; we use the simulated culture to perform handwritten digit recognition and rigorously evaluate its performance; we also show that it is possible to find improved simulation parameters for the specific task, which can guide the creation of real cultures.


title: The outburst of a 60 min AM CVn reveals peculiar color evolution: implications for outbursts in long period double white dwarfs

abstract: We report on multi-wavelength observations during quiescence and of the first detected outburst of the ~60 min orbital period AM CVn SDSS J113732+405458. Using X-ray and UV observations we determined an upper limit duration of the event of about one year. The amplitude of the outburst was remarkably small, of around one magnitude in r and 0.5 magnitudes in g. We have also investigated the color variations of SDSS J113732+405458 and other long period AM CVns in outbursts and identified a track on the color-magnitude diagram that is not compatible with the predictions of the disk instability model, suggesting that some outbursts in long period AM CVns are caused by enhanced mass-transfer. To our knowledge, these are the first studies of the color evolution in AM CVns. During quiescence we measured an X-ray luminosity for SDSS J113732+405458 of ~3E29 erg/s in the 0.5-10 keV band. This indicates a very low accretion rate, in agreement with the disk instability model for long period systems. However, such a model predicts stable disks at somewhat long periods. The discovery of this system outburst, along with similarities to the long period system SDSS J080710+485259 with a comparably long, weak outburst, indicates that these enhanced mass-transfer events may be more common in long period AM CVns. A larger sample would be needed to determine empirically at what period, if any, the disk instability stops functioning entirely. Finally, we identified an infrared excess in the quiescence spectrum attributable to the donor. This makes SDSS J113732+405458 the second AM CVn to have a directly detected donor.


title: Proceduray -- A light-weight engine for procedural primitive ray tracing

abstract: We introduce Proceduray, an engine for real-time ray tracing of procedural geometry. Its motivation is the current lack of mid-level abstraction tools for scenes with primitives involving intersection shaders. Those scenes impose strict engine design choices since they need flexibility in the shader table setup. Proceduray aims at providing a fair tradeoff between that flexibility and productivity. It also aims to be didactic. Shader table behavior can be confusing because parameters for indexing come from different parts of a system, involving both host and device code. This is different in essence from ray tracing triangle meshes (which must use a built-in intersection shader for all objects) or rendering with the traditional graphics or compute pipelines. Additionals goals of the project include fomenting deeper discussions about DirectX RayTracing (DXR) host code and providing a good starting point for developers trying to deal with procedural geometry using DXR.


title: Diffusion processes in germanium and silicon films grown on Si$_3$N$_4$ substrates

abstract: In this article, the results of investigation of processes occurring during the molecular-beam deposition of germanium layers on Si$_3$N$_4$ dielectric substrates within a wide range of the Ge film growth temperatures (30 to 600{\deg}C) are presented. The intensity of the IR absorption bands related to the vibrations of the N$-$H and Si$-$N bonds are established to decrease with the increase of the Ge deposition temperature. It appears that this phenomenon cannot be explained only as a thermally activated process. Simultaneously, the peak corresponding to the Ge$-$N vibration bonds emerges in the X-ray photoelectron spectra. We suppose that the deposition of germanium layers on Si3N4 dielectric substrates containing hydrogen atoms causes the diffusion of hydrogen atoms from the dielectric layer into the growing film. The experimental results may be interpreted in terms of a model, according to which the migration of hydrogen atoms from the Si$_3$N$_4$ layer into the growing germanium film is due to the difference in chemical potentials of hydrogen atoms in the dielectric layer and the germanium film. This process initiates the diffusion of germanium atoms in the opposite direction, into the Si$_3$N$_4$ layer, where they connect to free bonds of nitrogen atoms arising due to the escape of hydrogen atoms. The analogous processes occur during the deposition of silicon layers on Si$_3$N$_4$ substrates.


title: High-energy cosmic particles

abstract: A review of the status of the knowledge in the field of High-energy cosmic particles is presented. The spectrum, arrival direction distribution and composition measurements are summarized, together with some implications for the understanding of the cosmic ray origin and their propagation. Special emphasis is put in the ultra-high energy range, corresponding to particles of extragalactic origin.


title: When Machine Learning Meets Quantum Computers: A Case Study

abstract: Along with the development of AI democratization, the machine learning approach, in particular neural networks, has been applied to wide-range applications. In different application scenarios, the neural network will be accelerated on the tailored computing platform. The acceleration of neural networks on classical computing platforms, such as CPU, GPU, FPGA, ASIC, has been widely studied; however, when the scale of the application consistently grows up, the memory bottleneck becomes obvious, widely known as memory-wall. In response to such a challenge, advanced quantum computing, which can represent 2^N states with N quantum bits (qubits), is regarded as a promising solution. It is imminent to know how to design the quantum circuit for accelerating neural networks. Most recently, there are initial works studying how to map neural networks to actual quantum processors. To better understand the state-of-the-art design and inspire new design methodology, this paper carries out a case study to demonstrate an end-to-end implementation. On the neural network side, we employ the multilayer perceptron to complete image classification tasks using the standard and widely used MNIST dataset. On the quantum computing side, we target IBM Quantum processors, which can be programmed and simulated by using IBM Qiskit. This work targets the acceleration of the inference phase of a trained neural network on the quantum processor. Along with the case study, we will demonstrate the typical procedure for mapping neural networks to quantum circuits.


title: Geometric frustration on the trillium lattice in a magnetic metal-organic framework

abstract: In the dense metal-organic framework Na[Mn(HCOO)$_3$], Mn$^{2+}$ ions ($S=\frac{5}{2}$) occupy the nodes of a `trillium' hyperkagome net. We show that this material exhibits a variety of behaviour characteristic of geometric frustration: the N\'eel transition is suppressed well below the characteristic magnetic interaction strength; short-range magnetic order persists far above the N\'eel temperature; and the magnetic susceptibility exhibits a pseudo-plateau at $\frac{1}{3}$-saturation magnetisation. We demonstrate that a simple nearest-neighbour Heisenberg antiferromagnet model accounts quantitatively for each observation, and hence Na[Mn(HCOO)$_3$] is the first experimental realisation of this model on the trillium net. We develop a mapping between this trillium model and that on the two-dimensional Shastry-Sutherland lattice, and demonstrate how both link geometric frustration within the classical spin liquid regime to a strong magnetocaloric response at low fields.

categories: cond-mat.str-el cond-mat.mtrl-sci cond-mat.other cond-mat.stat-mech


title: The Inductive Graph Dimension from The Minimum Edge Clique Cover

abstract: In this paper we prove that the inductively defined graph dimension has a simple additive property under the join operation. The dimension of the join of two simple graphs is one plus the sum of the dimensions of the component graphs: $\mathrm{dim}\, (G_1+ G_2) = 1 +\mathrm{dim}\, G_1+ \mathrm{dim}\, G_2$. We use this formula to derive an expression for the inductive dimension of an arbitrary finite simple graph from its minimum edge clique cover. A corollary of the formula is that any arbitrary finite simple graph whose maximal cliques are all of order $N$ has dimension $N-1$. We finish by finding lower and upper bounds on the inductive dimension of a simple graph in terms of its clique number.


title: On the distribution of scrambled $(0,m,s)$-nets over unanchored boxes

abstract: We introduce a new quality measure to assess randomized low-discrepancy point sets of finite size $n$. This new quality measure, which we call "pairwise sampling dependence index", is based on the concept of negative dependence. A negative value for this index implies that the corresponding point set integrates the indicator function of any unanchored box with smaller variance than the Monte Carlo method. We show that scrambled $(0,m,s)-$nets have a negative pairwise sampling dependence index. We also illustrate through an example that randomizing via a digital shift instead of scrambling may yield a positive pairwise sampling dependence index.


title: Exact matrix model for generalized Ising model

abstract: The Ising model was generalized to a system of cells interacting exclusively by presence of shared spins. Within the cells there are interactions of any complexity, the simplest intracell interactions come down to the Ising model. The system may be not only one-dimensional but also two-dimensional, three-dimensional, etc. The purpose of the paper is to develop an approach to constructing the exact matrix model for any considered system in the simplest way. Without this, it is almost impossible to analyze complex systems. After trying a lot of ways, the approach has been simplified to suit first year students. Using the approach, the exact matrix model for a two-dimensional generalized Ising model was constructed. The following approaches have already been developed to obtaining and analyzing partition function and various thermodynamic functions. They will be considered in following papers.


title: Phase separation and multibody effects in three-dimensional active Brownian particles

abstract: Simulation studies of the phase diagram of repulsive active Brownian particles in three dimensions reveal that the region of motility-induced phase separation between a high and low density phase is enclosed by a region of gas-crystal phase separation. Near-critical loci and structural crossovers can additionally be identified in analogy with simple fluids. Motivated by the striking similarity to the behaviour of equilibrium fluids with short-ranged pair-wise attractions, we show that a direct mapping to pair potentials in the dilute limit implies interactions that are insufficiently attractive to engender phase separation. Instead, this is driven by the emergence of multi-body effects associated with particle caging that occurs at sufficiently high number density. We quantify these effects via information-theoretical measures of $n$-body effective interactions extracted from the configurational structure.


title: Learning Complex 3D Human Self-Contact

abstract: Monocular estimation of three dimensional human self-contact is fundamental for detailed scene analysis including body language understanding and behaviour modeling. Existing 3d reconstruction methods do not focus on body regions in self-contact and consequently recover configurations that are either far from each other or self-intersecting, when they should just touch. This leads to perceptually incorrect estimates and limits impact in those very fine-grained analysis domains where detailed 3d models are expected to play an important role. To address such challenges we detect self-contact and design 3d losses to explicitly enforce it. Specifically, we develop a model for Self-Contact Prediction (SCP), that estimates the body surface signature of self-contact, leveraging the localization of self-contact in the image, during both training and inference. We collect two large datasets to support learning and evaluation: (1) HumanSC3D, an accurate 3d motion capture repository containing $1,032$ sequences with $5,058$ contact events and $1,246,487$ ground truth 3d poses synchronized with images collected from multiple views, and (2) FlickrSC3D, a repository of $3,969$ images, containing $25,297$ surface-to-surface correspondences with annotated image spatial support. We also illustrate how more expressive 3d reconstructions can be recovered under self-contact signature constraints and present monocular detection of face-touch as one of the multiple applications made possible by more accurate self-contact models.


title: Gravitational Edge Modes, Coadjoint Orbits, and Hydrodynamics

abstract: The phase space of general relativity in a finite subregion is characterized by edge modes localized at the codimension-2 boundary, transforming under an infinite-dimensional group of symmetries. The quantization of this symmetry algebra is conjectured to be an important aspect of quantum gravity. As a step towards quantization, we derive a complete classification of the positive-area coadjoint orbits of this group for boundaries that are topologically a 2-sphere. This classification parallels Wigner's famous classification of representations of the Poincar\'e group since both groups have the structure of a semidirect product. We find that the total area is a Casimir of the algebra, analogous to mass in the Poincar\'e group. A further infinite family of Casimirs can be constructed from the curvature of the normal bundle of the boundary surface. These arise as invariants of the little group, which is the group of area-preserving diffeomorphisms, and are the analogues of spin. Additionally, we show that the symmetry group of hydrodynamics appears as a reduction of the corner symmetries of general relativity. Coadjoint orbits of both groups are classified by the same set of invariants, and, in the case of the hydrodynamical group, the invariants are interpreted as the generalized enstrophies of the fluid.


title: Basis properties of Fucik eigenfunctions

abstract: We establish sufficient assumptions on sequences of Fucik eigenvalues of the one-dimensional Laplacian which guarantee that the corresponding Fucik eigenfunctions form a Riesz basis in $L^2(0,\pi)$.


title: Upper and Lower Bounds on the Performance of Kernel PCA

abstract: Principal Component Analysis (PCA) is a popular method for dimension reduction and has attracted an unfailing interest for decades. More recently, kernel PCA (KPCA) has emerged as an extension of PCA but, despite its use in practice, a sound theoretical understanding of KPCA is missing. We contribute several lower and upper bounds on the efficiency of KPCA, involving the empirical eigenvalues of the kernel Gram matrix and new quantities involving a notion of variance. These bounds show how much information is captured by KPCA on average and contribute a better theoretical understanding of its efficiency. We demonstrate that fast convergence rates are achievable for a widely used class of kernels and we highlight the importance of some desirable properties of datasets to ensure KPCA efficiency.


title: Quantization and martingale couplings

abstract: Quantization provides a very natural way to preserve the convex order when approximating two ordered probability measures by two finitely supported ones. Indeed, when the convex order dominating original probability measure is compactly supported, it is smaller than any of its dual quantizations while the dominated original measure is greater than any of its stationary (and therefore any of its optimal) quadratic primal quantization. Moreover, the quantization errors then correspond to martingale couplings between each original probability measure and its quantization. This permits to prove that any martingale coupling between the original probability measures can be approximated by a martingale coupling between their quantizations in Wassertein distance with a rate given by the quantization errors but also in the much finer adapted Wassertein distance. As a consequence, while the stability of (Weak) Martingale Optimal Transport problems with respect to the marginal distributions has only been established in dimension $1$ so far, their value function computed numerically for the quantized marginals converges in any dimension to the value for the original probability measures as the numbers of quantization points go to $\infty$.


title: The first higher Stasheff-Tamari orders are quotients of the higher Bruhat orders

abstract: We prove the conjecture that the higher Tamari orders of Dimakis and M\"uller-Hoissen coincide with the first higher Stasheff--Tamari orders. To this end, we show that the higher Tamari orders may be conceived as the image of an order-preserving map from the higher Bruhat orders to the first higher Stasheff--Tamari orders. This map is defined by taking the first cross-section of a cubillage of a cyclic zonotope. We provide a new proof that this map is surjective and show further that the map is full, which entails the aforementioned conjecture. We explain how order-preserving maps which are surjective and full correspond to quotients of posets. Our results connect the first higher Stasheff--Tamari orders with the literature on the role of the higher Tamari orders in integrable systems.


title: The Cosmic-Ray Composition between 2 PeV and 2 EeV Observed with the TALE Detector in Monocular Mode

abstract: We report on a measurement of the cosmic ray composition by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector (FD). By making use of the Cherenkov light signal in addition to air fluorescence light from cosmic ray (CR) induced extensive air showers, the TALE FD can measure the properties of the cosmic rays with energies as low as $\sim 2$ PeV and exceeding 1 EeV. In this paper, we present results on the measurement of $X_{\rm max}$ distributions of showers observed over this energy range. Data collected over a period of $\sim 4$ years was analyzed for this study. The resulting $X_{\rm max}$ distributions are compared to the Monte Carlo (MC) simulated data distributions for primary cosmic rays with varying composition and a 4-component fit is performed. The comparison and fit are performed for energy bins, of width 0.1 or 0.2 in $\log_{10} (E/{\rm eV})$, spanning the full range of the measured energies. We also examine the mean $X_{\rm max}$ value as a function of energy for cosmic rays with energies greater than $10^{15.8}$ eV. Below $10^{17.3}$ eV, the slope of the mean $X_{\rm max}$ as a function of energy (the elongation rate) for the data is significantly smaller than that of all elements in the models, indicating that the composition is becoming heavier with energy in this energy range. This is consistent with a rigidity-dependent cutoff of events from galactic sources. Finally, an increase in the $X_{\rm max}$ elongation rate is observed at energies just above $10^{17}$ eV indicating another change in the cosmic rays composition.


title: Real-time observation of Cooper pair splitting showing strong non-local correlations

abstract: Controlled generation and detection of quantum entanglement between spatially separated particles constitute an essential prerequisite both for testing the foundations of quantum mechanics and for realizing future quantum technologies. Splitting of Cooper pairs from a superconductor provides entangled electrons at separate locations. However, experimentally accessing the individual split Cooper pairs constitutes a major unresolved issue as they mix together with electrons from competing processes. Here, we overcome this challenge with the first real-time observation of the splitting of individual Cooper pairs, enabling direct access to the time-resolved statistics of Cooper pair splitting. We determine the correlation statistics arising from two-electron processes and find a pronounced peak that is two orders of magnitude larger than the background. Our experiment thereby allows to unambiguously pinpoint and select split Cooper pairs with 99 % fidelity. These results open up an avenue for performing experiments that tap into the spin-entanglement of split Cooper pairs.


title: Imaging neutron capture cross sections: i-TED proof-of-concept and future prospects based on Machine-Learning techniques

abstract: i-TED is an innovative detection system which exploits Compton imaging techniques to achieve a superior signal-to-background ratio in ($n,\gamma$) cross-section measurements using time-of-flight technique. This work presents the first experimental validation of the i-TED apparatus for high-resolution time-of-flight experiments and demonstrates for the first time the concept proposed for background rejection. To this aim both $^{197}$Au($n,\gamma$) and $^{56}$Fe($n, \gamma$) reactions were measured at CERN n\_TOF using an i-TED demonstrator based on only three position-sensitive detectors. Two \cds detectors were also used to benchmark the performance of i-TED. The i-TED prototype built for this study shows a factor of $\sim$3 higher detection sensitivity than state-of-the-art \cds detectors in the $\sim$10~keV neutron energy range of astrophysical interest. This paper explores also the perspectives of further enhancement in performance attainable with the final i-TED array consisting of twenty position-sensitive detectors and new analysis methodologies based on Machine-Learning techniques.


title: Nonequilibrium Thermodynamics of Non-Ideal Chemical Reaction Networks

abstract: All current formulations of nonequilibrium thermodynamics of open chemical reaction networks rely on the assumption of non-interacting species. We develop a general theory which accounts for interactions between chemical species within a mean-field approach using activity coefficients. Thermodynamic consistency requires that rate equations do not obey to standard mass-action kinetics, but account for the interactions with concentration dependent kinetic constants. Many features of the ideal formulations are recovered. Crucially, the thermodynamic potential and the forces driving non-ideal chemical systems out of equilibrium are identified. Our theory is general and holds for any mean-field expression of the interactions leading to lower bounded free energies.


title: Identification of Metallic Objects using Spectral MPT Signatures: Object Characterisation and Invariants

abstract: The early detection of terrorist threats, such as guns and knives, through improved metal detection, has the potential to reduce the number of attacks and improve public safety and security. To achieve this, there is considerable potential to use the fields applied and measured by a metal detector to discriminate between different shapes and different metals since, hidden within the field perturbation, is object characterisation information. The magnetic polarizability tensor (MPT) offers an economical characterisation of metallic objects that can be computed for different threat and non-threat objects and has an established theoretical background, which shows that the induced voltage is a function of the hidden object's MPT coefficients. In this paper, we describe the additional characterisation information that measurements of the induced voltage over a range of frequencies offer compared to measurements at a single frequency. We call such object characterisations its MPT spectral signature. Then, we present a series of alternative rotational invariants for the purpose of classifying hidden objects using MPT spectral signatures. Finally, we include examples of computed MPT spectral signature characterisations of realistic threat and non-threat objects that can be used to train machine learning algorithms for classification purposes.


title: Progress in Understanding the Enigmatic Fast Radio Bursts

abstract: In less than a decade, fast radio bursts have gone from a single debated curiosity to a diverse extragalactic population with established host galaxies and energy scales. While a wide range of models remain viable, the central engines of FRBs are likely to involve energetic young magnetars, as confirmed by the recent discovery of a Galactic analog to these extragalactic bursts. Here we provide a brief introductory review of fast radio bursts, focusing on the rapid recent progress in observations of these enigmatic events, our understanding of their central engines, and their use as probes of the intergalactic medium. We caution against a rush to judgement on the mechanisms and classification of all FRBs: at this point, it remains plausible that there could be one dominant central engine, as well as the possibility that radio bursts are a generic feature produced by many different mechanisms. We also emphasize the importance of improved modeling of our Galaxy and Galactic halo, which otherwise impose systematic errors on every FRB line of sight. The future of science with fast radio bursts appears bright.


title: YouNiverse: Large-Scale Channel and Video Metadata from English-Speaking YouTube

abstract: YouTube plays a key role in entertaining and informing people around the globe. However, studying the platform is difficult due to the lack of randomly sampled data and of systematic ways to query the platform's colossal catalog. In this paper, we present YouNiverse, a large collection of channel and video metadata from English-language YouTube. YouNiverse comprises metadata from over 136k channels and 72.9M videos published between May 2005 and October 2019, as well as channel-level time-series data with weekly subscriber and view counts. Leveraging channel ranks from socialblade.com, an online service that provides information about YouTube, we are able to assess and enhance the representativeness of the sample of channels. Additionally, the dataset also contains a table specifying which videos a set of 449M anonymous users commented on. YouNiverse, publicly available at https://doi.org/10.5281/zenodo.4650046, will empower the community to do research with and about YouTube.


title: Tropical Normal Functions -- Higher Abel-Jacobi Invariants of Tropical cycles

abstract: We consider the variation of tropical Hodge structure (TVHS) associated to families of tropical varieties. The family of the tropical intermediate Jacobians of the associated tropical Hodge structure defines a bundle of tropical Jacobians, whose sections we call the tropical normal functions. We define formal sequential derivatives of these functions on the base with respect to the natural Gauss-Manin connection as the Hodge theoretic invariants detecting tropical cycles in the fibers. The associated invariants which are defined inductively are the higher Abel-Jacobi invariants in the tropical category. They naturally identify the tropical Bloch-Beilinson filtration on the tropical Chow group. We examine this construction on the moduli of tropical curves with marked points, in order to study the tropical tautological classes in the tautological ring of $\mathcal{M}_{g,n}^{\text{trop}}$. The expectation is the nontriviality of these cycles could be examined with less complexity in the tropical category. The construction is compatible with the tropicalization functor on the category of schemes, and the aforementioned procedure will also provide an alternative way to examine the relations in the tautological ring of $\mathcal{M}_{g,n}$ in the schemes category.


title: Evidence of a $J/\psi \varLambda$ structure and observation of excited $\varXi^-$ states in the $\varXi_b^-\to J/\psi \varLambda K^-$ decay

abstract: First evidence of a structure in the $J/\psi \varLambda$ invariant mass distribution is obtained from an amplitude analysis of $\varXi_b^- \to J/\psi \varLambda K^-$ decays. The observed structure is consistent with being due to a charmonium pentaquark with strangeness. Its mass and width are determined to be $4458.8\pm2.9\,^{+4.7}_{-1.1}$MeV and $17.3\pm6.5\,^{+8.0}_{-5.7}$MeV, where the quoted uncertainties are statistical and systematic, respectively. The structure is also consistent with being due to two resonances. In addition, the narrow excited $\varXi^-$ states, $\varXi(1690)^-$ and $\varXi(1820)^-$, are seen for the first time in a $\varXi_b^-$ decay, and their masses and widths are measured with improved precision. The analysis is performed using $pp$ collision data corresponding to a total integrated luminosity of 9fb$^{-1}$, collected with the LHCb experiment at centre-of-mass energies of $7$, $8$ and $13$TeV.


title: Breaking Space Inversion-Symmetry to Obtain Asymmetric Spin-Wave Excitation in Systems with Nonuniform Magnetic Exchange

abstract: We report on the consequences of non-uniform exchange in magnetic systems. The quantum mechanical exchange interaction between spins is responsible for the phenomenon of magnetic order, and is generally considered to be uniform across bulk magnetic systems. Partly inspired by the Dzyaloshinskii-Moriya interaction--also known as antisymmetric exchange--we use a linearly varying exchange interaction along a magnetic strip as a route to spatial inversion symmetry-breaking. We find that, in addition to asymmetric modes and localization, spatially-varying exchange can be used to design nonreciprocal magnetic signal excitation at frequencies that are tunable. Moreover, our work predicts nonreciprocity to occur across a vast range of frequencies up to hundreds of GHz. Such spin wave engineering is a key area of ongoing research in the fields of magnonics and spintronics, which are expected to enable the next generation of communication technology. Analogous nonreciprocity is expected to occur in other wave systems with gradient properties.


title: Stabilization method with Relativistic Configuration-interaction applied to two-electron resonances

abstract: We applied a relativistic configuration-interaction (CI) framework to the stabilization method as an approach for obtaining the autoionization resonance structure of heliumlike ions. In this method, the ion is confined within an impenetrable spherical cavity, the size of which determines the radial space available for electron wavefunctions and electron-electron interactions. By varying the size of the cavity, one can obtain the autoionization resonance position and width. The applicability of this method is tested on the resonances of He atom while comparing with benchmark data available in the literature. The present method is further applied to the determination of the resonance structure of heliumlike uranium ion, where a relativistic framework is mandatory. In the strong-confinement region, the present method can be useful to simulate the properties of an atom or ion under extreme pressure. An exemplary application of the present method to determine the structure of ions embedded in a dense plasma environment is briefly discussed.


title: Frequency-modulated comb in a VECSEL

abstract: Optical frequency combs based on mode-locked lasers have revolutionized many areas of science and technology, such as precision metrology, optical frequency synthesis or telecommunications. In recent years, a particular kind of frequency comb has been observed in edge-emitting semiconductor lasers where the phase difference between longitudinal laser modes is fixed but not zero. This results in a linearly chirped output in the time domain with nearly constant intensity. Here, by using coherent beatnote spectroscopy, we show that such a comb regime can also exist in vertical-external-cavity surface-emitting lasers (VECSELs), as evidenced for a specific part of the laser spectrum. Our findings may not only lead to a better understanding of the physics of frequency-modulated combs but also enable comb applications with high optical power per comb line and flexible emission wavelengths.


title: Hybrid Genetic Search for the CVRP: Open-Source Implementation and SWAP* Neighborhood

abstract: The vehicle routing problem is one of the most studied combinatorial optimization topics, due to its practical importance and methodological interest. Yet, despite extensive methodological progress, many recent studies are hampered by the limited access to simple and efficient open-source solution methods. Given the sophistication of current algorithms, reimplementation is becoming a difficult and time-consuming exercise that requires extensive care for details to be truly successful. Against this background, we use the opportunity of this short paper to introduce a simple -- open-source -- implementation of the hybrid genetic search (HGS) specialized to the capacitated vehicle routing problem (CVRP). This state-of-the-art algorithm uses the same general methodology as Vidal et al. (2012) but also includes additional methodological improvements and lessons learned over the past decade of research. In particular, it includes an additional neighborhood called SWAP* which consists in exchanging two customers between different routes without an insertion in place. As highlighted in our study, an efficient exploration of SWAP* moves significantly contributes to the performance of local searches. Moreover, as observed in experimental comparisons with other recent approaches on the classical instances of Uchoa et al. (2017), HGS still stands as a leading metaheuristic regarding solution quality, convergence speed, and conceptual simplicity.


title: Large Scale Features of Southwest Monsoon During 2015

abstract: During 2015, the southwest monsoon (SWM) rainfall over the country remained deficient with seasonal rainfall of about 86% of the long period average (Table 1.1). Last year, the seasonal rainfall deficiency over the country as a whole was 12% (www.imd.gov.in). Thus, this is a fourth episode of two consecutive years, with deficient monsoon, similar to 1904-05, 1965-66 and 1986-87 (www.imd.gov.in).


title: Future Climate Change Projections over the Indian Region

abstract: Assessments of impacts of climate change and future projections over the Indian region, have so far relied on a single regional climate model (RCM) - eg., the PRECIS RCM of the Hadley Centre, UK. While these assessments have provided inputs to various reports (e.g., INCCA 2010; NATCOMM2 2012), it is important to have an ensemble of climate projections drawn from multiple RCMs due to large uncertainties in regional-scale climate projections. Ensembles of multi-RCM projections driven under different perceivable socio-economic scenarios are required to capture the probable path of growth, and provide the behavior of future climate and impacts on various biophysical systems and economic sectors dependent on such systems. The Centre for Climate Change Research, Indian Institute of Tropical Meteorology (CCCR-IITM) has generated an ensemble of high resolution downscaled projections of regional climate and monsoon over South Asia until 2100 for the Intergovernmental Panel for Climate Change (IPCC)using a RCM (ICTP-RegCM4) at 50 km horizontal resolution, by driving the regional model with lateral and lower boundary conditions from multiple global atmosphere-ocean coupled models from the Coupled Model Intercomparison Project Phase 5 (CMIP5). The future projections are based on three Representation Concentration Pathway (RCP) scenarios (viz., RCP2.6, RCP4.5, RCP8.5) of the IPCC.


title: Regional Climate Change Datasets for South Asia

abstract: The Centre for Climate Change Research (CCCR;http://cccr.tropmet.res.in) at the Indian Institute of Tropical Meteorology (IITM; http://www.tropmet.res.in), Pune, launched in 2009 with the support of the Ministry of Earth Sciences (MoES), Government of India, focuses on the development of new climate modelling capabilities in India and South Asia to address issues concerning the science of climate change. CCCR-IITM has the mandate of developing an Earth System Model and to make the regional climate projections. An important achievement was made by developing an Earth System Model at IITM, which is an important step towards understanding global and regional climate response to long-term climate variability and climate change. CCCR-IITM has also generated an ensemble of high resolution dynamically downscaled future projections of regional climate over South Asia and Indian monsoon, which are found useful for impact assessment studies and for quantifying uncertainties in the regional projections. A brief overview of these core climate change modeling activities of CCCR-IITM was presented in an Interim Report on Climate Change over India (available at http://cccr.tropmet.res.in/home/reports.jsp)


title: aw_nas: A Modularized and Extensible NAS framework

abstract: Neural Architecture Search (NAS) has received extensive attention due to its capability to discover neural network architectures in an automated manner. aw_nas is an open-source Python framework implementing various NAS algorithms in a modularized manner. Currently, aw_nas can be used to reproduce the results of mainstream NAS algorithms of various types. Also, due to the modularized design, one can simply experiment with different NAS algorithms for various applications with awnas (e.g., classification, detection, text modeling, fault tolerance, adversarial robustness, hardware efficiency, and etc.). Codes and documentation are available at https://github.com/walkerning/aw_nas.


title: Reinforcement Learning for Unified Allocation and Patrolling in Signaling Games with Uncertainty

abstract: Green Security Games (GSGs) have been successfully used in the protection of valuable resources such as fisheries, forests and wildlife. While real-world deployment involves both resource allocation and subsequent coordinated patrolling with communication and real-time, uncertain information, previous game models do not fully address both of these stages simultaneously. Furthermore, adopting existing solution strategies is difficult since they do not scale well for larger, more complex variants of the game models. We therefore first propose a novel GSG model that combines defender allocation, patrolling, real-time drone notification to human patrollers, and drones sending warning signals to attackers. The model further incorporates uncertainty for real-time decision-making within a team of drones and human patrollers. Second, we present CombSGPO, a novel and scalable algorithm based on reinforcement learning, to compute a defender strategy for this game model. CombSGPO performs policy search over a multi-dimensional, discrete action space to compute an allocation strategy that is best suited to a best-response patrolling strategy for the defender, learnt by training a multi-agent Deep Q-Network. We show via experiments that CombSGPO converges to better strategies and is more scalable than comparable approaches. Third, we provide a detailed analysis of the coordination and signaling behavior learnt by CombSGPO, showing group formation between defender resources and patrolling formations based on signaling and notifications between resources. Importantly, we find that strategic signaling emerges in the final learnt strategy. Finally, we perform experiments to evaluate these strategies under different levels of uncertainty.


title: Deep Learning and the Global Workspace Theory

abstract: Recent advances in deep learning have allowed Artificial Intelligence (AI) to reach near human-level performance in many sensory, perceptual, linguistic or cognitive tasks. There is a growing need, however, for novel, brain-inspired cognitive architectures. The Global Workspace theory refers to a large-scale system integrating and distributing information among networks of specialized modules to create higher-level forms of cognition and awareness. We argue that the time is ripe to consider explicit implementations of this theory using deep learning techniques. We propose a roadmap based on unsupervised neural translation between multiple latent spaces (neural networks trained for distinct tasks, on distinct sensory inputs and/or modalities) to create a unique, amodal global latent workspace (GLW). Potential functional advantages of GLW are reviewed, along with neuroscientific implications.


title: Analytical solutions of the cylindrical bending problem for the relaxed micromorphic continuum and other generalized continua (including full derivations)

abstract: We consider the cylindrical bending problem for an infinite plate as modelled with a family of generalized continuum models, including the micromorphic approach. The models allow to describe length scale effects in the sense that thinner specimens are comparatively stiffer. We provide the analytical solution for each case and exhibit the predicted bending stiffness. The relaxed micromorphic continuum shows bounded bending stiffness for arbitrary thin specimens, while classical micromorphic continuum or gradient elasticity as well as Cosserat models [35] exhibit unphysical unbounded bending stiffness for arbitrary thin specimens. This finding highlights the advantage of using the relaxed micromorphic model, which has a definite limit stiffness for small samples and which aids in identifying the relevant material parameters.


title: Emulation as an Accurate Alternative to Interpolation in Sampling Radiative Transfer Codes

abstract: Computationally expensive Radiative Transfer Models (RTMs) are widely used} to realistically reproduce the light interaction with the Earth surface and atmosphere. Because these models take long processing time, the common practice is to first generate a sparse look-up table (LUT) and then make use of interpolation methods to sample the multi-dimensional LUT input variable space. However, the question arise whether common interpolation methods perform most accurate. As an alternative to interpolation, this work proposes to use emulation, i.e., approximating the RTM output by means of statistical learning. Two experiments were conducted to assess the accuracy in delivering spectral outputs using interpolation and emulation: (1) at canopy level, using PROSAIL; and (2) at top-of-atmosphere level, using MODTRAN. Various interpolation (nearest-neighbour, inverse distance weighting, piece-wice linear) and emulation (Gaussian process regression (GPR), kernel ridge regression, neural networks) methods were evaluated against a dense reference LUT. In all experiments, the emulation methods clearly produced more accurate output spectra than classical interpolation methods. GPR emulation performed up to ten times more accurately than the best performing interpolation method, and this with a speed that is competitive with the faster interpolation methods. It is concluded that emulation can function as a fast and more accurate alternative to commonly used interpolation methods for reconstructing RTM spectral data.

categories: physics.ao-ph cs.LG physics.comp-ph stat.AP


title: A deep network approach to multitemporal cloud detection

abstract: We present a deep learning model with temporal memory to detect clouds in image time series acquired by the Seviri imager mounted on the Meteosat Second Generation (MSG) satellite. The model provides pixel-level cloud maps with related confidence and propagates information in time via a recurrent neural network structure. With a single model, we are able to outline clouds along all year and during day and night with high accuracy.


title: Deep Learning for Climate Model Output Statistics

abstract: Climate models are an important tool for the assessment of prospective climate change effects but they suffer from systematic and representation errors, especially for precipitation. Model output statistics (MOS) reduce these errors by fitting the model output to observational data with machine learning. In this work, we explore the feasibility and potential of deep learning with convolutional neural networks (CNNs) for MOS. We propose the CNN architecture ConvMOS specifically designed for reducing errors in climate model outputs and apply it to the climate model REMO. Our results show a considerable reduction of errors and mostly improved performance compared to three commonly used MOS approaches.


title: Transfer Learning with Convolutional Networks for Atmospheric Parameter Retrieval

abstract: The Infrared Atmospheric Sounding Interferometer (IASI) on board the MetOp satellite series provides important measurements for Numerical Weather Prediction (NWP). Retrieving accurate atmospheric parameters from the raw data provided by IASI is a large challenge, but necessary in order to use the data in NWP models. Statistical models performance is compromised because of the extremely high spectral dimensionality and the high number of variables to be predicted simultaneously across the atmospheric column. All this poses a challenge for selecting and studying optimal models and processing schemes. Earlier work has shown non-linear models such as kernel methods and neural networks perform well on this task, but both schemes are computationally heavy on large quantities of data. Kernel methods do not scale well with the number of training data, and neural networks require setting critical hyperparameters. In this work we follow an alternative pathway: we study transfer learning in convolutional neural nets (CNN s) to alleviate the retraining cost by departing from proxy solutions (either features or networks) obtained from previously trained models for related variables. We show how features extracted from the IASI data by a CNN trained to predict a physical variable can be used as inputs to another statistical method designed to predict a different physical variable at low altitude. In addition, the learned parameters can be transferred to another CNN model and obtain results equivalent to those obtained when using a CNN trained from scratch requiring only fine tuning.


title: Cloud detection machine learning algorithms for PROBA-V

abstract: This paper presents the development and implementation of a cloud detection algorithm for Proba-V. Accurate and automatic detection of clouds in satellite scenes is a key issue for a wide range of remote sensing applications. With no accurate cloud masking, undetected clouds are one of the most significant sources of error in both sea and land cover biophysical parameter retrieval. The objective of the algorithms presented in this paper is to detect clouds accurately providing a cloud flag per pixel. For this purpose, the method exploits the information of Proba-V using statistical machine learning techniques to identify the clouds present in Proba-V products. The effectiveness of the proposed method is successfully illustrated using a large number of real Proba-V images.


title: A global Canopy Water Content product from AVHRR/Metop

abstract: Spatially and temporally explicit canopy water content (CWC) data are important for monitoring vegetation status, and constitute essential information for studying ecosystem-climate interactions. Despite many efforts there is currently no operational CWC product available to users. In the context of the Satellite Application Facility for Land Surface Analysis (LSA-SAF), we have developed an algorithm to produce a global dataset of CWC based on data from the Advanced Very High Resolution Radiometer (AVHRR) sensor on board Meteorological Operational (MetOp) satellites forming the EUMETSAT Polar System (EPS). CWC reflects the water conditions at the leaf level and information related to canopy structure. An accuracy assessment of the EPS/AVHRR CWC indicated a close agreement with multi-temporal ground data from SMAPVEX16 in Canada and Dahra in Senegal. The present study further evaluates the consistency of the LSA-SAF product with respect to the Simplified Level 2 Product Prototype Processor (SL2P) product, and demonstrates its applicability at different spatio-temporal resolutions using optical data from MSI/Sentinel-2 and MODIS/Terra and Aqua. We conclude that the EPS/AVHRR CWC product is a promising tool for monitoring vegetation water status at regional and global scales.


title: Hadamard

abstract: This is the Hadamard package for Macaulay2 which computes the Hadamard product of projective subvarieties.


title: Micromechanical description of the compaction of soft pentagon assemblies

abstract: We analyze the isotropic compaction of assemblies composed of soft pentagons interacting through classical Coulomb friction via numerical simulations. The effect of the initial particle shape is discussed by comparing packings of pentagons with packings of soft circular particles. We characterize the evolution of the packing fraction, the elastic modulus, and the microstructure (particle rearrangement, connectivity, contact force and particle stress distributions) as a function of the applied stresses. Both systems behave similarly; the packing fraction increases and tends asymptotically to a maximum value $\phi_{max}$, where the bulk modulus diverges. At the microscopic scale we show that particle rearrangements occur even beyond the jammed state, the mean coordination increases as a square root of the packing fraction and, the force and stress distributions become more homogeneous as the packing fraction increases. Soft pentagons present larger particle rearrangements than circular ones, and such behavior decreases proportionally to the friction. Interestingly, the friction between particles also contributes to a better homogenization of the contact force network in both systems. From the expression of the granular stress tensor, we develop a model that describes the compaction behavior as a function of the applied pressure, the Young modulus and the initial shape of the particles. This model, settled on the joint evolution of the particle connectivity and the contact stress, provides outstanding predictions from the jamming point up to very high densities.


title: Trademark filings and patent application count time series are structurally near-identical and cointegrated: Implications for studies in innovation

abstract: Through time series analysis, this paper empirically explores, confirms and extends the trademark/patent inter-relationship as proposed in the normative intellectual-property (IP)-oriented Innovation Agenda view of the science and technology (S&T) firm. Beyond simple correlation, it is shown that trademark-filing (Trademarks) and patent-application counts (Patents) have similar (if not, identical) structural attributes (including similar distribution characteristics and seasonal variation, cross-wavelet synchronicity/coherency (short-term cross-periodicity) and structural breaks) and are cointegrated (integration order of 1) over a period of approximately 40 years (given the monthly observations). The existence of cointegration strongly suggests a "long-run" equilibrium between the two indices; that is, there is (are) exogenous force(s) restraining the two indices from diverging from one another. Structural breakpoints in the chrono-dynamics of the indices supports the existence of potentially similar exogeneous forces(s), as the break dates are simultaneous/near-simultaneous (Trademarks: 1987, 1993, 1999, 2005, 2011; Patents: 1988, 1994, 2000, and 2011). A discussion of potential triggers (affecting both time series) causing these breaks, and the concept of equilibrium in the context of these proxy measures are presented. The cointegration order and structural co-movements resemble other macro-economic variables, stoking the opportunity of using econometrics approaches to further analyze these data. As a corollary, this work further supports the inclusion of trademark analysis in innovation studies. Lastly, the data and corresponding analysis tools (R program) are presented as Supplementary Materials for reproducibility and convenience to conduct future work for interested readers.


title: Using colloidal deposition to mobilize immiscible fluids from porous media

abstract: Colloidal particles hold promise for mobilizing and removing trapped immiscible fluids from porous media, with implications for key energy and water applications. Most studies focus on accomplishing this goal using particles that can localize at the immiscible fluid interface. Therefore, researchers typically seek to optimize the surface activity of particles, as well as their ability to freely move through a pore space with minimal deposition onto the surrounding solid matrix. Here, we demonstrate that deposition can, surprisingly, promote mobilization of a trapped fluid from a porous medium without requiring any surface activity. Using confocal microscopy, we directly visualize both colloidal particles and trapped immiscible fluid within a transparent, three-dimensional (3D) porous medium. We find that as non-surface active particles deposit on the solid matrix, increasing amounts of trapped fluid become mobilized. We unravel the underlying physics by analyzing the extent of deposition, as well as the geometry of trapped fluid droplets, at the pore scale: deposition increases the viscous stresses on trapped droplets, overcoming the influence of capillarity that keeps them trapped. Given an initial distribution of trapped fluid, this analysis enables us to predict the extent of fluid mobilized through colloidal deposition. Taken together, our work reveals a new way by which colloids can be harnessed to mobilize trapped fluid from a porous medium.

categories: cond-mat.soft cond-mat.mes-hall cond-mat.mtrl-sci physics.flu-dyn