title: Non Projected Calabi-Yau Supermanifolds over $\mathbb{P}^2$

abstract: We start a systematic study of non-projected supermanifolds, concentrating on supermanifolds with fermionic dimension 2 and with the reduced manifold a complex projective space. We show that all the non-projected supermanifolds of dimension $2|2$ over $\mathbb{P}^2$ are completely characterised by a non-zero 1-form $\omega$ and by a locally free sheaf $\mathcal{F}$ of rank $0|2$, satisfying $Sym^2 \mathcal{F} \cong K_{\mathbb{P}^2}$. Denoting such supermanifolds with $\mathbb{P}^{2}_\omega(\mathcal{F})$, we show that all of them are Calabi-Yau supermanifolds and, when $\omega \neq 0$, they are non-projective, that is they cannot be embedded into any projective superspace $\mathbb{P}^{n|m}$. Instead, we show that every non-projected supermanifolds over $\mathbb{P}^2$ admits an embedding into a super Grassmannian. By contrast, we give an example of a supermanifold $\mathbb P^{2}_\omega(\mathcal F)$ that cannot be embedded in any of the $\Pi$-projective superspaces $\mathbb P^{n}_{\Pi}$ introduced by Manin and Deligne. However, we also show that when $\mathcal F$ is the cotangent bundle over $\mathbb{P}^2$, then the non-projected $\mathbb{P}^2_\omega(\mathcal F)$ and the $\Pi$-projective plane $\mathbb P^{2}_{\Pi}$ do coincide.


title: Witnessing a Large-scale Slipping Magnetic Reconnection along a Dimming Channel during a Solar Flare

abstract: We report the intriguing large-scale dynamic phenomena associated with the M6.5 flare~(SOL2015-06-22T18:23) in NOAA active region 12371, observed by RHESSI, Fermi, and the Atmospheric Image Assembly (AIA) and Magnetic Imager (HMI) on the Solar Dynamic Observatory (SDO). The most interesting feature of this event is a third ribbon (R3) arising in the decay phase, propagating along a dimming channel (seen in EUV passbands) towards a neighboring sunspot. The propagation of R3 occurs in the presence of hard X-ray footpoint emission, and is broadly visible at temperatures from 0.6~MK to over 10~MK through the Differential Emission Measure (DEM) analysis. The coronal loops then undergo an apparent slipping motion following the same path of R3, after a $\sim$ 80~min delay. To understand the underlying physics, we investigate the magnetic configuration and the thermal structure of the flaring region. Our results are in favor of a slipping-type reconnection followed by the thermodynamic evolution of coronal loops. In comparison with those previously reported slipping reconnection events, this one proceeds across a particularly long distance ($\sim$60 Mm) over a long period of time ($\sim$50 min), and shows two clearly distinguished phases: the propagation of the footpoint brightening driven by nonthermal particle injection and the apparent slippage of loops governed by plasma heating and subsequent cooling.


title: On the rationality problem for quadric bundles

abstract: We classify all positive integers n and r such that (stably) non-rational complex r-fold quadric bundles over rational n-folds exist. We show in particular that for any n and r, a wide class of smooth r-fold quadric bundles over projective n-space are not stably rational if r lies in the interval from $2^{n-1}-1$ to $2^{n}-2$. In our proofs we introduce a generalization of the specialization method of Voisin and Colliot-Th\'el\`ene--Pirutka which avoids universally $CH_0$-trivial resolutions of singularities.


title: Characterization of multivariate Bernoulli distributions with given margins

abstract: We express each Fr\'echet class of multivariate Bernoulli distributions with given margins as the convex hull of a set of densities, which belong to the same Fr\'echet class. This characterisation allows us to establish whether a given correlation matrix is compatible with the assigned margins and, if it is, to easily construct one of the corresponding joint densities. % Such %representation is based on a polynomial expression of the distributions of a Fr\'echet class. We reduce the problem of finding a density belonging to a Fr\'echet class and with given correlation matrix to the solution of a linear system of equations. Our methodology also provides the bounds that each correlation must satisfy to be compatible with the assigned margins. An algorithm and its use in some examples is shown.


title: Quadric surface bundles over surfaces and stable rationality

abstract: We prove a general specialization theorem which implies stable irrationality for a wide class of quadric surface bundles over rational surfaces. As an application, we solve with the exception of two cases, the stable rationality problem for any very general complex projective quadric surface bundle over the projective plane, given by a symmetric matrix of homogeneous polynomials. Both exceptions degenerate over a plane sextic curve and the corresponding double cover is a K3 surface.


title: Supergeometry of $\Pi$-Projective Spaces

abstract: In this paper we prove that $\Pi$-projective spaces $\mathbb{P}^n_\Pi$ arise naturally in supergeometry upon considering a non-projected thickening of $\mathbb{P}^n$ related to the cotangent sheaf $\Omega^1_{\mathbb{P}^n}$. In particular, we prove that for $n \geq 2$ the $\Pi$-projective space $\mathbb{P}^n_\Pi$ can be constructed as the non-projected supermanifold determined by three elements $(\mathbb{P}^n, \Omega^1_{\mathbb{P}^n}, \lambda)$, where $\mathbb{P}^n$ is the ordinary complex projective space, $\Omega^1_{\mathbb{P}^n}$ is its cotangent sheaf and $\lambda $ is a non-zero complex number, representative of the fundamental obstruction class $\omega \in H^1 (\mathcal{T}_{\mathbb{P}^n} \otimes \bigwedge^2 \Omega^1_{\mathbb{P}^n}) \cong \mathbb{C}.$ Likewise, in the case $n=1$ the $\Pi$-projective line $\mathbb{P}^1_\Pi$ is the split supermanifold determined by the pair $(\mathbb{P}^1, \Omega^1_{\mathbb{P}^1} \cong \mathcal{O}_{\mathbb{P}^1} (-2)).$ Moreover we show that in any dimension $\Pi$-projective spaces are Calabi-Yau supermanifolds. To conclude, we offer pieces of evidence that, more in general, also $\Pi$-Grassmannians can be constructed the same way using the cotangent sheaf of their underlying reduced Grassmannians, provided that also higher, possibly fermionic, obstruction classes are taken into account. This suggests that this unexpected connection with the cotangent sheaf is characteristic of $\Pi$-geometry.


title: Cosmological Asymptotics in Higher-Order Gravity Theories

abstract: We study the early-time behavior of isotropic and homogeneous solutions in vacuum as well as radiation-filled cosmological models in the full, effective, four dimensional gravity theory with higher derivatives. We use asymptotic methods to analyze all possible ways of approach to the initial singularity of such universes. In order to do so, we construct autonomous dynamical systems that describe the evolution of these models, and decompose the associated vector fields. We prove that, at early times, all flat vacua as well as general curved ones are globally attracted by the "universal" square root scaling solution. Open vacua, on the other hand show in both, future and past directions a dominant asymptotic approach to horizon-free, Milne states that emerge from initial data sets of smaller dimension. Closed universes exhibit more complex logarithmic singularities. Our results on asymptotic stability show a possible relation to cyclic and ekpyrotic cosmologies at the passage through the singularity. In the case of radiation-filled universes of the same class we show the essential uniqueness and stability of the resulting asymptotic scheme, once more dominated by $t^{1/2}$, in all cases except perhaps that of the conformally invariant Bach-Weyl gravity. In all cases, we construct a formal series representation valid near the initial singularity of the general solution of these models and prove that curvature as well as radiation play a subdominant role in the dominating form. A discussion is also made on the implications of these results for the generic initial state of the theory.


title: An Integrated Quadratic Reconstruction for Finite Volume Schemes to Scalar Conservation Laws in Multiple Dimensions

abstract: We proposed a piecewise quadratic reconstruction method in multiple dimensions, which is in an integrated style, for finite volume schemes to scalar conservation laws. This integrated quadratic reconstruction is parameter-free and applicable on flexible grids. We show that the finite volume schemes with the new reconstruction satisfy a local maximum principle with properly setup on time steplength. Numerical examples are presented to show that the proposed scheme attains a third-order accuracy for smooth solutions in both 2D and 3D cases. It is indicated by numerical results that the local maximum principle is helpful to prevent overshoots in numerical solutions.


title: The Geometry of Nodal Sets and Outlier Detection

abstract: Let $(M,g)$ be a compact manifold and let $-\Delta \phi_k = \lambda_k \phi_k$ be the sequence of Laplacian eigenfunctions. We present a curious new phenomenon which, so far, we only managed to understand in a few highly specialized cases: the family of functions $f_N:M \rightarrow \mathbb{R}_{\geq 0}$ $$ f_N(x) = \sum_{k \leq N}{ \frac{1}{\sqrt{\lambda_k}} \frac{|\phi_k(x)|}{\|\phi_k\|_{L^{\infty}(M)}}}$$ seems strangely suited for the detection of anomalous points on the manifold. It may be heuristically interpreted as the sum over distances to the nearest nodal line and potentially hints at a new phenomenon in spectral geometry. We give rigorous statements on the unit square $[0,1]^2$ (where minima localize in $\mathbb{Q}^2$) and on Paley graphs (where $f_N$ recovers the geometry of quadratic residues of the underlying finite field $\mathbb{F}_p$). Numerical examples show that the phenomenon seems to arise on fairly generic manifolds.

categories: math.SP math-ph math.AP math.FA math.MP stat.ML


title: Stochastic Integration and Stochastic PDEs Driven by Jumps on the Dual of a Nuclear Space

abstract: We develop a novel theory of weak and strong stochastic integration for cylindrical martingale-valued measures taking values in the dual of a nuclear space. This is applied to develop a theory of SPDEs with rather general coefficients. In particular, we can then study SPDEs driven by general L\'{e}vy processes in this context.


title: On the Catalysis of the Electroweak Vacuum Decay by Black Holes at High Temperature

abstract: We study the effect of primordial black holes on the classical rate of nucleation of AdS regions within the standard electroweak vacuum at high temperature. We find that the energy barrier for transitions to the new vacuum, which determines the exponential suppression of the nucleation rate, can be reduced significantly, or even eliminated completely, in the black-hole background if the Standard Model Higgs is coupled to gravity through the renormalizable term $\xi {\cal R} h^2$.


title: Synchronization and separation in the Johnson schemes

abstract: Recently Peter Keevash solved asymptotically the existence question for Steiner systems by showing that $S(t,k,n)$ exists whenever the necessary divisibility conditions on the parameters are satisfied and $n$ is sufficiently large in terms of $k$ and $t$. The purpose of this paper is to make a conjecture which if true would be a significant extension of Keevash's theorem, and to give some theoretical and computational evidence for the conjecture. We phrase the conjecture in terms of the notions (which we define here) of synchronization and separation for association schemes. These definitions are based on those for permutation groups which grow out of the theory of synchronization in finite automata. In this theory, two classes of permutation groups (called \emph{synchronizing} and \emph{separating}) lying between primitive and $2$-homogeneous are defined. A big open question is how the permutation group induced by $S_n$ on $k$-subsets of $\{1,\ldots,n\}$ fits in this hierarchy; our conjecture would give a solution to this problem for $n$ large in terms of $k$. We prove the conjecture in the case $k=4$: our result asserts that $S_n$ acting on $4$-sets is separating for $n\ge10$ (it fails to be synchronizing for $n=9$).


title: Detection of a Substantial Molecular Gas Reservoir in a brightest cluster galaxy at z = 1.7

abstract: We report the detection of CO(2-1) emission coincident with the brightest cluster galaxy (BCG) of the high-redshift galaxy cluster SpARCS1049+56, with the Redshift Search Receiver (RSR) on the Large Millimetre Telescope (LMT). We confirm a spectroscopic redshift for the gas of z = 1.7091+/-0.0004, which is consistent with the systemic redshift of the cluster galaxies of z = 1.709. The line is well-fit by a single component Gaussian with a RSR resolution-corrected FWHM of 569+/-63 km/s. We see no evidence for multiple velocity components in the gas, as might be expected from the multiple image components seen in near-infrared imaging with the Hubble Space Telescope. We measure the integrated flux of the line to be 3.6+/-0.3 Jy km/s and, using alpha_CO = 0.8 Msun (K km s^-1 pc^2)^-1 we estimate a total molecular gas mass of 1.1+/-0.1x10^11 Msun and a M_H2/M_star ~ 0.4. This is the largest gas reservoir detected in a BCG above z > 1 to date. Given the infrared-estimated star formation rate of 860+/-130 Msun/yr, this corresponds to a gas depletion timescale of ~0.1Gyr. We discuss several possible mechanisms for depositing such a large gas reservoir to the cluster center -- e.g., a cooling flow, a major galaxy-galaxy merger or the stripping of gas from several galaxies -- but conclude that these LMT data are not sufficient to differentiate between them.


title: Symmetric cohomology of groups

abstract: We investigate the relationship between the symmetric, exterior and classical cohomologies of groups. The first two theories were introduced respectively by Staic and Zarelua. We show in particular, that there is a map from exterior cohomology to symmetric cohomology which is a split monomorphism in general and an isomorphism in many cases, but not always. We introduce two spectral sequences which help to explain the realtionship between these cohomology groups. As a sample application we obtain that symmetric and classical cohomologies are isomorphic for torsion free groups.


title: Chemical complexity induced by efficient ice evaporation in the Barnard 5 molecular cloud

abstract: Cold gas-phase water has recently been detected in a cold dark cloud, Barnard 5 located in the Perseus complex, by targeting methanol peaks as signposts for ice mantle evaporation. Observed morphology and abundances of methanol and water are consistent with a transient non-thermal evaporation process only affecting the outermost ice mantle layers, possibly triggering a more complex chemistry. We present the detection of the Complex Organic Molecules (COMs) acetaldehyde and methyl formate as well as formic acid and ketene, and the tentative detection of di-methyl ether towards the methanol hotspot of Barnard 5 located between two dense cores using the single dish OSO 20m, IRAM 30m, and NRO 45m telescopes. The high energy cis- conformer of formic acid is detected, suggesting that formic acid is mostly formed at the surface of interstellar grains and then evaporated. The detection of multiple transitions for each species allows us to constrain their abundances through LTE and non-LTE methods. All the considered COMs show similar abundances between $\sim 1$ and $\sim 10$ % relative to methanol depending on the assumed excitation temperature. The non-detection of glycolaldehyde, an isomer of methyl formate, with a [glycolaldehyde]/[methyl formate] abundance ratio lower than 6 %, favours gas phase formation pathways triggered by methanol evaporation. According to their excitation temperatures derived in massive hot cores, formic acid, ketene, and acetaldehyde have been designated as "lukewarm" COMs whereas methyl formate and di-methyl ether were defined as "warm" species. Comparison with previous observations of other types of sources confirms that "lukewarm" and "warm" COMs show similar abundances in low-density cold gas whereas the "warm" COMs tend to be more abundant than the "lukewarm" species in warm protostellar cores.


title: Renyi entropies and C_T for higher derivative free scalars and spinors on even spheres

abstract: General expressions for the R\'enyi entropies and central charges for higher derivative free spinors and scalars on even spheres are obtained using a direct spectral method on a compact lune division of the sphere. Formulae and numbers are rapidly obtained for any dimension and order of derivative. The relation between the conformal anomaly and the hyperbolic free energy is briefly explored using standard expansions. A field theoretic derivation of the central charge formula for higher derivative scalars in any (even) dimension, given by Osborn and Stergiou and by Gliozzi {\it et al}, is thereby provided. The corresponding extension to spinors is made. Generalised Bernoulli polynomials play an important technical role.


title: Powerful outflows in the central parsecs of the low-luminosity Active Galactic Nucleus NGC 1386

abstract: Low-luminosity Active Galactic Nuclei, i.e. L_bol/L_edd ~ 10^-6 - 10^-3, constitute the bulk population of Active Galactic Nuclei (AGNs). Powerful jets, common in these objects, are a crucial source of feedback energy driving mass outflows into the host galaxy and the intergalactic medium. This paper reports the first direct measurement of powerful mass outflows traced by the forbidden high ionization gas in the low luminosity AGN NGC1386 at scales of a few parsecs from the central engine. The high angular resolution of the data allows us to directly measure the location, morphology and kinematic of the outflow. This the form of two symmetrical expanding hot gas shells moving in opposite directions along the line of sight. The co-spatiality of the gas shells with radio emission seen at the same parsec scales and with X-rays indicates that this is a shock-driven outflow induced by an incipient core-jet. With a minimum number of assumptions, we derive a mass outflow rate of 11 solar masses/yr, comparable to those of powerful AGN. The result has strong implications in the global accounting of feedback mass and energy driven by a low-luminosity AGN into the medium and the corresponding galaxy evolution.


title: Intersections of three quadrics in $\mathbb{P}^7$

abstract: We study rationality properties of smooth complete intersections of three quadrics in $\mathbb{P}^7$. We exhibit a smooth family of such intersections with both rational and non-rational fibers.


title: Sparse Tableau Formulation for Optimal Power Flow Applications

abstract: Typical formulations of the optimal power flow (OPF) problem rely on what is termed the "bus-branch" model, with network electrical behavior summarized in the Ybus admittance matrix. From a circuit perspective, this admittance representation restricts network elements to be voltage controlled and limitations of the Ybus have long been recognized. A fixed Ybus is unable to represent an ideal circuit breaker, and more subtle limitations appear in transformer modeling. In power systems parlance, more detailed approaches to overcome these limitations are termed "node-breaker" representations, but these are often cumbersome, and are not widely utilized in OPF. This paper develops a general network representation adapted to the needs of OPF, based on the Sparse Tableau Formulation (STF) with following advantages for OPF: (i) conceptual clarity in formulating constraints, allowing a comprehensive set of network electrical variables; (ii) improved fidelity in capturing physical behavior and engineering limits; (iii) added flexibility in optimization solution, in that elimination of intermediate variables is left to the optimization algorithm. The STF is then applied to OPF numerical case studies which demonstrate that the STF shows little or no penalty in computational speed compared to classic OPF representations, and sometimes provides considerable advantage in computational speed.


title: Miscible-Immiscible Transition and Nonequilibrium Scaling in Two-Component Driven Open Condensate Wires

abstract: We investigate the steady state phase diagram of two-component driven open condensates in one dimension. We identify a miscible-immiscible transition which is predominantly driven by gapped density fluctuations and occurs upon increasing the inter-component inelastic coupling. Below the transition in the miscible phase, we find the effective long wavelength dynamics to be described by a two-component Kardar-Parisi-Zhang (KPZ) equation that belongs to the nonequilibrium universality class of the one-dimensional single-component KPZ equation at generic choices of parameters. Our results are relevant for different experimental realizations for two-component driven open condensates in exciton-polariton systems.


title: Self-consistent modelling of our Galaxy with Gaia data

abstract: Galaxy models are fundamental to exploiting surveys of our Galaxy. There is now a significant body of work on axisymmetric models. A model can be defined by giving the DF of each major class of stars and of dark matter. Then the self-consistent gravitational potential is determined. Other modelling techniques are briefly considered before an overview of some early work on non-axisymmetric models.


title: Effective Chern-Simons actions of particles coupled to 3D gravity

abstract: Point particles in 3D gravity are known to behave as topological defects, while gravitational field can be expressed as the Chern-Simons theory of the appropriate local isometry group of spacetime. In the case of the Poincar\'e group, integrating out the gravitational degrees of freedom it is possible to obtain the effective action for particle dynamics. We review the known results, both for single and multiple particles, and attempt to extend this approach to the (anti-)de Sitter group, using the factorizations of isometry groups into the double product of the Lorentz group and AN(2) group. On the other hand, for the de Sitter group one can also perform a contraction to the semidirect product of AN(2) and the translation group. The corresponding effective action curiously describes a Carrollian particle with the AN(2) momentum space. We derive this contraction in a more rigorous manner and further explore its properties, including a generalization to the multiparticle case.


title: Optimal Design Method of MIMO Antenna Directivities and Corresponding Current Distributions by Using Spherical Mode Expansion

abstract: This paper proposes a new methodology to design optimal antennas for MIMO (Multi-Input Multi-Output) communication systems by using spherical mode expansion. Given spatial channel properties of a MIMO channel, such as the angular profile at both sides, the optimal MIMO antennas should provide the largest channel capacity with a constraint of the limited implementation space (volume). In designing a conventional MIMO antenna, first the antenna structure (current distribution) is determined, second antenna directivity is calculated based on the current distribution, and thirdly MIMO channel capacity is calculated by using given angular profiles and obtained antenna directivity. This process is repeated by adjusting the antenna structure until the performance satisfies a predefined threshold. To the contrary, this paper solves the optimization problem analytically and finally gives near optimal antenna structure (current distribution) without any greedy search. In the proposed process, first the optimal directivity of MIMO antennas is derived by applying spherical mode expansion to the angular profiles, and second a far-near field conversion is applied on the derived optimal directivity to achieve near optimal current distributions on a limited surface. The effectiveness of the proposed design methodology is validated via numerical calculation of MIMO channel capacity as in the conventional design method while giving near optimal current distribution with constraint of an antenna structure derived from proposed methodology.


title: A note on the Prodi-Serrin conditions for the regularity of a weak solution to the Navier-Stokes equations

abstract: The paper is concerned with the regularity of weak solutions to the Navier-Stokes equations. The aim is to investigate on a relaxed Prodi-Serrin condition in order to obtain regularity for t > 0. The most interesting aspect of the result is that no compatibility condition is required to the initial data $v_0\in J^2(\OO) J2({\Omega})$.


title: Observation of the Black Widow B1957+20 millisecond pulsar binary system with the MAGIC telescopes

abstract: B1957+20 is a millisecond pulsar located in a black widow type compact binary system with a low mass stellar companion. The interaction of the pulsar wind with the companion star wind and/or the interstellar plasma is expected to create plausible conditions for acceleration of electrons to TeV energies and subsequent production of very high energy {\gamma} rays in the inverse Compton process. We performed extensive observations with the MAGIC telescopes of B1957+20. We interpret results in the framework of a few different models, namely emission from the vicinity of the millisecond pulsar, the interaction of the pulsar and stellar companion wind region, or bow shock nebula. No significant steady very high energy {\gamma}-ray emission was found. We derived a 95% confidence level upper limit of 3.0 x 10 -12 cm -2 s -1 on the average {\gamma}-ray emission from the binary system above 200 GeV. The upper limits obtained with MAGIC constrain, for the first time, different models of the high-energy emission in B1957+20. In particular, in the inner mixed wind nebula model with mono-energetic injection of electrons, the acceleration efficiency of electrons is constrained to be below ~(2-10)% of the pulsar spin down power. For the pulsar emission, the obtained upper limits for each emission peak are well above the exponential cut-off fits to the Fermi-LAT data, extrapolated to energies above 50 GeV. The MAGIC upper limits can rule out a simple power-law tail extension through the sub-TeV energy range for the main peak seen at radio frequencies.


title: Rare $\Lambda_b\to n l^+l^-$ decays in the relativistic quark-diquark picture

abstract: The form factors of the rare $\Lambda_b\to n l^+l^-$ decays are calculated in the framework of the relativistic quark-diquark picture of baryons with the consistent account of the relativistic effects. Their momentum transfer squared dependence is determined explicitly in the whole accessible kinematical range. The decay branching fractions, forward-backward asymmetries and the fractions of longitudinally polarized dileptons are determined. The branching fraction of the rare $\Lambda_b\to n\mu^+\mu^-$ decay are found to be ${\rm Br}(\Lambda_b\to n\mu^+\mu^-)=(3.75\pm0.38)\times10^{-8}$ and thus could be measured at the LHC. Prediction for the branching fraction of the rare radiative $\Lambda_b\to n \gamma$ decay is also given.


title: Automatic Response Assessment in Regions of Language Cortex in Epilepsy Patients Using ECoG-based Functional Mapping and Machine Learning

abstract: Accurate localization of brain regions responsible for language and cognitive functions in Epilepsy patients should be carefully determined prior to surgery. Electrocorticography (ECoG)-based Real Time Functional Mapping (RTFM) has been shown to be a safer alternative to the electrical cortical stimulation mapping (ESM), which is currently the clinical/gold standard. Conventional methods for analyzing RTFM signals are based on statistical comparison of signal power at certain frequency bands. Compared to gold standard (ESM), they have limited accuracies when assessing channel responses. In this study, we address the accuracy limitation of the current RTFM signal estimation methods by analyzing the full frequency spectrum of the signal and replacing signal power estimation methods with machine learning algorithms, specifically random forest (RF), as a proof of concept. We train RF with power spectral density of the time-series RTFM signal in supervised learning framework where ground truth labels are obtained from the ESM. Results obtained from RTFM of six adult patients in a strictly controlled experimental setup reveal the state of the art detection accuracy of $\approx 78\%$ for the language comprehension task, an improvement of $23\%$ over the conventional RTFM estimation method. To the best of our knowledge, this is the first study exploring the use of machine learning approaches for determining RTFM signal characteristics, and using the whole-frequency band for better region localization. Our results demonstrate the feasibility of machine learning based RTFM signal analysis method over the full spectrum to be a clinical routine in the near future.


title: Accurate thermal conductivities from optimally short molecular dynamics simulations

abstract: The evaluation of transport coefficients in extended systems, such as thermal conductivity or shear viscosity, is known to require impractically long simulations, thus calling for a paradigm shift that would allow to deploy state-of-the-art quantum simulation methods. We introduce a new method to compute these coefficients from optimally short molecular dynamics simulations, based on the Green-Kubo theory of linear response and the cepstral analysis of time series. Information from the full sample power spectrum of the relevant current for a single and relatively short trajectory is leveraged to evaluate and optimally reduce the noise affecting its zero-frequency value, whose expectation is proportional to the corresponding conductivity. Our method is unbiased and consistent, in that both the resulting bias and statistical error can be made arbitrarily small in the long-time limit. A simple data-analysis protocol is proposed and validated with the calculation of thermal conductivities in the paradigmatic cases of elemental and molecular fluids (liquid Ar and H$_2$O) and of crystalline and glassy solids (MgO and a-SiO$_2$). We find that simulation times of one to a few hundred picoseconds are sufficient in these systems to achieve an accuracy of the order of 10% on the estimated thermal conductivities.


title: Neuro-RAM Unit with Applications to Similarity Testing and Compression in Spiking Neural Networks

abstract: We study distributed algorithms implemented in a simplified biologically inspired model for stochastic spiking neural networks. We focus on tradeoffs between computation time and network complexity, along with the role of randomness in efficient neural computation. It is widely accepted that neural computation is inherently stochastic. In recent work, we explored how this stochasticity could be leveraged to solve the `winner-take-all' leader election task. Here, we focus on using randomness in neural algorithms for similarity testing and compression. In the most basic setting, given two $n$-length patterns of firing neurons, we wish to distinguish if the patterns are equal or $\epsilon$-far from equal. Randomization allows us to solve this task with a very compact network, using $O \left (\frac{\sqrt{n}\log n}{\epsilon}\right)$ auxiliary neurons, which is sublinear in the input size. At the heart of our solution is the design of a $t$-round neural random access memory, or indexing network, which we call a neuro-RAM. This module can be implemented with $O(n/t)$ auxiliary neurons and is useful in many applications beyond similarity testing. Using a VC dimension-based argument, we show that the tradeoff between runtime and network size in our neuro-RAM is nearly optimal. Our result has several implications -- since our neuro-RAM can be implemented with deterministic threshold gates, it shows that, in contrast to similarity testing, randomness does not provide significant computational advantages for this problem. It also establishes a separation between feedforward networks whose gates spike with sigmoidal probability functions, and well-studied deterministic sigmoidal networks, whose gates output real number sigmoidal values, and which can implement a neuro-RAM much more efficiently.


title: Sparse Stochastic Bandits

abstract: In the classical multi-armed bandit problem, d arms are available to the decision maker who pulls them sequentially in order to maximize his cumulative reward. Guarantees can be obtained on a relative quantity called regret, which scales linearly with d (or with sqrt(d) in the minimax sense). We here consider the sparse case of this classical problem in the sense that only a small number of arms, namely s < d, have a positive expected reward. We are able to leverage this additional assumption to provide an algorithm whose regret scales with s instead of d. Moreover, we prove that this algorithm is optimal by providing a matching lower bound - at least for a wide and pertinent range of parameters that we determine - and by evaluating its performance on simulated data.


title: Time-Varying Formation Controllers for Unmanned Aerial Vehicles Using Deep Reinforcement Learning

abstract: We consider the problem of designing scalable and portable controllers for unmanned aerial vehicles (UAVs) to reach time-varying formations as quickly as possible. This brief confirms that deep reinforcement learning can be used in a multi-agent fashion to drive UAVs to reach any formation while taking into account optimality and portability. We use a deep neural network to estimate how good a state is, so the agent can choose actions accordingly. The system is tested with different non-high-dimensional sensory inputs without any change in the neural network architecture, algorithm or hyperparameters, just with additional training.


title: Distinguishing Spin-Aligned and Isotropic Black Hole Populations With Gravitational Waves

abstract: The first direct detections of gravitational waves from merging binary black holes open a unique window into the binary black hole formation environment. One promising environmental signature is the angular distribution of the black hole spins; systems formed through dynamical interactions among already-compact objects are expected to have isotropic spin orientations whereas binaries formed from pairs of stars born together are more likely to have spins preferentially aligned with the binary orbital angular momentum. We consider existing gravitational wave measurements of the binary effective spin, the best-measured combination of spin parameters, in the four likely binary black hole detections GW150914, LVT151012, GW151226, and GW170104. If binary black hole spin magnitudes extend to high values we show that the data exhibit a $2.4\sigma$ ($0.015$ odds ratio) preference for an isotropic angular distribution over an aligned one. By considering the effect of 10 additional detections, we show that such an augmented data set would enable in most cases a preference stronger than $5\sigma$ ($2.9 \times 10^{-7}$ odds ratio). The existing preference for either an isotropic spin distribution or low spin magnitudes for the observed systems will be confirmed (or overturned) confidently in the near future.


title: Li-intercalated Graphene on SiC(0001): an STM study

abstract: We present a systematical study via scanning tunneling microscopy (STM) and low-energy electron diffraction (LEED) on the effect of the exposure of Lithium (Li) on graphene on silicon carbide (SiC). We have investigated Li deposition both on epitaxial monolayer graphene and on buffer layer surfaces on the Si-face of SiC. At room temperature, Li immediately intercalates at the interface between the SiC substrate and the buffer layer and transforms the buffer layer into a quasi-free-standing graphene. This conclusion is substantiated by LEED and STM evidence. We show that intercalation occurs through the SiC step sites or graphene defects. We obtain a good quantitative agreement between the number of Li atoms deposited and the number of available Si bonds at the surface of the SiC crystal. Through STM analysis, we are able to determine the interlayer distance induced by Li-intercalation at the interface between the SiC substrate and the buffer layer.


title: Strongly aperiodic subshifts of finite type on hyperbolic groups

abstract: We prove that a hyperbolic group admits a strongly aperiodic subshift of finite type if and only if it has at most one end.


title: The Fierz convergence criterion: a controlled approach to strongly-interacting systems with small embedded clusters

abstract: We present an embedded-cluster method, based on the TRILEX formalism. It turns the Fierz ambiguity, inherent to approaches based on a bosonic decoupling of local fermionic interactions, into a convergence criterion. It is based on the approximation of the three-leg vertex by a coarse-grained vertex computed by solving a self-consistently determined multi-site effective impurity model. The computed self-energies are, by construction, continuous functions of momentum. We show that, in three interaction and doping regimes of parameters of the two-dimensional Hubbard model, self-energies obtained with clusters of size four only are very close to numerically exact benchmark results. We show that the Fierz parameter, which parametrizes the freedom in the Hubbard-Stratonovich decoupling, can be used as a quality control parameter. By contrast, the $GW$+extended dynamical mean field theory approximation with four cluster sites is shown to yield good results only in the weak-coupling regime and for a particular decoupling. Finally, we show that the vertex has spatially nonlocal components only at low Matsubara frequencies.


title: Mendelian Randomization when Many Instruments are Invalid: Hierarchical Empirical Bayes Estimation

abstract: Estimating the causal effect of an exposure on an outcome is an important task in many economical and biological studies. Mendelian randomization, in particular, uses genetic variants as instruments to estimate causal effects in epidemiological studies. However, conventional instrumental variable methods rely on some untestable assumptions, which may be violated in real problems. In this paper, we adopt a Bayesian framework and build hierarchical models to incorporate invalid effects of instruments. We introduce an empirical Bayes estimator for which some of the instruments are invalid by utilizing a Gaussian mixture prior. Theoretical performance and algorithm implementations are provided and illustrated. The reliable performance of the proposed method is demonstrated in various simulation settings and on real datasets concerning the causal effects of HDL cholesterol and LDL cholesterol on type 2 diabetes.


title: Nilpotent gelfand pairs and Schwartz extensions of spherical transforms via quotient pairs

abstract: It has been shown that for several nilpotent Gelfand pairs (N,K) (i.e., with N a nilpotent Lie group, K a compact group of automorphisms of N and the algebra L^1(N)^K commutative) the spherical transform establishes a 1-to-1 correspondence between the space S(N)^K of K-invariant Schwartz functions on N and the space S({\Sigma}) of functions on the Gelfand spectrum {\Sigma} of L^1(N)^K which extend to Schwartz functions on Rd, once {\Sigma} is suitably embedded in Rd. We call this property (S). We present here a general bootstrapping method which allows to establish property (S) to new nilpotent pairs (N,K), once the same property is known for a class of quotient pairs of (N,K) and a K-invariant form of Hadamard formula holds on N. We finally show how our method can be recursively applied to prove property (S) for a significant class of nilpotent Gelfand pairs.


title: Reversed Dickson polynomials of the $(k+1)$-th kind over finite fields, II

abstract: Let $p$ be an odd prime. In this paper, we study the permutation behaviour of the reversed Dickson polynomials of the $(k+1)$-th kind $D_{n,k}(1,x)$ when $n=p^{l_1}+3$, $n=p^{l_1}+p^{l_2}+p^{l_3}$, and $n=p^{l_1}+p^{l_2}+p^{l_3}+p^{l_4}$, where $l_1, l_2$, $l_3$, and $l_4$ are non-negative integers. A generalization to $n=p^{l_1}+p^{l_2}+\cdots +p^{l_i}$ is also shown. We find some conditions under which $D_{n,k}(1,x)$ is not a permutation polynomial over finite fields for certain values of $n$ and $k$. We also present a generalization of a recent result regarding $D_{p^l-1,1}(1,x)$ and present some algebraic and arithmetic properties of $D_{n,k}(1,x)$.


title: Border-Crossing Model for the Diffusive Coarsening of Two-Dimensional and Quasi-Two-Dimensional Wet Foams

abstract: For dry foams, the transport of gas from small high-pressure bubbles to large low-pressure bubbles is dominated by diffusion across the thin soap films separating neighboring bubbles. For wetter foams, the film areas become smaller as the Plateau borders and vertices inflate with liquid. So-called "border-blocking" models can explain some features of wet-foam coarsening based on the presumption that the inflated borders totally block the gas flux, however, this approximation dramatically fails in the wet/unjamming limit where the bubbles become close-packed spheres and coarsening proceeds even though there are no films. Here, we account for the ever-present border-crossing flux by a new length scale defined by the average gradient of gas concentration inside the borders. We compute that it is proportional to the geometric average of film and border thicknesses, and we verify this scaling by numerical solution of the diffusion equation. We similarly consider transport across inflated vertices and surface Plateau borders in quasi-2d foams. And we show how the $dA/dt=K_0(n-6)$ von~Neumann law is modified by the appearance of terms that depend on bubble size and shape as well as the concentration gradient length scales. Finally, we use the modified von~Neumann law to compute the growth rate of the average bubble area, which is not constant.


title: Jump Type Stochastic Differential Equations with Non-Lipschitz Coefficients: Non Confluence, Feller and Strong Feller Properties, and Exponential Ergodicity

abstract: This paper considers multidimensional jump type stochastic differential equations with super linear growth and non-Lipschitz coefficients. After establishing a sufficient condition for nonexplosion, this paper presents sufficient non-Lipschitz conditions for pathwise uniqueness. The non confluence property for solutions is investigated. Feller and strong Feller properties under non-Lipschitz conditions are investigated via the coupling method. Sufficient conditions for irreducibility and exponential ergodicity are derived. As applications, this paper also studies multidimensional stochastic differential equations driven by L\'evy processes and presents a Feynman-Kac formula for L\'evy type operators.


title: Multi-Observation Elicitation

abstract: We study loss functions that measure the accuracy of a prediction based on multiple data points simultaneously. To our knowledge, such loss functions have not been studied before in the area of property elicitation or in machine learning more broadly. As compared to traditional loss functions that take only a single data point, these multi-observation loss functions can in some cases drastically reduce the dimensionality of the hypothesis required. In elicitation, this corresponds to requiring many fewer reports; in empirical risk minimization, it corresponds to algorithms on a hypothesis space of much smaller dimension. We explore some examples of the tradeoff between dimensionality and number of observations, give some geometric characterizations and intuition for relating loss functions and the properties that they elicit, and discuss some implications for both elicitation and machine-learning contexts.


title: Evaporation-triggered microdroplet nucleation and the four life phases of an evaporating Ouzo drop

abstract: Evaporating liquid droplets are omnipresent in nature and technology, such as in inkjet printing, coating, deposition of materials, medical diagnostics, agriculture, food industry, cosmetics, or spills of liquids. While the evaporation of pure liquids, liquids with dispersed particles, or even liquid mixtures has intensively been studied over the last two decades, the evaporation of ternary mixtures of liquids with different volatilities and mutual solubilities has not yet been explored. Here we show that the evaporation of such ternary mixtures can trigger a phase transition and the nucleation of microdroplets of one of the components of the mixture. As model system we pick a sessile Ouzo droplet (as known from daily life - a transparent mixture of water, ethanol, and anise oil) and reveal and theoretically explain its four life phases: In phase I, the spherical cap-shaped droplet remains transparent, while the more volatile ethanol is evaporating, preferentially at the rim of the drop due to the singularity there. This leads to a local ethanol concentration reduction and correspondingly to oil droplet nucleation there. This is the beginning of phase II, in which oil microdroplets quickly nucleate in the whole drop, leading to its milky color which typifies the so-called 'Ouzo-effect'. Once all ethanol has evaporated, the drop, which now has a characteristic non-spherical-cap shape, has become clear again, with a water drop sitting on an oil-ring (phase III), finalizing the phase inversion. Finally, in phase IV, also all water has evaporated, leaving behind a tiny spherical cap-shaped oil drop.


title: ToPs: Ensemble Learning with Trees of Predictors

abstract: We present a new approach to ensemble learning. Our approach constructs a tree of subsets of the feature space and associates a predictor (predictive model) - determined by training one of a given family of base learners on an endogenously determined training set - to each node of the tree; we call the resulting object a tree of predictors. The (locally) optimal tree of predictors is derived recursively; each step involves jointly optimizing the split of the terminal nodes of the previous tree and the choice of learner and training set (hence predictor) for each set in the split. The feature vector of a new instance determines a unique path through the optimal tree of predictors; the final prediction aggregates the predictions of the predictors along this path. We derive loss bounds for the final predictor in terms of the Rademacher complexity of the base learners. We report the results of a number of experiments on a variety of datasets, showing that our approach provides statistically significant improvements over state-of-the-art machine learning algorithms, including various ensemble learning methods. Our approach works because it allows us to endogenously create more complex learners - when needed - and endogenously match both the learner and the training set to the characteristics of the dataset while still avoiding over-fitting.


title: Multi-Skyrmions on $AdS_2 \times S_2$, Rational maps and Popcorn Transitions

abstract: By combining two different techniques to construct multi-soliton solutions of the (3+1)-dimensional Skyrme model, the generalized hedgehog and the rational map ansatz, we find multi-Skyrmion configurations in $AdS_{2}\times S_{2}$. We construct Skyrmionic multi-layered configurations such that the total Baryon charge is the product of the number of kinks along the radial $AdS_{2}$ direction and the degree of the rational map. We show that, for fixed total Baryon charge, as one increases the charge density on $\partial\left( AdS_{2}\times S_{2}\right) $, it becomes increasingly convenient energetically to have configurations with more peaks in the radial $AdS_{2}$ direction but a lower degree of the rational map. This has a direct relation with the so-called holographic popcorn transitions in which, when the charge density is high, multi-layered configurations with low charge on each layer are favored over configurations with few layers but with higher charge on each layer. The case in which the geometry is $M_{2}\times S_{2}$ can also be analyzed.


title: Generalized magnetic mirrors

abstract: We propose generalized magnetic mirrors that can be achieved by excitations of sole electric resonances. Conventional approaches to obtain magnetic mirrors rely heavily on exciting the fundamental magnetic dipoles, whereas here we reveal that besides magnetic resonances, electric resonances of higher orders can be also employed to obtain highly efficient magnetic mirrors. Based on the electromagnetic duality, it is also shown that electric mirrors can be achieved by exciting magnetic resonances. We provide direct demonstrations of the generalized mirrors proposed in a simple system of one-dimensional periodic array of all-dielectric wires, which may shed new light to many advanced fields of photonics related to resonant multipolar excitations and interferences.


title: Language Generation with Recurrent Generative Adversarial Networks without Pre-training

abstract: Generative Adversarial Networks (GANs) have shown great promise recently in image generation. Training GANs for language generation has proven to be more difficult, because of the non-differentiable nature of generating text with recurrent neural networks. Consequently, past work has either resorted to pre-training with maximum-likelihood or used convolutional networks for generation. In this work, we show that recurrent neural networks can be trained to generate text with GANs from scratch using curriculum learning, by slowly teaching the model to generate sequences of increasing and variable length. We empirically show that our approach vastly improves the quality of generated sequences compared to a convolutional baseline.


title: Electronic and energetic properties of Ge(110) pentagons

abstract: The electronic and energetic properties of the elementary building block, i.e. a five-membered atom ring (pentagon), of the Ge(110) surface was studied by scanning tunneling microscopy and spectroscopy at room temperature. The Ge(110) surface is composed of three types of domains: two ordered domains ((16x2) and c(8x10)) and a disordered domain. The elementary building block of all three domains is a pentagon. Scanning tunneling spectra recorded on the (16x2), c(8x10) and disordered domains are very similar and reveal three well-defined electronic states. Two electronic states are located 1.1 eV and 0.3 eV below the Fermi level respectively, whereas the third electronic state is located 0.4 eV above the Fermi level. The electronic states at -0.3 eV and 0.4 eV can be ascribed to the pentagons, whilst we tentatively assigned the electronic state at -1.1 eV to a Ge-Ge back bond or trough state. In addition, we have analyzed the straight [1-12] oriented step edges. From the kink density and kink-kink distance distributions we extracted the nearest neighbor interaction energy between the pentagons, which exhibit a strong preference to occur in twins, as well as the strain relaxation energy along the pentagon-twin chains.


title: The scaling structure of the global road network

abstract: Because of increasing global urbanization and its immediate consequences, including changes in patterns of food demand, circulation and land use, the next century will witness a major increase in the extent of paved roads built worldwide. To model the effects of this increase, it is crucial to understand whether possible self-organized patterns are inherent in the global road network structure. Here, we use the largest updated database comprising all major roads on Earth, together with global urban and cropland inventories, to suggest that road length distributions within croplands are indistinguishable from urban ones, once rescaled to account for the difference in mean road length. Such similarity extends to road length distributions within urban or agricultural domains of given area. We find two distinct regimes for the scaling of the mean road length with the associated area, holding in general at small and at large values of the latter. In suitably large urban and cropland domains, we find that mean and total road lengths increase linearly with their domain area, differently from earlier suggestions. Scaling regimes suggest that simple and universal mechanisms regulate urban and cropland road expansion at the global scale. As such, our findings bear implications for global road infrastructure growth based on land-use change and for planning policies sustaining urban expansions.


title: $k_T$-factorization approach to the Higgs boson production in $ZZ^*\to 4\ell$ channel at the LHC

abstract: We calculated a differential cross section of the Higgs boson production in the $h\to ZZ^*\to 4\ell$ decay channel within the framework of $k_T$-factorization. Results are obtained using an off-shell matrix element for the $g^*g^*\to h\to ZZ^*\to 4\ell$ process together with Ciafaloni-Catani-Fiorani-Marchesini (CCFM) evolution equations for an unintegrated gluon distribution function. We have presented a comparison of our results with the latest experimental measurements at $\sqrt{S}$ = 8 TeV and $\sqrt{S}$ = 13 TeV from the ATLAS and CMS collaborations at the LHC. In addition to this, we have compared our results with the results from the collinear factorization formalism calculated up to next-to-next-to-leading order plus next-to-next-to-leading logarithm (NNLO + NNLL) accuracy obtained using the HRes code for the Higgs boson production in the gluon-gluon fusion process. Our estimates are consistently close to NNLO + NNLL results obtained using a collinear factorization formalism and are also in agreement with experimental measurements.


title: Sign-changing self-similar solutions of the nonlinear heat equation with positive initial value

abstract: We consider the nonlinear heat equation $u_t - \Delta u = |u|^\alpha u$ on ${\mathbb R}^N$, where $\alpha >0$ and $N\ge 1$. We prove that in the range $0 < \alpha <\frac {4} {N-2}$, for every $\mu >0$, there exist infinitely many sign-changing, self-similar solutions to the Cauchy problem with initial value $u_0 (x)= \mu |x|^{-\frac {2} {\alpha }}$. The construction is based on the analysis of the related inverted profile equation. In particular, we construct (sign-changing) self-similar solutions for positive initial values for which it is known that there does not exist any local, nonnegative solution.


title: High efficiency photonic storage of single photons in cold atoms

abstract: An efficient quantum storage is highly desired for quantum information processing. As indicated by certain applications, a universal quantum storage is required to have a storage efficiency above 50% to beat the no-cloning limit. Although significant progress has been achieved in improving various quantum storage, the best storage efficiency of single photons is still below this criteria. By integrating a highly controllable single photon source with an optimized quantum storage, here we demonstrate an optical storage of single photons with storage efficiency of 65% in a cold atomic ensemble based on electromagnetically induced transparency. Meanwhile, the nonclassical characteristics of our storage are verified through the well-maintained nonclassical and single photon nature of the retrieved single photons.


title: The Classification of Branched Willmore Spheres in the $3$-Sphere and the $4$-Sphere

abstract: We extend the classification of Robert Bryant of Willmore spheres in $S^3$ to variational branched Willmore spheres $S^3$ and show that they are inverse stereographic projections of complete minimal surfaces with finite total curvature in $\mathbb{R}^3$ and vanishing flux. We also obtain a classification of variational branched Willmore spheres in $S^4$, generalising a theorem of Seb\'{a}stian Montiel. As a result of our asymptotic analysis at branch points, we obtain an improved $C^{1,1}$ regularity of the unit normal of variational branched Willmore surfaces in arbitrary codimension. We also prove that the width of Willmore sphere min-max procedures in dimension $3$ and $4$, such as the sphere eversion, is an integer multiple of $4\pi$.


title: NullHop: A Flexible Convolutional Neural Network Accelerator Based on Sparse Representations of Feature Maps

abstract: Convolutional neural networks (CNNs) have become the dominant neural network architecture for solving many state-of-the-art (SOA) visual processing tasks. Even though Graphical Processing Units (GPUs) are most often used in training and deploying CNNs, their power efficiency is less than 10 GOp/s/W for single-frame runtime inference. We propose a flexible and efficient CNN accelerator architecture called NullHop that implements SOA CNNs useful for low-power and low-latency application scenarios. NullHop exploits the sparsity of neuron activations in CNNs to accelerate the computation and reduce memory requirements. The flexible architecture allows high utilization of available computing resources across kernel sizes ranging from 1x1 to 7x7. NullHop can process up to 128 input and 128 output feature maps per layer in a single pass. We implemented the proposed architecture on a Xilinx Zynq FPGA platform and present results showing how our implementation reduces external memory transfers and compute time in five different CNNs ranging from small ones up to the widely known large VGG16 and VGG19 CNNs. Post-synthesis simulations using Mentor Modelsim in a 28nm process with a clock frequency of 500 MHz show that the VGG19 network achieves over 450 GOp/s. By exploiting sparsity, NullHop achieves an efficiency of 368%, maintains over 98% utilization of the MAC units, and achieves a power efficiency of over 3TOp/s/W in a core area of 6.3mm$^2$. As further proof of NullHop's usability, we interfaced its FPGA implementation with a neuromorphic event camera for real time interactive demonstrations.


title: Towards a Flow- and Path-Sensitive Information Flow Analysis: Technical Report

abstract: This paper investigates a flow- and path-sensitive static information flow analysis. Compared with security type systems with fixed labels, it has been shown that flow-sensitive type systems accept more secure programs. We show that an information flow analysis with fixed labels can be both flow- and path-sensitive. The novel analysis has two major components: 1) a general-purpose program transformation that removes false dataflow dependencies in a program that confuse a fixed-label type system, and 2) a fixed-label type system that allows security types to depend on path conditions. We formally prove that the proposed analysis enforces a rigorous security property: noninterference. Moreover, we show that the analysis is strictly more precise than a classic flow-sensitive type system, and it allows sound control of information flow in the presence of mutable variables without resorting to run-time mechanisms.


title: Milne-like spacetimes and their role in Cosmology

abstract: Milne-like spacetimes are a class of FLRW models which admit $C^0$ spacetime extensions through the big bang. The boundary of a Milne-like spacetime can be identified with a null cone in the extension. We find that the comoving observers all emanate from a single point in the extension. This suggests that something physical may have happened there. Next we ask what role these spacetimes have in modern cosmology. We find that Milne-like spacetimes are inconsistent with classical cosmology theory in which a radiation-dominated era follows immediately after the big bang. However, Milne-like spacetimes are consistent with inflationary theory. In fact the very nature of Milne-like spacetimes solves the horizon problem. We note that this has already been observed for the Milne universe in [2] and [25]. Next we show, assuming a finite initial condition on the energy density, that Milne-like spacetimes naturally exhibit a slow-roll inflationary era. This makes Milne-like spacetimes good models for inflationary theory. We also show that Milne-like spacetimes still admit $C^0$ extensions under well-behaved perturbations.


title: On the spectra of cardinalities of branches of Kurepa trees

abstract: We are interested in the possible sets of cardinalities of branches of Kurepa trees in models of $ZFC$ $+$ $CH$. In this paper we present a sufficient condition (for sets of cardinals) to be consistently the set of cardinalities of branches of Kurepa trees.


title: Decay Estimates for 1-D Parabolic PDEs with Boundary Disturbances

abstract: In this work decay estimates are derived for the solutions of 1-D linear parabolic PDEs with disturbances at both boundaries and distributed disturbances. The decay estimates are given in the L2 and H1 norms of the solution and discontinuous disturbances are allowed. Although an eigenfunction expansion for the solution is exploited for the proof of the decay estimates, the estimates do not require knowledge of the eigenvalues and the eigenfunctions of the corresponding Sturm-Liouville operator. Examples show that the obtained results can be applied for the stability analysis of parabolic PDEs with nonlocal terms.


title: Effect of sea quarks on single-spin asymmetries $A^{W^{\pm}}_{N}$ in transversely polarized pp collisions at RHIC

abstract: We calculate the single-spin asymmetries $A^{W^{\pm}}_{N}$ of $W^{\pm}$ bosons produced in transversely polarized pp collisions with the valence part of the up (u) and down (d) quark Sivers functions treated by an available parametrization and the light-cone quark spectator-diquark model respectively, while the sea part Sivers functions of u and d quarks treated as parametrization. Comparing our results with those from experimental data at RHIC, we find that the Sivers functions of sea quarks play an important role in the determination of the shapes of $A^{W^{\pm}}_{N}$. It is shown that $A^{W^{-}}_{N}$ is sensitive to u sea Sivers function, while $A^{W^{+}}_{N}$ to d sea Sivers function intuitively. The results show that the contributions of u and d sea Sivers functions are rather sizable and of the same sign, and their signs agree with that of d valence quarks and are opposite to that of u valence quarks.


title: Boundary regularity for mean curvature flows of higher codimension

abstract: In this paper, we derive global bounds for the H\"older norm of the gradient of solutions of graphic mean curvature flow with boundary of arbitrary codimension.


title: Measuring the meteoroid environments of the planets with meteor detectors on Earth

abstract: We describe how meteors recorded at the Earth can be used to partly reconstruct the meteoroid environments of the planets if a large sample (i.e. millions of orbits at a minimum) is available. The process involves selecting from the Earth-based sample those meteors which passed near a planet's orbit prior to arriving at the Earth, and so carry information about the planetary meteoroid environment. Indeed this process can be extended to any location in the Solar System, though some regions of space are better sampled than others. From such a reconstruction performed with data from the Canadian Meteor Orbit Radar (CMOR), we reveal that Mars has apex, helion, anti-helion and toroidal sporadic sources, much as Earth does. Such reconstructions, albeit partial, have the potential to provide a wealth of detail about planetary meteoroid environments, and to allow for the ground-truthing of model meteoroid populations without in situ sampling.


title: A fast direct solver for boundary value problems on locally perturbed geometries

abstract: Many applications involve solving several boundary value problems on geometries that are local perturbations of an original geometry. The boundary integral equation for a problem on a locally perturbed geometry can be expressed as a low rank update to the original system. A fast direct solver for the new linear system is presented in this paper. The solution technique utilizes a precomputed fast direct solver for the original geometry to efficiently create the low rank factorization of the update matrix and to accelerate the application of the Sherman-Morrison formula. The method is ideally suited for problems where the local perturbation is the same but its placement on the boundary changes and problems where the local perturbation is a refined discretization on the same geometry. Numerical results illustrate that for fixed local perturbation the method is three times faster than building a new fast direct solver from scratch.


title: Higher-twist corrections to gluon TMD factorization

abstract: We calculate power corrections to TMD factorization for particle production by gluon-gluon fusion in hadron-hadron collisions.


title: The numerical initial boundary value problem for the generalized conformal field equations

abstract: In this paper we study a numerical implementation for the initial boundary value formulation for the generalized conformal field equations. We propose a formulation which is well suited for the study of the long-time behaviour of perturbed exact solutions such as a Schwarzschild or even a Kerr black hole. We describe the derivation of the implemented equations which we give in terms of the space-spinor formalism. We discuss the conformal Gauss gauge, and a slight generalization thereof which seems to be particularly useful in the presence of boundaries. We discuss the structure of the equations at the boundary and propose a method for imposing boundary conditions which allow the correct number of degrees of freedom to be freely specified while still preserving the constraints. We show that this implementation yields a numerically well-posed system by testing it on a simple case of gravitational perturbations of Minkowski space-time and subsequently with gravitational perturbations of Schwarzschild space-time.


title: A method for the online construction of the set of states of a Markov Decision Process using Answer Set Programming

abstract: Non-stationary domains, that change in unpredicted ways, are a challenge for agents searching for optimal policies in sequential decision-making problems. This paper presents a combination of Markov Decision Processes (MDP) with Answer Set Programming (ASP), named {\em Online ASP for MDP} (oASP(MDP)), which is a method capable of constructing the set of domain states while the agent interacts with a changing environment. oASP(MDP) updates previously obtained policies, learnt by means of Reinforcement Learning (RL), using rules that represent the domain changes observed by the agent. These rules represent a set of domain constraints that are processed as ASP programs reducing the search space. Results show that oASP(MDP) is capable of finding solutions for problems in non-stationary domains without interfering with the action-value function approximation process.


title: Learning Whenever Learning is Possible: Universal Learning under General Stochastic Processes

abstract: This work initiates a general study of learning and generalization without the i.i.d. assumption, starting from first principles. While the traditional approach to statistical learning theory typically relies on standard assumptions from probability theory (e.g., i.i.d. or stationary ergodic), in this work we are interested in developing a theory of learning based only on the most fundamental and necessary assumptions implicit in the requirements of the learning problem itself. We specifically study universally consistent function learning, where the objective is to obtain low long-run average loss for any target function, when the data follow a given stochastic process. We are then interested in the question of whether there exist learning rules guaranteed to be universally consistent given only the assumption that universally consistent learning is possible for the given data process. The reasoning that motivates this criterion emanates from a kind of optimist's decision theory, and so we refer to such learning rules as being optimistically universal. We study this question in three natural learning settings: inductive, self-adaptive, and online. Remarkably, as our strongest positive result, we find that optimistically universal learning rules do indeed exist in the self-adaptive learning setting. Establishing this fact requires us to develop new approaches to the design of learning algorithms. Along the way, we also identify concise characterizations of the family of processes under which universally consistent learning is possible in the inductive and self-adaptive settings. We additionally pose a number of enticing open problems, particularly for the online learning setting.


title: Asymptotically normal estimators for Zipf's law

abstract: Zipf's law states that sequential frequencies of words in a text correspond to a power function. Its probabilistic model is an infinite urn scheme with asymptotically power distribution. The exponent of this distribution must be estimated. We use the number of different words in a text and similar statistics to construct asymptotically normal estimators of the exponent.


title: Ipanema-\beta : tools and examples for HEP analysis on GPU

abstract: We present here a set of examples, classes and tools which can be used for statistical analysis in Graphics Processing Units (GPU). This includes binned and unbinned maximum likelihood fits, pseudo-experiment generation, convolutions, Markov Chain Monte Carlo method implementations, and limit setting techniques.


title: Large flares (M1 - X7) in solar cycle 24

abstract: The large (X-ray class > M1) and very large (X-ray class >X1) flares (according to the observations of GOES-15 and Preliminary data from Current Catalog of Flare Events) in solar activity cycle 24 were analyzed. The monthly average values of optical Flare Index for 2010 - 2016 were calculated. The values of the total energy of the flare E (J m-2) in the 0.1 -- 0.8 nm range at the level of the earth's atmosphere were estimated. The energy spectrum (the dependence of the number of flares with the full energy E from the value of this full energy) for 115 flares of M5 - X7 classes was built. The comparative study of monthly averaged values of several indices of solar activity in current cycle 24: the relative sunspot numbers (SSN), the 10.7 cm radio flux, the radiation flux in the Lyman-alpha line, the solar constant (TSI) and the Flare Index (FI) was made.


title: Atomic and Electronic Structure of Si Dangling Bonds in Quasi-Free-Standing Monolayer Graphene

abstract: Si dangling bonds without H termination at the interface of quasi-free standing monolayer graphene (QFMLG) are known scattering centers that can severely affect carrier mobility. In this report, we study the atomic and electronic structure of Si dangling bonds in QFMLG using low-temperature scanning tunneling microscopy/spectroscopy (STM/STS), atomic force microscopy (AFM), and density functional theory (DFT) calculations. Two types of defects with different contrast were observed on a flat terrace by STM and AFM. Their STM contrast varies with bias voltage. In STS, they showed characteristic peaks at different energies, 1.1 and 1.4 eV. Comparison with DFT calculations indicates that they correspond to clusters of 3 and 4 Si dangling bonds, respectively. The relevance of these results for the optimization of graphene synthesis is discussed.


title: Cladistical analysis of the Jovian and Saturnian satellite systems

abstract: Jupiter and Saturn each have complex systems of satellites and rings. These satellites can be classified into dynamical groups, implying similar formation scenarios. Recently, a larger number of additional irregular satellites have been discovered around both gas giants that have yet to be classified. The aim of this paper is to examine the relationships between the satellites and rings of the gas giants, using an analytical technique called cladistics. Cladistics is traditionally used to examine relationships between living organisms, the `tree of life'. In this work, we perform the first cladistical study of objects in a planetary science context. Our method uses the orbital, physical and compositional characteristics of satellites to classify the objects in the Jovian and Saturnian systems. We find that the major relationships between the satellites in the two systems, such as families, as presented in previous studies, are broadly preserved. In addition, based on our analysis of the Jovian system, we identify a new retrograde irregular family, the Iocaste family, and suggest that the Phoebe family of the Saturnian system can be further divided into two subfamilies. We also propose that the Saturnian irregular families be renamed, to be consistent with the convention used in Jovian families. Using cladistics, we are also able to assign the new unclassified irregular satellites into families. Taken together, the results of this study demonstrate the potential use of the cladistical technique in the investigation of relationships between orbital bodies.


title: A Partial Entropic Stabilization of Lattice Boltzmann MHD

abstract: The entropic lattice Boltzmann algorithm of Karlin et. al. is partially extended to magnetohydrodynamics, based on the Dellar model of introducing a vector distribution for the magnetic field. This entropic ansatz is now applied only to the scalar particle distribution function so as to permit the many problems entailing magnetic field reversal. A 9-bit lattice is employed for both particle and magnetic distributions for our two dimensional simulations. The entropic ansatz is benchmarked against our earlier multiple relaxation lattice-Boltzmann model for the Kelvin-Helmholtz instability in a magnetized jet. Other two dimensional simulations are performed and compared to results determined by more standard direct algorithms: in particular the switch over between the Kelvin-Helmholtz/tearing mode instability of Chen et. al., and the generalized Orszag-Tang vortex model of Biskamp-Welter. Very good results are achieved.


title: The validity of 21 cm spin temperature as a kinetic temperature indicator in atomic and molecular gas

abstract: The gas kinetic temperature ($T_K$) of various interstellar environments is often inferred from observations that can deduce level populations of atoms, ions, or molecules using spectral line observations; H I 21 cm is perhaps the most widely used with a long history. Usually the H I 21 cm line is assumed to be in thermal equilibrium and the populations are given by the Boltzmann distribution. A variety of processes, many involving Lyman alpha ($Ly\alpha$), can affect the 21 cm line. Here we show how this is treated in the spectral simulation code Cloudy, and present numerical simulations of environments where this temperature indicator is used, with a detailed treatment of the physical processes that determine level populations within $H^0$. We discuss situations where this temperature indicator traces $T_K$, cases where they fail, as well as the effects of $Ly\alpha$ pumping on the 21 cm spin temperature. We also show that the $Ly\alpha$ excitation temperature rarely traces the gas kinetic temperature.


title: Double Sparsity Kernel Learning with Automatic Variable Selection and Data Extraction

abstract: Learning with Reproducing Kernel Hilbert Spaces (RKHS) has been widely used in many scientific disciplines. Because a RKHS can be very flexible, it is common to impose a regularization term in the optimization to prevent overfitting. Standard RKHS learning employs the squared norm penalty of the learning function. Despite its success, many challenges remain. In particular, one cannot directly use the squared norm penalty for variable selection or data extraction. Therefore, when there exists noise predictors, or the underlying function has a sparse representation in the dual space, the performance of standard RKHS learning can be suboptimal. In the literature,work has been proposed on how to perform variable selection in RKHS learning, and a data sparsity constraint was considered for data extraction. However, how to learn in a RKHS with both variable selection and data extraction simultaneously remains unclear. In this paper, we propose a unified RKHS learning method, namely, DOuble Sparsity Kernel (DOSK) learning, to overcome this challenge. An efficient algorithm is provided to solve the corresponding optimization problem. We prove that under certain conditions, our new method can asymptotically achieve variable selection consistency. Simulated and real data results demonstrate that DOSK is highly competitive among existing approaches for RKHS learning.


title: A simple neural network module for relational reasoning

abstract: Relational reasoning is a central component of generally intelligent behavior, but has proven difficult for neural networks to learn. In this paper we describe how to use Relation Networks (RNs) as a simple plug-and-play module to solve problems that fundamentally hinge on relational reasoning. We tested RN-augmented networks on three tasks: visual question answering using a challenging dataset called CLEVR, on which we achieve state-of-the-art, super-human performance; text-based question answering using the bAbI suite of tasks; and complex reasoning about dynamic physical systems. Then, using a curated dataset called Sort-of-CLEVR we show that powerful convolutional networks do not have a general capacity to solve relational questions, but can gain this capacity when augmented with RNs. Our work shows how a deep learning architecture equipped with an RN module can implicitly discover and learn to reason about entities and their relations.


title: On the correspondence between thermodynamics and inference

abstract: We expand upon a natural analogy between Bayesian statistics and statistical physics in which sample size corresponds to inverse temperature. This analogy motivates the definition of two novel statistical quantities: a learning capacity and a Gibbs entropy. The analysis of the learning capacity, corresponding to the heat capacity in thermal physics, leads to new insight into the mechanism of learning and explains why some models have anomalously-high learning performance. We explore the properties of the learning capacity in a number of examples, including a sloppy model. Next, we propose that the Gibbs entropy provides a natural device for counting distinguishable distributions in the context of Bayesian inference. We use this device to define a generalized principle of indifference (GPI) in which every distinguishable model is assigned equal a priori probability. This principle results in a new solution to a long-standing problem in Bayesian inference: the definition of an objective or uninformative prior. A key characteristic of this new approach is that it can be applied to analyses where the model dimension is unknown and circumvents the automatic rejection of higher-dimensional models in Bayesian inference.


title: Analytical Solutions of Transient Drift-Diffusion in P-N Junction Pixel Sensors

abstract: Radiation detection in applications ranging from high energy physics to medical imaging rely on solid state detectors, often hybrid pixel detectors with (1) reverse biased p-n junction pixel sensors and (2) readout ASICs, attached by flip-chip-bonding. Transient signals characteristics are important in, e.g., matching ASIC and sensor design, modeling and optimizing detector parameters and describing timing and charge sharing properties. Currently analytical forms of transient signals are available for only a few limited cases (e.g., drift or diffusion) or for the steady state (which is not relevant for high energy radiation detection). Tools are available for (relatively slow) numerical evaluation of the transient charge transport. We present here the first analytical solutions of partial differential equations describing drift-diffusion-recombination charge transport in planar p-n junction sensors in a variety of conditions: (1) undepleted, (2) fully depleted, (3) taking into account the gradual velocity saturation, and (4) overdepleted. We deduce the Green's functions which can be applied to any detection problem through simple convolution with the initial conditions. We compare the analytical solutions with Monte Carlo simulations and industry standard simulations (Synopsys Sentaurus), demonstrating good agreement. Using the analytical equations enables fast modeling of the influence of various detector parameters on tracking, imaging and timing performance, describing performance and enabling optimizations for different applications. Finally, we illustrate this model with applications in 3D+T (x,y,z,time) photon tracking and 4D+T (x,y,theta,phi,time) relativistic charged particle tracking.


title: Toward a refining of the topological theory of phase transitions

abstract: The topological theory of phase transitions was proposed on the basis of different arguments, the most important of which are: a direct evidence of the relation between topology and phase transitions for some exactly solvable models; an explicit relation between entropy and topological invariants of certain submanifolds of configuration space, and, finally, two theorems stating that, for a wide class of physical systems, phase transitions should necessarily stem from topological changes of some submanifolds of configuration space. It has been recently shown that the $2D$ lattice $\phi^4$-model provides a counterexample that falsifies the mentioned theorems. On the basis of a numerical investigation, the present work indicates the way to overcome this difficulty: in spite of the absence of critical points of the potential in correspondence of the transition energy, the phase transition of this model stems from an asymptotic ($N\to\infty$) change of topology of the energy level sets.


title: Finite Groups with a Trivial Chermak-Delgado Subgroup

abstract: The Chermak-Delgado lattice of a finite group is a modular, self-dual sublattice of the lattice of subgroups of $G$. The least element of the Chermak-Delgado lattice of $G$ is known as the Chermak-Delgado subgroup of $G$. This paper concerns groups with a trivial Chermak-Delgado subgroup. We prove that if the Chermak-Delgado lattice of such a group is lattice isomorphic to a Cartesian product of lattices, then the group splits as a direct product, with the Chermak-Delgado lattice of each direct factor being lattice isomorphic to one of the lattices in the Cartesian product. We establish many properties of such groups and properties of subgroups in the Chermak-Delgado lattice. We define a CD-minimal group to be an indecomposable group with a trivial Chermak-Delgado subgroup. We establish lattice theoretic properties of Chermak-Delgado lattices of CD-minimal groups. We prove an extension theorem for CD-minimal groups, and use the theorem to produce twelve examples of CD-minimal groups, each having different CD lattices. Curiously, quasi-antichain $p$-group lattices play a major role in the author's constructions.


title: Strategic Equilibria in Queues with Dynamic Service Rate and Full Information

abstract: We consider the problem of customer equilibrium behavior of a single server Markovian queue with dynamic control of the service rate. Customers arrive according a Poisson procedure and the system administrator makes a service rate choice between a low and a high value according to a $T$-threshold dynamic service policy, where the decision for switching to the higher service rate is made when the number of customers exceeds T without any additional cost. We assume that customers are identical and they are making join decisions regarding the maximization of their expected net benefit, receiving a fixed reward for service completion and incurring a waiting cost. In addition, we consider the observable case of the model where customers are fully informed on the service policy and the queue length upon arrival.


title: Visual Interaction Networks

abstract: From just a glance, humans can make rich predictions about the future state of a wide range of physical systems. On the other hand, modern approaches from engineering, robotics, and graphics are often restricted to narrow domains and require direct measurements of the underlying states. We introduce the Visual Interaction Network, a general-purpose model for learning the dynamics of a physical system from raw visual observations. Our model consists of a perceptual front-end based on convolutional neural networks and a dynamics predictor based on interaction networks. Through joint training, the perceptual front-end learns to parse a dynamic visual scene into a set of factored latent object representations. The dynamics predictor learns to roll these states forward in time by computing their interactions and dynamics, producing a predicted physical trajectory of arbitrary length. We found that from just six input video frames the Visual Interaction Network can generate accurate future trajectories of hundreds of time steps on a wide range of physical systems. Our model can also be applied to scenes with invisible objects, inferring their future states from their effects on the visible objects, and can implicitly infer the unknown mass of objects. Our results demonstrate that the perceptual module and the object-based dynamics predictor module can induce factored latent representations that support accurate dynamical predictions. This work opens new opportunities for model-based decision-making and planning from raw sensory observations in complex physical environments.


title: Positron cooling and annihilation in noble gases

abstract: Understanding the dynamics of positron cooling in gases, including the fraction of positrons surviving to thermalisation, is critical for accurate interpretation of positron lifetime spectra, for the development of efficient positron cooling in traps and accumulators, and for a cryogenically cooled, ultra-high-energy-resolution, trap-based positron beam. Here, positron cooling and annihilation in noble gases is simulated using accurate scattering and annihilation cross sections calculated \emph{ab initio} with many-body theory. It is shown that a strikingly small fraction of positrons survive to thermalisation: $\sim$0.1 in He, $\sim$0 in Ne (due to cooling effectively stalling in the relatively deep momentum-transfer cross-section minimum), $\sim$0.15 in Ar, $\sim$0.05 in Kr and $\sim$0.01 in Xe. For Xe, the time-varying annihilation rate $\bar{Z}_{\rm eff}(\tau)$ is shown to be highly sensitive to the depletion of the distribution due to annihilation, conclusively explaining the long-standing discrepancy between gas-cell and trap-based measurements in Xe. The \emph{ab initio} calculations enable the first simultaneous probing of the energy dependence of the the scattering cross section and annihilation rate. Overall, the use of the accurate atomic data gives $\bar{Z}_{\rm eff}(\tau)$ in close agreement with experiment for all noble gases except Ne, the experiment for which is proffered to have suffered from incomplete knowledge of the fraction of positrons surviving to thermalisation and/or the presence of impurities.


title: Hamiltonian Monte Carlo Methods for Subset Simulation in Reliability Analysis

abstract: This paper studies a non-random-walk Markov Chain Monte Carlo method, namely the Hamiltonian Monte Carlo (HMC) method in the context of Subset Simulation used for structural reliability analysis. The HMC method relies on a deterministic mechanism inspired by Hamiltonian dynamics to propose samples following a target probability distribution. The method alleviates the random walk behavior to achieve a more effective and consistent exploration of the probability space compared to standard Gibbs or Metropolis-Hastings techniques. After a brief review of the basic concepts of the HMC method and its computational details, two algorithms are proposed to facilitate the application of the HMC method to Subset Simulation in structural reliability analysis. Next, the behavior of the two HMC algorithms is illustrated using simple probability distribution models. Finally, the accuracy and efficiency of Subset Simulation employing the two HMC algorithms are tested using various reliability examples. The supporting source code and data are available for download at (the URL that will become available once the paper is accepted).


title: Extraction of conformal data in critical quantum spin chains using the Koo-Saleur formula

abstract: We study the emergence of two-dimensional conformal symmetry in critical quantum spin chains on the finite circle. Our goal is to characterize the conformal field theory (CFT) describing the universality class of the corresponding quantum phase transition. As a means to this end, we propose and demonstrate automated procedures which, using only the lattice Hamiltonian $H = \sum_j h_j$ as an input, systematically identify the low-energy eigenstates corresponding to Virasoro primary and quasiprimary operators, and assign the remaining low-energy eigenstates to conformal towers. The energies and momenta of the primary operator states are needed to determine the primary operator scaling dimensions and conformal spins -- an essential part of the conformal data that specifies the CFT. Our techniques use the action, on the low-energy eigenstates of $H$, of the Fourier modes $H_n$ of the Hamiltonian density $h_j$. The $H_n$ were introduced as lattice representations of the Virasoro generators by Koo and Saleur [Nucl. Phys. B 426, 459 (1994)]. In this paper we demonstrate that these operators can be used to extract conformal data in a nonintegrable quantum spin chain.


title: Exploring the determinants of Bitcoin's price: an application of Bayesian Structural Time Series

abstract: Currently, there is no consensus on the real properties of Bitcoin. The discussion comprises its use as a speculative or safe haven assets, while other authors argue that the augmented attractiveness could end accomplishing money's functions that economic theory demands. This paper explores the association between Bitcoin's market price and a set of internal and external factors using Bayesian Structural Time Series Approach. I aim to contribute to the discussion by differentiating among several attractiveness sources and employing a method that provides a more flexible analytic framework that decompose each of the components of the time series, apply variable selection, include information on previous studies, and dynamically examine the behavior of the explanatory variables, all in a transparent and tractable setting. The results show that the Bitcoin price is negatively associated with a neutral investor's sentiment, gold's price and Yuan to USD exchange rate, while positively related to stock market index, USD to Euro exchange rate and variated signs among the different countries' search trends. Hence, I find that Bitcoin has mixed properties since still seems to act as a speculative, safe haven and a potential a capital flights instrument.


title: Some basic properties of bounded solutions of parabolic equations with p-Laplacian diffusion

abstract: We provide a detailed (and fully rigorous) derivation of several fundamental properties of bounded weak solutions to initial-value problems for general conservative 2nd-order parabolic equations with p-Laplacian diffusion and (arbitrary) bounded and integrable initial data.


title: Observation of photonic Fermi arcs in topological metamaterials

abstract: The discovery of topological phases introduces new perspectives and platforms for various interesting physics originally investigated in quantum context and then, on an equal footing, in classic wave systems, such as photonics, acoustics and mechanics. As a characteristic feature, nontrivial Fermi arcs, connecting between topologically distinct Fermi surfaces, play vital roles in the clarification of Dirac and Weyl semimetals and have been observed in quantum materials very recently. However, in classical systems, no direct experimental observation of Fermi arcs in the momentum space has been reported so far. Here, using near-field scanning measurements, we show the observation of a photonic Fermi arcs connecting topologically distinct bulk states in chiral hyperbolic metamaterials. To verify the topological nature of this system, we further observe backscattering-immune propagation of nontrivial surface wave across a three-dimension physical step. Our results demonstrate metamaterial approach towards topological photonics and offer a deeper understanding of topological phases in three-dimension classical systems.


title: Coherent inflationary dynamics for Bose-Einstein condensates crossing a quantum critical point

abstract: Quantum phase transitions, transitions between many-body ground states, are of extensive interest in research ranging from condensed matter physics to cosmology. Key features of the phase transitions include a stage with rapidly growing new order, called inflation in cosmology, followed by the formation of topological defects. How inflation is initiated and evolves into topological defects remains a hot debate topic. Ultracold atomic gas offers a pristine and tunable platform to investigate quantum critical dynamics. We report the observation of coherent inflationary dynamics across a quantum critical point in driven Bose-Einstein condensates. The inflation manifests in the exponential growth of density waves and populations in well-resolved momentum states. After the inflation stage, extended coherent dynamics is evident in both real and momentum space. We present an intuitive description of the quantum critical dynamics in our system and demonstrate the essential role of phase fluctuations in the formation of topological defects.


title: The extension of holomorphic functions on a non-pluriharmonic locus

abstract: Let $n \geq 4$ and let $\Omega$ be a bounded hyperconvex domain in $\mathbb{C}^{n}$. Let $\varphi$ be a negative exhaustive smooth plurisubharmonic function on $\Omega$. We show that any holomorphic function defined on a connected open neighborhood of the support of $(i\partial \overline{\partial}\varphi)^{n-3}$ can be extended to the holomorphic function on $\Omega$.


title: The Capacity of Private Information Retrieval from Byzantine and Colluding Databases

abstract: We consider the problem of single-round private information retrieval (PIR) from $N$ replicated databases. We consider the case when $B$ databases are outdated (unsynchronized), or even worse, adversarial (Byzantine), and therefore, can return incorrect answers. In the PIR problem with Byzantine databases (BPIR), a user wishes to retrieve a specific message from a set of $M$ messages with zero-error, irrespective of the actions performed by the Byzantine databases. We consider the $T$-privacy constraint in this paper, where any $T$ databases can collude, and exchange the queries submitted by the user. We derive the information-theoretic capacity of this problem, which is the maximum number of \emph{correct symbols} that can be retrieved privately (under the $T$-privacy constraint) for every symbol of the downloaded data. We determine the exact BPIR capacity to be $C=\frac{N-2B}{N}\cdot\frac{1-\frac{T}{N-2B}}{1-(\frac{T}{N-2B})^M}$, if $2B+T < N$. This capacity expression shows that the effect of Byzantine databases on the retrieval rate is equivalent to removing $2B$ databases from the system, with a penalty factor of $\frac{N-2B}{N}$, which signifies that even though the number of databases needed for PIR is effectively $N-2B$, the user still needs to access the entire $N$ databases. The result shows that for the unsynchronized PIR problem, if the user does not have any knowledge about the fraction of the messages that are mis-synchronized, the single-round capacity is the same as the BPIR capacity. Our achievable scheme extends the optimal achievable scheme for the robust PIR (RPIR) problem to correct the \emph{errors} introduced by the Byzantine databases as opposed to \emph{erasures} in the RPIR problem. Our converse proof uses the idea of the cut-set bound in the network coding problem against adversarial nodes.


title: Types of Cognition and its Implications for future High-Level Cognitive Machines

abstract: This work summarizes part of current knowledge on High-level Cognitive process and its relation with biological hardware. Thus, it is possible to identify some paradoxes which could impact the development of future technologies and artificial intelligence: we may make a High-level Cognitive Machine, sacrificing the principal attribute of a machine, its accuracy.


title: Berezinskii-Kosterlitz-Thouless transition and criticality of an elliptic deformation of the sine-Gordon model

abstract: We introduce and study the properties of a periodic model interpolating between the sine-- and the sinh--Gordon theories in $1+1$ dimensions. This model shows the peculiarities, due to the preservation of the functional form of their potential across RG flows, of the two limiting cases: the sine-Gordon, not having conventional order/magnetization at finite temperature, but exhibiting Berezinskii-Kosterlitz-Thouless (BKT) transition; and the sinh-Gordon, not having a phase transition, but being integrable. The considered interpolation, which we term as {\em sn-Gordon} model, is performed with potentials written in terms of Jacobi functions. The critical properties of the sn-Gordon theory are discussed by a renormalization-group approach. The critical points, except the sinh-Gordon one, are found to be of BKT type. Explicit expressions for the critical coupling as a function of the elliptic modulus are given.


title: Batched Large-scale Bayesian Optimization in High-dimensional Spaces

abstract: Bayesian optimization (BO) has become an effective approach for black-box function optimization problems when function evaluations are expensive and the optimum can be achieved within a relatively small number of queries. However, many cases, such as the ones with high-dimensional inputs, may require a much larger number of observations for optimization. Despite an abundance of observations thanks to parallel experiments, current BO techniques have been limited to merely a few thousand observations. In this paper, we propose ensemble Bayesian optimization (EBO) to address three current challenges in BO simultaneously: (1) large-scale observations; (2) high dimensional input spaces; and (3) selections of batch queries that balance quality and diversity. The key idea of EBO is to operate on an ensemble of additive Gaussian process models, each of which possesses a randomized strategy to divide and conquer. We show unprecedented, previously impossible results of scaling up BO to tens of thousands of observations within minutes of computation.